24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 whiteman 的 10 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

whiteman

铜虫 (正式写手)

应助: 0 (幼儿园)
金币: 1672.9
散金: 3101
红花: 1
帖子: 903
在线: 1355.3小时
虫号: 470154
注册: 2007-11-29
专业: 功能陶瓷

[交流] 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函数求极小值

回复此楼

» 猜你喜欢

1楼 2010-01-07 19:32:41

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

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9800.3
散金: 89
红花: 19
帖子: 3609
在线: 377小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

★
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 ]