bohua1208 发表于 2011-4-5 22:37

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

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

无水1324 发表于 2011-4-6 08:33

不可以！！！

bohua1208 发表于 2011-4-6 08:51

回复 2 # 无水1324 的帖子

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

meiyongyuandeze 发表于 2011-4-6 09:33

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

bohua1208 发表于 2011-4-6 10:02

回复 4 # meiyongyuandeze 的帖子

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

meiyongyuandeze 发表于 2011-4-6 10:07

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

meiyongyuandeze 发表于 2011-4-6 10:37

回复 5 # bohua1208 的帖子

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

bohua1208 发表于 2011-4-6 11:07

回复 7 # meiyongyuandeze 的帖子

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

bohua1208 发表于 2011-4-6 11:09

回复 5 # bohua1208 的帖子

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

meiyongyuandeze 发表于 2011-4-6 11:18

回复 9 # bohua1208 的帖子

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

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

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

clear all
global c
zhang=[];
M=;
counter=1;
for counter=1:length(M)
       c=M(counter);
    i=2;
xmax=0;
xmaxold=0;
frmdata=[];
error=0;
tspan=;
var=1;
y0=;
=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=;
            end
                xmaxold=xmax;
      end
      i=i+1;
end
r= length(frmdata)-20:length(frmdata);
fradata1(1,r)=frmdata(r,1);
zhang=;
end
plot(M,zhang,'k.','markersize',1);
      xlabel(sprintf('c'));
ylabel(sprintf('x'))

bohua1208 发表于 2011-4-6 15:26

回复 10 # meiyongyuandeze 的帖子

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

meiyongyuandeze 发表于 2011-4-6 15:29

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

a139 发表于 2012-5-9 10:56

路过，支持一下
页: [1]
查看完整版本: 双频激励分岔图用什么方法画啊