实验六 信号与系统复频域分析

《信号与系统》实验指导书

函数。

function dplxy(k,r,A,B)

%The function to draw the frequency response of discrete system p=roots(A); %求极点 q=roots(B);

figure(1) ljdt(A,B)

w=0:l*pi/k:r*pi; y=exp(i*w);

N=length(p); M=length(q);

yp=ones(N,1)*y;

圆向量

yq=ones(M,1)*y; 圆向量

vp=yp-p*ones(1,r*k+1); 向量

vq=yq-q*ones(1,r*k+1); 向量 Ai=abs(vp);

的模 Bj=abs(vq);

%求零点

%画零极点图

%定义单位圆上的k个频率等分点 %求极点个数 %求零点个数

%定义行数为极点个数的单位

%定义行数为零点个数的单位

%定义极点到单位圆上各点的

%定义零点到单位圆上各点的

%求出极点到单位圆上各点的向量

%求出零点到单位圆上各点的向量

9

《信号与系统》实验指导书

的模

Ci=angle(vp); 的相角 Dj=angle(vq); 的相角

fai=sum(Dj,1)-sum(Ci,1); H=prod(Bj,1)./prod(Ai,1); figure(2) plot(w,H);

%绘制幅频特性曲线 %求系统相频响应 %求系统幅频响应

%求出零点到单位圆上各点的向量

%求出极点到单位圆上各点的向量

title('离散系统幅频特性曲线') xlabel('角频率') ylabel('幅度') figure(3) plot(w,fai)

title('离散系统的相频特性曲线') xlabel('角频率') ylabel('相位')

5/4(1?z?1)已知系统函数H(z)?，画频率响应和零极点图。

1?1/4z?1A=[1 -1/4]; B=[5/4 -5/4];

dplxy(500,2,A,B) %绘制系统2π频率范围内500个频

10

《信号与系统》实验指导书

率点的幅频和相频特性曲线及零极点图

三、上机实验内容

1．验证实验原理中所述的相关程序； 2．求信号f(t)?te?3tu(t)的拉普拉斯变换

s3?5s2?9s?73．求函数F(s)?的反变换 2s?3s?24．已知连续系统的系统函数如下，试用MATLAB绘制系统的

零极点图，并根据零极点图判断系统的稳定性

5．系统函数是

1?5z?5z?1?2s2?s?2H(s)=3 23s?5s?4s?6?z?3 求频率响应。

11