| 查看: 400 | 回复: 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函数求极小值 |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有17人回复
-
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
2楼2010-01-07 22:01:08
★ ★ ★
小木虫(金币+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
30