| 4 | 1/1 | 返回列表 |
| 查看: 726 | 回復: 3 | ||
alittlehug新蟲 (初入文壇)
|
[求助]
1stopt求解變參數(shù)非線性方程組結果只有一組解 已有2人參與
|
|
用1stopt5.0運行的時候,非線性方程只有一組解,實際上應該有多組解 Constant A=2,gamma1=0.1,gamma2=0.2,k=500,w0=15,eps=0.1; LoopConstant omega=[10:1:20] ; Parameters x1,x2,x3,x4; Function -x2*omega+gamma1*x1+w0^2*x2/omega+3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3=A; x1*omega+x2*gamma1-w0^2*x1/omega-3*k*x3*(x3^2+x4^2)/(4*omega^3)+gamma2*x4=0; x1*gamma1-x4*omega+w0^2*x2/omega+(1+1/eps)*(3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3)=A; x2*gamma1+x3*omega-w0^2*x1/omega-(1+1/eps)*(3*k*x3*(x3^2+x4^2)/(4*omega^3)-gamma2*x4)=0; 結果顯示: Constant A=2,gamma1=0.1,gamma2=0.2,k=500,w0=15,eps=0.1; LoopConstant omega=[10:1:20] ; Parameters x1,x2,x3,x4; Function -x2*omega+gamma1*x1+w0^2*x2/omega+3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3=A; x1*omega+x2*gamma1-w0^2*x1/omega-3*k*x3*(x3^2+x4^2)/(4*omega^3)+gamma2*x4=0; x1*gamma1-x4*omega+w0^2*x2/omega+(1+1/eps)*(3*k*x4*(x3^2+x4^2)/(4*omega^3)+gamma2*x3)=A; x2*gamma1+x3*omega-w0^2*x1/omega-(1+1/eps)*(3*k*x3*(x3^2+x4^2)/(4*omega^3)-gamma2*x4)=0; 這里我只放了omega=20的時候方程組的解。 |
鐵桿木蟲 (職業(yè)作家)
|
版本問題吧,有11組解: 循環(huán)常數(shù) omega 目標函數(shù)值 x1 x2 x3 x4 10 4.58635664500889E-26 0.00380294459721743 0.160294854647957 -0.0277043774026065 0.156228877092279 11 7.91721570819433E-27 0.00669308907034421 0.211859317209918 -0.0315743501340042 0.208493086819324 12 1.34137963790563E-24 0.0132338824083539 0.296344281124706 -0.0365113493760948 0.294971542948296 13 1.01936620072096E-23 0.0328237831973817 0.461292243502208 -0.040457396825942 0.468686530302756 14 9.15535540666567E-23 0.138909893810943 0.895184889404847 -0.00267827329179549 0.989325928727004 15 1.27077263082487E-20 0.746988968174331 1.13840265715959 0.822882483116559 2.42193671971171 16 1.52179093252339E-21 1.20677974111929 0.376059697385337 2.69161830257015 2.00602076273742 17 3.86982718524109E-24 0.042338001880337 -0.52754140284278 0.105065051557478 -0.52181906319924 18 1.78381813142587E-25 0.0196379178791081 -0.36165925952281 0.0595504858703478 -0.356709103497778 19 6.20181994642092E-26 0.011584651009512 -0.278277662864219 0.0405814648006044 -0.274614864703752 20 1.31099134650033E-27 0.00775740966522293 -0.227852327097322 0.0303034605957675 -0.225094128630944 |
銅蟲 (小有名氣)

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +6 | xiiiia 2026-03-04 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
[考研] 材料328求調(diào)劑 +7 | 一個蘿卜02 2026-03-03 | 7/350 |
|
|
[考研] 歡迎采礦、地質、巖土、計算機、人工智能等專業(yè)的同學報考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
[考研] 320材料與化工,求調(diào)劑 +4 | 鶴遨予卿 2026-03-04 | 5/250 |
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
[考研] 0703化學調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
[考研] 264求調(diào)劑 +3 | thext 2026-03-03 | 3/150 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
[考研] 中國林科院林化所(南京)2026年招收化學/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
[基金申請] 請問大家,研究風險與應對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
[考研] 26考研報考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 11408,學碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
[考研] 085600材料工程一志愿中科大總分312求調(diào)劑 +9 | 吃宵夜1 2026-02-28 | 11/550 |
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|