| 查看: 269 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
whiteman铜虫 (正式写手)
|
[交流]
【求助】200金币求如何在matlab里实现循环求极小值？
|
||
|
如M函数： function f=c(p); z=p(1); T=377;u1=0.012;u2=u1; a1=(4.273e+5)*(T-377.0); a11=-6.36e+8; a111=2.81e+9; a1111=1.74e+10; Q11=0.12;Q12=-0.053;Q44=0.052;s11=4.6e-12;s12=-1.1e-12;s44=11.1e-12; aa3=a1-(u1+u2)*Q12/(s11+s12); aa33=a11+Q12^2/(s11+s12); f=aa3*z^2+aa33*z^4+a111*z^6+a1111*z^8； 我现在只会在m函数中改变u1然后在命令行输入[x,f]=fminbnd('c',0,1)求不同u1下的极小值。怎样才能在Matlab上自动求u1=-0.015，-0.014，。。。0.015函数求极小值 |
回复此楼» 猜你喜欢
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
要不要辞职读博？
已经有6人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+2,VIP+0):谢谢应助！ 1-8 09:26
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+2,VIP+0):谢谢应助！ 1-8 09:26
|
看来起楼主的u1可能的取值也就31个啊，再循环一下也就循环了 事先声明，我这是个笨办法。 fminbnd要求单参数的函数，但可以接受anonymous function，帮助文档里面就有一个例子： x = fminbnd(@(x) sin(x*x),x1,x2); 楼主这里可以这样，不详细写了，大概可以看出个意思 for u1=-0.015:0.001:0.015 [x, f] = fminbun(@(x) 含x和u1的函数表达式, 0, 1); % display result disp([u1, x, f]); end |
3楼2010-01-08 00:28:07