ssmile37

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 14.4
帖子: 87
在線: 36.8小時
蟲號: 1416175
注冊: 2011-09-25
性別: MM
專業(yè): 熱力學(xué)

[求助] 1stopt運(yùn)行程序，總是出現(xiàn)字符函數(shù)錯誤的提示，請高手幫助

不知是哪里出現(xiàn)了錯誤，點(diǎn)擊運(yùn)行就會出現(xiàn)錯誤提示，現(xiàn)將程序和數(shù)據(jù)發(fā)上來，請大家?guī)兔Σ檎蚁洛e誤
Title 2;
Parameters a,a1,b,b1,c,c1,d,d1,e,e1;
Variable x1,x2,x3,x4,n2,x5,T;
conststr tao145=a+a1*(T-298.15);
conststr tao514=b+b1*(T-298.15);
conststr tao135=c+c1*(T-298.15);
conststr tao513=d+d1*(T-298.15);
conststr tao2414=e+e1*(T-298.15);
conststr tao1424=-e-e1*(T-298.15);
conststr tao4232=-1.0268-9.008*10^-3*(T-298.15);
conststr tao3242=1.0268+9.008*10^-3*(T-298.15);
conststr tao235=-5.7508+5.37*10^-3*(T-298.15);
conststr tao523=11.0427-8.159*10^-3*(T-298.15);
conststr tao245=-3.9126;
conststr tao524=8.2457;
conststr g15=(x3*exp(-0.2*tao135)+2*x4*exp(-0.2*tao145))/(x3+2*x4);
conststr g25=(x3*exp(-0.2*tao235)+2*x4*exp(-0.2*tao245))/(x3+2*x4);
conststr g35=(x1*exp(-0.2*tao135)+x2*exp(-0.2*tao235))/(x1+x2);
conststr g45=(x1*exp(-0.2*tao145)+x2*exp(-0.2*tao245))/(x1+x2);
conststr tao15=-ln(g15)/0.2;
conststr tao25=-ln(g25)/0.2;
conststr tao35=-ln(g35)/0.2;
conststr tao45=-ln(g45)/0.2;
conststr tao5131=tao15-tao135+tao513;
conststr tao5141=tao15-tao145+tao514;
conststr tao5424=tao45-tao245+tao524;
conststr tao5313=tao35-tao135+tao513;
conststr tao5414=tao45-tao145+tao514;
conststr tao5232=tao25-tao235+tao523;
conststr tao5242=tao25-tao245+tao524;
conststr tao5414=tao45-tao145+tao514;
conststr tao5424=tao45-tao245+tao524;
conststr tao5323=tao35-tao235+tao523;
conststr lnr11=x3/(x3+2*x4)*x5*exp(-0.2*tao5131)*tao5131/(x3+x5*exp(-0.2*tao5131))+2*x4/(2*x4+x3)*x5*exp(-0.2*tao5141)/(2*x4+x5*exp(-0.2*tao5141))+x5*g15/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5)*(tao15-(x1*g15*tao15+x2*g25*tao25+x3*g35*tao35+2*x4*g45*tao45)/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5))-x1/(x1+x2)*x3/(x1+x5*exp(-0.2*tao5313))*x5*exp(-0.2*tao5313)*tao5313/(x1+x5*exp(-0.2*tao5313))-x1/(x1+x2)*2*x4/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414))*(x2*tao2414*exp(-0.2*tao2414)+x5*tao5414*exp(-0.2*tao5414))/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414))+x2/(x1+x2)*2*x4*exp(-0.2*tao1424)/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424))*(tao1424-(x1*tao1424*exp(-0.2*tao1424)+x5*tao5424*exp(-0.2*tao5424))/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424)));
conststr lnr21=x3/(x3+2*x4)*(2*x4*tao4232*exp(-0.2*tao4232)+x5*tao5232*exp(-0.2*tao5232))/(x3+2*x4*exp(-0.2*tao4232)+x5*exp(-0.2*tao5232))+2*x4/(x3+2*x4)*(x3*tao3242*exp(-0.2*tao3242)+x5*tao5242*exp(-0.2*tao5242))/(x3*exp(-0.2*tao3242)+2*x4+x5*exp(-0.2*tao5242))+x5*g25/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5)*(tao25-(x1*tao15*g15+x2*g25*tao25+x3*g35*tao35+2*x4*g45*tao45)/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5))-x2/(x1+x2)*x3/(x2+x5*exp(-0.2*tao5323))*x5*tao5323*exp(-0.2*tao5323)/(x2+x5*exp(-0.2*tao5323))+x1/(x1+x2)*2*x4*exp(-0.2*tao2414)/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414))*(tao2414-(x2*tao2414*exp(-0.2*tao2414)+x5*tao5414*exp(-0.2*tao5414))/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414)))-x2/(x1+x2)*2*x4/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424))*(x1*tao1424*exp(-0.2*tao1424)+x5*tao5424*exp(-0.2*tao5424))/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424));
conststr A=-61.44534*exp(T/273.15-1)+2.864468*(exp(T/273.15-1))^2+183.5379*ln(T/273.15)-0.6820223*(T-273.15)+0.0007875695*(T^2-(273.15)^2)+58.95788*(273.15/T);
conststr I=0.5*(x1+x2+x3+4*x4);
conststr M=(n2*98+x5*18)/(n2+x5);
conststr lnr12=-sqrt(1000/M)*A*(2/14.9*ln(1+14.9*sqrt(I))+(sqrt(I)-2*I^1.5)/(1+14.9*sqrt(I)));
conststr lnr22=-sqrt(1000/M)*A*(2/14.9*ln(1+14.9*sqrt(I))+(sqrt(I)-2*I^1.5)/(1+14.9*sqrt(I)));
conststr r1=exp(lnr11+lnr12);
conststr r2=exp(lnr21+lnr22);
Function 0.189=r1*x1/(r2*x2);
Data;
0.01073 0.003444 0.00477 0.01417 0.0009472 0.9798 289.10
0.01166 0.002988 0.004278 0.01465 0.0009463 0.9789 292.65
0.01249 0.002595 0.003825 0.01508 0.0009455 0.9781 296.55
0.01459 0.001654 0.002626 0.01624 0.0009435 0.9760 307.05
0.01621 0.0009806 0.001647 0.01719 0.0009419 0.9744 317.75
0.01712 0.0006286 0.001079 0.01774 0.0009411 0.9735 323.05
0.01804 0.0002739 0.00049 0.01831 0.0009402 0.9726 328.95
0.01924 5.184E-13 9.564E-13 0.01888 0.0009390 0.9714 336.65

