| 查看: 972 | 回复: 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就没法跑起来了，求高手指点！ |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有239人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 呼叫版主，在线紧急求助，关于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
20