| 查看: 1558 | 回复: 1 | |||
[求助]
常微分方程的参数拟合 已有1人参与
|
|
参数拟合一直出错，但是找不出错误。哪位能帮忙给解决一下？先谢谢了！ function Par_SI clear clc ID=[1 3 25 72 222 282 256 233 189 123 70 25 11 4]; X0=[762,1];%初始条件 lb=[0,0];%约束条件，参数值的下界 ub=[0.1,1];%约束条件，参数值的下界 par1guess=[0.01 0.1];%参数值的初始 options=optimset('Display','final','MaxIter',2000,'MaxFunEvals',2000);%求最优解的优化条件 [par1,fval]=fmincon(LSmin,par1guess,[],[],[],[],lb,ub,[],options,ID,X0); %利用优化算法得到的参数拟合数据 [T1,X1]=ode45(SImodel,[0 14],X0,[],par1);%利用估计参数重新运行模型 figure(1) plot(T1,X1(:,1),'k-'); hold on plot(T1,X1(:,2),'k-',0:14,ID,'ko') %定义最小二乘的目标函数 function d=LSmin(par1,ID,X0) [T,x]=ode45(SImodel,[0:1:14],X0,[],par1) d=norm(x(:,2)'-ID)^2; %定义微分方程 function dy=SImodel(t,y,par1) dy=zeros(2,1); beta=par1(1); gamma=par1(2); dy=[-beta*y(1)*y(2);beta*y(1)*y(2)-gamma*y(2)]; |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有168人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教一下各位大神怎么求解 Runge–Kutta method的微分方程的各个系数
已经有12人回复
- 1stopt软件拟合非线性方程，计算结果全为0，这是为什么
已经有3人回复
- 求大神帮你拟合一组数据~万分感激。。。
已经有6人回复
- 使用matlab最优化方法拟合获得多个动力学参数中的问题
已经有4人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 用1stopt软件做微分方程组的参数估计
已经有13人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 通过实验数据拟合，求解公式中的参数
已经有14人回复
- 紧急求助，利用Matlab对实验数据进行拟合求解参数。
已经有27人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 微分方程组参数拟合的问题
已经有11人回复
- 关于matlab微分、及非线性拟合问题
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 回归和拟合有啥区别？
已经有13人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- 求解四阶微分方程组
已经有14人回复
- matlab-常微分方程参数估计
已经有12人回复
- lsqnonlin函数拟合微分方程组参数拟合问题
已经有10人回复
- 【求助】多元方程组求解及拟合
已经有9人回复
- 重金加人民币求拟合wilson方程和NRTL方程参数的程序
已经有17人回复
- 【求助】常微分方程组的解析解
已经有3人回复
- 【求助】最小二乘法求参数
已经有6人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-08-16 16:59:14
20