数字信号实验1

1.课本第6页程序（图1）

程序：n=-5:5;

x=sin(pi*n/5);

subplot(3,2,1);stem(n,x,'.');line([-5,6],[0

,0])

axis([-5,6,-1.2,1.2]);xlabel('n');ylabel('y(

n)')

2.课本第14页程序（图2）

程序：xn=[1,1,1,1];hn=[1,1,1,1];

yn=conv(xn,hn);n=0:length(yn)-1; subplot(3,2,1);stem(n,yn,'.')

axis([0,7,0,4.2]);xlabel('n');ylabel('y(n)')

3.课本第21页程序（图3）

程序：a=0.8;ys=1;

xn=[1,zeros(1,30)]; B=1;A=[1,-a]; xi=filtic(B,A,ys); yn=filter(B,A,xn,xi); n=0:length(yn)-1;

subplot(3,2,1);stem(n,yn,'.') xlabel('n');ylabel('y(n)')

4.实验一第（2）题

①输入序列为x（n）=R8（n）（图4.1） 程序：xn=[1,ones(1,7)];

B=[0.05,0.05];A=[1,-0.9]; yn=filter(B,A,xn,xi); n=0:length(yn)-1;

subplot(3,2,1);stem(n,yn,'.') xlabel('n');ylabel('y(n)')

②输入序列为x（n）=u（n）（图4.2）

程序：xn=[1,ones(1,30)];

B=[0.05,0.05];A=[1,-0.9]; yn=filter(B,A,xn,xi);

n=0:length(yn)-1; subplot(3,2,1);stem(n,yn,'.') xlabel('n');ylabel('y(n)') axis([0,32,0,1.2])

图1

图2

图3

图4.1

图4.2

数字信号处理实验一 班级：06普本电信2班 学号：060102004 姓名：洪霞

③系统的单位脉冲响应（图4.3） 程序：xn=[1,zeros(1,30)];

B=[0.05,0.05];A=[1,-0.9]; yn=filter(B,A,xn,xi); n=0:length(yn)-1; subplot(3,2,1);stem(n,yn,'.') xlabel('n');ylabel('y(n)') axis([0,32,0,1])

5.实验一第(3)题

①系统的单位脉冲响应为h1(n) (图5.1) 程序：xn=[1,ones(1,7)];

hn=[1,ones(1,9)];

yn=conv(xn,hn);n=0:length(yn)-1;

subplot(3,2,1);stem(n,yn,'.') xlabel('n');ylabel('y(n)') axis([0,18,0,10])

②系统的单位脉冲响应为h2(n) (图5.2) 程序：xn=[1,ones(1,7)];

hn=[1,2.5,2.5,1];

yn=conv(xn,hn);n=0:length(yn)-1; subplot(3,2,1);stem(n,yn,'.') xlabel('n');ylabel('y(n)') axis([0,11,0,8])

6.实验一第(4)题

①输入信号为u(n) (图6.1) 程序：xn=[1,ones(1,30)];

ys=[0,0];b0=1/100.49;

B=[b0,-b0];A=[1,-1.8237,0.9801]; xi=filtic(B,A,ys);

yn=filter(B,A,xn,xi);n=0:length(yn)-1; subplot(3,2,1);stem(n,yn,'.') xlabel('n');ylabel('y(n)') axis([0,31,-0.04,0.04])

图4.3

图5.1

图5.2

图6.1

数字信号处理实验一 班级：06普本电信2班 学号：060102004 姓名：洪霞

②输入信号为x(n)=sin(0.014n)+sin(0.4n) 如图6.2 程序：n=0:40;

xn=sin(0.014*n)+sin(0.4*n); ys=[0,0];b0=1/100.49;

B=[b0,-b0];A=[1,-1.8237,0.9801]; xi=filtic(B,A,ys);

yn=filter(B,A,xn,xi);ny=0:length(yn)-1; subplot(3,2,1);stem(ny,yn,'.') xlabel('n');ylabel('y(n)') axis([0,42,-0.2,0.2])

图6.2