dm13505.m

the code of the book come form the book of marlix laboratory with the book of singal

%dm13505
%绘制离散系统幅频特性曲线
num1=[1 -1];                   %定义系统（a）分子多项式系数向量
num1=num1*3/4;                 %定义系统（a）分母多项式系数向量
den1=[1 -1/2];
[H1,w1]=freqz(num1,den1);      %计算系统（a）频率响应样值
H1m=abs(H1);                   %计算幅频响应样值
w1=w1/pi;                      %归一化角频率
plot(w1,H1m)                   %绘制系统幅频特性曲线
title('幅频特性曲线')
xlabel('\theta/\pi');
ylabel('Magnitude')
%绘制系统（b）幅频特性曲线
num2=[1 0 -1];
num2=num2/4;
den2=[1 -3/4 1/2];
[H2,w2]=freqz(num2,den2);
H2m=abs(H2);
w2=w2/pi;                      
plot(w2,H2m)
title('幅频特性曲线')
xlabel('\theta/\pi');
ylabel('Magnitude')
%绘制系统（c）幅频特性曲线
num3=[1 -1 1];
num3=num3*3/4;
den3=[1 -3/4 1/2];
[H3,w3]=freqz(num3,den3);
H3m=abs(H3);
w3=w3/pi;                      
plot(w3,H3m)
title('幅频特性曲线')
xlabel('\theta/\pi');
ylabel('Magnitude')
%绘制系统（d）幅频特性曲线
num4=[1 zeros(1,9) 1];
den4=[1 zeros(1,10)];
num4=num4/2;
[H4,w4]=freqz(num4,den4);
H4m=abs(H4);
w4=w4/pi;                      
plot(w4,H4m)
title('幅频特性曲线')
xlabel('\theta/\pi');
ylabel('Magnitude')