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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab中fmincon的fun是矩阵运算怎么办？
61/1返回列表
查看: 2459  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

goodjackzj

新虫 (小有名气)

[求助] matlab中fmincon的fun是矩阵运算怎么办？ 已有1人参与

想用优化算法求解一个矩阵运算式子的最小值，上网看了有关fmincon的资料，发现fmincon中的fun式子一般都是能写出x1，x2表达式的，比如
function y = fun1(x)
y = x(1)-0.25*x(2)+0.6;
如果我的fun是一个矩阵运算如：
min （h2-G2*phi2)'*W2*(h2-G2*phi2)
subject to: phi2>0
其中h2,G2,W2已知，phi2是一个3*1的向量，无法写出包含phi2（i）的具体式子，请问应该怎么办？
​小弟初学优化算法，还请各位大虾多多指教~~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-30 09:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小米2011

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
goodjackzj: 金币+5 2014-01-10 18:24:54
xzhdty: 金币+1, 谢谢参与 2014-01-10 19:52:33
矩阵运算也可以的。fmincon可以求 scalar, vector, or matrix的解。具体操作时无需写成 phi2(i)的形式，就将phi2作为一个vector就行：
function y = fun1(phi2)
y=（h2-G2*phi2)'*W2*(h2-G2*phi2);
赞一下(1人) 回复此楼
2楼2014-01-02 02:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goodjackzj

新虫 (小有名气)
引用回帖:
2楼: Originally posted by 小米2011 at 2014-01-02 02:37:28
矩阵运算也可以的。fmincon可以求 scalar, vector, or matrix的解。具体操作时无需写成 phi2(i)的形式，就将phi2作为一个vector就行：
function y = fun1(phi2)
y=（h2-G2*phi2)'*W2*(h2-G2*phi2);

请问我这样写对吗？
f = @(x)(h2-G2*x)'*W2*(h2-G2*x);
A = (-1)*eye(3);
B=[0,0,0]';
x0=[1,1,1]';
phi2=fmincon(f,x0,A,B);
初始点x0随便设的，我想求phi2的全局最小点，请问请应该怎么做呢？
赞一下 回复此楼
3楼2014-01-05 23:06:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小米2011

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by goodjackzj at 2014-01-05 16:06:06
请问我这样写对吗？
f = @(x)(h2-G2*x)'*W2*(h2-G2*x);
A = (-1)*eye(3);
B=';
x0=';
phi2=fmincon(f,x0,A,B);
初始点x0随便设的，我想求phi2的全局最小点，请问请应该怎么做呢？...

就从你这些语句来看，我觉得你写的是对的。
赞一下 回复此楼
4楼2014-01-06 23:18:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goodjackzj

新虫 (小有名气)
引用回帖:
4楼: Originally posted by 小米2011 at 2014-01-06 23:18:48
就从你这些语句来看，我觉得你写的是对的。...

已经解决了，是我问题没有描述清楚。
function mf=min(phi2)
h2=[];   %你求出的值
G2=[];   %你求出的值
W2=[];   %你求出的值
mf=(h2-G2*phi2)'*W2*(h2-G2*phi2);

end

fun=@min;
A=[];
b=[];
Aeq=[];
beq=[];
lb=[0 0 0];
ub=[];
nonlcon=[];
x0=phi2;   %初始值
options=optimset('Algorithm','interior-point');

[phi2_min,fval]=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options);
赞一下 回复此楼
5楼2014-01-10 18:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onlymindyou

铁杆木虫 (著名写手)
引用回帖:
5楼: Originally posted by goodjackzj at 2014-01-10 18:24:06
已经解决了，是我问题没有描述清楚。
function mf=min(phi2)
h2=[];   %你求出的值
G2=[];   %你求出的值
W2=[];   %你求出的值
mf=(h2-G2*phi2)'*W2*(h2-G2*phi2);

end

fun=@min;
A=[];
b=[];
Aeq= ...

我遇到的也是变量是矩阵。请问您这个有function的是直接粘进去运行么？为啥我的报错啊，Function definitions are not permitted in this context.
赞一下 回复此楼
6楼2016-02-02 21:51:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 goodjackzj 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 求调剂 +5 一定有学上- 2026-03-12 5/250 2026-03-13 18:31 by ms629
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
关闭
请填处理意见
关闭