版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

goodjackzj

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 188.3
帖子: 53
在线: 10.2小时
虫号: 2214639
注册: 2012-12-30
专业: 信号理论与信号处理

[求助] 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）的具体式子，请问应该怎么办？
小弟初学优化算法，还请各位大虾多多指教~~

回复此楼

» 猜你喜欢

读博已经有5人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有13人回复
博士申请都是内定的吗？已经有6人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有5人回复
博士读完未来一定会好吗已经有29人回复
投稿精细化工已经有4人回复
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！已经有4人回复
导师想让我从独立一作变成了共一第一已经有9人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复

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

matlab中如何建立符号矩阵？已经有12人回复
我用matlab编写了一个风力机叶片的优化设计程序，有点问题已经有13人回复
MATLAB中自定义含参量的矩阵已经有5人回复
matlab简单矩阵运算问题请教已经有7人回复
matlab数值计算时迭代项分母为零的时候怎么办？已经有15人回复
matlab中带一行均为0的矩阵，如何输入？已经有6人回复
求教matlab线性优化求最小值的问题已经有7人回复
.dat文件在matlab中直接应用吗，以矩阵形式已经有12人回复
matlab中如何随机不重复选取矩阵中的元素已经有4人回复
matlab中如何通过for语句得到已有矩阵的部分数据已经有3人回复
菜鸟求助~matlab怎么样让一个1xn的向量中的每个元素是个矩阵？已经有11人回复
matlab中建立稀疏矩阵组并操作已经有4人回复
matlab fmincon优化函数怎样知道每次得到的迭代点已经有3人回复
matlab中系数是矩阵的积分已经有5人回复
【求助】matlab矩阵中寻找距离最近的行向量已经有4人回复
小弟新学matlab，想借用 fmincon求个函数的最小值，求助错误该法，感谢！已经有4人回复

1楼2013-12-30 09:27:22

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

小米2011

新虫 (初入文坛)

应助: 7 (幼儿园)
金币: 205.4
帖子: 27
在线: 69.5小时
虫号: 1236969
注册: 2011-03-18
性别: MM
专业: 信息处理方法与技术

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 188.3
帖子: 53
在线: 10.2小时
虫号: 2214639
注册: 2012-12-30
专业: 信号理论与信号处理

引用回帖:

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

新虫 (初入文坛)

应助: 7 (幼儿园)
金币: 205.4
帖子: 27
在线: 69.5小时
虫号: 1236969
注册: 2011-03-18
性别: MM
专业: 信息处理方法与技术

引用回帖:

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

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 188.3
帖子: 53
在线: 10.2小时
虫号: 2214639
注册: 2012-12-30
专业: 信号理论与信号处理

引用回帖:

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

铁杆木虫 (著名写手)

应助: 0 (幼儿园)
金币: 7188.1
红花: 1
帖子: 1186
在线: 234.1小时
虫号: 1644760
注册: 2012-02-26
专业: 化学动力学

引用回帖:

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 的主题更新

返回列表