应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt提示函数字符串有误！真心求教！
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 1951  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

947683166

新虫 (初入文坛)

[求助] 1stopt提示函数字符串有误！真心求教！ 已有2人参与

我用的是1.5版本的，遇到这个问题好久没解决
这个函数表达式我排查了，y=(temp-20)^m
在这里有问题，会提示函数字符串有误，我就不明白了这么简单的个式子怎么就能有误，而且y=(temp+20)^m可以，y=(20-temp)^m也可以运行
是不是软件本身的问题，真心求教

具体的如下，数据只是一部分，太多了

Title "Type your title here";
Parameters p1,p2,p3,n,m;
Variable strain,temp,rate,y;
Function y=(p1+p2*strain^n)*(1+p3*ln(rate/0.1))*(1-(temp-20)^m/1580^m);
Data;
0.00018337        550        0.1        55.34468
0.00049922        550        0.1        59.02854
0.00045379        550        0.1        60.79586
0.00109        550        0.1        66.21667
0.00121        550        0.1        68.24785
0.00166        550        0.1        71.17947
0.00176        550        0.1        73.48994
0.0023        550        0.1        79.47316
0.0028        550        0.1        83.24459
0.00312        550        0.1        86.55612
0.00356        550        0.1        90.90448
0.00414        550        0.1        96.946
0.0042        550        0.1        99.26967
0.00505        550        0.1        104.85363
0.00545        550        0.1        109.29731
0.00573        550        0.1        113.73305
0.00643        550        0.1        118.39214
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-04-16 17:27:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
947683166: 金币+2, ★★★很有帮助 2014-04-17 23:04:09
问题在^m那里
改成这样
Title "Type your title here";
Parameters p1,p2,p3,n,m;
Variable strain,temp,rate,y;
Function y=(p1+p2*strain^n)*(1+p3*ln(rate/0.1))*(1-(temp-20)/(m^(-1))/1580^m);
Data;
0.00018337        550        0.1        55.34468
0.00049922        550        0.1        59.02854
0.00045379        550        0.1        60.79586
0.00109        550        0.1        66.21667
0.00121        550        0.1        68.24785
0.00166        550        0.1        71.17947
0.00176        550        0.1        73.48994
0.0023        550        0.1        79.47316
0.0028        550        0.1        83.24459
0.00312        550        0.1        86.55612
0.00356        550        0.1        90.90448
0.00414        550        0.1        96.946
0.0042        550        0.1        99.26967
0.00505        550        0.1        104.85363
0.00545        550        0.1        109.29731
0.00573        550        0.1        113.73305
0.00643        550        0.1        118.39214

相关系数(R): 0.999027475750151
相关系数之平方(R^2): 0.998055897303719
决定系数(DC): 0.998055897303719
卡方系数(Chi-Square): 0.077807687887097
F统计(F-Statistic): 7700.6417861533

参数        最佳估算
----------        -------------
p1                 53.638622036976
p2                 7016.04664269701
p3                 0
n                  0.926716508704312
m                  14.7890352501302

====== 结果输出 =====

No        实测值y        计算值y
1        55.34468        56.0555132
2        59.02854        59.7528994
3        60.79586        59.2354858
4        66.21667        66.2460558
5        68.24785        67.5273193
6        71.17947        72.2560954
7        73.48994        73.2931937
8        79.47316        78.8247902
9        83.24459        83.8612114
10        86.55612        87.0492133
11        90.90448        91.3941683
12        96.946        97.0623738
13        99.26967        97.6452760
14        104.85363        105.8415271
15        109.29731        109.6625735
16        113.73305        112.3250053
17        118.39214        118.9404581
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-04-16 19:07:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

947683166

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-04-16 19:07:04
问题在^m那里
改成这样
Title "Type your title here";
Parameters p1,p2,p3,n,m;
Variable strain,temp,rate,y;
Function y=(p1+p2*strain^n)*(1+p3*ln(rate/0.1))*(1-(temp-20)/(m^(-1))/1580^m); ...

非常感谢你的回帖，我看了一下这里不明白

(temp-20)/(m^(-1))/1580^m  这个连除是（temp-20）除以m的-1次方就是(temp-20)*m再除1580^m
我的那个式子是(temp-20)的m次方除以1580的m次方
不太明白，还是我理解有误，请教
赞一下 回复此楼
3楼2014-04-16 22:19:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
用6.0版运行1#的代码没任何问题，只是楼主的公式有过拟合现象，目标函数值稳定唯一，但参数组值不唯一。要么简化方程，要么给参数加上范围约束：
1：
均方差(RMSE):0.857868362167362
残差平方和(SSE):12.5109481557311
相关系数(R): 0.999027475750151
相关系数之平方(R^2): 0.998055897303719
确定系数(DC): 0.998055897303719
卡方系数(Chi-Square): 0.0778076878084914
F统计(F-Statistic): 1543.12835723077

参数        最佳估算
----------        -------------
p1        74.1260859093751
p2        9695.85060001215
p3        719.448785319653
n        0.926716502902412
m        1.17728704826398

2：
均方差(RMSE):0.857868362167359
残差平方和(SSE):12.510948155731
相关系数(R): 0.999027475750151
相关系数之平方(R^2): 0.998055897303719
确定系数(DC): 0.998055897303719
卡方系数(Chi-Square): 0.0778076879103115
F统计(F-Statistic): 1543.12835723078

