matlab 快速傅立叶变换 FFT 实验报告

《matlab 快速傅立叶变换 FFT 实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、南京邮电大学 实验报告 实验名称__熟悉MATLAB环境_____________ 快速傅里叶变换（FFT）及其应用 IIR数字滤波器的设计 FIR数字滤波器的设计 课程名称数字信号处理A 班级学号B08020221________ 姓名___孟祥熙_____ 开课时间2009/2010学年，第一学期 实验一 一、观察采样引起的混叠。（来源：p57习题1-3）设模拟信号为EMBEDEquation.3，t的单位为毫秒(ms)。设采样频率为3kHz，确定与EMBEDEquation.3混叠的采样重建信号EMBEDEquation.3。 画出EMBEDEquation.3和EMBEDEqu

2、ation.3在EMBEDEquation.3范围内的连续波形。（因数字计算机无法真正画出连续波形，可用较密的离散点的连线来近似。） 分别用EMBEDEquation.3和EMBEDEquation.3在两信号波形上标记出3kHz采样点。两信号波形是否相同？采样后的两序列是否相同？ 程序与分析： %problem1 %============= %作用与汇编中的；以及C++编程中的//作用一样，起注释作用 Clear %清屏 %estimatex(t)andxa(t)withamuchhighersamplingfreq.'fs1' %因数字计算机无法真正画出连续波形，可用较密的离散点的连

3、线来近似。 %因此，为画出x(n)与xa(t)的波形，取fs1=50KHZ。 time_period=6;%unit:ms fs1=50;%unit:kHz T1=1/fs1;%unit:ms n1=0:fix(time_period/T1); %用两个冒号组成等增量语句，其格式为：t=[初值：增量：终值]。 %如：t=[0：0.02：0.08]，t=00.020.040.060.08 %当增量为1时，这个增量值可以略去，因而该语句只有一个冒号。 %如：k=1：6，k=123456 %fix向0舍入为整数，舍尾法 x=cos(5*pi*n1*T1)+4*sin(2*pi*n1*T1).*s

4、in(3*pi*n1*T1); xa=cos(pi*n1*T1); %obtainx(nT)andxa(nT)withgivensamplingfreq.'fs' fs=3; T=1/fs; n=0:fix(time_period/T); x_sample=cos(5*pi*n*T)+4*sin(2*pi*n*T).*sin(3*pi*n*T); xa_sample=cos(pi*n*T); figure,plot(n1*T1,x,'r',n1*T1,xa,'b',n*T,x_sample,'ro'),%figure指令可以打开一个新的图形绘制窗口，在下一条这样的指令出现之前，%所有的绘图

5、指令都会在这个窗口上绘制。%plot函数可以同时画多组数值图像，即以（X坐标数组1，Y坐标数组1，X %坐标数组2，Y坐标数组2，……）的形式书写。每组坐标书写后，可用一个‘’%表示前面的一组坐标用怎样的形式绘制：线型，点型和颜色。如‘r’就表示 %用绿色的线绘制，‘b’蓝色。‘ro’红色圆圈等等。 holdon,stem(n*T,xa_sample,'b:x') %hold保持当前图形窗的内容即：在画完一张图后用hold命令保持住，继续在%当前窗画图。 %stem离散序列绘图 %b:x表示蓝色虚线x号 legend('x(t)','xa(t)','x(nT)','xa(nT)'),xla

6、bel('t(ms)') %legend这个指令可以在图像的右上角绘出一个图例，表示每条曲线代表什么， %名称需编程者指定。 %xlabelx轴标注，ylabel：y轴标注，title：三维坐标标注。 上面程序结果如下图所示： 波形分析： x(t)与xa(t)两波形并不相同，但采样后，波形一致。 二、判别离散时间系统的时不变性。（来源：p105例3.2.2）设输入序列为EMBEDEquation.3，系统EMBEDEquation.3实现对EMBEDEquation.3的抽取。设EMBEDEquation.3。取延迟量D（例如D＝30）。记EMBEDEquation.3，画出

7、EMBEDEquation.3、EMBEDEquation.3的序列波形。 编程求出系统对EMBEDEquation.3的响应EMBEDEquation.3以及对EMBEDEquation.3的响应EMBEDEquation.3 画出EMBEDEquation.3、EMBEDEquation.3的波形。该系统是否为时不变的？ 程序与分析： %============= %problem2 %============= c