内容发布更新时间 : 2024/12/24 2:15:15星期一 下面是文章的全部内容请认真阅读。

数字信号处理实验报告

姓名：蔡由卫 学号：0702100202 班级：通信1102

实验一：系统响应及系统稳定

1.给定一个低通滤波器的差分方程为y(n)=0.05x(n)+0.05x(n-1)+0.9y(n-1)，输入信号 x1(n)=R8(n), x2(n)=u(n)

（1）分别求出x1(n)=R8(n)和x2(n)=u(n)的系统响应，并画出其波形。 （2）求出系统的单位脉冲响应，画出其波形。 clc clear all close all ys=0

x1=[ones(1,7) zeros(1,24)]; x2=[ones(1,30)]; x3=[1,zeros(1,29)];

B=[0.05,0.05];A=[1,-0.9];%差分方程系数

xi=filtic(B,A,ys); %由初始条件计算等效初始条件的输入序列xi y1=filter(B,A,x1,xi);%调用filter解差分方程，求系统输出信号y1(n) y2=filter(B,A,x2,xi); %调用filter解差分方程，求系统输出信号y2(n) y3=filter(B,A,x3,xi); %调用filter解差分方程，求系统输出信号y3(n) n=0:29;

subplot(3,1,1);stem(y1,'.');xlabel('n');ylabel('y_1(n)');%画y1(n)的波形 subplot(3,1,2);stem(y2,'.');xlabel('n');ylabel('y_2(n)'); %画y2(n)的波形 subplot(3,1,3);stem(y3,'.');xlabel('n');ylabel('h(n)'); %画h(n)的波形

. . .

2.给定系统的单位脉冲响应为 h1 (n)=R10(n)

h2(n)=??n? +2.5?（n-1+2.5?(n-2)+?(n-3) ）用线性卷积法求x1（n）=R8（n）分别对系统h1（n）和h2（n）的输出响应并画出波形。

clc clear all close all h1n=[ones(1,9)]; h2n=[1 2.5 2.5 1]; xn=[ones(1,7)];

y1n=conv(xn,h1n); %y1n为xn与h1n的线性卷积 y2n=conv(xn,h2n); %y2n为xn与h2n的线性卷积 subplot(2,1,1);stem(0:length(y1n)-1,y1n,'.');xlabel('n');ylabel('y_1(n)');title('y_1(n)=x(n)*h_1(n)'); %画y1n的波形

subplot(2,1,2);stem(0:length(y2n)-1,y2n,'.');xlabel('n');ylabel('y_2(n)');title('y_2(n)=x(n)*h_2(n)'); %画y2n的波形

. . .

y1(n)=x(n)*h1(n)86y(n)142002468101214ny2(n)=x(n)*h2(n)86y(n)242001234n56789

3.给定一谐振器的差分方程为

y(n)=1.8237y(n-1)-0.9801y(n-2)+b0x(n)-b0x(n-2) 令b0=1/100.49，谐振器的谐振频率为0.4rad。

（1）用实验方法检查系统是否稳定。输入信号为u（n）时，画出系统输出波形。 （2）给定输入信号为 x（n）=sin（0.014n）+sin（0.4n）

求出系统的输出响应，并画出其波形。

clc close all clear all

un=ones(1,256)；%输入un n=0:255;

xn=sin(0.014*n)+sin(0.4*n); %输入xn A=[1 -1.8237 0.9801];%查分方程系数 B=[1/100.49 -1/100.49];%查分方程系数

y1n=filter(B,A,un); %调用filter解差分方程，求系统输出信号y1(n)

y2n=filter(B,A,xn); %调用filter解差分方程，求系统输出信号y2(n)

subplot(2,1,1);stem(0:length(y1n)-1,y1n,'.');xlabel('n');ylabel('y_1(n)');title('输入u

(n)输出波形');%画输出波形y1n

. . .