matlab求出过冲的最大幅值,频谱分析为什么最大幅值处频率为0？

这是我的程序

function [wa]=Wa0(x0,y0)

syms t tao

[m,~]=size(x0);  %求xx的维数

xx=0:0.001:13;  %定义绘图坐标轴

t=0:0.001:13;

plot(x0,y0,'*');

hold on;

fs=2000;N=20000;   %采样频率和数据点数

n=0:N-1;

t=n/fs;   %时间序列

%x=y0; %信号

W1=detrend(y0); % 消除时间序列中的线性趋势项

W1=detrend(W1,'constant'); % 消除时间序列中的均值

subplot(2,2,1),plot(x0,W1,'*');

yy=fftshift(fft(W1,N));    %对信号进行快速Fourier变换

mag=abs(yy);     %求得Fourier变换后的振幅

mag=mag/(N/2);   %换算成实际的幅度

mag(1)=mag(1)/2;

f=(-N/2:N/2-1)*fs/N;%进行对应的频率转换

subplot(2,2,2),plot(f,mag);   %绘出随频率变化的振幅

xlabel('频率/Hz');

ylabel('振幅');grid on;

subplot(2,2,3),plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅

xlabel('频率/Hz');

ylabel('振幅');

title('振幅-频率曲线图');grid on;

Pyy=[1:N/2];

for i=1:N/2

Pyy(i)=phase(yy(i)); %计算相位

Pyy(i)=Pyy(i)*180/pi; %换算为角度

end;

subplot(2,2,4),plot(f(1:N/2),Pyy(1:N/2));   %显示相位图

title('相位-频率曲线图');

xlabel('相位/度');

ylabel('振幅');grid on;

1.jpg

(169.4 KB, 下载次数: 6)

2015-10-14 11:13 上传

2015-10-14 11:12 上传

点击文件名下载附件

2.81 MB, 下载次数: 40