| 查看: 961 | 回复: 1 | ||
[求助]
求助修改该matlab拟合带参数的时滞微分方程程序
|
|
clear all clc tspan = [2 4 6 8 10 12 15 20 25 30 35 40 50 60 80 100 120 140 160 180]; x0=[300; 410]; k0 = [1 1 1 1 1 1 1 1]*1e-2; lb = [0 0 0 0 0 0 0 0 0 0]; ub = [1 1 1 1 1 1 1 1 1 1]*1e5; data=... [ 2 238 410 4 224 408 6 211 315 8 202 294 10 198 243 12 190 322 15 177 282 20 160 247 25 151 253 30 135 241 40 123 209 50 111 173 60 99 121 80 85 118 100 74 47 120 88 44 140 80 60 160 77 63 180 83 58 ]; yexp = data(:,2:3); [b,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,yexp); function f = ObjFunc(b,tspan,yexp) % 目标函数 sol=dde23(@ddefun,[20],@ddehist,tspan); % 20是时滞 f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2))]; function dydt=ddefun(t,y,Z) % dde函数 ylag1=Z(:,1); dGdt=-b(1)*y(1)-(b(4)*b(8)*y(1)*y(2))/(y(1)+b(8))+b(7); dIdt=-b(2)*y(2)+b(6)*ylag1(1) dydt=[dGdt; dIdt]; function S=ddehist(t) % 初始值 S=[300;410]; 这段程序是我根据论坛上类似的ode45的帖子中的程序段改编的，但是不知道为什么函数换成了dde23就没法跑起来了，求高手指点！ |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有165人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 怎样求满足这些条件的平面曲线方程？
已经有10人回复
- 求助一个最小二乘法matlab拟合程序。
已经有7人回复
- 高分求助，用MATLAB模拟一曲线方程，已知曲线上散点
已经有3人回复
- 求助，matlab自定义微分方程拟合实验数据来求方程中的参数
已经有4人回复
- 紧急求助，利用Matlab对实验数据进行拟合求解参数。
已经有27人回复
- 微分方程组参数拟合的问题（多参数动力系统参数估计）
已经有11人回复
- 微分方程组参数拟合的问题
已经有11人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- （十万火急）切换退化时滞系统（带约束条件的微分方程）的matlab作图程序
已经有4人回复
- matlab数据拟合求助！！希望高手帮忙写一下matlab的程序 谢谢了！！
已经有15人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- 偏微分方程数值解有基本理论和一些例子（化工类）MATLAB源程序
已经有296人回复
- lsqnonlin函数拟合微分方程组参数拟合问题
已经有10人回复
- 急切求助用matlab曲线拟合色散方程
已经有10人回复
- 求助matlab隐函数参数拟合
已经有16人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】Zero-order，First-order和Higuchi扩散方程的中文意思？
已经有3人回复
2楼2015-06-12 16:27:59
投票:
