 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 742 | 回復(fù): 3 | ||
中將木蟲 (小有名氣)
|
[求助]
求高手幫我分析下我這個matlab求解非線性微分方程組解的情況
|
|
我這個是一組非線性方程組，可是每次給不同的初值都達不到我要求結(jié)果，不知道怎么辦，求大神指教，不知道是方法不對還是怎么的，我希望結(jié)果輸出的都是正數(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)); ]; 下面是調(diào)用 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 |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有92人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab如何求解一個非線性微分方程組
已經(jīng)有8人回復(fù)
- 求助，用matlab解微分方程組，希望高手能給予指點
已經(jīng)有17人回復(fù)
- matlab解偏微分方程求助
已經(jīng)有12人回復(fù)
- 求高手解答Matlab求解微分方程的初始值和邊界值問題
已經(jīng)有5人回復(fù)
- 二階非線性微分方程求解
已經(jīng)有14人回復(fù)
- 請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。
已經(jīng)有17人回復(fù)
- 求Matlab高手解決線性方程組的迭代求解問題
已經(jīng)有17人回復(fù)
- 請教用matlab求解一個非線性偏微分方程組的數(shù)值解
已經(jīng)有8人回復(fù)
- 高手進來用matlab求解常微分方程，十分感謝！
已經(jīng)有17人回復(fù)
- matlab有限差分法求解復(fù)雜的偏微分方程組
已經(jīng)有19人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復(fù)
- Matlab求解偏微分方程組
已經(jīng)有22人回復(fù)
- matlab求解非線性方程組
已經(jīng)有16人回復(fù)
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復(fù)
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- matlab解微分方程組
已經(jīng)有15人回復(fù)
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復(fù)
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復(fù)
中將
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1958.7
- 散金: 86
- 紅花: 1
- 帖子: 221
- 在線: 76小時
- 蟲號: 886327
- 注冊: 2009-10-28
- 性別: GG
- 專業(yè): 分離過程
2樓2013-06-22 16:52:34
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-06-24 17:16:47
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-06-24 17:16:47
|
與初值選取不當(dāng)有個吧，試試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
中將
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1958.7
- 散金: 86
- 紅花: 1
- 帖子: 221
- 在線: 76小時
- 蟲號: 886327
- 注冊: 2009-10-28
- 性別: GG
- 專業(yè): 分離過程
4樓2013-06-27 12:24:25
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 中將 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856，材料與化工321分求調(diào)劑 +4 | 大饞小子 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +7 | lin0039 2026-03-26 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +15 | 顧九笙要謙虛 2026-03-24 | 21/1050 |
|
|
|
[考研] 269專碩求調(diào)劑 +9 | 金恩貝 2026-03-21 | 9/450 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 321求調(diào)劑 +5 | 材料cailiao 2026-03-21 | 5/250 |
|
|
|
[考研] 08開頭275求調(diào)劑 +3 | 拉誰不重要 2026-03-26 | 3/150 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 機械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 11/550 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
|
[考研] 生物學(xué)一志愿985，分數(shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分數(shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|