参数        最佳估算
----------        -------------
p1        60.9175220909949
p2        7968.14241441522
p3        -10.9271853024159
n        0.926716510072309
m        1.94501055312701

3：
均方差(RMSE):0.85786836216736
残差平方和(SSE):12.510948155731
相关系数(R): 0.999027475750151
相关系数之平方(R^2): 0.998055897303719
确定系数(DC): 0.998055897303719
卡方系数(Chi-Square): 0.0778076878560482
F统计(F-Statistic): 1543.12835723078

参数        最佳估算
----------        -------------
p1        -75.3934868066866
p2        -9861.62944305136
p3        2796.62330171067
n        0.926716505606552
m        -0.491933401387309
赞一下(1人) 回复此楼
4楼2014-04-17 10:20:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

947683166

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2014-04-17 10:20:12
用6.0版运行1#的代码没任何问题，只是楼主的公式有过拟合现象，目标函数值稳定唯一，但参数组值不唯一。要么简化方程，要么给参数加上范围约束：
1：
均方差(RMSE):0.857868362167362
残差平方和(SSE):12.510948 ...

谢谢你的回复，6.0的版本就是厉害啊，数据有4000多行，这里给出的少，可能不是过拟合吧，式子没法简化这就是原公式数据是实验值，既然你能跑通，没出现字符串有误，说明1.5的这个版本太低，有缺陷，哎，高版本都是收费版
赞一下 回复此楼
5楼2014-04-17 14:18:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
3楼: Originally posted by 947683166 at 2014-04-16 22:19:54
非常感谢你的回帖，我看了一下这里不明白

(temp-20)/(m^(-1))/1580^m  这个连除是（temp-20）除以m的-1次方就是(temp-20)*m再除1580^m
我的那个式子是(temp-20)的m次方除以1580的m次方
不太明白，还是我理解有 ...

不好意思，我把原方程看错了。看来确实软件版本的问题。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2014-04-17 14:35:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★
947683166: 金币+3, ★★★很有帮助 2014-04-17 23:03:46
个人觉得你的公式完全可以化简：
原公式中：
y=(p1+p2*strain^n)*(1+p3*ln(rate/0.1))*(1-(temp-20)/(m^(-1))/1580^m);
因为rate和temp都是常数（550和0.1），则下面部分：
(1+p3*ln(rate/0.1))*(1-(temp-20)/(m^(-1))/1580^m)
等同于：p3*m，再代入原公式并展开:
y=p1*p3*m+p2*p3*m*strain^n
该公式有等同于：y=p1+p2*strain^n
只剩3个参数：p1, p2 和n，效果是完全一样的，参数也唯一。
赞一下(1人) 回复此楼
7楼2014-04-17 14:42:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

947683166

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by 月只蓝 at 2014-04-17 14:35:51
不好意思，我把原方程看错了。看来确实软件版本的问题。...

谢谢你的回帖和支持。
回复此楼
8楼2014-04-17 22:46:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

947683166

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by dingd at 2014-04-17 14:42:07
个人觉得你的公式完全可以化简：
原公式中：
y=(p1+p2*strain^n)*(1+p3*ln(rate/0.1))*(1-(temp-20)/(m^(-1))/1580^m);
因为rate和temp都是常数（550和0.1），则下面部分：
(1+p3*ln(rate/0.1))*(1-(temp-20)/( ...

公式是这个公式，你写那个是回帖的写错的一个，
y=(p1+p2*strain^n)*(1+p3*ln(rate/0.1))*(1-(temp-20)^m/1580^m)
实话说，简化过程没太明白，帮个忙，帮我跑跑吧，我这个估计是不行了，版本低，用mathlab会不会数据太多了，不适用。
一直想传上来，没刷出来传送文件的界面，问题结束了不知还能给解决不，谢谢。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 程序.wps
    • 2014-04-19 19:57:04, 3.13 M
9楼2014-04-19 00:04:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 947683166 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料0856 英一数二 323 求调剂 +11 袁sy 2026-04-01 11/550 2026-04-05 11:30 by 风雨无晴
[考研] 348求调剂 +5 wukira 2026-04-04 5/250 2026-04-05 11:05 by xiayan13521
[考研] 一志愿西北农林畜牧专硕336分求调剂 +3 5ourr 2026-04-03 3/150 2026-04-05 10:40 by JOKER0401
[考研] 341求调剂 +3 学无止境，冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 考研调剂 +11 小sun要好运 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 土木304求调剂 +4 兔突突突， 2026-03-31 4/200 2026-04-04 13:34 by 1753564080
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊！
[考研] 326分求调剂 +3 于是乎呢 2026-04-01 5/250 2026-04-03 14:23 by 于是乎呢
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 286求调剂 +7 Faune 2026-03-30 7/350 2026-04-03 10:14 by linyelide
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 一志愿厦门大学材料工程专硕354找调剂！！！ +8 贝呗钡钡 2026-03-30 8/400 2026-04-03 09:41 by hypershenger
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 一志愿9初试366 本双非求调剂 +4 运气来得若有似� 2026-04-02 4/200 2026-04-02 09:56 by guanxin1001
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[考研] 环境工程调剂 +9 hyzzzzzzz. 2026-04-01 9/450 2026-04-01 14:20 by salamander`
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
信息提示
关闭
请填处理意见
关闭