数字信号处理第三版高西全实验.pdfVIP

数字信号处理高西全版课后上机实验

1.实验一：系统响应及系统稳定性

clc

clear

x1=[ones(1,8),zeros(1,20)];

subplot(321)

stem(x1,.)

title(u8(n));xlabel(n);ylabel(u(n));

x2=ones(1,30);

subplot(322)

stem(x2,.)

title(u(n));xlabel(n);ylabel(u(n));

B=[0.05,0.05];

A=[1,-0.9];

y1=filter(B,A,x1);

subplot(323)

stem(y1,.)

title(y(n));xlabel(n);ylabel(y(n));

y2=filter(B,A,x2);

subplot(324)

stem(y2,.)

title(y(n));xlabel(n);ylabel(y(n));

x3=[1,zeros(1,30)];

subplot(325)

stem(x3,.)

title(单位脉冲);xlabel(n);ylabel(u(n));

y3=filter(B,A,x3);

subplot(326)

stem(y3,.)

title(单位脉冲响应);xlabel(n);ylabel(y(n));

图1实验(1)波形

clc

clear

x1=[ones(1,8),zeros(1,20)];

subplot(321)

stem(x1,.)

title(u8(n));xlabel(n);ylabel(u(n));

h1=[ones(1,10),zeros(1,20)];

subplot(322)

stem(h1,.)

title(h1(n));xlabel(n);ylabel(h(n));

h2=[1,2.5,2.5,1,zeros(1,20)];

subplot(323)

stem(h2,.)

title(h2(n));xlabel(n);ylabel(h(n));

y1=conv(x1,h1)

subplot(324)

stem(y1,.)

axis([0,30,0,10])

title(y1(n));xlabel(n);ylabel(y(n));

y2=conv(x1,h2)

subplot(325)

stem(y2,.)

axis([0,30,0,10])

title(y2(n));xlabel(n);ylabel(y(n));

图2实验(2)波形

clc

clear

clf

x1=ones(1,50);

subplot(221)

stem(x1,.)

axis([0,50,0,1])

title(u(n));xlabel(n);ylabel(u(n));

A=[1,-1.8237,0.9801];

B=[1/100.49,0,-1/100.49];

y1=filter(B,A,x1);

subplot(223)

stem(y1,.)

axis([0,50,-0.05,0.05])

title(y(n));xlabel(n);ylabel(y(n));

n=0:100;

x2=sin(0.014*n)+sin(0.4*n);

subplot(222)

stem(x2,.)

axis([0,100,-1,2])

title(x(n));xlabel(n);ylabel(x(n));

y2=filter(B,A,x2);

subplot(224)

stem(y2,.)

axis([0,100,-1,1])

title(y(n));xlabel(n);ylabel(y(n));

图3实验(3)波形

2.实验二：时域采样与频域采样

A=444.128;

a=50*sqrt(2)*pi;

w=50*sqrt(2)*pi;

fs1=1000;

N1=0.064*fs1;

n=0:N1-1;

xn1=A.*exp(-a*n/fs1).*sin(w*n/fs1).*(n=0);

subplot(321)