《通信电子线路》顾宝良 编著附录Matlab程序A.1信号采样预备知识A1.2要求频率分辨率为1KHz以下

 

 

%《通信电子线路》顾宝良 编著
% A.1 信号采样预备知识
% A.1.2 要求频率分辨率为1KHz以下
clear
clf;
clc;

usefftshift = 1;% 经过fftshift移位，画图得到图1.2a。
usefftshift = 0;% 未经过fftshift移位，画图得到图1.2b。



f0 = 2e3;     %信号频率
t0 = 1/f0;    %信号周期

fs = 1e5;     %采样频率
ts = 1/fs;    %采样周期

df=1e3;

t=[-t0/2:ts:t0/2];
N = length(t);

xn = cos(2*pi*f0*t);

subplot(2,1,1);
stem(t*1000,xn,'fill');
grid on;
title('采样时间序列x(n)');
xlabel('T(ms)');
ylabel('Amplitute');

[X,x,df1] = fftseq(xn,ts,df);
X1=X/fs;
f1 = [0:df1:df1*(length(x)-1)]-fs/2;

subplot(2,1,2);
if usefftshift ==1
    plot(f1/1000,abs(fftshift(X1))); %使用fftshift
    grid on;
    title('经补零增长后的FFT变换（经过fftshift排列移位）');
    xlabel('f(kHz)');
    ylabel('Amplitute');
else

    plot(f1/1000,abs(X1));
    grid on;
    title('经补零增长后的FFT变换（未经过fftshift排列移位）');
    xlabel('f(kHz)');
    ylabel('Amplitute');

end