plotspec.m

来自「利用MATLAB实现正交滤波器组的设计。」· M 代码 · 共 15 行

function plotspec(x,Ts)

N=length(x);                               % 信号长度
t=Ts*(1:N);                                % 定义一个时间向量 
ssf=(-N/2:N/2-1)/(Ts*N);                   % 频率向量
fx=fft(x(1:N));                            % 做FFT变换
fxs=fftshift(fx);                           
subplot(2,1,1), plot(t,x)                  % 花时域波形
xlabel('seconds'); ylabel('amplitude')     
subplot(2,1,2), plot(ssf,abs(fxs))         %花频率波形
xlabel('frequency'); ylabel('magnitude')