example1_20a.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 19 行

%例1-20程序，example1_20a
%模拟信号
Dt=0.00005;t=-0.005:Dt:0.005;xa=exp(-1000*abs(t));
%离散时间信号
Ts=0.0002;n=-25:1:25;
x=exp(-1000*abs(n*Ts));
%离散傅立叶变换
K=500;k=0:1:K;w=pi*k/K;
X=x*exp(-j*n'*w);X=real(X);
w=[-fliplr(w),w(2:K+1)];
X=[fliplr(X),X(2:K+1)];
subplot(1,1,1)
subplot(2,1,1);plot(t*1000,xa);
xlabel('t毫秒');ylabel('x1(n)');
title('离散信号');hold on
stem(n*Ts*1000,x);gtext('Ts=0.2毫秒');hold off
subplot(2,1,2);plot(w/pi,X);
xlabel('频率（单位：π）');ylabel('X1(w)');
title('离散时间傅立叶变换');