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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】fsolve求解问题
31/1返回列表
查看: 323  |  回复: 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
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
关闭
请填处理意见
关闭