| 查看: 875 | 回复: 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就没法跑起来了，求高手指点！ |
回复此楼» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有140人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有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