ch3example11prg2.m

清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB和SIMULINK代码

% ch3example11prg2.m
w_m=40;     % 截断频率
T=pi/w_m;   % 采样间隔
L=5;
t=0: T :L ; % 时域截断
x_t=exp(-t).*(t>0); %信号序列
N=length(x_t);      %序列长度(点数)
%-----------------------------------
X_k=fft(x_t);                % FFT计算
%-----------------------------------
w0=2*pi/(N*T);               % 离散频率间隔
kw=2*pi/(N*T) .*[0: N-1];    % 离散频率样点
X_kw=T.*X_k;                 % 乘以T得到连续傅利叶变换频谱的样值
%-----------------------------------
subplot(2,1,1); 
plot(kw, abs(X_kw),'.','MarkerSize', 14); % 作出数值计算的幅度谱点
axis([-40 90 -0.1 1.1]);xlabel('freq (rad/s)');
hold on;                                  % 保持前面作图曲线不被擦除
w=-40: 0.1 : 40;
X_w=1./(1+j*w);                           % 理论计算频谱表达式
plot(w, abs(X_w));                        % 作图对比
subplot(2,1,2); 
plot(kw, angle(X_kw),'.','MarkerSize', 14); % 作出数值计算的相位谱点
hold on;
plot(w,  angle(F_w));                       % 频域相位谱
axis([-40 90 -pi pi]);xlabel('freq (rad/s)');