| 查看: 301 | 回复: 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函数求极小值 |
回复此楼» 猜你喜欢
071000生物学求调剂，初试成绩343
已经有7人回复
-
0703化学调剂，求导师收
已经有8人回复
-
0856材料化工调剂 总分330
已经有11人回复
-
本科新能源科学与工程，一志愿华理能动285求调剂
已经有9人回复
-
332求调剂
已经有4人回复
-
085600，专业课化工原理，321分求调剂
已经有3人回复
-
材料工程294求调剂
已经有3人回复
-
295材料工程专硕求调剂
已经有7人回复
-
305求调剂
已经有5人回复
-
食品工程专硕一志愿中海洋309求调剂
已经有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
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51005.2
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
2楼2010-01-07 22:01:08
