版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

fugitivesk

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 68.2
帖子: 42
在線: 36小時(shí)
蟲(chóng)號(hào): 22326762
注冊(cè): 2020-05-25
性別: GG
專(zhuān)業(yè): 電磁場(chǎng)與波

[求助] 求助！三個(gè)未知數(shù) 四個(gè)方程用1stopt(Anto2Fit)約束優(yōu)化已有2人參與

計(jì)算小白，新來(lái)的沒(méi)啥金幣

用的Anto2Fit 5.5試用版（明明三個(gè)parameters為啥說(shuō)超出四個(gè)了）
方程都比較長(zhǎng)
不知道為什么三個(gè)未知數(shù)會(huì)出現(xiàn)四個(gè)結(jié)果
第四個(gè)結(jié)果好像是我給的常數(shù)變量或者變參量（偷梁換柱？）
貼上mathematica里寫(xiě)的表達(dá)式（這個(gè)軟件里運(yùn)行過(guò)三個(gè)方程三個(gè)變量的求解，方程本身沒(méi)問(wèn)題）

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;

//方程我用了兩種寫(xiě)法第二種太長(zhǎng)了

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);

求助！三個(gè)未知數(shù) 四個(gè)方程用1stopt(Anto2Fit)約束優(yōu)化

22.JPG@月只藍(lán)

回復(fù)此樓

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生已經(jīng)有4人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有117人回復(fù)
計(jì)算三階力常數(shù)用thirdorder擴(kuò)胞時(shí)參數(shù)如何選擇？臨近原子數(shù)怎么取值已經(jīng)有7人回復(fù)
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有2人回復(fù)

1樓 2020-06-24 18:48:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

shikang999

新蟲(chóng) (小有名氣)

計(jì)算強(qiáng)帖: 1
應(yīng)助: 27 (小學(xué)生)
金幣: 1454.7
紅花: 7
帖子: 182
在線: 206.4小時(shí)
蟲(chóng)號(hào): 2201981
注冊(cè): 2012-12-23
性別: GG
專(zhuān)業(yè): 固體力學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
獨(dú)孤神宇: 金幣+5, 鼓勵(lì)交流 2020-06-28 20:41:18
fugitivesk: 金幣+7, ★★★很有幫助 2020-06-28 21:14:32

引用回帖:

3樓: Originally posted by fugitivesk at 2020-06-28 19:11:15
請(qǐng)問(wèn)求解思路是什么，用matlab還是mathemetica求解的嗎？...

1. 主要使用自己寫(xiě)的MathSword軟件進(jìn)行求解。（PS：附件可以下載MathSword，在打開(kāi)程序界面最上端選擇【高級(jí)優(yōu)化】，在優(yōu)化模塊右鍵選擇上傳的附件 solve.vb文件，然后點(diǎn)擊求解，就能看到你問(wèn)題的答案了。）

2. 至于你說(shuō)的方程用什么約束比較合適，不太明白你的意思。因?yàn)橐粋€(gè)方程本身就是一個(gè)約束，有約束就可以構(gòu)建誤差，一般的優(yōu)化器內(nèi)部會(huì)根據(jù)你構(gòu)建的1個(gè)或多個(gè)誤差進(jìn)行優(yōu)化。

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : MathSword.2020.06.27.rar

2020-06-28 20:08:56, 2.97 M

附件 2 : solve.vb

2020-06-28 20:09:09, 6.46 K

我還是那個(gè)我，過(guò)去如此，未來(lái)亦如此！

5樓2020-06-28 20:26:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 11 個(gè)回答

shikang999

新蟲(chóng) (小有名氣)

計(jì)算強(qiáng)帖: 1
應(yīng)助: 27 (小學(xué)生)
金幣: 1454.7
紅花: 7
帖子: 182
在線: 206.4小時(shí)
蟲(chóng)號(hào): 2201981
注冊(cè): 2012-12-23
性別: GG
專(zhuān)業(yè): 固體力學(xué)

【答案】應(yīng)助回帖

★
fugitivesk: 金幣+1, ★有幫助 2020-06-28 19:26:45

