求助幅值谱图

sssssxxxxx921

clc
clear all
x0=[0;1.0;1 ];
[t,y]=ode45('duffing',[0:0.01:100],x0);
%频谱图
N=50;           %采样点数
x=y(9950:10000,1);  %生成信号（Y轴方向的位移）
Y=fft(x,N);
Pyy = Y.* conj(Y) /N;
df=fs/(N-1);          %频率分辨率Hz
Y(1)=[];
f=(1:N-1)*df;    %频域序列
figure
% plot(f,Pyy);
% figure
plot(f(1:fix(N/2)),Pyy(1:fix(N/2)));
xlabel('Frequency(HZ)');
ylabel('Magnitude');
title('功率谱图');
hold on

sssssxxxxx921

上边的是个功率谱图   功率谱和幅值谱我在上边已经说了
功率谱在很多书上都会讲到  我就不多说了，表达能力有限，不过提醒的一点是功率谱和功率谱密度是不一样的，但是一般人很少区分，它们也不是一种固定的做法，不过相同的一点是都经过傅立叶变换后平方（也就是表达功率的概念吧）
对于幅值谱与功率谱的区别是它不进行平方

秋月

谢谢指点!
采样频率应该取多少啊?
画出的图怎么会是那样啊?
和书上的一点都不一样啊!
是那有问题啊???

sssssxxxxx921

clc
clear all
x0=[0;1.0;1 ];
[t,y]=ode45('duffing',[0:0.01:100],x0);
%频谱图
N=1024;           %采样点数
x=y(10000-1024:10000,1);  %生成信号（Y轴方向的位移）
Y=fft(x,N);
Y(1)=[];
Pyy = Y.* conj(Y) /N;  %或者改变为Pyy = Y.* conj(Y) ;  
df=100;          %频率分辨率Hz

f=(1:N/2)*df/N;    %频域序列   
figure
% plot(f,Pyy);
% figure
plot(f,Pyy(1:N/2));    %幅值谱改变为plot(f,Y(1:N/2));   
xlabel('Frequency(HZ)');
ylabel('Magnitude');
title('功率谱图');

sssssxxxxx921

我这都没有运行，不过这个程序应该没问题了  呵呵

秋月

非常感谢!
sssssxxxxx921谢谢你的指点!
还望多指点!