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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

stone924

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 19.5
帖子: 2
在線: 57分鐘
蟲號: 1345084
注冊: 2011-07-13
專業(yè): 礦床學(xué)

[求助] 1stOpt函數(shù)中含有間接變量怎么辦

剛接觸1stOpt就遇到了很大的困難：由于函數(shù)表達(dá)式過于復(fù)雜，不得不在其中引入間接變量（如代碼中的F11,F12等等），而運(yùn)行結(jié)果卻并不是設(shè)定的parameter。請高手指點(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ù)此樓

» 猜你喜歡

284求調(diào)劑已經(jīng)有10人回復(fù)
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑已經(jīng)有4人回復(fù)
07化學(xué)280分求調(diào)劑已經(jīng)有4人回復(fù)
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有12人回復(fù)
求材料，環(huán)境專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有7人回復(fù)
一志愿吉大化學(xué)322求調(diào)劑已經(jīng)有4人回復(fù)
環(huán)境學(xué)碩288求調(diào)劑已經(jīng)有8人回復(fù)
341求調(diào)劑(一志愿湖南大學(xué)070300) 已經(jīng)有6人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

1stopt為何結(jié)果是計(jì)算值全為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擬合藍(lán)寶石襯底的GaN薄膜紅外反射光譜已經(jīng)有6人回復(fù)
1stopt解方程組已經(jīng)有12人回復(fù)
用1stopt軟件做微分方程組的參數(shù)估計(jì) 已經(jīng)有13人回復(fù)
用1stOpt進(jìn)行復(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ù)
運(yùn)算過程有復(fù)數(shù)出現(xiàn)，求高版本1stOpt 已經(jīng)有6人回復(fù)
1stopt運(yùn)行程序，總是出現(xiàn)字符函數(shù)錯誤的提示，請高手幫助已經(jīng)有5人回復(fù)
大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時候有沒有遇到過這種情況？已經(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ù)

1樓 2014-07-31 14:24:18

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

feixiaolin

榮譽(yù)版主 (文壇精英)

專家經(jīng)驗(yàn): +518
應(yīng)助: 942 (博后)
貴賓: 1.275
金幣: 3430
散金: 58785
紅花: 532
沙發(fā): 11
帖子: 24215
在線: 2601.8小時
蟲號: 2139575
注冊: 2012-11-21
專業(yè): 光學(xué)信息獲取與處理
管轄: 數(shù)學(xué)

ConstStr 定義中間變量

贊一下

回復(fù)此樓

2樓2014-07-31 16:20:25

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

stone924

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 19.5
帖子: 2
在線: 57分鐘
蟲號: 1345084
注冊: 2011-07-13
專業(yè): 礦床學(xué)

引用回帖:

2樓: Originally posted by feixiaolin at 2014-07-31 16:20:25
ConstStr 定義中間變量

能詳細(xì)一點(diǎn)嗎，太謝謝了

贊一下

回復(fù)此樓

3樓2014-08-07 21:53:42

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

feixiaolin

榮譽(yù)版主 (文壇精英)

專家經(jīng)驗(yàn): +518
應(yīng)助: 942 (博后)
貴賓: 1.275
金幣: 3430
散金: 58785
紅花: 532
沙發(fā): 11
帖子: 24215
在線: 2601.8小時
蟲號: 2139575
注冊: 2012-11-21
專業(yè): 光學(xué)信息獲取與處理
管轄: 數(shù)學(xué)

引用回帖:

3樓: Originally posted by stone924 at 2014-08-07 21:53:42
能詳細(xì)一點(diǎn)嗎，太謝謝了...

ConstStr F11=6*u*Ia1/(2*pi*(3*K+4*u));
代替
Function F11=6*u*Ia1/(2*pi*(3*K+4*u));

贊一下

回復(fù)此樓

4樓2014-08-07 21:56:15

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 stone924 的主題更新

返回列表

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

[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 熱情沙漠
[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 269求調(diào)劑 +4	我想讀研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 303求調(diào)劑 +4	元夕元 2026-03-20	4/200	2026-03-23 19:00 by macy2011
[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 北科281學(xué)碩材料求調(diào)劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 求調(diào)劑院校信息 +6	CX 330 2026-03-21	6/300	2026-03-22 15:25 by 無懈可擊111
[考研] 318求調(diào)劑 +4	plum李子 2026-03-21	7/350	2026-03-22 14:17 by ColorlessPI
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 085600材料與化工306 +4	z1z2z3879 2026-03-21	4/200	2026-03-21 23:44 by ms629
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑 +11	吃吃吃才有意義 2026-03-19	11/550	2026-03-21 18:23 by 學(xué)員8dgXkO
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-21 00:39 by JourneyLucky
[考研] 288求調(diào)劑 +16	于海海海海 2026-03-19	16/800	2026-03-20 22:28 by JourneyLucky
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

stone924

[求助] 1stOpt函數(shù)中含有間接變量怎么辦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

feixiaolin

stone924

feixiaolin

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频