【答案】应助回帖

» 猜你喜欢

申请26博士已经有5人回复
职称评审没过，求安慰已经有22人回复
垃圾破二本职称评审标准已经有15人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有20人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
聘U V热熔胶研究人员已经有10人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复
三无产品还有机会吗已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

Origin自定义模拟与matlab中cftool工具箱拟合结果比较已经有4人回复
通过Matlab软件自行计算了一个自由曲面拟合成非球面多项式或者zernike多项式的问题已经有3人回复
关于使用origin数据拟合的问题请教。已经有6人回复
【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助已经有16人回复

我就是我

1楼2015-11-25 19:35:56

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +1

两组值在哪儿？公式最好先给出文本格式

2楼2015-11-25 21:00:47

liuliping3

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 60.8
散金: 50
帖子: 54
在线: 36.6小时
虫号: 2125053
注册: 2012-11-13
性别: MM
专业: 分离过程

引用回帖:

2楼: Originally posted by dingd at 2015-11-25 21:00:47
两组值在哪儿？公式最好先给出文本格式

x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24];
y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0];

y=exp(k*P*Y*M*t/S)/(exp(k*P*Y*M*t/S)+exp(k*ρ*R*T*x/u)-1);

P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29

我就是我

3楼2015-11-25 21:56:49

liuliping3

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 60.8
散金: 50
帖子: 54
在线: 36.6小时
虫号: 2125053
注册: 2012-11-13
性别: MM
专业: 分离过程

引用回帖:

2楼: Originally posted by dingd at 2015-11-25 21:00:47
两组值在哪儿？公式最好先给出文本格式

x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24];
y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0];

y= exp(k*P*Y*M*t/S)/( exp(k*P*Y*M*t/S)+exp(k*ρ*R*T*x/u)-1);

P=101325, Y=0.0112, M=0.032, t=4/60, S=9.08, ρ=2857.14， R=8.314，T=973，u=815.29

我就是我

4楼2015-11-25 21:59:44

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

不是两组x, y数值吗？对应不同的P, Y, M, t, S, ρ， R，T，u值

如果只有这一组，1stOpt试试：

CODE:

Constant P=101325, Y1=0.0112, M=0.032, t1=4/60, S=9.08, p1=2857.14, R=8.314,T=973,u=815.29;

Function y= exp(k*P*Y1*M*t1/S)/( exp(k*P*Y1*M*t1/S)+exp(k*p1*R*T*x/u)-1);

Data;

x=[0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 8.7 10.4 12.1 13.8 15.5 17.2 18.9 20.6 22.3 24];

y=[0.96564 0.96454 0.95738 0.9543 0.82423 0.6315 0.51949 0.36035 0.2554 0.20826 0.14967 0.06784 0.05 0.01123 0.00749 0.000220264 0 0 0 0 0 0 0 0];

均方差(RMSE):0.13879469945569
残差平方和(SSE):0.462335246327889
相关系数(R): 0.958172545145573
相关系数之平方(R^2): 0.918094626270745
确定系数(DC): 0.860178141633635
卡方系数(Chi-Square): 1.00096945974885

参数最佳估算
---------- -------------
k 8.31057566523805E-6

赞一下(1人)

5楼2015-11-25 22:46:38

liuliping3

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 60.8
散金: 50
帖子: 54
在线: 36.6小时
虫号: 2125053
注册: 2012-11-13
性别: MM
专业: 分离过程

引用回帖:

5楼: Originally posted by dingd at 2015-11-25 22:46:38
不是两组x, y数值吗？对应不同的P, Y, M, t, S, ρ， R，T，u值

如果只有这一组，1stOpt试试：

Constant P=101325, Y1=0.0112, M=0.032, t1=4/60, S=9.08, p1=2857.14, R=8.314,T=973,u=815.29;
Function y= ...

谢谢你的回复。可能我表达的不清楚，我现在差不多有22组x,y值，P，R，p1是始终不变的，其余的都变，我本来想的是让大家帮我编一个只有一组x,y的程序，然后计算的时候我自己改变数。我没学过MATLAB，不懂，所以是可以将全部数据都上传吗？另外1stOpt也没用过，没有相关的软件。

我就是我

6楼2015-11-26 09:09:52

liuliping3

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 60.8
散金: 50
帖子: 54
在线: 36.6小时
虫号: 2125053
注册: 2012-11-13
性别: MM
专业: 分离过程

引用回帖:

6楼: Originally posted by liuliping3 at 2015-11-26 09:09:52
谢谢你的回复。可能我表达的不清楚，我现在差不多有22组x,y值，P，R，p1是始终不变的，其余的都变，我本来想的是让大家帮我编一个只有一组x,y的程序，然后计算的时候我自己改变数。我没学过MATLAB，不懂，所以是可 ...

再麻烦您一下，那个1stOpt软件拟合的时候可以画出原图和拟合图吗？另外您有1stOpt的程序没？

我就是我

7楼2015-11-26 09:23:17