| 查看: 1526 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
水中的梦88木虫 (正式写手)
|
[求助]
微分方程拟合求参数问题已有1人参与
|
|||
|
最近需要做微分方程拟合求参数的问题，但是以前并没有用过Matlab，只是学过C语言，研究了一段时间了。 照葫芦画瓢，但是始终都弄不出来。  求各位大神帮忙，需要代码。我可以再根据代码继续研究。跪求啊！ 下面是我自己弄的。再这个的基础上希望能够算出相关系数。 format long clear all clc tspan = [0 7 14 28 56 90 180]; x0 = 0.353542884; k0 = 1e-3; lb = 1e-5; ub = 1; yexp = [0.353542884 0.355002172 0.372553224 0.377598226 0.39054644 0.403335059 0.420130491]; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,x0,yexp); ci = nlparci(k,residual,jacobian); fprintf('\n\n使用函数lsqnonlin()估计得到的参数值为:\n') fprintf('\tk1 = %.6e\n',k) fprintf(' The sum of the squares is: %.6e\n\n',resnorm) function f = ObjFunc(k,tspan,x0,yexp) % 目标函数 [t Xsim] = ode45(@KineticsEqs,tspan,x0,[],k); ysim = Xsim(2:end); f = ysim-yexp; function dAdt = KineticsEqs(t,A,k) % ODE模型方程 dAdt=k*3*(1-A)^(2/3)/(2-2*(1-A)^(1/3)); |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有20人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有20人回复
-
垃圾破二本职称评审标准
已经有12人回复
-
EST投稿状态问题
已经有7人回复
-
谈谈两天一夜的“延安行”
已经有15人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
水中的梦88
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 2378.9
- 散金: 26
- 红花: 5
- 帖子: 323
- 在线: 195.7小时
- 虫号: 1308847
- 注册: 2011-05-28
- 性别: MM
- 专业: 结构工程
11楼2015-10-20 17:39:31
水中的梦88
木虫 (正式写手)
- 应助: 52 (初中生)
- 金币: 2378.9
- 散金: 26
- 红花: 5
- 帖子: 323
- 在线: 195.7小时
- 虫号: 1308847
- 注册: 2011-05-28
- 性别: MM
- 专业: 结构工程
2楼2015-10-07 09:15:50
水中的梦883楼2015-10-07 09:39:17
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68121.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +1
|
用1stopt计算，代码和结果如下： 代码： Variables t,y; OdeFunction y'=k*3*(1-y)^(2/3)/(2-2*(1-y)^(1/3)); Data; 0 0.3535 7.0000 0.3550 14.0000 0.3726 28.0000 0.3776 56.0000 0.3905 90.0000 0.4033 180.0000 0.4201 结果： Root of Mean Square Error (RMSE): 0.00878112685716906 Residual Sum of Squares (SSE): 0.000539757322171871 Correlation Coef. (R): 0.990208373887058 R-Square: 0.980512623716052 Determination Coef. (DC): 0.974279082009332 Parameters Best Estimate -------------------- -------------k 6.18772301033453E-5 ====== Output Results ===== File: Data file - 1 No Obs. y Cal. y 1 0.355 0.357062297934808 2 0.3726 0.360570784579331 3 0.3776 0.367433567235337 4 0.3905 0.380587817889355 5 0.4033 0.3956451569869 6 0.4201 0.431527943596942 |
» 本帖已获得的红花（最新10朵）
4楼2015-10-07 10:14:39