实验一 数字信号处理 实验报告.pdf
古之立大事者,不惟有超世之才,亦必有坚忍不拔之志。——苏轼
07电信35号林烟泉
1.已知系统的差分方程如下式:
y1(n)=0.9y1(n-1)+x(n)
程序编写如下:
(1)输入信号x(n)=R10(n),初始条件y1(-1)=1,试用递推法求解输出y1(n);
a=0.9;ys=1;%设差分方程系数a=0.9,初始状态:y(-1)=1
xn=ones(1,10);%矩型序列R10(n)=u(n)-u(n-10),定义其宽度为0~9
n=1:35;%设差分方程系数a=0.9,初始状态:y(-1)=1
xn=sign(sign(10-n)+1);
B=1;A=[1,-a];%差分方程系数
xi=filtic(B,A,ys);%由初始条件计算等效初始条件输入序列xi
yn=filter(B,A,xn,xi);%调用filter解差分方程,求系统输出y(n)
n=0:length(yn)-1;
subplot(2,1,1);stem(n,yn,linewidth,2);axis([-5,15,0,8]);gridon
title(图(a)y1(n)=0.9y1(n-1)+x(n)初始条件y1(-1)=1);
xlabel(n);ylabel(y(n))
(2)输入信号x(n)=R10(n),初始条件y1(-1)=0,试用递推法求解输出y1(n)。
a=0.9;ys=0;%设差分方程系数a=0.9,初始状态:y(-1)=1
xn=ones(1,10);%矩型序列R10(n)=u(n)-u(n-10)
B=1;A=[1,-a];%差分方程系数
xi=filtic(B,A,ys);%由初始条件计算等效初始条件输入序列xi
yn=filter(B,A,xn,xi);%调用filter解差分方程,求系统输出y(n)
n=0:length(yn)-1;
subplot(2,1,2);stem(n,yn,linewidth,2);axis([-5,15,0,8]);gridon
title(图(b)y1(n)=0.9y1(n-1)+x(n)初始条件y1(-1)=0);
xlabel(n);ylabel(y(n))图形输出如下:
07电信35号林烟泉1
古之立大事者,不惟有超世之才,亦必有坚忍不拔之志。——苏轼
07电信35号林烟泉
图(a)y1(n)=0.9y1(n-1)+x(n)初始条件y1(-1)=1
8
6
)
n4
(
y
2
0