回復(fù)此樓

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有250人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有62人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)

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

1stOpt的函數(shù)作圖不靠譜嗎？已經(jīng)有4人回復(fù)
用origin或1stopt擬合一條隱函數(shù)方程曲線已經(jīng)有4人回復(fù)

只要相信，總會有奇跡

1樓 2012-10-10 22:02:44

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

dingd

鐵桿木蟲 (職業(yè)作家)

計算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 歡迎交流 2012-10-11 09:39:15
ssmile37: 金幣+5, ★有幫助, 謝謝參與 2012-10-12 20:51:08

其它的沒細(xì)看，但這句：
Function 0.189=r1*x1/(r2*x2);

對擬合而言有問題吧？等式左邊應(yīng)該是因變量出現(xiàn)才對。

贊一下

回復(fù)此樓

2樓2012-10-11 08:46:51

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

ssmile37

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 14.4
帖子: 87
在線: 36.8小時
蟲號: 1416175
注冊: 2011-09-25
性別: MM
專業(yè): 熱力學(xué)

引用回帖:

2樓: Originally posted by dingd at 2012-10-11 08:46:51
其它的沒細(xì)看，但這句：
Function 0.189=r1*x1/(r2*x2);

對擬合而言有問題吧？等式左邊應(yīng)該是因變量出現(xiàn)才對。

嗯，因變量在左邊也有同樣的問題

贊一下

回復(fù)此樓

只要相信，總會有奇跡

3樓2012-10-11 11:35:09

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

dingd

鐵桿木蟲 (職業(yè)作家)

計算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

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

★ ★ ★ ★ ★
ssmile37: 金幣+5, ★有幫助 2012-10-12 20:50:43

把正確的擬合公式貼出來，指明哪個是自變量、哪個是因變量。

贊一下

回復(fù)此樓

4樓2012-10-11 15:06:32

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

ssmile37

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 14.4
帖子: 87
在線: 36.8小時
蟲號: 1416175
注冊: 2011-09-25
性別: MM
專業(yè): 熱力學(xué)

引用回帖:

4樓: Originally posted by dingd at 2012-10-11 15:06:32
把正確的擬合公式貼出來，指明哪個是自變量、哪個是因變量。

已經(jīng)解決了，謝謝

回復(fù)此樓

只要相信，總會有奇跡

5樓2012-10-12 20:50:33

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

wdhua2008

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 24.7
散金: 6
紅花: 1
帖子: 30
在線: 22.6小時
蟲號: 1620061
注冊: 2012-02-16
專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域

引用回帖:

5樓: Originally posted by ssmile37 at 2012-10-12 20:50:33
已經(jīng)解決了，謝謝...

請問是怎么解決的，我也遇到了相同的問題

贊一下

回復(fù)此樓

路漫漫其修遠(yuǎn)兮，上下而求索！

6樓2014-02-05 11:05:13

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

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

[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +18	rare12345 2026-03-18	18/900	2026-03-19 22:36 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3	@taotao 2026-03-19	3/150	2026-03-19 14:07 by peike
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 材料專碩英一數(shù)二306 +5	z1z2z3879 2026-03-18	5/250	2026-03-19 07:43 by BruceLiu320
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無懈可擊111
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實(shí)用型專利一項省級立項 +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen