example9_1.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 53 行

%例9-1;example9_1
%(1)
h=[-4,1,-1,-2,5,6,5,-2,-1,1,-4];
M=length(h);n=0:M-1
[Hr,w,f,L]=hr_type(h);
%画图
figure(1)
subplot(2,2,1);stem(n,h);
m1=[0,10];m2=[0,0];
line(m1,m2);
title('脉冲响应偶对称，N=奇数')
subplot(2,2,3);plot(w/pi,Hr);
xlabel('pi');
title('类型1频率振幅响应')
%(2)
h=[-4,1,-1,-2,5,6,6,5,-2,-1,1,-4];
M=length(h);n=0:M-1
[Hr,w,f,L]=hr_type(h);
%画图
subplot(2,2,2);stem(n,h);
axis([0 12 -4 6]);
m1=[0,12];m2=[0,0];
line(m1,m2);
title('脉冲响应偶对称，N=偶数')
subplot(2,2,4);plot(w/pi,Hr);
axis([0 2 -30 30]);xlabel('pi');
title('类型2频率振幅响应')
%(3)
h=[-4,1,-1,-2,5,0,-5,2,1,-1,4];
M=length(h);n=0:M-1
[Hr,w,f,L]=hr_type(h);
%画图
figure(2)
subplot(2,2,1);stem(n,h);
m1=[0,10];m2=[0,0];
line(m1,m2);
title('脉冲响应奇对称，N=奇数')
subplot(2,2,3);plot(w/pi,Hr);
axis([0 2 -30 30]);xlabel('pi');
title('类型3频率振幅响应')
%(4)
h=[-4,1,-1,-2,5,6,-6,-5,2,1,-1,4];
M=length(h);n=0:M-1
[Hr,w,f,L]=hr_type(h);
%画图
subplot(2,2,2);stem(n,h);
axis([0 12 -4 6]);
m1=[0,12];m2=[0,0];
line(m1,m2);
title('脉冲响应奇对称，N=偶数')
subplot(2,2,4);plot(w/pi,Hr);
axis([0 2 -10 30]);xlabel('pi');
title('类型4频率振幅响应')