 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 743 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
中將木蟲 (小有名氣)
|
[求助]
求高手幫我分析下我這個(gè)matlab求解非線性微分方程組解的情況
|
||
|
我這個(gè)是一組非線性方程組，可是每次給不同的初值都達(dá)不到我要求結(jié)果，不知道怎么辦，求大神指教，不知道是方法不對(duì)還是怎么的，我希望結(jié)果輸出的都是正數(shù)在0和1之間，而且最后的total要求最好小于1。。。。。。。。這樣的非線性方程組怎么求好呢，有什么好的方法呢，求教求教 這個(gè)是我函數(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ù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有300人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收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)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab如何求解一個(gè)非線性微分方程組
已經(jīng)有8人回復(fù)
- 求助，用matlab解微分方程組，希望高手能給予指點(diǎn)
已經(jīng)有17人回復(fù)
- matlab解偏微分方程求助
已經(jīng)有12人回復(fù)
- 求高手解答Matlab求解微分方程的初始值和邊界值問(wèn)題
已經(jīng)有5人回復(fù)
- 二階非線性微分方程求解
已經(jīng)有14人回復(fù)
- 請(qǐng)問(wèn)這個(gè)4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。
已經(jīng)有17人回復(fù)
- 求Matlab高手解決線性方程組的迭代求解問(wèn)題
已經(jīng)有17人回復(fù)
- 請(qǐng)教用matlab求解一個(gè)非線性偏微分方程組的數(shù)值解
已經(jīng)有8人回復(fù)
- 高手進(jìn)來(lái)用matlab求解常微分方程，十分感謝！
已經(jīng)有17人回復(fù)
- matlab有限差分法求解復(fù)雜的偏微分方程組
已經(jīng)有19人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(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ù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題
已經(jīng)有21人回復(fù)
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 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)有個(gè)吧，試試1stOpt，不用初值，不過(guò)本題似乎有多組解： 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小時(shí)
- 蟲號(hào): 886327
- 注冊(cè): 2009-10-28
- 性別: GG
- 專業(yè): 分離過(guò)程
2樓2013-06-22 16:52:34
中將
木蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1958.7
- 散金: 86
- 紅花: 1
- 帖子: 221
- 在線: 76小時(shí)
- 蟲號(hào): 886327
- 注冊(cè): 2009-10-28
- 性別: GG
- 專業(yè): 分離過(guò)程
4樓2013-06-27 12:24:25
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程求調(diào)劑 +8 | 多久上課 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +3 | 手機(jī)用戶 2026-03-26 | 3/150 |
|
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +8 | 路癡小琪 2026-03-27 | 8/400 |
|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
|
[考研] 22408 359分調(diào)劑 +3 | Qshers 2026-03-27 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +4 | 星野? 2026-03-26 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +7 | 安憶靈 2026-03-22 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +10 | 呼呼？~+123456 2026-03-24 | 10/500 |
|
|
|
[考研] 316求調(diào)劑 +5 | 江辭666 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
|
[考研] 尋找調(diào)劑 +5 | 倔強(qiáng)芒? 2026-03-21 | 8/400 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研]
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|