幅频响应正反问题？

tudexin

非线性系统里在画幅频响应曲线的时候，有正着算和反着算问题，我正算和反算后好像图像是一样的，也不知道是不是程序的问题，请论坛里的前辈看看我的程序对不对？我感觉我的程序可能有问题，按照这个程序，正算和反算就应该是一样的好像

下面贴上我的程序，请大家帮忙看看
正着算的时候；
function zuidazhi
clear;
options=odeset('RelTol',1e-6,'AbsTol',1e-6);
for f=0.1:0.01:10  
      f
tt=1/f;
y0=[0,0,0,0];
for n=1:1500
[t,y]=ode45('qiujie411',[0,tt],y0,options,f);
y0=y(end,:);
end
for n=1:100
    [t,y]=ode45('qiujie411',[0,tt],y0,options,f);
    y0=y(end,:);
    yd=abs(y(:,1));
    yn(n)=max(yd);
end
hold on
plot(f,yn,'k.')
end
xlabel('频率/Hz','Fontsize',16);
ylabel('Xb/m','Fontsize',16');
反着算的时候：
function zuidazhi
clear;
options=odeset('RelTol',1e-6,'AbsTol',1e-6);
for f=10:-0.01:0.1
    f
tt=1/f;
y0=[0,0,0,0];
for n=1:1500
[t,y]=ode45('qiujie411',[0,tt],y0,options,f);
y0=y(end,:);
end
for n=1:100
    [t,y]=ode45('qiujie411',[0,tt],y0,options,f);
    y0=y(end,:);
    yd=abs(y(:,1));
    yn(n)=max(yd);
end
hold on
plot(f,yn,'k.')
end
xlabel('频率/Hz','Fontsize',16);
ylabel('Xb/m','Fontsize',16');
请教前辈们，我这程序可有问题？这种程序图画出来时一样的。。。。。。。

tudexin

哪位前辈懂啊，给我指导指导啊！

octopussheng

你这里说的正算和反算，是怎样的物理意义？

tudexin

octopussheng 发表于 2011-1-5 21:33

                               
登录/注册后可看大图

你这里说的正算和反算，是怎样的物理意义？

正算指的就是走0.1到10 递增的画图
反算就是走10到0.1递减的画图
这2者的程序是不是就是循环语句里的一个是for f=0.1:0.01:1，一个是for 10:-0.01:0.1这么简单的区别？感觉对于同样的f值的话 递增或者递减这样的话 不就一样了吗？因为都是稳定周期后取点

凌绝顶

tudexin 发表于 2011-1-4 14:33

                               
登录/注册后可看大图

非线性系统里在画幅频响应曲线的时候，有正着算和反着算问题，我正算和反算后好像图像是一样的，也不知道是 ...

你的y0取错了，除在计算第一个频率时的响应时根据自己需要自己取定初始值外，剩下的频率的初始值应该这样选择：以上一个频率计算出的一组稳态解作为下一个频率的初始值。
根据你的程序，你每次计算都是从零初始值开始计算的，显然这样你的两种情况画出来的是一样的结果，因为你每个频率之间并没有建立任何联系，你的正画反画无非只是打点的顺序不一样而已。

hsfy919

回复 4 # tudexin 的帖子

你是想表现出升速和降速的滞后效应吧，要把上一次循环的稳态解作为下一次循环的初值

tudexin

凌绝顶 发表于 2011-1-5 23:49

                               
登录/注册后可看大图

你的y0取错了，除在计算第一个频率时的响应时根据自己需要自己取定初始值外，剩下的频率的初始值应该这样 ...

哦 了解了

tudexin

hsfy919 发表于 2011-1-6 09:49

                               
登录/注册后可看大图

回复 4 # tudexin 的帖子

你是想表现出升速和降速的滞后效应吧，要把上一次循环的稳态解作为下一次循环的初 ...

恩 我按照新思路再来算一次