24小时热门版块排行榜    

查看: 1916  |  回复: 35

rzsyh

铁虫 (小有名气)

[求助] 求解优化方程 已有1人参与

求解一优化方程:x=argmin||xA-B||22+d||deltx||22
22为二范数的平方,deltx表示x的梯度,A和B各位图像矩阵,d为一参数,求解满足xA=B,且x最平滑,的解3。
新虫,钱很少。哪位看到的网友给以指点??跪谢。
回复此楼

» 收录本帖的淘帖专辑推荐

五味杂陈

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

catbin

银虫 (小有名气)

你这里x是向量还是矩阵?若x是矩阵,你确定第一项的范数是矩阵2范数而不是Frobenius范数?后者才有解析解。
静水流深
3楼2014-10-11 06:47:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

http://www.mathworks.cn/matlabce ... variation-denoising
程序里有使用Kronecker内积实现算子的代码
19楼2014-10-16 14:57:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

引用回帖:
19楼: Originally posted by yongcailiu at 2014-10-16 14:57:33
http://www.mathworks.cn/matlabcentral/fileexchange/36278-split-bregman-method-for-total-variation-denoising
程序里有使用Kronecker内积实现算子的代码

找到公式了,详见图片
求解优化方程
(W{7G}RY[ZI)VAV$S1~NALN.png

23楼2014-10-22 15:38:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

生医小工程

银虫 (初入文坛)

公式看着麻烦,你还是截图上传的好些。不过我也不会,我是来支持的。
2楼2014-10-10 20:07:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rzsyh

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by catbin at 2014-10-11 06:47:45
你这里x是向量还是矩阵?若x是矩阵,你确定第一项的范数是矩阵2范数而不是Frobenius范数?后者才有解析解。

x是矩阵吧,第一项的范数应该是矩阵2范数呀。
求解优化方程-1
无标题.png

4楼2014-10-11 09:18:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rzsyh

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by 生医小工程 at 2014-10-10 20:07:16
公式看着麻烦,你还是截图上传的好些。不过我也不会,我是来支持的。

原方程如图:
求解优化方程-2
无标题.png

5楼2014-10-11 09:19:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rzsyh

铁虫 (小有名气)

会者不难,难者不会。
人只说用共轭梯度就解出来了,会的虫友们,教教俺呗。
6楼2014-10-11 09:49:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
6楼: Originally posted by rzsyh at 2014-10-11 09:49:23
会者不难,难者不会。
人只说用共轭梯度就解出来了,会的虫友们,教教俺呗。

查看最优函数的first order condition, 将其对\beta求导令其等于0,问题转化为线性方程组的求解,然后用cg求解。
7楼2014-10-11 10:28:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rzsyh

铁虫 (小有名气)

引用回帖:
7楼: Originally posted by yongcailiu at 2014-10-11 10:28:24
查看最优函数的first order condition, 将其对\beta求导令其等于0,问题转化为线性方程组的求解,然后用cg求解。...

能再详细一点吗?我不懂优化啊
8楼2014-10-11 11:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongcailiu

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
rzsyh: 金币+5, ★★★很有帮助, 非常感谢 2014-10-13 09:10:53
引用回帖:
8楼: Originally posted by rzsyh at 2014-10-11 11:28:16
能再详细一点吗?我不懂优化啊...

求解极值问题

$\min_\beta \|\beta \ln H-\ln J\|_2^2+\lambda \|\nabla \beta\|_2^2$

对$\beta$求导,并令其等于零,即

$2\ln H (\beta \ln H-\ln J)+2\lambda \nabla^{*}(\nabla \beta) = 0$

消去系数,并利用关系$\nabla^{*}\nabla=-\triangle$,得到

$\ln H (\beta \ln H-\ln J)-\lambda \triangle \beta = 0$

接着就是求解关于$\beta$的方程组了。推导过程可参见《Total Variation Denoising using Split
Bregman》一文中(46)-(47)式,具体符号的离散计算参见原文。

小木虫里贴图我不会,上面是latex里的整理内容。
9楼2014-10-11 15:26:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

catbin

银虫 (小有名气)

引用回帖:
9楼: Originally posted by yongcailiu at 2014-10-11 15:26:25
求解极值问题

$\min_\beta \|\beta \ln H-\ln J\|_2^2+\lambda \|\nabla \beta\|_2^2$

对$\beta$求导,并令其等于零,即

$2\ln H (\beta \ln H-\ln J)+2\lambda \nabla^{*}(\nabla \beta) = 0$

消去系 ...

Euler-Lagrange Equation里面第一项括号外面的\ln H少了个转置。

之前楼主楼主回我的帖说,变量都为矩阵,第一项的范数是矩阵2范数。矩阵2范数的变分不是你给出的这个,要复杂一些需要用到奇异值分解来定义。但我怀疑楼主给的模型是个普通的带正则项的最小二乘问题,要么变量都本质上是向量,或者第一个范数是Frobenius范数。
静水流深
10楼2014-10-12 10:51:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rzsyh 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 20:50 by xiangy672
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见