这个空间需要大家的积极维护，希望大家能将一些好的程序贴上，以期与大家共享！
管理 给我的留言
-
qingzheng25 留言于2008-05-11 20:12:41
- 李老师：
我抄你的程序Matlab code for poltting the Poicanre sections of autonomous chaotic systems 出错Input argument "x" is undefined. 为什么 帮我改一下，谢谢
-
lilly_chen 留言于2008-04-29 15:15:37
- 李老师：
最近想画duffing方程中Lyapunov指数随驱动力幅值r变化的图形，从网上找了计算Lyapunov指数的程序经修改后总得不到预期的效果。如果您有这样的程序，能够发给我一份，谢谢了
-
linzx270 留言于2008-04-24 20:52:11
- 您好:
我在用matlab画时滞系统图遇到了一个问题想请教您:
m文件如下:function dydt = dde_chenf(t,y,Z,a,b,c)
ylag = Z(:,1);
dydt = [ -a*y(1)+a*y(2)-ylag(1)
(c-a)*y(1)+c*y(2)-ylag(2)-y(1)*y(3)
(-b)*y(3)+ylag(3)+y(1)*y(2) ];
程序:
a=35;
b=3;
c=28;
tau=2;
tspan = [0,25];
opts = ddeset('RelTol',1e-5,'AbsTol',1e-8);
history = [-3;-4;2];
t=linspace(0,50,length(sol.y(1,4000:end)));
sol = dde23(@dde_chenf,tau,history,tspan,opts,a,b,c)
plot(怎么写?)
问题就是怎么画x(t-tau)-x(t)
先谢谢了 您可以给我回信!
linzx270@nenu.edu.cn
-
yibanlixue 留言于2008-04-21 10:57:39
- 李老师：
您好，导师要求画最大Lyapunov指数随参数变化图，一些文献上是用wolf，但它是针对离散时间序列的，请问你有针对转子系统的程序吗？多谢了，9922633@163.com
-
wlm_shooker 留言于2008-04-14 22:17:13
- 李老师：
最近想画Lyapunov指数随参数变化图，从网上找了计算Lyapunov指数的程序经修改后总得不到预期的效果----和已有文献差别很大，如果您有这样的程序，能够发给我一份，学习学习，先谢了。我和您联系过。王立明
邮箱：wlm_shooker@yahoo.com.cn
-
a1234 留言于2007-10-11 09:52:39
- ★国内外公开发行的正规期刊 发表论文请提前联系
本刊是国家新闻出版总署批准在国内外公开发行的综合性科技经济类期刊，半月刊，双刊号，国际标准大16开本。
内设：建筑工程、经济论谈、创新教育、信息科学、医药与保健等栏目。
投稿邮箱: hljkjxx@sina.com
QQ：63413536
咨询电话：13349312686
联系人：王兴杰 主任
来稿要求：作者必须提供论文的题目、作者、单位、邮编、所在省、市、摘要、关键词
本刊对采用的稿件有文字上的删改权，不同意删改者，请在来稿中申明。（来稿在3000字以内为佳）
-
- ★国内外公开发行的正规期刊 发表论文请提前联系
本刊是国家新闻出版总署批准在国内外公开发行的综合性科技经济类期刊，半月刊，双刊号，国际标准大16开本。
内设：建筑工程、经济论谈、创新教育、信息科学、医药与保健等栏目。
投稿邮箱: hljkjxx@sina.com
QQ：63413536
咨询电话：13349312686
联系人：王兴杰 主任
来稿要求：作者必须提供论文的题目、作者、单位、邮编、所在省、市、摘要、关键词
本刊对采用的稿件有文字上的删改权，不同意删改者，请在来稿中申明。（来稿在3000字以内为佳）
-
xueyongzhou 留言于2007-08-28 09:58:09
- 最近我在研究系统分支.
能帮我画个分支图吗?
function dy = rigid(t,y)
global r a2 beta c c1 c2 k k1 k2
dy=zeros(3,1); % a column vector
dy(1)=r.*y(1).*(1-y(1)/k-y(2)/k)-beta.*y(1).*y(2);
dy(2)=beta.*y(1).*y(2)-c.*y(2)-c1.*y(3).*y(2)/(y(2)+k1);
dy(3)=(a2-c2.*y(3)/(y(2)+k2)).*y(3);
clc;clear all
global r a2 beta c c1 c2 k k1 k2
r=5;a2=1.2;beta=8.50;c=1.5;c1=2;c2=1;k=3;k1=0.6;k2=0.5;
options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-4]);
[t,y]=ode45('epidemic',[0 100],[0.5 1 0.8],options);
figure(1);
plot(t,y(:,1),'b',t,y(:,2),'r');
figure(2);
plot(t,y(:,3),'g');
figure(3);
plot3(y(:,1),y(:,2),y(:,3),'r');
先谢谢了.
阿勇
