| 查看: 237 | 回复: 3 | |||
我要飞铁虫 (正式写手)
|
[求助]
为何出错？
|
 请大家看看这个程序为何出错？如何修改？谢谢了！function [T, Y]=FOChen(parameters, orders, TSim, Y0) % format long e % Numerical Solution of the Fractional-Order Chen's System % % D^q1 x(t) = a(y(t)-x(t)) % D^q2 y(t) = dx(t) - x(t)z(t) + cy(t) % D^q3 z(t) = x(t)y(t) - bz(t) % % function [T, Y] = FOChen(parameters, orders, TSim, Y0) % % Input: parameters - model parameters [a, b, c, d] % orders - derivatives orders [q1, q2, q3] % TSim - simulation time (0 - TSim) in sec % Y0 - initial conditions [Y0(1), Y0(2), Y0(3)] % % Output: T - simulation time (0 : Tstep : TSim) % Y - solution of the system (x=Y(1), y=Y(2), z=Y(3)) % % Author: (c) Ivo Petras (ivo.petras@tuke.sk), 2010. % % time step: h=0.005; % number of calculated mesh points: n=100/h; %orders of derivatives, respectively: q1=0.99; q2=0.96; q3=0.92; % constants of Chen's system: a=10; b=1; c=1; d=5; % binomial coefficients calculation: cp1=1; cp2=1; cp3=1; for j=1:n c1(j)=(1-(1+q1)/j)*cp1; c2(j)=(1-(1+q2)/j)*cp2; c3(j)=(1-(1+q3)/j)*cp3; cp1=c1(j); cp2=c2(j); cp3=c3(j); end % initial conditions setting: x(1)=0.5; y(1)=0.6; z(1)=0.1; % calculation of phase portraits /numerical solution/: for i=2:n x(i)=(a*(y(i-1)-x(i-1)))*h^q1 - memo(x, c1, i); y(i)=(-d*x(i)-x(i)*z(i-1)+c*y(i-1))*h^q2 - memo(y, c2, i); z(i)=(x(i)*y(i)-b*z(i-1))*h^q3 - memo(z, c3, i); end for j=1:n Y(j,1)=x(j); Y(j,2)=y(j); Y(j,3)=z(j); end T=h:h:TSim; % [ Last edited by 我要飞 on 2011-11-18 at 10:00 ] |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有129人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
2楼2011-11-18 10:18:17
我要飞
铁虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1515.3
- 散金: 144
- 红花: 5
- 帖子: 458
- 在线: 591.2小时
- 虫号: 671362
- 注册: 2008-12-09
- 专业: 计算数学与科学工程计算
3楼2011-11-18 10:25:27
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
4楼2011-11-18 10:30:05
