利用MATLAB实现连续信号的采样与重构仿真课程设计 2
- 格式：doc
- 大小：526.00 KB
- 文档页数：15
华北水利水电大学之答禄夫天创作
创作时间：贰零贰壹年柒月贰叁拾日
课程设计
课程名称：连续信号的采样与重构
专业班级：通信工程目 录
1、摘要1
2、正文2
2.1、设计目的2
、设计原理
(1)、连续时间信号2
(2)、采样定理3
(3)、信号重构5
、信号采样与恢复的程序5
（1）设计连续信号6
(2)设计连续信号的频谱7
（3）设计采样信号8
（4）设计采样信号的频谱图9
（5）设计低通滤波器10
（6）恢复原信号12
3、总结与致谢13
4、参考文献14
本次课程设计应用MATLAB实现连续信号的采样与重构仿真，了解MATLAB软件，学习应用MATLAB软件的仿真技术。它主创作时间：贰零贰壹年柒月贰叁拾日
创作时间：贰零贰壹年柒月贰叁拾日 要偏重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。初步掌握线性系统的设计方法，培养独立工作能力。
加深理解采样与重构的概念，掌握利用MATLAB分析系统频率响应的方法和掌握利用MATLAB实现连续信号采取与重构的方法。计算在临界采样、过采样、欠采样三种分歧条件下重构信号的误差，并由此总结采样频率对信号重构误差的影响。
要做到以下基本要求：
1. 掌握利用MATLAB分析系统频率响应的方法，增加对仿真软件MATLAB的感性认识，学会该软件的操纵和使用方法。
2. 掌握利用MATLAB实现连续信号采取与重构的方法，加深理解采样与重构的概念。
3 . 初步掌握线性系统的设计方法，培养独立工作能力。
4. 学习MATLAB中信号暗示的基本方法及绘图函数的调用，实现对经常使用连续时间信号的可视化暗示，加深对各种电信号的理解。
5. 加深理解采样对信号的时域和频域特性的影响；验证信号与系统的基本概念、基本理论，掌握信号与系统的分析方法。
6. 加深对采样定理的理解和掌握，以及对信号恢复的需要性；掌握对连续信号在时域的采样与重构的方法。
2.1 设计目的与要求 创作时间：贰零贰壹年柒月贰叁拾日
创作时间：贰零贰壹年柒月贰叁拾日 对连续信号进行采样，在满足采样定理和不满足采取定理两种情况下对连续信号和采样信号进行FFT频谱分析。
2.2 设计原理
（1） 连续时间信号
连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点以外，信号都有确定的值与之对应。严格来说，MATLAB其实不克不及处理连续信号，而是用等时间间隔点的样值来近似暗示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。
在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值来暗示，而且可以用这些样本值把信号完全恢复过来。这样，抽样定理为连续时间信号与离散时间信号的相互转换提供了理论依据。通过观察采样信号的频谱，发现它只是原信号频谱的线性重复搬移，只要给它乘以一个门函数，就可以在频域恢复原信号的频谱，在时域是否也能恢复原信号时，利用频域时域的对称关系，得到了信号。
(2)采样定理
模拟信号经过 (A/D) 变换转换为数字信号的过程称为采样，信号采样后其频谱发生了周期延拓，每隔一个采样频率
fs，重复出现一次。为包管采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。 时域采样定理从采样信号恢复原信号必须满足两个条件： 创作时间：贰零贰壹年柒月贰叁拾日
创作时间：贰零贰壹年柒月贰叁拾日 a、必须是带限信号，其频谱函数在＞各处为零；（对信号的要求，即只有带限信号才干适用采样定理。）
b、 取样频率不克不及过低，必须＞2（或＞2）。（对取样频率的要求，即取样频率要足够大，采得的样值要足够多，才干恢复原信号。）
如图1所示，给出了信号采样原理图
图1 信号采样原理图
由图1可见，)()()(ttftfsTs，其中，冲激采样信号)(tsT的表达式为： nsTnTtts)()(
其傅立叶变换为nssn)(，其中ssT2。设)(jF，)(jFs分别为)(tf，)(tfs的傅立叶变换，由傅立叶变换的频域卷积定理，可得
nssnsssnjFTnjFjF)]([1)(*)(21)(（1）
若设)(tf是带限信号，带宽为m，)(tf经过采样后的频谱创作时间：贰零贰壹年柒月贰叁拾日
创作时间：贰零贰壹年柒月贰叁拾日 )(jFs就是将)(jF在频率轴上搬移至,,,,,02nsss处（幅度为原频谱的sT1倍）。因此，当ms2时，频谱不发生混叠；而当ms2时，频谱发生混叠。
一个理想采样器可以看成是一个载波为理想单位脉冲序列)(tT的幅值调制器，即理想采样器的输出信号)(*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。设采样频率ms2，则由式（1）知)(jFs是以s为周期的谱线。现选取一个频率特性ccsTjH0)(（其中截止频率c满足2scm）的理想低通滤波器与)(jFs相乘，得到的频谱即为原信号的频谱)(jF。
2.3 信号采样与恢复的程序
现在以正弦函数为例，进行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;
发生的图形如下：创作时间：贰零贰壹年柒月贰叁拾日
创作时间：贰零贰壹年柒月贰叁拾日
（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|');
波形如下：
（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; 波形如下：
（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时，波形图如下：
（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;
波形如下：
