p32_18.m

eda实例DSP上机程序

%%%%%%%%%%%%%%%%%% 数字信号处理第一章习题P32第18小题程序%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%思路：若xn=单位脉冲信号，则输出yn=hn（即整个系统的单位脉冲响应）%%%%%%%%%%%%%%%%%%%%%%%%
%%%需要调用conv和filter函数
clc
clear
h1n=[1,1/2,1/4,1/8,1/16,1/32];
h2n=ones(1,6);
h3n=[1/4,1/2,1/4,zeros(1,97)];
%%%%%计算vn=h1n*h2n+h3n
h12n0=conv(h1n,h2n);
h12n=[h12n0,zeros(1,89)];
vn=h12n+h3n;
%%%调用filter计算hn等于T4对vn的响应(即整个系统的hn)
B4=[1,1];
A4=[1,-0.9,0.81];
hn=filter(B4,A4,vn);
%%%以下为绘图部分
n=0:length(hn)-1;
subplot(1,1,1);
stem(n,hn,'o');
%stem(n,hn,'*');
%stem(n,hn,'.');
xlabel('n');
ylabel('h(n)');