声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: LinSarah2010

[分形与混沌] 请教无水等高手混沌问题。。。

[复制链接]
发表于 2010-11-8 10:17 | 显示全部楼层
本帖最后由 zdeming 于 2010-11-8 10:18 编辑

LZ的问题不知道解决没有?
觉得,程序没什么问题,也可以写成:
function dy=myode11(t,x)
%---------------------------------------------
M=2;C=2;K=2e5;A=14500;w=2.4;e=0.004;
t1=t*w/pi;t2=fix(t1/2)*2;
if t1>t2&&t1<t2+1
    mf=A*sin(w*t);
else
    mf=0;
end
%---------------------------------------------
if x(1)>e
    mydd=x(1)-e;
elseif x(1)<-e
    mydd=x(1)+e;
elseif x(1)>-e&&x(1)<e
    mydd=0;
end
dy=zeros(2,1);
dy(1)=x(2);
dy(2)=(mf-C*x(2)-K*mydd)/M;
%---------------------------------------------

%---------------------------------------------
%  这个是做poincare的程序:
clear all;clc;
%options = odeset('RelTol',1e-7);
[t,x] = ode45(@myode11,[0:2*pi/10:800*2*pi],[0,0]);
i=3000:10:8000;
plot(x(i,1),x(i,2),'.')
%plot(x(:,1),x(:,2));
%plot(t,x(:,1));
%---------------------------------------------

但是:1)模型以及参数不知是否合理;2)量纲化了吗。
看不出什么结果来。
回复 支持 反对
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 17:47 , Processed in 0.266051 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表