利用MATLAB绘制信号时域波形和信号的频谱

[x,Fs] = wavread('test1');%读取音频数据
x = x(:,1);
x = x';
N = length(x);%求取抽样点数
t = (0:N-1)/Fs;%显示实际时间
y = fft(x);%对信号进行傅里叶变换
f = Fs/N*(0:round(N/2)-1);%显示实际频点的一半
subplot(211);
plot(t,x,'g');%绘制时域波形
axis([0 max(t) -1 1]);
xlabel('Time / (s)');ylabel('Amplitude');
title('信号的波形');
grid;
subplot(212);
plot(f,abs(y(1:round(N/2))));
xlabel('Frequency / (s)');ylabel('Amplitude');
title('信号的频谱');
grid;
%% 以实际频点显示或者显示频率的一半
f = Fs/N*(0:N-1);
f = Fs/N*(0:round(N/2)-1);
%%以实际频点显示，显示的是Fs的一半
[h,w] = freqz(a,b,512,Fs);
%% w介于0--pi之间，除以pi，归一化到0--1，乘以奈奎斯特频率，实现显示频率的一半
[h,w] = freqz(a,b,512);
f = Fs*w/(2*pi);
%% 显示实际的时间，抽样点数除以抽样频率
t = (0:N-1)/Fs;
%% 信号频率，采样频率，关系

T/deltat = Fs/fs = T * Fs = N %抽样点数等于信号周期T乘以抽样频率，即每秒抽样多少次，乘以时间