实验报告模板实验一
- 格式：docx
- 大小：53.38 KB
- 文档页数：5
百度文库-让每个人平等地提升自我
11
实验一离散时间信号的MATLAB实现
课程名称：数字信号处理 实验日期： ________________
班级： 数媒1403 姓名： 李欣 _______________ 学号： 2
一、 实验目的
掌握:(1)序列的生成；(2)离散序列的基本运算与波形变换的 MATLA实现;
(3)线性时不变系统求解的MATLA实现;(4)序列的傅里叶变换与信号重构
MATLA实现
二、 实验要求
按照实验材料内容，利用 MATLA编程实现相关任务。
三、 实验步骤
(以实验材料第四部分“序列的傅里叶变换与信号重构 MATLA实现”为例，
给出实验步骤描述、相关代码等。)
例111亦=严切「求其傅立叶超J期:
箜：严裕说，在MATLAB申不燮毛巧mbdi匚工具箱杲不能分析模拟值号的，佢是当采椿时 简闾陽充分小的时候，可产生平滑的團賂 当时間足鴨长，可显示出所有的模型，也就是可 以近似的分析.此例中.心⑵为応WOCOIE的需限信号，因比取
1 r = 5 * 10 << ------------- = 25 *10 ! 2 * 20IM
程序如下：
Dt=;
t=: Dt:;
xa=exp(-1000*abs(t));
%模拟信号
Wmax=2*pi*2000;
K=500;
k=0:1:K;
W=k*Wmax/K;
%连续时间傅立叶变换
Xa=xa*exp(-1i*t'*W)*Dt;
Xa=real(Xa);
W=[-fliplr(W),W(2:501)];
Xa=[fliplr(Xa),Xa(2:501)];
figure(1)
subplot(2,1,1); 百度文库-让每个人平等地提升自我
12
plot(t*1000,xa,'.');百度文库-让每个人平等地提升自我
22
xlabel('t in msec'); ylabel('xa(t)'); gtext('模拟信号'); subplot(2,1,2);
plot(W/(2*pi*1000),Xa*1000,'.'); xlabel('Freque nee in KHz');
ylabel('Xa(jw)*1000');
gtext('连续时间傅立叶变换');
程序运行结果:
例1T2 1-11申的耳⑴说匪耒舞扳幸对麺城特性的影响，分别取采样频率为
^=?OOQHz^ f=lOOOBz.養出 曲裟二
程序如下：
Dt=;
t=:Dt:; %模拟信号
xa=exp(-1000*abs(t));
Ts=;
n=-25:1:25; %离散时间信号
x=exp(-1000*abs( n*Ts));
K=500;
k=0:1:K;
w=pi*k/K; %离散时间傅立叶变换
X=x*exp(-1i* n'*w);
X=real(X);
w=[-fliplr(w),w(2:501)];
X=[fliplr(X),X(2:501)];
figure
subplot(2,1,1); plot(t*1000,xa,T);百度文库-让每个人平等地提升自我
33
ylabel('xa1(t)');
gtext('离散信号');
hold on ;
stem( n*Ts*1000,x);
hold off ;
subplot(2,1,2);
plot(w/pi,X,'.');
ylabel('X1(jw)');
gtext('离散时间傅立叶变换');
gtext('Ts=')
程序运行结果:
例1 13对例1小中产生的离散字列血⑹和紐来月sine函数进等内擂重构“ 程序如下：
Ts1=;
Fs1=1Ts1;
n1= -25:1:25;
n Ts仁n 1*Ts1; %离散时间信号
x1= exp(-1000*abs( nTs1));
Ts2=;
Fs2=1T s2; n2=-5:1:5;
n Ts2=n 2*Ts2;
x2=exp(-1000*abs( nTs2)); \ /
Dt=;
t=:Dt:; %模拟信号重构
xa1=x1*si nc(Fs1*(o nes(le ngth( nTs1),1)*t -n Ts1'*o nes(1,le ngth(t)))); xa2=x2*si nc(Fs2*(o nes(le ngth( nTs2),1)*t -n Ts2'*o nes(1,le ngth(t))));百度文库-让每个人平等地提升自我
44
subplot(2,1,1);
plot(t*1000,xa1,'.');
ylabel('xa1(t)');
title('从x1(n)重构模拟信号xa(t)');
hold on;
stem( n1*Ts1*1000,x1);
hold off;
subplot(2,1,2);plot(t*1000,xa2,'.');
ylabel('xa2(t)');
title('从x2(n)重构模拟信号xa(t)');
hold on;
stem( n2*Ts2*1000,x2);
hold off
四、实验体会
(1) 掌握了 MATLA的基本操作以及学会了序列的生成；
(2) 学会了用MATLA实现离散序列的基本运算与波形变换
(3) 学会了用MATLA实现线性时不变系统求解；
(4) 学会了用MATLA实现序列的傅里叶变换与信号重构 程序运行结果:
-5 4 -J -1 -I 0 I J J 4 5
U.x^nj l：r,