按你第一種寫(xiě)法, 如果最后四個(gè)方程是等于0的約束的話，求得一個(gè)如下的解(其中四個(gè)方程的誤差平方和為2.07045657911154E-05）（注意，下面不是1stopt的結(jié)果，所以不保證和1stopt的結(jié)果一致）
k1 = 48.5939462992421, k2 = -0.161774554371271, dz = 0.588775267394145

贊一下

回復(fù)此樓

我還是那個(gè)我，過(guò)去如此，未來(lái)亦如此！

2樓2020-06-28 15:36:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fugitivesk

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 68.2
帖子: 42
在線: 36小時(shí)
蟲(chóng)號(hào): 22326762
注冊(cè): 2020-05-25
性別: GG
專(zhuān)業(yè): 電磁場(chǎng)與波

引用回帖:

2樓: Originally posted by shikang999 at 2020-06-28 15:36:14
按你第一種寫(xiě)法, 如果最后四個(gè)方程是等于0的約束的話，求得一個(gè)如下的解(其中四個(gè)方程的誤差平方和為2.07045657911154E-05）（注意，下面不是1stopt的結(jié)果，所以不保證和1stopt的結(jié)果一致）
k1 = 48.5939462992421 ...

請(qǐng)問(wèn)求解思路是什么，用matlab還是mathemetica求解的嗎？

贊一下

回復(fù)此樓

3樓2020-06-28 19:11:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fugitivesk

鐵蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 68.2
帖子: 42
在線: 36小時(shí)
蟲(chóng)號(hào): 22326762
注冊(cè): 2020-05-25
性別: GG
專(zhuān)業(yè): 電磁場(chǎng)與波

引用回帖:

我原先就只是四選三作為方程組求解，感覺(jué)不能保證誤差。四個(gè)方程的約束用什么方法比較合適呢?

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

4樓2020-06-28 19:24:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 11 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 282分材料專(zhuān)業(yè)求調(diào)劑院校 +10	楓橋ZL 2026-03-09	13/650	2026-03-10 00:47 by wdwdoo
[考研] 294 英二數(shù)二物化求調(diào)劑 +6	米飯團(tuán)不好吃 2026-03-09	6/300	2026-03-09 23:55 by barlinike
[考研] 295復(fù)試調(diào)劑 +4	簡(jiǎn)木ChuFront 2026-03-09	4/200	2026-03-09 21:49 by wangjihu
[考研] 調(diào)劑 +4	調(diào)劑的考研學(xué)生 2026-03-09	4/200	2026-03-09 20:22 by wsxw
[考研] 0832食品科學(xué)與工程293調(diào)劑 +3	東東不東 2026-03-07	3/150	2026-03-09 15:58 by ls_13456
[考研] 求調(diào)劑，數(shù)一英一274分 +4	小菲會(huì)努力 2026-03-08	4/200	2026-03-09 12:40 by 一定上岸喲_
[考研] 0856求調(diào)劑 +3	squirtle11 2026-03-07	3/150	2026-03-09 09:54 by @颯颯颯颯
[考研] 297求調(diào)劑 +3	胡達(dá)靈 2026-03-05	5/250	2026-03-08 23:13 by 清風(fēng)月
[考研] 求調(diào)劑，一志愿江南大學(xué)，食品科學(xué)與工程，總分，320 +3	yyyyyukino 2026-03-07	3/150	2026-03-08 23:07 by 清風(fēng)月
[考研] 安徽農(nóng)業(yè)大學(xué)材料與化學(xué)學(xué)院0856材料招收調(diào)劑 +3	akakk47 2026-03-05	9/450	2026-03-08 21:03 by deng.z.1314
[考研] 一志愿中科大能動(dòng)297求調(diào)劑，本科川大 +4	邵11 2026-03-03	4/200	2026-03-08 20:16 by xiongyaxuan
[考研] 化學(xué)工程求調(diào)劑 +12	化工人999 2026-03-04	12/600	2026-03-08 11:17 by 我的船我的海
[考研] 308求調(diào)劑 +7	倘若起風(fēng)了呢 2026-03-05	9/450	2026-03-06 22:59 by 倘若起風(fēng)了呢
[考研] 085602 293分求調(diào)劑 +3	SivanNano. 2026-03-05	3/150	2026-03-05 19:41 by caszguilin
[考研] 紡織、生物、化學(xué)、材料等專(zhuān)業(yè) +3	Eember. 2026-03-05	7/350	2026-03-05 16:35 by Eember.
[考研] 化學(xué) 0703求調(diào)劑總分293 一志愿211 +9	土土小蟲(chóng) 2026-03-03	11/550	2026-03-05 09:26 by vfdbghn
[考研] 307求調(diào)劑 +4	xuyuan111 2026-03-04	5/250	2026-03-04 20:39 by kakakapanpan
[考研] 266求調(diào)劑 +7	哇塞王帥 2026-03-03	7/350	2026-03-04 20:30 by 一切OK
[考研] 成績(jī)276，專(zhuān)業(yè)代碼0856求調(diào)劑 +10	小陳朵 2026-03-03	10/500	2026-03-04 20:16 by 1124861332
[考研] 298求調(diào)劑 +3	人間唯你是清歡 2026-03-03	4/200	2026-03-03 18:09 by ms629