|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
小弟近来很疑惑,作出的频谱和功率谱在幅值上差距很大,不知道是什么,肯定高手解惑
功率谱
下面是程序:
y0=[0;0];Omega=0.8883;T1=2*pi/Omega;T=round(T1);tt=T1/100;max=500*T1;
[t,x]=ode45('vdp',[0:tt:max],y0);
y=x(:,1);dy=x(:,2);
h=x(round(((max/tt)+1)/2):end,1);nfft=2^nextpow2(length(h));hh=2*(abs(fft(h,nfft)))/length(h);fs=2*pi/tt;f=(0:(nfft-1))/nfft*fs;%求频率分布
figure;plot(f,hh);xlim([0,6]);xlabel('频率');ylabel('幅值');title('频谱');
h=x(round(((max/tt)+1)/2):end,1);nfft=2^nextpow2(length(h));hh=2*(abs(fft(h,nfft)).^2)/length(h);fs=2*pi/tt;f=(0:(nfft-1))/nfft*fs;%求频率分布
figure;plot(f,hh);xlim([0,6]);xlabel('频率');ylabel('幅值');title('功率谱')
function fy=vdp(t,x);
a=0.2;b=0.888523;c=0.394737;Omega=0.8883;g=0.2;
fy=[x(2);-a*x(2)-b^2*x(1)-3*c*x(1)^2-c*x(1)^3-g*cos(Omega*t)*x(1)-g*cos(Omega*t)];
end
|
|