应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 图像处理 » 求解优化方程
361/4返回列表1234下一页
查看: 1972  |  回复: 35
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

rzsyh

铁虫 (小有名气)

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

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

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

五味杂陈

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-10-10 15:52:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

catbin

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

yongcailiu

金虫 (小有名气)
http://www.mathworks.cn/matlabce ... variation-denoising
程序里有使用Kronecker内积实现算子的代码
赞一下(1人) 回复此楼
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
赞一下(1人) 回复此楼
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 的主题更新
361/4返回列表1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯，如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