| 查看: 890 | 回复: 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就没法跑起来了，求高手指点！ |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有266人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 呼叫版主，在线紧急求助，关于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
