基于matlab的信号抽样系统设计与仿真——汪智宽

第 1 页

课程设计报告

题 目： 基于matlab的信号采样系统的设计

与仿真

课 程： 电子课程设计

学生姓名： 汪智宽

学生学号： 1214030230

专 业： 通信工程

班 级： 2班

指导教师： 王千春

电子工程学院制

2015年5月

第 2 页 基于matlab信号抽样系统的设计与仿真

学生：汪智宽

指导教师：王千春

电子工程学院: 通信工程

1 课程设计的任务与要求

1.1课程设计的任务

（1） 掌握利用MATLAB分析系统频率响应的方法，增加对仿真软件MATLAB的感性认识，学会该软件的操作和使用方法。

（2）掌握利用MATLAB实现连续信号采用与重构的方法，加深理解采样与重构的概念。

（3） 初步掌握线性系统的设计方法，培养独立工作能力。

（4） 学习MATLAB中信号表示的基本方法及绘图函数的调用，实现对常用连续时间信号的可视化表示，加深对各种电信号的理解。

（5） 加深理解采样对信号的时域和频域特性的影响；验证信号与系统的基本概念、基本理论，掌握信号与系统的分析方法。

（6） 加深对采样定理的理解和掌握，以及对信号恢复的必要性；掌握对连续信号在时域的采样与重构的方法。

1.2课程设计的要求

对连续信号进行采样，在满足采样定理和不满足采用定理两种情况下对连续信号和采样信号进行FFT频谱分析。

2 仿真系统的设计原理

(1)连续信号

连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。

在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值来表示，并且可以用这些样本值把信号完全恢复过来。这样，抽样定理为连续时间信号与离散时间信号的相互转换提供了理论依据。通过观察采样信号的频谱，发现它只 第 3 页 是原信号频谱的线性重复搬移，只要给它乘以一个门函数，就可以在频域恢复原信号的频谱，在时域是否也能恢复原信号时，利用频域时域的对称关系，得到了信号。

(2)采样定理

模拟信号经过 (A/D) 变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率 fs，重复出现一次。为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。

时域采样定理从采样信号恢复原信号必需满足两个条件：

a、必须是带限信号，其频谱函数在 ＞ 各处为零；（对信号的要求，即只有带限信号才能适用采样定理。）

b、 取样频率不能过低，必须 ＞2 （或 ＞2）。（对取样频率的要求，即取样频率要足够大，采得的样值要足够多，才能恢复原信号。）

如图1所示，给出了信号采样原理图

图1 信号采样原理图3

见，，其中，冲激采样信号)(tsT的表达为： 由图1可nsTnTtts)()(其傅立叶变换为nssn)(，其中ssT2。设)(jF，)(jFs分别为)(tf，)(tfs的傅立叶变换，由傅立叶变换的频域卷积定理，可得

nssnsssnjFTnjFjF)]([1)(*)(21)(

若设)(tf是带限信号，带宽为m， )(tf经过采样后的频谱)(jFs就是将)(jF在频率轴上搬移至,,,,,02nsss处（幅度为原频谱的sT1倍）。因此，当ms2时，频谱不发生混叠；而当ms2时，频谱发生混叠。

一个理想采样器可以看成是一个载波为理想单位脉冲序列)(tT的幅值调制器，)()()(ttftfsTs 第 4 页 即理想采样器的输出信号)(*te，是连续输入信号)(te调制在载波)(tT上的结果，如图2所示。

图2 信号的采样

用数学表达式描述上述调制过程，则有 )()()(*tteteT

理想单位脉冲序列)(tT可以表示为 0)()(nTnTtt

其中)(nTt是出现在时刻nTt，强度为1的单位脉冲。由于

的数值仅在采样瞬时才有意义，同时，假设

00)(tte

所以)(*te又可表示为 *0()()()netenTtnT

（3）信号重构

设信号)(tf被采样后形成的采样信号为)(tfs，信号的重构是指由)(tfs经过内插处理后，恢复出原来信号)(tf的过程,又称为信号恢复。

若设)(tf是带限信号，带宽为m，经采样后的频谱为)(jFs。设采样频率ms2，则由式（9）知)(jFs是以s为周期的谱线。现选取一个频率特性ccsTjH0)(（其中截止频率c满足2scm）的理想低通滤波器与)(jFs相乘，得到的频谱即为原信号的频谱)(jF。

3利用matlab进行仿真

3.1matlab软件介绍

MATLAB（Matrix Laboratory）是1984年美国Math Works公司产品，Matlab的)(te 第 5 页 推出得到了各个领域专家学者的广泛关注，并越来越多的应用到我们的学习生活中来，是目前通信工程上最广泛应用的软件之一。Matrix Laboratory意为“矩阵实验室”，最初的MATLAB只是一个数学计算工具。但现在的MATLAB已经远不仅仅是一个“矩阵实验室”，它已经成为一个集概念设计、算法开发、建模仿真，实时实现于一体的集成环境，它拥有许多衍生子集工具。

MATLAB的基本数据单位是矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完全相同的事情简捷得多.在新的版本中也加入了对C,FORTRAN,c++ ,JAVA的支持.可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。

3.2仿真程序及仿真结果

现在以正弦函数为例，进行MATLAB仿真实验。

（1）设计连续信号。

先制作一个程序，使之产生一个正弦连续信号。所用程序如下所示：

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;

产生的图形如下: 第 6 页

图3连续时间信号波形

（2）设计连续信号的频谱。

设计一频谱程序，使其产生频谱波形图。程序如下：

n=0:511; %长度

N=512; %设采样点的N值

Xk=abs(fft(x,N));

figure(2); %频域波形

plot(n,Xk);

axis([0 N 1.1*min(Xk) 1.1*max(Xk)]);

%可用axis函数来调整图轴的范围

xlabel('时域频谱波形图');

ylabel('|Xk|');

波形如下： 第 7 页

图4时域频谱波形

（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;

波形如下： 第 8 页

图5采样信号波形

（4）设计采样信号的频谱图。

设计出该采样信号的频谱程序，程序如下：

figure(4);plot(w,abs([X(1:256)])); %频谱图的绘制

xlabel('Hz');ylabel('频率响应幅度');

title('采样频谱波形图'); %命名

grid;

波形如下： 第 9 页

图6采样频谱波形

（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('频率响应幅度');

title('低通滤波器波形图'); %命名

grid;

当采样频率f=350Hz时，波形图如下： 第 10 页

图7低通滤波器波形

当采样频率f=500Hz时，波形图如下：

图8低通滤波器波形

4.总结与致谢

经过此次MATLAB课程设计我学到了很多知识和学习方法。仅凭我在信号与系统实验课上所学的那点知识显然是不够的。所以为了做好这次的课程设计，我查阅了很多资料，并上网搜索了许多与此有关的知识，这个过程中我也学会了好多，通过浏览MATLAB论坛，看到了好多有关学习MATLAB的方法。

在这次设计中，同样也学到了对信号的采样定理的应用，以及信号的重构，并通过观察MATLAB所生成的频谱图，进一步了解了有关信号的采样与重构，对信号的采样程度进行比较其误差，了解不同采样程度的重构信号和原信号所产生的差异。

同时，感谢老师的热心帮助下，使我的能够顺利的完成课程设计。同时，感谢各位同学的帮助，使我更好的完成课程设计。