appl6_3.m

MATLAB7.x数字信号处理 光盘内容

%Appl6_3
load iir_signaldata.txt; %调入数据文件
Xt=iir_signaldata; %得到原始信号序列
Fs=50; %采样率为50Hz
subplot(2,1,1); %时域坐标方框图
N=length(Xt); n=0:N-1; t=n/Fs;
plot(t,Xt); %画出时域中的原始波形图
xlabel('时间/s'),title('时间域'); %时域图的X轴坐标标识
ylabel('振幅'); %时域图的Y轴坐标标识
Xf=fft(Xt); %对原始波形数据进行快速Fourier变换
subplot(2,1,2); %频域坐标方框图
mag=abs(Xf); %得到信号序列Fourier变换后的振幅
semilogy([0:N-1]/(N/Fs),mag*2/N); %画出频域中的FFT波形图(y轴为对数坐标)
xlabel('频率/Hz'),title('频率域'); %频域图的X轴坐标标识
ylabel('振幅'); %频域图的Y轴坐标标识
xlim([0 Fs/2]); %频率轴只画出Nyquist频率之前的频率成分