p_74_32.m

eda实例DSP上机程序

%%%%%%%%%%%%%%%%%% 数字信号处理第二章习题P74第32小题程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc
clear
A=[1,-1.6,0.9425];
B1=1;
B2=[1,-0.3];
B3=[1,-0.8];
B4=[1,-1.6,0.8];
b1=[1,0,0];
b2=[1,-0.3,0];
b3=[1,-0.8,0];
b4=[1,-1.6,0.8];
p=roots(A);   %%%%%计算H1(z),H2(z),H3(z),H4(z)的零点
z1=roots(b1);    %%%%%计算H1(z)的零点
z2=roots(b2);    %%%%%计算H2(z)的零点
z3=roots(b3);    %%%%%计算H3(z)的零点
z4=roots(b4);    %%%%%计算H4(z)的零点
[h1n,n]=impz(B1,A,100);    %%%%计算单位脉冲响应h1(n)的100个样值
[h2n,n]=impz(B2,A,100);    %%%%计算单位脉冲响应h2(n)的100个样值
[h3n,n]=impz(B3,A,100);    %%%%计算单位脉冲响应h3(n)的100个样值
[h4n,n]=impz(B4,A,100);    %%%%计算单位脉冲响应h4(n)的100个样值
%%%%%以下是绘图部分
figure(1);
subplot(2,2,1);
zplane(B1,A);     %%%%绘制H1(z)的零极点分布图
subplot(2,2,2);
stem(n,h1n,'.');    %%%%绘制h1(n)的波形图
line([0,100],[0,0]);
xlabel('n');
ylabel('h1(n)');
subplot(2,2,3);
zplane(B2,A);       %%%%绘制H2(z)的零极点分布图
subplot(2,2,4);
stem(n,h2n,'.');     %%%%绘制h2(n)的波形图
line([0,100],[0,0]);
xlabel('n');
ylabel('h2(n)');
figure(2);
subplot(2,2,1);
zplane(B3,A);       %%%%绘制H3(z)的零极点分布图
subplot(2,2,2);
stem(n,h3n,'.');     %%%%绘制h3(n)的波形图
line([0,100],[0,0]);
xlabel('n');
ylabel('h3(n)');
subplot(2,2,3);
zplane(B4,A);       %%%%绘制H4(z)的零极点分布图
subplot(2,2,4);
stem(n,h4n,'.');     %%%%绘制h4(n)的波形图
line([0,100],[0,0]);
xlabel('n');
ylabel('h4(n)');