数字信号处理实验报告110.docxVIP

中南民族大学South-Central University for Nationalities实 验 报 告 2017-2018学年 第 1 学期 课程名称： 数字信号处理 学院： 生物医学工程学院 专业： 生物医学工程 年级： 2015级 姓名： 张慧 学号： 201521121032 实验一：系统响应及系统稳定性实验目的（1）掌握求系统响应的方法（2）掌握时域离散系统的时域特性（3）分析、观察及检查系统的稳定性实验代码close all; clear all A=[1,-0.9];B=[0.05,0.05];x1n=[1 1 1 1 1 1 1 1 zeros(1,50)];x2n=ones(1,128);hn=impz(B,A,58);%subplot (2,2,1);y=h(n);tstem(hn,y);subplot(2,2,1);y=h(n);tstem(hn,y);title((a)系统单位脉冲响应h(h))y1n=filter(B,A,x1n);subplot(2,2,2);y=y1(n);tstem(y1n,y);title((b)系统对R8（n）的响应yl(n))y2n=filter(B,A,x2n);subplot(2,2,4);y=y2(n);tstem(y2n,y);title((c)系统对u(n)的响应y2(n))x1n=[1 1 1 1 1 1 1 1 ];h1n=[ones(1,10) zeros(1,10)];h2n=[1 2.5 2.5 1 zeros(1,10)];y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);figure(2)subplot (2,2,1);y=h1(n);tstem(h1n,y);title((d)系统单位脉冲响应h1（n）)subplot(2,2,2);y=y21(n);tstem(y21n,y);title((e)h1(n)与R8(n)的卷积y21(n))subplot(2,2,3);y=h2(n);tstem(h2n,y);title((f)系统单位脉冲响应h2（n）)subplot(2,2,4);y=y22(n);tstem(y22n,y);title((g)h2(n)与R8(n)的卷积y22(n))un=ones(1,256);n=0:255;xsin=sin(0.014*n)+sin(0.4*n);A=[1,-1.8237,0.9801];B=[1/100.49,0,-1/100.49];y31n=filter(B,A,un);y32n=filter(B,A,xsin);figure(3)subplot(2,1,1);y=y31(n);tstem(y31n,y);title((h)谐振器对u(n)的响应y31(n))subplot(2,1,2);y=y32(n);tstem(y32n,y);title((i)谐振器对正弦信号的响应y32(n))实验结果 实验二：时域采样与频域采样一、实验目的时域采样理论与频域采样时数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何采样频率才能使采样后的信号不丢失信息；要求掌握频域采样会引起时域周期化的概念，以及频域采样定理及其对频域采样点数选择的指导作用。二、实验代码1、时域采样理论的验证程序Tp=64/1000;Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp(-alph*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M);yn=xa(nT);subplot(3,2,1);tstem(xnt,yn);box on;title((a)Fs=1000Hz);k=0:M-1;fk=k/Tp;subplot(3,2,2); plot(fk,abs(Xk));title((a)T*FT[xa(nT)],Fs=1000Hz);xlabel(f(Hz));ylabel(幅度);axis([0,Fs,0,1.2*max(abs(Xk))])%==============================%Fs=300Hz和Fs=200Hz的程序同上Tp=64/1000;Fs=300;T=1/Fs;M1=Tp*Fs;n=0:M1-1;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;x1nt=A*exp(-alph*n*T).*sin(omega*n*T);X1k=T*fft(x1nt,M1);y1n=xa(nT);subplot(3,2,3);tstem(x1nt,y1n)