原创 实验二：离散系统的频域分析

实验二：离散系统的频域分析

1： 目的：了解函数
 （1）：离散时间傅立叶变换:freqz()
 （2）：离散傅立叶变换：fft（），ifft（）
 （3）：z变换zplane（），tf2zp（），zp2sos（）
2：函数功能以及源码：

%默认长为1，re=dsin(f,num)
function re="dsin"(f,num)
t=0:1/num:1;
w=2*pi*f;
re=sin(pi*2*f*t);
freqz():
已知序列p0,p1,p2,p3……………pm。
求出其离散时间傅立叶变换
Init;
numt=1000;
numw=1024;%(2的幂)
w=0:2*pi/(numw-1):4*pi;%由0到4*pi=2倍的0到2*pi，则重复两次
sin1=dsin(20,numt);
h=freqz(sin1,1,w);
hh=abs(h);
plot(hh);
ylabel('freqz');
结果：2-1

说明：y的值大的在20附近……。
 （2）：fft()
 已知序列in；求其离散傅立叶变换
 Init
 numt=1000;
sin2=dsin(20,numt);
h2=fft(sin2,1000);%把2*pi分为1000份去打印，与2-1图的一般相同，数越大越逼真
hh2=abs(h2);
plot(hh2);
ylabel('fft');
结果：2-2-1：取1000点，即把2*pi分为1000份

如果取100点：2-2-2

 (3):
1：zplane（）画出零级点图
已知：G(Z)=f(z)/g(z);中f(z)与g(z)的多项式系数
画出复数零级点图
Eg:
zplane([1 1 2 3],[9 8 0 7]);

结果：2-3-1:

2：tf2zp（）求出零级点值
已知：G(Z)=f(z)/g(z);中f(z)与g(z)的多项式系数，
求出零级点值
Eg:
[z,p,k]=tf2zp([1 1 2 3],[9 8 0 7]);%(k可以省略)
其中：z是零点，p为极点，k是增益常数
3：zp2sos????????????????????????????????????
已知：系统的零级点
求系统的G(Z)=f(z)/g(z)形式的f(z)与g(z)的多项式系数。
Eg：
cx=zp2sos(z,p,k)% k是增益常数, k可以省略.