| 查看: 867 | 回复: 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就没法跑起来了，求高手指点！ |
回复此楼» 猜你喜欢
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有233人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有1人回复
-
投稿chemical physical letters不送审？
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 呼叫版主，在线紧急求助，关于matlab中微分方程组参数拟合得问题！
已经有12人回复
- 怎样求满足这些条件的平面曲线方程？
已经有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