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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】fsolve求解问题
31/1返回列表
查看: 322  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

gg0520

铜虫 (初入文坛)

[交流] 【求助】fsolve求解问题

function xFsovle5  % Pn+M--->Dn+P1
clear all;
clc;
format short g;
global Ka Ki Kp Ktm Ktal Kth Kd Qti0 Qt0 Qp0 QH0 Mti Mt Mp MH T V Qv Qs
Ka=0.0013;Ki=108.85;Kp=108.85;Ktm=0.012;Ktal=0.12;Kth=4.8;Kd=0.001;  %速率常数
Qti0=3000;Qt0=10000;Qp0=817500000;QH0=4000000;                     Mti=193.9;                                                           Mt=114;                                                             Mp=42;                                                               MH=2;                                                                T=3600;                                                            
V=90;                                                               
Qv=4.608206;                                                         
Qs=0.00855230;                                                      
x0=[6.05e-006  1.1413e+003 120.5079 2.8374 0.0049237 3.6821 5.5035e+003 18.1943 1.3606e+004 2.0337e+007 0.0043327]'; [x,FVAL,EXITFLAG,output]=fsolve(@NonlinEqs,x0)     
Mn=42*(x(6)+x(9))/(x(5)+x(8))
Mw=42*(x(7)+x(10))/(x(6)+x(9))
PDI=Mw/Mn
%-------------------------------------
function f =NonlinEqs(x)
global Ka Ki Kp Ktm Ktal Kth Kd Qti0 Qt0 Qp0 QH0 Mti Mt Mp MH T V Qv Qs
f(1)=Ka*x(11)+(Kth*x(3)^0.5+Ktal*x(4))*x(5)-Ki*x(1)*x(2)-Kd*x(1)+(0-x(1)*Qs)/V; %
f(2)=-Ki*x(1)*x(2)-(Kp*x(2)+Ktm*x(2))*x(5)+(Qp0/(Mp*T)-x(2)*Qv)/V;
f(3)=-Kth*x(5)*x(3)^0.5+(QH0/(MH*T)-x(3)*Qv)/V;
f(4)=-Ktal*x(4)*x(5)+(Qt0/(Mt*T)-x(4)*Qs)/V;
f(5)=Ki*x(1)*x(2)-(Kth*x(3)^0.5+Ktal*x(4)+Kd)*x(5)+(0-x(5)*Qs)/V;
f(6)=Ki*x(1)*x(2)+Kp*x(2)*x(5)-(Kth*x(3)^0.5+Ktm*x(2)+Ktal*x(4)+Kd)*x(6)+Ktm*x(2)*x(5)+(0-x(6)*Qs)/V;
f(7)=Ki*x(1)*x(2)+2*Kp*x(2)*x(6)+Kp*x(2)*x(5)-(Kth*x(3)^0.5+Ktm*x(2)+Ktal*x(4)+Kd)*x(7)+Ktm*x(2)*x(5)+(0-x(7)*Qs)/V;
f(8)=(Kth*x(3)^0.5+Ktm*x(2)+Ktal*x(4)+Kd)*x(5)+(0-x(8)*Qs)/V;
f(9)=(Kth*x(3)^0.5+Ktm*x(2)+Ktal*x(4)+Kd)*x(6)+(0-x(9)*Qs)/V;
f(10)=(Kth*x(3)^0.5+Ktm*x(2)+Ktal*x(4)+Kd)*x(7)+(0-x(10)*Qs)/V;
f(11)=-Ka*x(11)+(1.35/T-x(11)*Qs)/V;
-------------------------------------------------------------------------------
结果为：
Maximum number of function evaluations reached:
increase options.MaxFunEvals.

x =

    0.0033832
       1141.3
       11.212
       2.7032
     0.002532
      0.66931
       5503.5
       18.194
        13606
  2.0337e+007
    0.0043327


FVAL =

  Columns 1 through 7

      -420.24      -733.23        5.558  -0.00080748       420.24       714.72       1409.1

  Columns 8 through 11

     0.074467       18.849  1.6369e+005 -1.8776e-006


EXITFLAG =

     0


output =

       iterations: 91
        funcCount: 1104
        algorithm: 'trust-region dogleg'
    firstorderopt: 6.4012e+008
          message: [1x78 char]


Mn =

        31405


Mw =

        62792


PDI =

       1.9994
感觉程序就根本没有运行一样，初值基本不变，特别是最后几个，无论你怎么改变他们的值，都不会变，而是相应的改变前面的值，请问这是为什么~谢谢~
回复此楼

» 猜你喜欢
1楼 2009-06-26 09:46:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mininurse

金虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
kuhailangyu(金币+3,VIP+0):欢迎常来！ 6-26 14:33
1）这个解并没有收敛，所以exitflag=0，应该增加options.MaxFunEvals
2）这么多个参数在非线性方程中作优化，本身就是挺难获得最优解的，建议改进算法，考虑其它办法，而不是只是依赖于这个函数进行求解
赞一下(11人) 回复此楼
2楼2009-06-26 14:04:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你是在做参数估值吧，用lsqnonlin  (L-M法)，lsqunicon（单纯形法）求解
赞一下(1人) 回复此楼
3楼2009-06-28 09:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gg0520 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否申诉？ +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 中青基了要发朋友圈吗？ +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