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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 拟合求助：1stopt正版sharedmodel计算求助，我自己写好代码，帮求4组
51/1返回列表
查看: 1419  |  回复: 10
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

962321702

银虫 (小有名气)

[求助] 拟合求助：1stopt正版sharedmodel计算求助，我自己写好代码，帮求4组已有1人参与

代码：
parameters  k12,k13;  
constant k9=9.92471*10^(-15), k10=4.36826*10^(-7), k11=4.63956*10^(-11), h=3251.622, m=3.285;     
variable a,p;  
parvariable t=[0.000001,0.000000000001];        
sharedmodel;
function a=(p/h+k10*p/(h*t)+k10*k11*p/(h*t^2)+k10*m*(p/h)/(k12*k13*t+k10*p/h+k13*t^2))/m;
               t=(1+k13*m*t/(k12*k13*t+k10*p/h+k13*t^2))*(k10*p/(h*t)+k10*m*(p/h)/(k12*k13*t+k10*p/h+k13*t^2)+2*k10*k11*p/(h*t^2)+k9/t);
data;
0.638        34.5
0.667        46.5
0.684        60
0.713       78

对应其他三组数据：
k9          2.86026*10^(-14)       9.38162*10^(-14)      2.49148*10^(-13)
k10        5.01604*10^(-7)         5.23156*10^(-7)        4.79116*10^(-7)
k11        5.99747*10^(-11)       7.23603*10^(-11)      7.57401*10^(-11)
h           4078.309                     5344.079                    6791.551
m          3.285                           3.285                          3.285
a,p      
//0.598        34.5
0.612        46.5
0.636        60
0.66        78
  
//0.538        34.5
0.562        46.5
0.583        60
0.598        78

//0.465        34.5
0.495        46.5
0.517        60
0.533        78@月只蓝@beefly
回复此楼

» 猜你喜欢
1楼2018-10-31 14:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)
引用回帖:
7楼: Originally posted by 962321702 at 2018-11-01 20:50:30
我之前参数没有给范围，所有拟合的结果不好，能重新帮我算一下吗？
第一组数据：
Parameters  k12,k13;  //拟合有机胺质子化常数，氨基甲酸盐稳定常数
Constant k9=9.92471E-15, k10=4.36826E-07, k11=4.63956E ...

已经发给你了。
回复此楼
数值计算
8楼2018-11-01 20:56:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

1028942887
顶一下

发自小木虫Android客户端
2楼2018-10-31 15:25:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
962321702: 金币+200, ★★★★★最佳答案 2018-11-01 22:34:27
模型公式: a = (p/3251.622+4.36826E-7*p/(3251.622*t)+4.36826E-7*4.63956E-11*p/(3251.622*t^2)+4.36826E-7*3.285*(p/3251.622)/(k12*k13*t+4.36826E-7*p/3251.622+k13*t^2))/3.285
          a = (p/3251.622+4.36826E-7*p/(3251.622*t)+4.36826E-7*4.63956E-11*p/(3251.622*t^2)+4.36826E-7*3.285*(p/3251.622)/(238.258176720815*3.69008634786073e-5*t+4.36826E-7*p/3251.622+3.69008634786073e-5*t^2))/3.285
模型公式: t = (1+k13*3.285*t/(k12*k13*t+4.36826E-7*p/3251.622+k13*t^2))*(4.36826E-7*p/(3251.622*t)+4.36826E-7*3.285*(p/3251.622)/(k12*k13*t+4.36826E-7*p/3251.622+k13*t^2)+2*4.36826E-7*4.63956E-11*p/(3251.622*t^2)+9.92471E-15/t)
          t = (1+3.69008634786073e-5*3.285*t/(238.258176720815*3.69008634786073e-5*t+4.36826E-7*p/3251.622+3.69008634786073e-5*t^2))*(4.36826E-7*p/(3251.622*t)+4.36826E-7*3.285*(p/3251.622)/(238.258176720815*3.69008634786073e-5*t+4.36826E-7*p/3251.622+3.69008634786073e-5*t^2)+2*4.36826E-7*4.63956E-11*p/(3251.622*t^2)+9.92471E-15/t)

迭代数: 24
计算用时(时:分:秒:微秒): 00:00:01:647
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0
残差平方和(SSR): 0
相关系数(R): 1
相关系数之平方(R^2): 1
修正R平方(Adj. R^2): -100
确定系数(DC): 1

参数                  最佳估算
--------------------        -------------
k12        238.258176720815
k13        3.69008634786073E-5
t0        3.09325866287863E-7
t1        3.70106253183998E-7
t2        4.45708163439644E-7
t3        5.1208433542262E-7

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

文件: 数据文件-1
No        目标 a        计算 a
1        0.638        0.638
2        0.667        0.667
3        0.684        0.684
4        0.713        0.713
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

962321702 962321702
数值计算
3楼2018-10-31 15:33:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

模型公式: a = (p/4078.309+5.01604E-7*p/(4078.309*t)+5.01604E-7*5.99747E-11*p/(4078.309*t^2)+5.01604E-7*3.285*(p/4078.309)/(k12*k13*t+5.01604E-7*p/4078.309+k13*t^2))/3.285
          a = (p/4078.309+5.01604E-7*p/(4078.309*t)+5.01604E-7*5.99747E-11*p/(4078.309*t^2)+5.01604E-7*3.285*(p/4078.309)/(3.55859929088445e-10*6027.61803316927*t+5.01604E-7*p/4078.309+6027.61803316927*t^2))/3.285
模型公式: t = (1+k13*3.285*t/(k12*k13*t+5.01604E-7*p/4078.309+k13*t^2))*(5.01604E-7*p/(4078.309*t)+5.01604E-7*3.285*(p/4078.309)/(k12*k13*t+5.01604E-7*p/4078.309+k13*t^2)+2*5.01604E-7*5.99747E-11*p/(4078.309*t^2)+2.86026E-14/t)
          t = (1+6027.61803316927*3.285*t/(3.55859929088445e-10*6027.61803316927*t+5.01604E-7*p/4078.309+6027.61803316927*t^2))*(5.01604E-7*p/(4078.309*t)+5.01604E-7*3.285*(p/4078.309)/(3.55859929088445e-10*6027.61803316927*t+5.01604E-7*p/4078.309+6027.61803316927*t^2)+2*5.01604E-7*5.99747E-11*p/(4078.309*t^2)+2.86026E-14/t)

迭代数: 24
计算用时(时:分:秒:微秒): 00:00:01:616
优化算法: 通用全局优化算法(UGO1)
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0
残差平方和(SSR): 0
相关系数(R): 1
相关系数之平方(R^2): 1
修正R平方(Adj. R^2): -100
确定系数(DC): 1

参数                  最佳估算
--------------------        -------------
k12        3.55859929088445E-10
k13        6027.61803316927
t0        6.94110112812778E-7
t1        7.847386865608E-7
t2        8.49821236885364E-7
t3        9.23503826465456E-7

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

文件: 数据文件-1
No        目标 a        计算 a
1        0.598        0.598
2        0.612        0.612
3        0.636        0.636
4        0.66        0.66
赞一下 回复此楼
数值计算
4楼2018-10-31 15:33:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