 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 907 | 回復(fù): 3 | |||
stone924新蟲 (初入文壇)
|
[求助]
1stOpt函數(shù)中含有間接變量怎么辦
|
剛接觸1stOpt就遇到了很大的困難：由于函數(shù)表達(dá)式過于復(fù)雜，不得不在其中引入間接變量（如代碼中的F11,F12等等），而運(yùn)行結(jié)果卻并不是設(shè)定的parameter。請(qǐng)高手指點(diǎn)迷津，不勝感激 結(jié)果如下： 參數(shù) 最佳估算 ---------- ------------- x1 0 x2 0 x3 0 rate1 0 rate2 0 rate3 0 r1 0 r2 0 r3 0 rs 0 代碼如下： Constant K=23.05,u=21.58,c1=4.6948,cs=0.00001,phi=0.1412,pi=3.1415926; Parameters x1,a1,x2,a2,x3,a3; Variable p,c; Function x1+x3+x3=phi; Ia1=2*pi*a1*(arccos(a1)-a1*(1-a1^2)^0.5)/(1-a1^2)^1.5; F11=6*u*Ia1/(2*pi*(3*K+4*u)); F12=6*u*(3*Ia1-4*pi)/(4*pi*(3*K+4*u)); F13=a1^2*(3-9*Ia1/(4*pi))*(6*K+2*u)/(2*(1-a1^2)*(3*K+4*u))+(6*u*Ia)/(8*pi*(3*K+4*u)); F14=0.5*((3-9*Ia1/(4*pi))*(6*K+2*u)/(2*(1-a1^2)*(3*K+4*u))-3*u*(1-Ia1/pi)/(3*K+4*u)); rate1=-p/(F11-F12*F13/(F13+F14))/K; //if rate1 <-1 // rate1=-1 //end Ia2=2*pi*a2*(arccos(a2)-a2*(1-a2^2)^0.5)/(1-a2^2)^1.5; F21=6*u*Ia2/(2*pi*(3*K+4*u)); F22=6*u*(3*Ia2-4*pi)/(4*pi*(3*K+4*u)); F23=a2^2*(3-9*Ia2/(4*pi))*(6*K+2*u)/(2*(1-a2^2)*(3*K+4*u))+(6*u*Ia)/(8*pi*(3*K+4*u)); F24=0.5*((3-9*Ia2/(4*pi))*(6*K+2*u)/(2*(1-a2^2)*(3*K+4*u))-3*u*(1-Ia2/pi)/(3*K+4*u)); rate2=-p/(F21-F22*F23/(F23+F24))/K; Ia3=2*pi*a3*(arccos(a3)-a3*(1-a3^2)^0.5)/(1-a3^2)^1.5; F31=6*u*Ia3/(2*pi*(3*K+4*u)); F32=6*u*(3*Ia3-4*pi)/(4*pi*(3*K+4*u)); F33=a3^2*(3-9*Ia3/(4*pi))*(6*K+2*u)/(2*(1-a3^2)*(3*K+4*u))+(6*u*Ia)/(8*pi*(3*K+4*u)); F34=0.5*((3-9*Ia3/(4*pi))*(6*K+2*u)/(2*(1-a3^2)*(3*K+4*u))-3*u*(1-Ia3/pi)/(3*K+4*u)); rate3=-p/(F31-F32*F33/(F33+F34))/K; L1z=1/(1-(a1*(1+rate1))^2)-(a1*(1+rate1))*arccos((a1*(1+rate1)))/(1-(a1*(1+rate1))^2)^1.5; L1x=(1-L1z)/2; L1y=L1x; R1=1/9*(1/(L1x*c1+(1-L1x)*c)+1/(L1y*c1+(1-L1y)*c)+1/(L1z*c1+(1-L1z)*c)); L2z=1/(1-(a2*(1+rate2))^2)-(a2*(1+rate2))*arccos((a2*(1+rate2)))/(1-(a2*(1+rate2))^2)^1.5; L2x=(1-L2z)/2; L2y=L2x; R2=1/9*(1/(L2x*c1+(1-L2x)*c)+1/(L2y*c1+(1-L2y)*c)+1/(L2z*c1+(1-L2z)*c)); L3z=1/(1-(a3*(1+rate3))^2)-(a3*(1+rate3))*arccos((a3*(1+rate3)))/(1-(a3*(1+rate3))^2)^1.5; L3x=(1-L3z)/2; L3y=L3x; R3=1/9*(1/(L3x*c1+(1-L3x)*c)+1/(L3y*c1+(1-L3y)*c)+1/(L3z*c1+(1-L3z)*c)); Rs=1/9*(1/(1/3*cs+(1-1/3)*c)+1/(1/3*cs+(1-1/3)*c)+1/(1/3*cs+(1-1/3)*c)); (1-(x1*(1+rate1)+x2*(1+rate2)+x3*(1+rate3)))*(cs-c)*Rs+x1*(1+rate1)*(c1-c)*R1+x2*(1+rate2)*(c1-c)*R2+x3*(1+rate3)*(c1-c)*R3=0; data; 0.008 0.1325 0.115 0.1309 0.020 0.1299 0.026 0.1290 0.040 0.1277 0.060 0.1269 |
回復(fù)此樓» 猜你喜歡
0856材料專碩353求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
293求調(diào)劑
已經(jīng)有4人回復(fù)
-
306求0703調(diào)劑一志愿華中師范
已經(jīng)有10人回復(fù)
-
【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】
已經(jīng)有7人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
材料與化工328分調(diào)劑
已經(jīng)有4人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有20人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271
已經(jīng)有6人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 1stopt為何結(jié)果是 計(jì)算值全為0
已經(jīng)有8人回復(fù)
- 1stopt自變量超多時(shí)是否有簡單語法？AND 自變量多 數(shù)據(jù)量大的話1stopt都會(huì)算的特慢么
已經(jīng)有7人回復(fù)
- 求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù)
已經(jīng)有10人回復(fù)
- 求高手給做一道題，用MATLAB或1stOpt，請(qǐng)把結(jié)果直接給出來，我的軟件就不行，謝謝
已經(jīng)有5人回復(fù)
- 1stopt提示函數(shù)字符串有誤！真心求教！
已經(jīng)有8人回復(fù)
- 求助--用1stopt擬合藍(lán)寶石襯底的GaN薄膜紅外反射光譜
已經(jīng)有6人回復(fù)
- 1stopt解方程組
已經(jīng)有12人回復(fù)
- 用1stopt軟件做微分方程組的參數(shù)估計(jì)
已經(jīng)有13人回復(fù)
- 用1stOpt進(jìn)行復(fù)數(shù)曲線擬合時(shí)，結(jié)果出錯(cuò)。急求教�。�！
已經(jīng)有5人回復(fù)
- 這個(gè)編號(hào)后為什么1stopt沒有反應(yīng)，不能執(zhí)行
已經(jīng)有5人回復(fù)
- 哪位朋友有注冊(cè)版1stOpt軟件？急需
已經(jīng)有5人回復(fù)
- 1stOPT1.5版每次仿真結(jié)果不一樣，求助！急！
已經(jīng)有3人回復(fù)
- 運(yùn)算過程有復(fù)數(shù)出現(xiàn)，求高版本1stOpt
已經(jīng)有6人回復(fù)
- 1stopt運(yùn)行程序，總是出現(xiàn)字符函數(shù)錯(cuò)誤的提示，請(qǐng)高手幫助
已經(jīng)有5人回復(fù)
- 大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒有遇到過這種情況？
已經(jīng)有8人回復(fù)
- 1stopt4階非線性常微分方程，幫小弟看一下，跪求。
已經(jīng)有7人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- 麻煩高人指點(diǎn)一下1stopt和spss的非線性回歸中weibull的兩參數(shù)的公式應(yīng)該怎么表達(dá)��？
已經(jīng)有9人回復(fù)
feixiaolin
榮譽(yù)版主 (文壇精英)
-
專家經(jīng)驗(yàn): +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3430
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時(shí)
- 蟲號(hào): 2139575
- 注冊(cè): 2012-11-21
- 專業(yè): 光學(xué)信息獲取與處理
- 管轄: 數(shù)學(xué)
2樓2014-07-31 16:20:25
stone924
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 19.5
- 帖子: 2
- 在線: 57分鐘
- 蟲號(hào): 1345084
- 注冊(cè): 2011-07-13
- 專業(yè): 礦床學(xué)
3樓2014-08-07 21:53:42
feixiaolin
榮譽(yù)版主 (文壇精英)
-
專家經(jīng)驗(yàn): +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3430
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時(shí)
- 蟲號(hào): 2139575
- 注冊(cè): 2012-11-21
- 專業(yè): 光學(xué)信息獲取與處理
- 管轄: 數(shù)學(xué)
4樓2014-08-07 21:56:15
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 stone924 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +10 | Zhao anqi 2026-03-22 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +5 | niko- 2026-03-22 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +10 | 梁梁校校 2026-03-17 | 10/500 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|