应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 200金币求如何在matlab里实现循环求极小值？
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 425  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 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函数求极小值
回复此楼

» 猜你喜欢
1楼 2010-01-07 19:32:41
已阅   回复此楼   关注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
赞一下(6人) 回复此楼
板凳要做十年冷文章不发一个字
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 ]
赞一下(11人) 回复此楼
板凳要做十年冷文章不发一个字
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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +20 Y-cap 2026-03-29 25/1250 2026-04-01 23:49 by 欣喜777
[考研] 274求调剂 +4 薛定谔的虎。 2026-04-01 4/200 2026-04-01 22:56 by wxiongid
[考研] 295材料工程专硕求调剂 +19 1428151015 2026-03-27 19/950 2026-04-01 22:34 by peike
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 290求调剂 +5 dfffsar 2026-03-29 5/250 2026-04-01 19:45 by 6781022
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 314求调剂 +6 1xiaojun23 2026-03-31 6/300 2026-04-01 11:28 by 今天星期八LQ
[考研] 326求调剂 +4 崽崽仔 2026-03-31 4/200 2026-04-01 09:58 by 我的船我的海
[考研] 一志愿南师大0703化学 275求调剂 +5 Ripcord上岸 2026-03-27 5/250 2026-03-31 19:52 by mg1014
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 085602化工求调剂（331分） +8 111@127 2026-03-30 8/400 2026-03-30 21:23 by 研究僧导导
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 21:07 by peike
[考研] 304材料求调剂 +5 钟llll 2026-03-26 5/250 2026-03-30 19:37 by 源_2020
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +5 贝呗钡钡 2026-03-30 5/250 2026-03-30 18:16 by 无际的草原
[考研] 一志愿双一流机械285分求调剂 +4 幸运的三木 2026-03-29 5/250 2026-03-29 14:49 by Miko19
[考研] 数一英一271专硕（085401）求调剂，可跨 +7 前行必有光 2026-03-28 8/400 2026-03-28 23:22 by 小木虫tim
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 266求调剂 +11 阳阳哇塞 2026-03-27 12/600 2026-03-27 17:56 by yu221
[考研] 调剂求收留 +7 果然有我 2026-03-26 7/350 2026-03-27 00:26 by wxiongid
信息提示
关闭
请填处理意见
关闭