利用MATLAB实现连续信号的采样与重构仿真课程设计

TOC\o"1-5"\h\z\o"CurrentDocument"、摘要 12正文 2设计目的 2设计原理 2（1）、MTLAB简介 2连续时间信号 2、 采样定理 3、 信号重构 5\o"CurrentDocument"2.3 、 信号采样和恢复的程序 5计连续信号 6..设计连续信号的频谱 7设计采样信号 错误！未定义书签。设计采样信号的频谱图 9设计低通滤波器 10恢复原信号 12、总结和致谢 错误！未定义书签。15、参考文献15摘要本次课程设计使用MATLA殴现连续信号的采样和重构仿真，了解MATLA歆件，学习使用MATLA歆件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。初步掌握线性系统的设计方法，培养独立工作能力。加深理解采样和重构的概念，掌握利用MATLA分析系统频率响应的方法和掌握利用MATLA寞现连续信号采用和重构的方法。计算在临界采样、过采样、欠采样三种不同条件下重构信号的误差，并由此总结采样频率对信号重构误差的影响。要做到以下基本要求：.掌握利用MATLA盼析系统频率响应的方法，增加对仿真软件MATLAB勺感性认识，学会该软件的操作和使用方法。.掌握利用MATLA寞现连续信号采用和重构的方法，加深理解采样和重构的概念。.初步掌握线性系统的设计方法，培养独立工作能力。.学习MATLA计信号表示的基本方法及绘图函数的调用，实现对常用连续时间信号的可视化表示，加深对各种电信号的理解。.加深理解采样对信号的时域和频域特性的影响； 验证信号和系统的基本概念、基本理论，掌握信号和系统的分析方法。.加深对采样定理的理解和掌握， 以及对信号恢复的必要性； 掌握对连续信号在时域的采样和重构的方法。正文设计目的和要求对连续信号进行采样，在满足采样定理和不满足采用定理两种情况下对连续信号和采样信号进行FFT频谱分析。设计原理(1)、MTLABJ介MATLA(BMatrixLaboratory )是 1984年美国 MathWorks公司产品，Matlab的推出得到了各个领域专家学者的广泛关注，并越来越多的使用到我们的学习生活中来， 是目前通信工程上最广泛使用的软件之一。MatrixLaboratory意为“矩阵实验室”，最初的MATLAB是一个数学计算工具。但现在的MATLAB经远不仅仅是一个“矩阵实验室”,它已经成为一个集概念设计、算法开发、建模仿真，实时实现于一体的集成环境，它拥有许多衍生子集工具。MATLAB的 基 本 数 据 单 位 是 矩 阵 ， 它的指令表达式和数TOC\o"1-5"\h\z学 ,工程 中 常 用 的 形 式 十 分 相 似 , 故 用 MATLAB来解算 问 题要 比 用 C, FORTRAN 等 语 言 完 全 相 同 的 事情简 捷 得 多 . 在 新的 版 本中 也 加 入 了 对 C,FORTRAN, c++,JAVA 的 支 持 . 可 以直 接 调 用 ,用 户 也 可 以 将 自 己 编 写的实用 程 序 导 入 到MATLAB函数库 中 方 便 自 己 以 后 调 用 。(2) 连续时间信号连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值和之对应。严格来说，MATLA冲不能处理连续信号，而是用等时间问隔点的样值来近似表示连续信号。当取样时间问隔足够小时，这些离散的样值就能较好地近似连续信号。在一定条件下，一个连续时间信号完全可以用该信号在等时间问隔上的瞬时值来表示，并且可以用这些样本值把信号完全恢复过来。 这样,抽样定理为连续时间信号和离散时间信号的相互转换提供了理论依据。通过观察采样信号的频谱，发现它只是原信号频谱的线性重复搬移，只要给它乘以一个门函数，就可以在频域恢复原信号的频谱，在时域是否也能恢复原信号时，利用频域时域的对称关系，得到了信号。（3）采样定理模拟信号经过（A/D）变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率fs,重复出现一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。时域采样定理从采样信号工。恢复原信号/Q）必需满足两个条件：a、共）必须是带限信号，其频谱函数在|同＞为各处为零；（对信号的要求，即只有带限信号才能适用采样定理。）b、取样频率不能过低，必须0＞2%（或工＞24）。（对取样频率的要求，即取样频率要足够大，采得的样值要足够多，才能恢复原信号。）如图1所示，给出了信号采样原理图今相乘心⑥今相乘心⑥%)图1信号采样原理图由图1可见，fs(t)=f(t)^Ts(t),其中，冲激采样信号6Ts(t)的cd表达式为： 不s(t)=£6(t—nTs)n=一二二其傅立叶变换为牝£飞(8-nss),其中。,=至。设F(jco),Fs(jco)1 Ts分别为f(t),fs(t)的傅立叶变换，由傅立叶变换的频域卷积定理，可得1 .二 1.二Fs(j)=一F(j)*'s,('—n's)=—F[j('—n's)]2二 n-.:： Tsn=若设f(t)是带限信号，带宽为8m,f(t)经过采样后的频谱Fs(j8)就是将F(j。)在频率轴上搬移至0,土缸,±0s,…，故ns,…处(幅度为原频谱的1/Ts倍)。因此，当赳之2Gm时，频谱不发生混叠；而当缸＜23时，频谱发生混叠。一个理想采样器可以看成是一个载波为理想单位脉冲序列 每1(t)的幅值调制器，即理想采样器的输出信号e*(t),是连续输入信号e(t)调制在载波6T(t)上的结果，如图2所示。图2信号的采样用数学表达式描述上述调制过程，则有qQ理想单位脉冲序列6T(t)可以表示为 6T(t)=(t—nT)n=0其中6(t—nT)是出现在时刻t=nT,强度为1的单位脉冲。由于的数值仅在采样瞬时才有意义，同时，假设e(t)=0 -t<0QO所以e(t)又可表小为 e(t) e(nT)6(t-nT)n=0(4)信号重构设信号f(t)被采样后形成的采样信号为 fs(t)，信号的重构是指由fs(t)经过内插处理后，恢复出原来信号f(t)的过程，又称为信号恢复。若设f(t)是带限信号，带宽为«m,经采样后的频谱为Fs(j«)o设采样频率牝22露，则由式(9)知Fs(j0)是以缸为周期的谱线。现选取一Ts ®<«c ©个频率特性H(jo)=」 (其中截止频率0c满足3£露£二)06A牝 2的理想低通滤波器和Fs(j^)相乘，得到的频谱即为原信号的频谱F(jo)02.3信号采样和恢复的程序现在以正弦函数为例，进行MATLABf真实验。设计连续信号。先制作一个程序，使之产生一个正弦连续信号。所用程序如下所示： f1=50;t=(1:50)/2000; % 时间轴步距x=sin(2*pi*t*f1);figure(1);plot(x);%绘制x(t)的图形图片号加底框xlabel('t');ylabel('x(t)');title('连续时间信号波形 '); %图片命名grid;形如产生

