双频激励分岔图用什么方法画啊

bohua1208

双频激励分岔图用什么方法画啊，能用闪频法吗，可以的话用哪个频率啊

无水1324

不可以！！！

bohua1208

回复 2 # 无水1324 的帖子

那应该用什么方法啊，请指点?

meiyongyuandeze

用一般的方法吧，比如找最大值方法及可以吧

bohua1208

回复 4 # meiyongyuandeze 的帖子

用poincare截面法，取x=y可以吗

meiyongyuandeze

庞加莱截面要求要与相轨线相横截，不能出现截面与轨线相切的情况。具体你的问题取X=Y是不是处处横截，我也不好下定论！不过你可以尝试下！

meiyongyuandeze

回复 5 # bohua1208 的帖子

你好，我刚才和我们以为做非线性的教授交流了一下！用闪频法做你的问题可能还是有些问题的，有两个外力作用时而且有非线性作用时，选取某个频率计算分岔图是很有可能会滤掉很多的信息，建议你用最大值方法做分岔图！

bohua1208

回复 7 # meiyongyuandeze 的帖子

好的，谢了，你能提供最大值方法做非自治系统分叉图的程序吗

bohua1208

回复 5 # bohua1208 的帖子

怎么判断是否处处横截呢？

meiyongyuandeze

回复 9 # bohua1208 的帖子

【再一次总结帖】怎么样才能绘制出最准确、最精确的分岔图？？
http://forum.vibunion.com/thread-101072-1-1.html

论坛上由很多这样的程序，上面的帖子有一个程序，你可根据相应的方程记性修改！

这是Rossler系统的分岔程序！
function u =rosser(t,x)
global  c
a=0.2; b=0.2;
x0=[0.5,0.4,0.3]';
u =[-x(2)-x(3), x(1)+a*x(2), b+x(3)*(x(1)-c)]';

clear all
global c
zhang=[];
M=[2.5:0.001:6.5];
counter=1;
for counter=1:length(M)
       c=M(counter);
    i=2;
xmax=0;
xmaxold=0;
frmdata=[];
error=0;
tspan=[0 300];
var=1;
y0=[0.5;0.3;0.2];
[t,x]=ode45(@rosser,tspan,y0);
while i < size(x,1)
        if x(i-1,var) < x(i,var) & x(i+1,var) <= x(i,var)
            xmax=x(i,var);
            if xmaxold ~= 0
                frmdata=[frmdata ; xmax xmaxold];
            end
                xmaxold=xmax;
        end
        i=i+1;
end
r= length(frmdata)-20:length(frmdata);
fradata1(1,r)=frmdata(r,1);
zhang=[zhang;fradata1];
end
plot(M,zhang,'k.','markersize',1);
        xlabel(sprintf('c'));
ylabel(sprintf('x'))

bohua1208

回复 10 # meiyongyuandeze 的帖子

这程序是自治系统的吧，还有运行你给的这个程序时一直busy，没结果啊。

meiyongyuandeze

你自己调试下，给你的程序是给做参考的，可能有些许问题，又不是针对你的问题专门写的程序！

a139

路过，支持一下