 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 913 | 回復(fù): 3 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
1stOpt函數(shù)中含有間接變量怎么辦
|
|||
剛接觸1stOpt就遇到了很大的困難：由于函數(shù)表達式過于復(fù)雜，不得不在其中引入間接變量（如代碼中的F11,F12等等），而運行結(jié)果卻并不是設(shè)定的parameter。請高手指點迷津，不勝感激 結(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ù)此樓» 猜你喜歡
材料化工340求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿廈門大學材料工程專碩354找調(diào)劑！��！
已經(jīng)有6人回復(fù)
-
271求調(diào)劑
已經(jīng)有20人回復(fù)
-
一志愿085600中科院寧波所276分求調(diào)劑
已經(jīng)有24人回復(fù)
-
一志愿北交大材料工程總分358
已經(jīng)有4人回復(fù)
-
0856材料與化工調(diào)劑，339
已經(jīng)有11人回復(fù)
-
一志愿鄭大材料工程290求調(diào)劑
已經(jīng)有19人回復(fù)
-
266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑
已經(jīng)有9人回復(fù)
-
372分材料與化工（085600）一志愿湖南大學求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿南師大0703化學 275求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 1stopt為何結(jié)果是 計算值全為0
已經(jīng)有8人回復(fù)
- 1stopt自變量超多時是否有簡單語法？AND 自變量多 數(shù)據(jù)量大的話1stopt都會算的特慢么
已經(jīng)有7人回復(fù)
- 求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù)
已經(jīng)有10人回復(fù)
- 求高手給做一道題，用MATLAB或1stOpt，請把結(jié)果直接給出來，我的軟件就不行，謝謝
已經(jīng)有5人回復(fù)
- 1stopt提示函數(shù)字符串有誤！真心求教！
已經(jīng)有8人回復(fù)
- 求助--用1stopt擬合藍寶石襯底的GaN薄膜紅外反射光譜
已經(jīng)有6人回復(fù)
- 1stopt解方程組
已經(jīng)有12人回復(fù)
- 用1stopt軟件做微分方程組的參數(shù)估計
已經(jīng)有13人回復(fù)
- 用1stOpt進行復(fù)數(shù)曲線擬合時，結(jié)果出錯。急求教�。�！
已經(jīng)有5人回復(fù)
- 這個編號后為什么1stopt沒有反應(yīng)，不能執(zhí)行
已經(jīng)有5人回復(fù)
- 哪位朋友有注冊版1stOpt軟件？急需
已經(jīng)有5人回復(fù)
- 1stOPT1.5版每次仿真結(jié)果不一樣，求助！急！
已經(jīng)有3人回復(fù)
- 運算過程有復(fù)數(shù)出現(xiàn)，求高版本1stOpt
已經(jīng)有6人回復(fù)
- 1stopt運行程序，總是出現(xiàn)字符函數(shù)錯誤的提示，請高手幫助
已經(jīng)有5人回復(fù)
- 大家在使用1stOpt進行數(shù)據(jù)擬合的時候有沒有遇到過這種情況？
已經(jīng)有8人回復(fù)
- 1stopt4階非線性常微分方程，幫小弟看一下，跪求。
已經(jīng)有7人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- 麻煩高人指點一下1stopt和spss的非線性回歸中weibull的兩參數(shù)的公式應(yīng)該怎么表達��？
已經(jīng)有9人回復(fù)
3樓2014-08-07 21:53:42
feixiaolin
榮譽版主 (文壇精英)
-
專家經(jīng)驗: +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3880
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時
- 蟲號: 2139575
- 注冊: 2012-11-21
- 專業(yè): 光學信息獲取與處理
- 管轄: 數(shù)學
2樓2014-07-31 16:20:25
feixiaolin
榮譽版主 (文壇精英)
-
專家經(jīng)驗: +518 - 應(yīng)助: 942 (博后)
- 貴賓: 1.275
- 金幣: 3880
- 散金: 58785
- 紅花: 532
- 沙發(fā): 11
- 帖子: 24215
- 在線: 2601.8小時
- 蟲號: 2139575
- 注冊: 2012-11-21
- 專業(yè): 光學信息獲取與處理
- 管轄: 數(shù)學
4樓2014-08-07 21:56:15
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料考研調(diào)劑 +9 | Gs大王 2026-04-02 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 324求調(diào)劑 +5 | 想上學求調(diào) 2026-04-01 | 6/300 |
|
|
|
[考研] 材料調(diào)劑 +12 | 一樣YWY 2026-04-01 | 12/600 |
|
|
|
[考研] 一志愿北京科技大學085601材料工程英一數(shù)二初試總分335求調(diào)劑 +7 | 雙馬尾痞老板2 2026-04-01 | 7/350 |
|
|
|
[考研] 江蘇科技大學招材料研究生 +4 | Su032713. 2026-04-01 | 5/250 |
|
|
|
[考研] 292求調(diào)劑 +17 | 木蟲er12138 2026-04-01 | 17/850 |
|
|
|
[考研] 食品學碩362求調(diào)劑 +3 | xuanxianxian 2026-04-01 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +7 | AZMK 2026-03-30 | 13/650 |
|
|
|
[考研] 調(diào)劑 +5 | 好好讀書。 2026-03-28 | 7/350 |
|
|
|
[論文投稿]
chinese chemical letters英文版投稿求助
 120+4
|
Yishengeryi 2026-03-30 | 5/250 |
|
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +9 | 吃的不少 2026-03-30 | 9/450 |
|
|
|
[考研] 求調(diào)劑，一志愿北林食品與營養(yǎng)095500，301分，已過六級，有科研經(jīng)歷 +4 | 快樂儲蓄罐 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
|
[考研] 一志愿中海洋材料357 +4 | 麥恩莉. 2026-03-30 | 4/200 |
|
|
|
[考研]
 304材料求調(diào)劑
+5
|
鐘llll 2026-03-26 | 5/250 |
|
|
|
[考研]
105500藥學求調(diào)劑，一志愿山東大學藥學，348分
+3
|
gr哈哈哈 2026-03-28 | 3/150 |
|
|
|
[考研] 一志愿中南大學化學0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
|
[考研] 本科新能源科學與工程，一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|