预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共14页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

/NUMPAGES14 目录 摘要…………………………………………………………………………………2 正文 一、设计目的与要求………………………………………………………………3 二、设计原理………………………………………………………………………4 三、设计容和步骤………………………………………………………………5 1.用MATLAB产生连续信号y=sin(t)和其对应的频谱…………………6 2.对连续信号y=sin(t)进行抽样并产生其频谱…………………………7 3.通过低通滤波恢复原连续信号…………………………………………9 四、总结……………………………………………………………………………12 五、致……………………………………………………………………………13 六、参考文献………………………………………………………………………14 摘要 数字信号处理是一门理论与实践紧密结合的课程。做大量的习题和上机实验,有助于进一步理解和巩固理论知识,还有助于提高分析和解决实际问题的能力。过去用其他算法语言,实验程序复杂,在有限的实验课时所做的实验容少。MATLAB强大的运算和图形显示功能,可使数字信号处理上机实验效率大大提高。特别是它的频谱分析和滤波器分析与设计功能很强,使数字信号处理工作变得十分简单、直观。 本实验设计的题目是:信号的采样与恢复、采样定理的仿真。通过产生一个连续时间信号并生成其频谱,然后对该连续信号抽样,并对采样后的频谱进行分析,最后通过设计低通滤波器滤出抽样所得频谱中多个周期中的一个周期频谱,并显示恢复后的时域连续信号。实验中,原连续信号的频谱由于无法实现真正的连续,所以通过扩大采样点的数目来代替,理论上当采样点数无穷多的时候即可实现连续,基于此尽可能增加采样点数并以此来产生连续信号的频谱。信号采样过程中,通过采样点的不同控制采样频率实现大于或小于二倍最高连续信号的频率,从而可以很好的验证采样定理。信号恢复,滤波器的参数需要很好的设置,以实现将抽样后的信号进行滤波恢复原连续信号。 一、设计目的与要求 1.设计目的和要求 1.掌握利用MATLAB在数字信号处理中的基本应用,并会对结果用所学知识进行分析。 2.对连续信号进行采样,在满足采样定理和不满足采用定理两种情况下对连续信号和采样信号进行FFT频谱分析。 3.从采样信号中恢复原信号,对不同采样频率下的恢复信号进行比较分析。 4.基本要求:每组一台电脑,电脑安装MATLAB6.5版本以上软件。 二、设计原理 本实验主要涉与采样定理的相关容以与低通滤波器恢复原连续信号的相关知识。 1.采样定理: 设连续信号属带限信号,最高截止频率为,如果采样角频率,那么让采样性信号通过一个增益为T、截止频率为的理想低通滤波器,可以唯一地恢复出原连续信号。否则,会造成采样信号中的频谱混叠现象,不可能无失真地恢复原连续信号。 对连续信号进行等间隔采样形成采样信号,对其进行傅里叶变换可以发现采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。 对模拟信号进行采样可以看做一个模拟信号通过一个电子开关S,设电子开关每隔周期T和上一次,每次和上的时间为τ,在电子开关的输出端得到采样信号x^a(t)。用公式表示如下: (2.2.1) 图1对模拟信号进行采样 2.信号的恢复: 可用传输函数的理想低通滤波器不失真地将原模拟信号恢复出来,只是一种理想恢复。 因为 (2.2.2) 理想低通滤波器的输入输出和, =*=(2.2.3) 三、设计容和步骤 1.用MATLAB产生连续信号y=sin(t)和其对应的频谱 %................时域连续信号和频谱...................................... x1=0:pi/10:(8*pi); w=linspace(0,8*pi,length(x1)); figure subplot(211) plot(x1,sin(x1));%原时域连续信号y=sin(t) xlabel('t');ylabel('x(t)'); title('原时域连续信号y=sin(t)'); grid sin1=sin(x1); n=0:(length(x1)-1); subplot(212) plot(w,fft1(w,sin1,n));%其对应频域信号Y=FFT(sin(t)) xlabel('w');ylabel('x(w)'); title('其对应频域信号Y=FT(sin(t))'); grid 其中要用到子函数fft1,程序代码如下: functionresult=fft1(w,ha