matlab z变换

clear all;
close all;
clc;
%  syms n k z;
% f=1/n;         %定义离散信号
%  F=ztrans(f)              %z变换
%  pretty(F);
%  fk=iztrans(F,k)
%  pretty(fk)
%
% b=[0 0 10 0];%分子的系数数组
% a=[1 -5 8 -4]; %分母的系数数组
% zplane(b,a)% 使用zplane函数绘制如下系统的零极点分布图
% [z,p,k]=tf2zp(b,a)%使用tf2zp求出系统函数的零、极点和增益运算结果
% [b,a] = zp2tf(z,p,k) %根据求出的零、极点和增益，然后自学使用zp2tf还原出 分子和分母的系数
% [sos,g]=tf2sos(b,a)%使用tf2sos将系统函数分解成一系列二阶子系统的级联形式
% [b,a]=sos2tf(sos,g)%根据求出的一系列二阶子系统，使用sos2tf还原出)(zH分子和分母的系数

b=[0 1 1 1]; %分子的系数数组
a=[1 2 2 2]; %分母的系数数组
n=(0:500)*pi/500; %在pi范围内取501个采样点
[h,w]=freqz(b,a,n);%求系统的频率响应
subplot(2,1,1),plot(n/pi,abs(h));
grid %作系统的幅度频响图
axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]);
ylabel('幅度');
subplot(2,1,2),plot(n/pi,angle(h));
grid %作系统的相位频响图
axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]);
ylabel('相位');
xlabel('以pi为单位的频率');

  33:  figure;plot(w,abs(h));

“Through these he has given us his very great and precious promises, so that through them you may participate in the divine nature, having escaped the corruption in the world caused by evil desires.”