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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 图像处理 » 求解优化方程
361/4返回列表1234下一页
查看: 1943  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +7 QUjhNVAcOSff 2026-09-13 7/350 2026-09-16 14:36 by i3KBRD19ERkG
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-16 06:51 by LH5NkK5Ud4bh
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:36 by LH5NkK5Ud4bh
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 06:11 by LH5NkK5Ud4bh
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 04:51 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-16 04:23 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-16 03:38 by DgNGHc3h5tPl
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:35 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:26 by zNcEhTcH7ihq
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-15 19:35 by CgyNCDVNhGVg
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:59 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-15 14:49 by VMDqgnlKimtZ
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
信息提示
关闭
请填处理意见
关闭