[工学]DSP实验报告.docVIP

实 验 报 告 实验名称 实验一 采样、系统性质及滤波 实验二_系统频率响应和样本处理算法实现 实验三 加窗和离散傅氏变换 实验四 数字滤波器设计 课程名称 数字信号处理A（双语） 班级学号________ _______ 姓 名______ ____________ 开课时间 实验一 采样、系统性质及滤波 实验目的和任务 （1）熟悉MATLAB的主要操作命令。 （2）学会简单的矩阵输入和数据读写。 （3）掌握简单的绘图命令。 （4）用MATLAB编程并学会创建函数。 （5）观察采样引起的混叠。 （6）判别离散时间系统的时不变性。 （7）卷积计算 实验内容与结果 A、观察采样引起的混叠 设模拟信号为，t的单位为毫秒(ms)。 设采样频率为3kHz，确定与混叠的采样重建信号。 画出和在范围内的连续波形。（因数字计算机无法真正画出连续波形，可用较密的离散点的连线来近似。） 分别用和在两信号波形上标记出3kHz采样点。两信号波形是否相同？采样后的两序列是否相同？ ※ 实验程序： clear % estimate x(t) and xa(t) with a much higher sampling freq. fs1 time_period=6; % unit: ms fs1=50; % unit: kHz T1=1/fs1; % unit: ms n1=0:fix(time_period/T1); 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), hold on, stem(n*T,xa_sample,b:x) legend(x(t),xa(t),x(nT),xa(nT)),xlabel(t(ms)) ※ 运行结果： B、判别离散时间系统的时不变性。 设输入序列为，系统实现对的抽取。 设。取延迟量D（例如D＝30）。记，画出、的序列波形。 编程求出系统对的响应以及对的响应 画出、的波形。 ※ 实验程序： clear % plot x(n) and x(n-D) D=30; N=500; n=1:N; x=sin(2*pi/100*n); for n=1:N+D, if (n-D)=0, xD(n)=0; else xD(n)=x(n-D); end end figure,subplot(2,1,1), plot(1:N,x,r:,1:length(xD),xD,b), legend(x(n),xD(n)),,xlabel(n) % plot y(n) and yD(n) for n=1:fix(N/2) y(n)=x(2*n); end for n=1:length(y)+D, if (n-D)=0, y_delay(n)=0; else y_delay(n)=y(n-D); end end for n=1:fix(length(xD)/2) yD(n)=xD(2*n); end subplot(2,1,2), plot(1:length(y),y,r:,1:length(y_delay),y_delay,r.:,1:length(yD),yD,b.), legend(y(n),y(n-D),yD(n)),,xlabel(n) axis([0 530 -1 1]) ※ 运行结果： C、利用卷积计算出输入信号通过FIR滤波器的输出，并观察输出信号的input-on暂态、input-off暂态和稳态部分。（计算卷积可用conv命令） 考虑下面两个滤波器，第一个的单位脉冲响应为，另一个的单位脉冲响应为；输入为周期方波，在一个周期内。 分别画出两个滤波器的输出的波形，并与书上p144例4.1.8的两幅图比较是否一致。 计算出图中稳态部分的值。 ※ 实验程序： clear h1=0.25*0.75.^(0