形如1(2)设计连续信号的频谱。1设计一频谱程序n=0:511; %N=512; %Xk=abs(fft(x,N));figure(2); %plot(n,Xk);使其产生频谱波形图。程序如下长度设采样点的使其产生频谱波形图。程序如下长度设采样点的N直频域波形胸用axis函数来调整图轴的范围xlabel('时域频谱波形图’)；ylabel('|Xk|');波形如下:波形如下:(3)设计采样信号。设计一采样程序，使之输出采样波形。程序如下：X=fft(x,512);w=(0:255)/256*500;T=4*t;x=sin(2*pi*T*f1);figure(3);stem(x); % 图形x(n)的绘制xlabel('n');ylabel('x(n)');title(' 采样信号波形图'); % 图形命名

grid；波形如下:(4)设计采样信号的频谱图。设计出该采样信号的频谱程序，程序如下：figure(4);plot(w,abs([X(1:256)]));% 频谱图的绘制xlabel('Hz');ylabel('频率响应幅度');title(' 采样频谱波形图'); % 命名grid;波形如下:

30256-O帽MW30256-O帽MW售网景柔样频谱波形图(5)设计低通滤波器。设计一低通滤波器，使之具有滤波作用。程序如下：[B,A]=butter(8,350/500);% 巴特沃斯低通滤波器的设计[H,w]=freqz(B,A,512,2000);figure(5);subplot(2,1,1);plot(w*2000/(2*pi),abs(H));%低通频谱图的绘制xlabel('Hz');ylabel('频率响应幅度xlabel('Hz');ylabel('频率响应幅度’);title（' 低通滤波器波形图’）； ％命名grid；当采样频率f=350Hz时，波形图如下— 低通滤波器波形图1-5（ ■ । p । ■ ।当采样频率f=350Hz时，波形图如下低通滤波器波形图15低通滤波器波形图15(6)恢复原信号。设计程序，对采样信号频谱进行滤波，弁输出该信号所恢复频谱信号和连续信号，程序如下：y=filter(B,A,x);figure(6);subplot(2,1,1);plot(y);xlabel('t');ylabel('x(t)');title('连续信号波形')；grid；Y=fft(y,512);w=(0:255)/256*500;subplot(2,1,2);

plot(w,abs([Y(1:256)]));额谱图的绘制xlabel('Hz');ylabel('频率响应幅度')；title('恢复后的频谱波形图')；grid；波形如下：㈣陛回誓招场波形如下：㈣陛回誓招场为采样f=150Hz时的图2030连虢信号液形专□16 10 15 20 25 30 35 40 46 50回复后的频谱诚悠图al2010□口L41III-512030连虢信号液形专□16 10 15 20 25 30 35 40 46 50回复后的频谱诚悠图al2010□口L41III-51形:总结和致谢经过此次MATLAB程设计我学到了很多知识和学习方法。仅凭我在信号和系统实验课上所学的那点知识显然是不够的。所以为了做好这次的课程设计，我查阅了很多资料，并上网搜索了许多和此有关的知识，这个过程中我也学会了好多，通过浏览MATLAB坛，看到了好多有关学习MATLAB方法。在这次设计中，同样也学到了对信号的采样定理的使用，以及信号的重构，并通过观察MATLA所生成的频谱图，进一步了解了有关信号的采样和重构，对信号的采样程度进行比较其误差，了解不同采样程度的重构信号和原信号所产生的差异。同时，感谢老师的热心帮助下，使我的能够顺利的完成课程设计。同时，感谢各位同学的帮助，使我更好的完成课程设计。参考资料．董长虹 .Matlab信号处理和使用 [M].北京：国防工业出版社，2005.01。.甘俊英.基于MATLAB信号和系统实验指导[M].北京：清华大学出版社， 2007.8。.吴大正.信号和线性系统分析 [M