| 查看: 2360 | 回復: 10 | ||
fugitivesk鐵蟲 (初入文壇)
|
[求助]
求助!三個未知數(shù) 四個方程用1stopt(Anto2Fit)約束優(yōu)化 已有2人參與
|
計算小白,新來的沒啥金幣![]() 用的Anto2Fit 5.5試用版(明明三個parameters為啥說超出四個了) 方程都比較長 不知道為什么三個未知數(shù)會出現(xiàn)四個結果 第四個結果好像是我給的常數(shù)變量或者變參量(偷梁換柱?) 貼上mathematica里寫的表達式(這個軟件里運行過三個方程三個變量的求解,方程本身沒問題) Constant cc = 3*10^8; la = cc/10^9; k0 = 2*pi/la; r1 =0.17532; r2 =0.21216; t1 =0.77139; t2 =-0.3507; Parameter k1[30,60],k2[-2,0],dz[0.5,0.7]; Function a=sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz)); b=cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz)); c=cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)); d=sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)); kk=2*(k1^2+k2^2); (r1*c)-(r2*d)+r1*(b*k1-a*k2)/(2*k0)+r1*(a*k0*k2+b*k0*k1)/kk-r2*(a*k1+b*k2)/(2*k0)-r2*(a*k0*k1-b*k0*k2)/kk-(b*k1-a*k2)/(2*k0)+(a*k0*k2+b*k0*k1)/kk; (r2*c)+(r1*d)+r1*(a*k1+b*k2)/(2*k0)+r1*(a*k0*k1-b*k0*k2)/kk+r2*(b*k1-a*k2)/(2*k0)+r2*(a*k0*k2+b*k0*k1)/kk-(a*k1+b*k2)/(2*k0)-(b*k0*k2-a*k0*k1)/kk; (t1*c)-(t2*d)+t1*(b*k1-a*k2)/(2*k0)+t1*(a*k0*k2+b*k0*k1)/kk-t2*(a*k1+b*k2)/(2*k0)-t2*(a*k0*k1-b*k0*k2)/kk+2; (t2*c)+(t1*d)+t1*(a*k1+b*k2)/(2*k0)+t1*(a*k0*k1-b*k0*k2)/kk+t2*(b*k1-a*k2)/(2*k0)+t2*(a*k0*k2+b*k0*k1)/kk; //方程我用了兩種寫法第二種太長了 ![]() Function (r1*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))-(r2*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+r1*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+r1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2)-r2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)-r2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)-(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2); (r2*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))+(r1*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+r1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)+r1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)+r2*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+r2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2)-(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)-(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1)/2*(k1^2+k2^2); (t1*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))-(t2*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+t1*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+t1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2)-t2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)-t2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)+2; (t2*cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)))+(t1*sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)))+t1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k1+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k2)/(2*k0)+t1*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k1-cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k2)/2*(k1^2+k2^2)+t2*(cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k1-sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k2)/(2*k0)+t2*(sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz))*k0*k2+cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz))*k0*k1)/2*(k1^2+k2^2); 22.JPG@月只藍 |
新蟲 (小有名氣)

鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
新蟲 (小有名氣)
|
1. 主要使用自己寫的MathSword軟件進行求解。(PS:附件可以下載MathSword,在打開程序界面最上端選擇【高級優(yōu)化】,在優(yōu)化模塊右鍵選擇上傳的附件 solve.vb文件,然后點擊求解,就能看到你問題的答案了。) 2. 至于你說的方程用什么約束比較合適,不太明白你的意思。因為一個方程本身就是一個約束,有約束就可以構建誤差,一般的優(yōu)化器內部會根據(jù)你構建的1個或多個誤差進行優(yōu)化。 |

版主 (知名作家)
|
代碼寫的有問題。 k1: 48.5939435307918 k2: -0.161774549680272 dz: 0.588775301298235 ------------------------------------------------ Constant cc = 3*10^8, r1 =0.17532, r2 =0.21216, t1 =0.77139, t2 =-0.3507, la = cc/10^9, k0 = 2*pi/la; Parameter k1=[30,60],k2=[-2,0],dz=[0.5,0.7]; ConstStr kk=2*(k1^2+k2^2), a=sin(k1*dz)*(exp(-k2*dz)+exp(k2*dz)), b=cos(k1*dz)*(exp(-k2*dz)-exp(k2*dz)), c=cos(k1*dz)*(exp(-k2*dz)+exp(k2*dz)), d=sin(k1*dz)*(exp(-k2*dz)-exp(k2*dz)); Function (r1*c)-(r2*d)+r1*(b*k1-a*k2)/(2*k0)+r1*(a*k0*k2+b*k0*k1)/kk-r2*(a*k1+b*k2)/(2*k0)-r2*(a*k0*k1-b*k0*k2)/kk-(b*k1-a*k2)/(2*k0)+(a*k0*k2+b*k0*k1)/kk=0; (r2*c)+(r1*d)+r1*(a*k1+b*k2)/(2*k0)+r1*(a*k0*k1-b*k0*k2)/kk+r2*(b*k1-a*k2)/(2*k0)+r2*(a*k0*k2+b*k0*k1)/kk-(a*k1+b*k2)/(2*k0)-(b*k0*k2-a*k0*k1)/kk=0; (t1*c)-(t2*d)+t1*(b*k1-a*k2)/(2*k0)+t1*(a*k0*k2+b*k0*k1)/kk-t2*(a*k1+b*k2)/(2*k0)-t2*(a*k0*k1-b*k0*k2)/kk+2=0; (t2*c)+(t1*d)+t1*(a*k1+b*k2)/(2*k0)+t1*(a*k0*k1-b*k0*k2)/kk+t2*(b*k1-a*k2)/(2*k0)+t2*(a*k0*k2+b*k0*k1)/kk=0; |

鐵蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
鐵蟲 (初入文壇)
鐵蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料復試調劑 +7 | 學材料的點 2026-03-01 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿985材料與化工 326分求調劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 一志愿西交化工專碩288專業(yè)課93求調劑 +5 | 好運好運接接 2026-03-04 | 5/250 |
|
|
[考研] 本科太原理工采礦工程,求調劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
[考研] 求調劑 +3 | 泡了個椒 2026-03-04 | 3/150 |
|
|
[考研]
085600 英一數(shù)二272求調劑
5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
[考研] 085600求調劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
[考研]
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|
|
[考研] 298求調劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè),申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 291求調劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 288求調劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 0856材料調劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 材料調劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
[考研] 275求調劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
[基金申請] 成果系統(tǒng)訪問量大,請一小時后再嘗試。---NSFC啥時候好哦,已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
[考研] 307求調劑 +4 | 73372112 2026-02-28 | 6/300 |
|