24小时热门版块排行榜    

查看: 1459  |  回复: 5

ssmile37

铁虫 (小有名气)

[求助] 1stopt运行程序,总是出现字符函数错误的提示,请高手帮助

不知是哪里出现了错误,点击运行就会出现错误提示,现将程序和数据发上来,请大家帮忙查找下错误
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
回复此楼
只要相信,总会有奇迹
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
csgt0: 金币+1, 欢迎交流 2012-10-11 09:39:15
ssmile37: 金币+5, 有帮助, 谢谢参与 2012-10-12 20:51:08
其它的没细看,但这句:
Function 0.189=r1*x1/(r2*x2);

对拟合而言有问题吧?等式左边应该是因变量出现才对。
2楼2012-10-11 08:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ssmile37

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2012-10-11 08:46:51
其它的没细看,但这句:
Function 0.189=r1*x1/(r2*x2);

对拟合而言有问题吧?等式左边应该是因变量出现才对。

嗯,因变量在左边也有同样的问题
只要相信,总会有奇迹
3楼2012-10-11 11:35:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
ssmile37: 金币+5, 有帮助 2012-10-12 20:50:43
把正确的拟合公式贴出来,指明哪个是自变量、哪个是因变量。
4楼2012-10-11 15:06:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ssmile37

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2012-10-11 15:06:32
把正确的拟合公式贴出来,指明哪个是自变量、哪个是因变量。

已经解决了,谢谢
只要相信,总会有奇迹
5楼2012-10-12 20:50:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wdhua2008

铁虫 (初入文坛)

引用回帖:
5楼: Originally posted by ssmile37 at 2012-10-12 20:50:33
已经解决了,谢谢...

请问是怎么解决的,我也遇到了相同的问题
路漫漫其修远兮,上下而求索!
6楼2014-02-05 11:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ssmile37 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 27/1350 2026-08-18 21:44 by Noways
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +3 瞬息宇宙 2026-08-15 3/150 2026-08-18 21:40 by xiaohuixin
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见