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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 信息科学 » 图像处理 » 求解优化方程
361/4返回列表1234下一页
查看: 1939  |  回复: 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.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 7/350 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中，邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗？ +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！ +4 SCSIOyxguo 2026-09-03 7/350 2026-09-05 10:39 by Honglianlian
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
关闭
请填处理意见
关闭