| 5 | 1/1 | 返回列表 |
| 查看: 728 | 回復: 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 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工304求調(diào)劑 +3 | 邱gl 2026-03-05 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑材料學碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
[基金申請] 沒有青基直接申請面上,感覺自己瘋了 +6 | kevin63t 2026-03-02 | 12/600 |
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 266求調(diào)劑 +7 | 哇塞王帥 2026-03-03 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個椒 2026-03-04 | 3/150 |
|
|
[考研] 能動297求調(diào)劑,本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
[考研] 理學,工學,農(nóng)學調(diào)劑,少走彎路,這里歡迎您! +8 | likeihood 2026-03-02 | 11/550 |
|
|
[考研] 085700資環(huán)求調(diào)劑,初始279,六級已過,英語能力強 +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
[考研] 材料學碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
[考博] 26申博 +4 | 想申博! 2026-02-26 | 6/300 |
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|