 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 744 | 回復: 3 | ||
中將木蟲 (小有名氣)
|
[求助]
求高手幫我分析下我這個matlab求解非線性微分方程組解的情況
|
|
我這個是一組非線性方程組，可是每次給不同的初值都達不到我要求結果，不知道怎么辦，求大神指教，不知道是方法不對還是怎么的，我希望結果輸出的都是正數(shù)在0和1之間，而且最后的total要求最好小于1。。。。。。。。這樣的非線性方程組怎么求好呢，有什么好的方法呢，求教求教 這個是我函數(shù)方程 function F=fun(x) ya0=0.05;ye0=0.2;yf0=0.75; Kp1=6.6090e9;Kp2=568.35;Kp3=2.4474e5;Kp5=6.6946e16;Kp6=3.3131e-3;Kp8=1.4006e13; F=[Kp1*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))-x(2)*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6)); Kp2*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))- x(2)*x(3); Kp3*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))^2-x(3)^3; Kp5*((1+2/3*x(2)-2*x(5)-2*x(4))*ya0-1/3*(3*x(1)+2*x(2)+x(3)+x(5)+x(4)+x(6)))*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6))^3-10e6*x(5)^3*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))^7; Kp6*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))*x(5)-x(4)*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6)); Kp8*x(3)*((1+2/3*x(2)-2*x(5)-2*x(4))*ye0-1/3*(3*x(1)+x(2)+2*x(3)-7*x(5)+5*x(6)-4*x(4)))-x(6)*((1+2/3*x(2)-2*x(5)-2*x(4))*yf0+x(1)-x(5)+x(6)); ]; 下面是調用 x0=[0.1;0.1;0.1;0.1;0.1;0.1]; options=optimset('Display','iter'); y=fsolve(@fun,x0,optimset('fsolve')) yA=((1+2/3*y(2)-2*y(5)-2*y(4))*0.1-1/3*(3*y(1)+2*y(2)+y(3)+y(5)+y(4)+y(6))) yE=((1+2/3*y(2)-2*y(5)-2*y(4))*0.5-1/3*(3*y(1)+y(2)+2*y(3)-7*y(5)+5*y(6)-4*y(4))) yF=((1+2/3*y(2)-2*y(5)-2*y(4))*0.4+y(1)-y(5)+y(6)) > total=y(1)+y(2)+y(3)+y(4)+y(5)+y(6)+yA+yE+yF |
回復此樓» 猜你喜歡
調劑
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有86人回復
-
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職
已經有0人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。
已經有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人
已經有0人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- matlab如何求解一個非線性微分方程組
已經有8人回復
- 求助，用matlab解微分方程組，希望高手能給予指點
已經有17人回復
- matlab解偏微分方程求助
已經有12人回復
- 求高手解答Matlab求解微分方程的初始值和邊界值問題
已經有5人回復
- 二階非線性微分方程求解
已經有14人回復
- 請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。
已經有17人回復
- 求Matlab高手解決線性方程組的迭代求解問題
已經有17人回復
- 請教用matlab求解一個非線性偏微分方程組的數(shù)值解
已經有8人回復
- 高手進來用matlab求解常微分方程，十分感謝！
已經有17人回復
- matlab有限差分法求解復雜的偏微分方程組
已經有19人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經有4人回復
- Matlab求解偏微分方程組
已經有22人回復
- matlab求解非線性方程組
已經有16人回復
- matlab數(shù)值求解邊界條件微分方程組
已經有7人回復
- 高金求助matlab解微分方程組
已經有12人回復
- matlab解微分方程組
已經有15人回復
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經有11人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經有21人回復
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經有22人回復
- 【求助】積分微分方程matlab求解
已經有6人回復
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經有13人回復
中將
木蟲 (小有名氣)
- 應助: 2 (幼兒園)
- 金幣: 1958.7
- 散金: 86
- 紅花: 1
- 帖子: 221
- 在線: 76小時
- 蟲號: 886327
- 注冊: 2009-10-28
- 性別: GG
- 專業(yè): 分離過程
2樓2013-06-22 16:52:34
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應助回帖
★
感謝參與，應助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-06-24 17:16:47
感謝參與，應助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-06-24 17:16:47
|
與初值選取不當有個吧，試試1stOpt，不用初值，不過本題似乎有多組解： 1： x1: 0.452159092308927 x2: -3.5576820815459E-17 x3: -1.6802621629643E-17 x4: 0.70583278260955 x5: 0.289090165134331 x6: 0.579315494431816 2： x1: 2.3247555351878 x2: -3.51606195970296 x3: 2.54582708579748E-8 x4: -1.05916256056018 x5: 1.08911300900033 x6: -0.182685684997559 |
3樓2013-06-24 09:06:15
中將
木蟲 (小有名氣)
- 應助: 2 (幼兒園)
- 金幣: 1958.7
- 散金: 86
- 紅花: 1
- 帖子: 221
- 在線: 76小時
- 蟲號: 886327
- 注冊: 2009-10-28
- 性別: GG
- 專業(yè): 分離過程
4樓2013-06-27 12:24:25
相關版塊跳轉
我要訂閱樓主 中將 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學求調劑，各位老師看看我�。。� +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 化學調劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 考研調劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 308求調劑 +6 | 墨墨漠 2026-03-27 | 6/300 |
|
|
|
[考研] 299求調劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
|
[考研] 085601 材料工程 313分 求調劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿鄭大085600，310分求調劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 292求調劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 316求調劑 +5 | 江辭666 2026-03-26 | 5/250 |
|
|
|
[考研] 求調劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
|
[考研]
 一志愿廈門大學化學學碩307求調劑
+8
|
y7czhao 2026-03-26 | 8/400 |
|
|
|
[考研] 340求調劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 資源與環(huán)境 調劑申請(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
|
[考研] 332求調劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 318求調劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
|
[考研] 0703化學調劑，求導師收 +7 | 天天好運來上岸?/a> 2026-03-24 | 7/350 |
|
|
|
[考研] 080500求調劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 285求調劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
|
[考研] 求調劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[考研] 330求調劑0854 +3 | assdll 2026-03-21 | 3/150 |
|