matlab_快速傅立葉变换_FFT_实验报告.docVIP

南京邮电大学 实 验 报 告 实验名称__熟悉MATLAB环境_____________ 快速傅里叶变换（FFT）及其应用 IIR数字滤波器的设计 FIR数字滤波器的设计 课程名称 数字信号处理A 班级学号 _______ 姓 名___孟祥熙_____ 开课时间 2009 /2010 学年， 第 一学期 实验一 一、观察采样引起的混叠。（来源：p57 习题1-3） 设模拟信号为EMBED Equation.3，t的单位为毫秒(ms)。 设采样频率为3kHz，确定与EMBED Equation.3混叠的采样重建信号EMBED Equation.3。 画出EMBED Equation.3和EMBED Equation.3在EMBED Equation.3范围内的连续波形。（因数字计算机无法真正画出连续波形，可用较密的离散点的连线来近似。） 分别用EMBED Equation.3和EMBED Equation.3在两信号波形上标记出3kHz采样点。 两信号波形是否相同？采样后的两序列是否相同？ 程序与分析： % problem 1 % ============= %作用与汇编中的；以及C++编程中的//作用一样，起注释作用 Clear %清屏 % estimate x(t) and xa(t) with a much higher sampling freq. fs1 %因数字计算机无法真正画出连续波形，可用较密的离散点的连线来近似。 %因此，为画出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=0 0.02 0.04 0.06 0.08 %当增量为1时，这个增量值可以略去，因而该语句只有一个冒号。 %如：k=1：6，k=1 2 3 4 5 6 %fix 向0舍入为整数，舍尾法 x=cos(5*pi*n1*T1)+4*sin(2*pi*n1*T1).*sin(3*pi*n1*T1); xa=cos(pi*n1*T1); % obtain x(nT) and xa(nT) with given sampling freq. 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指令可以打开一个新的图形绘制窗口，在下一条这样的指令出现之前， %所有的绘图指令都会在这个窗口上绘制。 %plot函数可以同时画多组数值图像，即以（X坐标数组1，Y坐标数组1，X %坐标数组2，Y坐标数组2，……）的形式书写。每组坐标书写后，可用一个‘’ %表示前面的一组坐标用怎样的形式绘制：线型，点型和颜色。如‘r’就表示 %用绿色的线绘制，‘b’蓝色。‘ro’红色圆圈等等。 hold on, stem(n*T,xa_sample,b:x) %hold保持当前图形窗的内容即：在画完一张图后用hold命令保持住，继续在%当前窗画图。 %stem离散序列绘图 %b:x表示蓝色虚线x号 legend(x(t),xa(t),x(nT),xa(nT)),xlabel(t(ms)) %legend这个指令可以在图像的右上角绘出一个图例，表示每条曲线代表什么， %名称需编程者指定。 %xlabel x轴标注，ylabel：y轴标注，title：三维坐标标注。 上面程序结果如下图所示： 波形分析： x(t)与xa(t)两波形并不相同，但采样后，波形一致。 二、判别离散时间系统的时不变性。（来源：p105 例3.2.2） 设输入序列为 EMBED Equation.3 ，系统 EMBED Equation.3 实现对 EMBED Equation.3 的抽取。 设 EMBED Equation.3 。取延迟量D（例如D＝30）。记 EMBED Equation.3 ，画出 EMBED Equation.3 、 EMBED Equation.3 的序列波形。