rectangular_ctft.m

经典《信号与系统》教程的matlab例程,对深入理解信号与系统相关概念有很大帮助

%==========================================================================
clear;close all
p=0.01;
t=-5:p:5;
%x=u(t+1)-u(t-1);%input('input the expression of signal x(t)：');
x=sin(4*t+eps)./(4*t+eps);
W1=6*pi;                           %确定信号的频率范围
N=512;           
k=-N:N;                              %对信号的角频率进行采样的点数为2N+1点
W=k*W1/N;                            %计算频率的采样间隔
%==========================================================================
%    下段程序计算信号的傅里叶变换，包括幅度谱和相位谱计算
%==========================================================================
X=x*exp(-j*t'*W)*p;                  %进行傅里叶数值积分运算（矩阵运算）
X1=abs(X);                           %求模——幅度谱
phi=57.3*angle(X);                     %求相位谱
%==========================================================================
%                            下段程序绘制图形
%==========================================================================
subplot(221);plot(t,x,'r');grid on;
xlabel('t');
title('The signal x(t)');axis([-5,5,-0.3*max(x),1.1*max(x)])
subplot(222);plot(W,X,'r');grid on;axis([-W1,W1,-0.3*max(X1),1.1*max(X1)])
title('The Fourier transform of x(t) ');xlabel('w');

%===========================================================================