声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1885|回复: 2

[稳定性与分岔] 做分岔图时,程序中这个语句是什么含义(含程序)

[复制链接]
发表于 2011-9-8 16:53 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 蓝云天 于 2011-9-8 21:31 编辑

各位大哥大姐,做分岔图的时候,ode45后面括号里面有一项的含义搞不清楚,已经用红色标出来了,哪位指教一下里面都是什么含义?r应该是参数,其他的数值是什么意思呢?谢谢了?
clear all
t0=[0 100];%积分时间
%bifurcation
for r=linspace(1,500,1000);  %r的变化精度
    [t,y]=ode45('Lorenz',t0,[1;1;1;16;r;4]);    [Xmax]=getmax(y(:,1));
    plot(r,Xmax,'b','markersize',1)
    hold on
    clear Xmax
end

function [Xmax] = getmax(y)
a=length(y);
j=1;
for i=(a-1)/2:a
    b=(y(i,1)-y(i-2,1))/2;
    c=(y(i,1)+y(i-2,1))/2-y(i-1,1);
    if y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)&c==0
        Xmax(j)=y(i-1,1);
        j=j+1;
    elseif y(i-2,1)<=y(i-1,1)&y(i-1,1)>=y(i,1)
        Xmax(j)=y(i-2,1)-b^2/(4*c);
        j=j+1;
    end
end
回复
分享到:

使用道具 举报

发表于 2011-9-8 20:07 | 显示全部楼层
三维的,初值[1,1,1],其它是参数赋值。对应着lorenz系统应该能看出来!!

评分

1

查看全部评分

 楼主| 发表于 2011-9-8 20:19 | 显示全部楼层
自己晕了,确实是,呵呵,刚才看了一下lorenz函数把系数也当状态向量处理了,我说怎么那么多初值。好久不用都忘了怎么回事了?谢谢lihaitao123 热心帮忙哈
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-17 15:30 , Processed in 0.057610 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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