 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2453 | 回復(fù): 4 | ||
[求助]
求助 （已經(jīng)-無量綱化）的非線性方程組（已經(jīng)仔細(xì)檢查，能算，但exitflag=0求高手
|
|
其實是后邊要用matcont做常微分方程組延拓分析的。但是那個延拓軟件必須有一個穩(wěn)態(tài)解（充分發(fā)展的，dx/dt=0時的解）作為初值，才能開始延拓。 代碼如下：已經(jīng)仔細(xì)檢查了一周了。能算，但是exitflag=0，估計初值太敏感。菜鳥求助啊。。。（%后邊是多余的項，不用管，直接貼代碼進(jìn)編輯器就行。。） 小女再次感謝！ steady.m X0=[0.278212173322273,0.912568632338739,0.085310950018308,0.912568731077048,0.087431259625095,0.087431260178122,0.499999721583934,0.000000000015842,0.087431268922953,0.053882265858466,0.087432024167792,0.087431276521867,0.500000000000000,0.087431268922956,0.087431268922956,0.500000000000000,0.087431273571452,0.500000000000001,0.078613189530066]; options=optimset('Display','iter','MaxFunEvals',100000,'MaxIter',800) %'MaxFunEvals',100000000); [x,fval,exitflag] =fsolve('steadys',X0,options) steadys.m function pp=steadys(x) pp(1)=(1-x(1))+0.00001*48.432471*(0.0132913*x(3)*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*1.6037e+08*2.71828^(-18.8930/x(1))*x(6)+5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(10))-2.782248*(x(1)-x(19)); pp(2)=1-x(2)-8.522*10^(-3)*1.3028e+12*2.71828^(-27.8955/x(1))*x(2)^2-5.6228*10^(-3)*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(3)-0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(7)-0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(10)-2.493604*10^(-4)*6.6735e+19*2.71828^(-45.6473/x(1))*x(9); %-y7*c4*x(2)*x(10) pp(3)=-x(3)+1.1722874*1.3028e+12*2.71828^(-27.8955/x(1))*x(2)^2-1.546946*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(3)-3.198*10^(-4)*0.01*x(3)*x(8); %-y10*c5*x(3)*x(10); pp(4)=1-x(4)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(4)-0.01923587*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(6); pp(5)=-x(5)+2.445938*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(3)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)-8.496798*17.6895*2.71828^(-2.8730/x(1))*x(5)*x(8)+3.2672*10^10*1.1663e+16*2.71828^(-36.9952/x(1))*x(9)-(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(5)*x(10); %+y16*c4*x(2)*x(10);+y18*c11 pp(6)=-x(6)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)-2.8425*10^10*1.1663e+16*2.71828^(-36.9952/x(1))*x(4)-8.496798*10^6*17.6895*2.71828^(-2.8730/x(1))*x(6)*x(8); pp(7)=-x(7)*2.44595*10^8*1.6037e+08*2.71828^(-18.8930/x(1))*x(2)*x(4)-5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(7)-(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(7)*x(10); %-y23*c5*x(3)*x(10);+y18*c11 pp(8)=-x(8)-5.871963*10^13*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(10)+5.05856*1.1663e+16*2.71828^(-36.9952/x(1))*x(13)+4.515847*10^5*1.1663e+16*2.71828^(-36.9952/x(1))*x(4)-134.9877*17.6895*2.71828^(-2.8730/x(1))*x(5)*x(8)+5.871963*10^13*1.1663e+16*2.71828^(-36.9952/x(1))*x(9)-8.033211*0.01*x(3)*x(8); pp(9)=-x(9)+1.699359668*10^(-2)*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(5)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(9)-0.216944*6.6735e+19*2.71828^(-45.6473/x(1))*x(9); pp(10)=-x(10)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)*65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(13)-x(10)*(8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)+4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1))*x(10)); %+y7*c4*x(2);+y18*c11 pp(11)=-x(11)+0.001*(0.0132913*x(3)*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)+5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(10))+65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(14)-x(11)*(8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)+(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)); %+y7*c4*x(2);+y32*c5*x(3);+y18*c11 pp(12)=-x(12)+0.00001*(0.0132913*x(3)*x(2)*x(7)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(5)+0.0132913*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*x(6)+5.781275*1.6306e+04*2.71828^(-9.6993/x(1))*x(2)*(200*x(11)+x(10))+65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(15)-x(12)*(150.3041*2.71828^(-2.8730/x(1))*x(8)+0.036761679*2.71828^(-2.0910/x(1)))*x(10)); %+y7*c4*x(2);+y32*c5*x(3);+y18*c11 pp(13)=-x(13)*8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(10)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(13); pp(14)=-x(14)+8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(11)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(14); pp(15)=-x(15)+8.496798*17.6895*2.71828^(-2.8730/x(1))*x(8)*x(12)-65.3447*1.1663e+16*2.71828^(-36.9952/x(1))*x(14); pp(16)=-x(16)*(0.5*4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)^2; %+y18*c11;+x(10)*(y7*c4*x(2)+y32*c5*x(3)); pp(17)=-x(17)+(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)*x(11); %+y18*c11;+x(11)*(y7*c4*x(2)+y32*c5*x(3)); pp(18)=-x(18)*(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)*x(12)+4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1))*x(11); %+y18*c11;+x(12)*(y7*c4*x(2)+y32*c5*x(3)) pp(19)=(0.727148704-12.7016)*x(19)+4.897673*(x(1)-x(19)); |
回復(fù)此樓» 猜你喜歡
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有101人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有26人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有34人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助數(shù)學(xué)大神 如何做圖像的相似性分析
已經(jīng)有5人回復(fù)
- 遼寧大學(xué)考博求助，遼寧大學(xué)考博英語近年是統(tǒng)一的還是針對學(xué)院的
已經(jīng)有6人回復(fù)
- pets5經(jīng)驗求助
已經(jīng)有4人回復(fù)
- 求助混凝土用環(huán)氧膩子，解決膩子滲透性問題，謝謝
已經(jīng)有7人回復(fù)
- 求助關(guān)于SPSS數(shù)據(jù)處理中的問題。。。。。高人呢。。。
已經(jīng)有8人回復(fù)
- 求助，納米復(fù)合磁性材料變形過程中的亞穩(wěn)相！
已經(jīng)有7人回復(fù)
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較
已經(jīng)有4人回復(fù)
- 請教Matlab高手解多元非線性方程組中的參數(shù)
已經(jīng)有6人回復(fù)
- 求助：origin8.0 自定義非線性擬合，參數(shù)擬合值與設(shè)定初始值一樣
已經(jīng)有18人回復(fù)
- 【求助】用mathematica 5.0求解一個非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 非線性方程組的迭代法（數(shù)值計算高手請進(jìn)）
已經(jīng)有7人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★
感謝參與，應(yīng)助指數(shù) +1
dbb627: 金幣+1, 歡迎交流 2012-05-11 21:58:10
感謝參與，應(yīng)助指數(shù) +1
dbb627: 金幣+1, 歡迎交流 2012-05-11 21:58:10
|
用1stOpt跑的結(jié)果： x1: 0.336178337939607 x2: 1.00000000001363 x3: 1.45731439404896E-15 x4: 1 x5: -1.80058046716552E-16 x6: -7.0862384417084E-17 x7: -0.213180808520247 x8: 0.00029442958151592 x9: -1.51795582475496E-15 x10: -4.95455914027804E-12 x11: -1.94482439925578E-14 x12: 3.65096737849878E-14 x13: 0.23526045284234 x14: 4.68097286349597E-16 x15: 1.79081692185712E-14 x16: 2.7607895322144 x17: -8.50846322024266E-17 x18: -0.0647548483684659 x19: 0.0975865006697487 不知道你的公式有誤問題，比如其中： pp(16)=-x(16)*(0.5*4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x(1)))*x(10)^2; 似乎不大對吧？ |
2樓2012-05-11 15:34:20
3樓2012-05-12 14:07:43
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應(yīng)助回帖
★ ★
cenwanglai: 金幣+2, 謝謝回帖~ 2012-05-12 17:26:32
cenwanglai: 金幣+2, 謝謝回帖~ 2012-05-12 17:26:32
|
如果都要求大于0（開始就應(yīng)該把這個條件說明）： x1: 0.336178337925227 x2: 0.999999999862475 x3: 2.27071941197048E-11 x4: 0.999999999999997 x5: 4.50305709438193E-12 x6: 1.22744048689865E-12 x7: 2.15050770012278 x8: 5.18658727927569E-14 x9: 4.62640607748565E-11 x10: 2.8267681811832E-11 x11: 2.22086150893169E-11 x12: 5.49036141029378E-12 x13: 0.695061141521546 x14: 2.9846697983009E-16 x15: 1.34476078944492E-14 x16: 0.00186288613726707 x17: 4.46468502401335E-11 x18: 0.00164117436817047 x19: 0.0975865006648326 自己帶回方程驗證，如果方程值均趨于0，可以說明至少在理論上是其中一組解。 |
4樓2012-05-12 14:17:44
|
您好！非常感謝！想請問您這個理論上能用1stopt1.5算出來嗎？為什么我在1stopt里報錯說“參數(shù)雖然被定義，但沒有在表達(dá)式中出現(xiàn)呢？” Parameter x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19; Function(1-x1)+0.00001*48.432471*(0.0132913*x3*x2*x7+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x5+0.0132913*16306*2.71828^(-9.6993/x1)*1.6037e+08*2.71828^(-18.8930/x1)*x6+5.781275*16306*2.71828^(-9.6993/x1)*x2*x10)-2.782248*(x1-x(19))=0 1-x2-8.522*10^(-3)*1.3028*10^12*2.71828^(-27.8955/x1)*x2^2-5.6228*10^(-3)*1.6037e+08*2.71828^(-18.8930/x1)*x2*x3-0.0132913*16306*2.71828^(-9.6993/x1)*x2*x7-0.0132913*16306*2.71828^(-9.6993/x1)*x2*x6-5.781275*16306*2.71828^(-9.6993/x1)*x2*x10-2.493604*10^(-4)*6.6735e+19*2.71828^(-45.6473/x1)*x9=0 x3+1.1722874*1.3028e+12*2.71828^(-27.8955/x1)*x2^2-1.546946*1.6037e+08*2.71828^(-18.8930/x1)*x2*x3-3.198*10^(-4)*0.01*x3*x8=0 1-x4-65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x4-0.01923587*17.6895*2.71828^(-2.8730/x1)*x8*x6=0 -x5+2.445938*1.6037e+08*2.71828^(-18.8930/x1)*x2*x3-5.781275*16306*2.71828^(-9.6993/x1)*x2*x5-8.496798*17.6895*2.71828^(-2.8730/x1)*x5*x8+3.2672*10^10*1.1663e+16*2.71828^(-36.9952/x1)*x9-(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1))*x5*x10=0 -x6-5.781275*16306*2.71828^(-9.6993/x1)*x2*x6-2.8425*10^10*1.1663e+16*2.71828^(-36.9952/x1)*x4-8.496798*10^6*17.6895*2.71828^(-2.8730/x1)*x6*x8=0 -x7*2.44595*10^8*1.6037e+08*2.71828^(-18.8930/x1)*x2*x4-5.781275*16306*2.71828^(-9.6993/x1)*x2*x7-(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1))*x7*x10=0 -x8-5.871963*10^13*17.6895*2.71828^(-2.8730/x1)*x8*x10+5.05856*1.1663e+16*2.71828^(-36.9952/x1)*x13+4.515847*10^5*1.1663e+16*2.71828^(-36.9952/x1)*x4-134.9877*17.6895*2.71828^(-2.8730/x1)*x5*x8+5.871963*10^13*1.1663e+16*2.71828^(-36.9952/x1)*x9-8.033211*0.01*x3*x8=0 -x9+1.699359668*10^(-2)*17.6895*2.71828^(-2.8730/x1)*x8*x5-65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x9-0.216944*6.6735e+19*2.71828^(-45.6473/x1)*x9=0 -x10+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x7+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x5+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x6*65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x13-x10*(8.496798*17.6895*2.71828^(-2.8730/x1)*x8+4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1)*x10)=0 -x11+0.001*(0.0132913*x3*x2*x7+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x5+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x6+5.781275*16306*2.71828^(-9.6993/x1)*x2*x10)+65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x14-x11*(8.496798*17.6895*2.71828^(-2.8730/x1)*x8+(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1))*x10)=0 -x12+0.00001*(0.0132913*x3*x2*x7+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x5+0.0132913*16306*2.71828^(-9.6993/x1)*x2*x6+5.781275*16306*2.71828^(-9.6993/x1)*x2*(200*x11+x10)+65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x15-x12*(150.3041*2.71828^(-2.8730/x1)*x8+0.036761679*2.71828^(-2.0910/x1))*x10)=0 -x13*8.496798*17.6895*2.71828^(-2.8730/x1)*x8*x10-65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x13=0 -x14+8.496798*17.6895*2.71828^(-2.8730/x1)*x8*x11-65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x14=0 -x15+8.496798*17.6895*2.71828^(-2.8730/x1)*x8*x12-65.3447*1.1663e+16*2.71828^(-36.9952/x1)*x14=0 -x16*(0.5*4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1))*x10^2=0 -x17+(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1))*x10*x11=0 -x18*(4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1))*x10*x12+4.54117*10^(-3)*8.0952*2.71828^(-2.0910/x1)*x11=0 (0.727148704-12.7016)*x(19)+4.897673*(x1-x19)=0 |
5樓2012-05-18 17:06:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 nanpse 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[基金申請]
有必要更換申報口嗎
 20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|