| 查看: 265 | 回复: 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人回复
-
什么时候开奖？
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有194人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
我要飞
铁虫 (正式写手)
- 应助: 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
- 专业: 污染控制化学
- 管辖: 计算模拟
2楼2011-11-18 10:18:17
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
