查看: 407  |  回复: 4
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 whiteman 的 10 个金币

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函数求极小值
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):恭喜抢沙发,给个红包
whiteman(金币+10,VIP+0):能不能写具体点,我是新手 。如果正确我可以追加金币200个 1-8 10:02
把M函数的参数改成2个,即f=c(p,u1)
u1=0.012就不要写进M函数了

在主程序中
for u1=-0.015:0.001:0.015
....
  [x,f]=fminbnd('c',0,1);
  disp([x f]);
.....
end
板凳要做十年冷文章不发一个字
2楼2010-01-07 22:21:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whiteman

铜虫 (正式写手)

能不能写具体点,我是新手
3楼2010-01-08 10:01:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫


formleaf(金币+1,VIP+0):知道兄弟不差钱,只是表示荣誉 1-8 13:48
哥们,你看我却却金币吗,呵呵


****************M函数*****************
function f=c(p,u1);
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;
****************************************************

************************出程序***********************

p=。。。

for u1=-0.015:0.001:0.015

  找到min f(p,u1)

end
****************************



思路是这样的,具体自己去调试呵呵

[ Last edited by bluesine on 2010-1-8 at 11:21 ]
板凳要做十年冷文章不发一个字
4楼2010-01-08 11:13:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whiteman

铜虫 (正式写手)

??? Input argument "u1" is undefined.

Error in ==> c at 4
T=450;u1;u2=-u1;

Error in ==> fminbnd at 176
fa = funfcn(ax,varargin{:});

不能运行
5楼2010-01-08 11:19:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whiteman 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见