MATLAB基础与应用 精品.pptVIP

二、周期信号通过系统的响应 [例2] 周期方波通过RC系统的响应。 * 三、离散系统频率响应的计算 计算频率响应的MATLAB函数 b 分子的系数； a 分母系数； w 抽样的频率点(至少2点)， w在0~2p之间。 h = freqz(b,a,w) * 三、离散系统频率响应的计算 b=[1]; a1=[1 -0.9]; a2=[1 0.9]; w=linspace(0,2*pi,512); h1=freqz(b,a1,w); h2=freqz(b,a2,w); plot(w/pi,abs(h1),w/pi,abs(h2),:); legend(\alpha=0.9,\alpha=-0.9); 解： [例3]画出 的幅度响应曲线。 * 三、离散系统频率响应的计算 [例3]画出 的幅度响应曲线。 * 利用MATLAB进行连续系统的复频域分析 一、部分分式展开的MATLAB实现 二、H(s)的零极点与系统特性的MATLAB计算 * 一、部分分式展开的MATLAB实现 [r,p,k]=residue(num,den) num,den分别为X(s)分子多项式和分母多项式的系数向量。 r为部分分式的系数；p为极点；k为多项式的系数。若为真分式，则k为空。 * 二、H(s)的零极点与系统特性的MATLAB计算 计算多项式根roots的函数可用于计算H(s)的零极点。 r=roots(D) %计算多项式D(s)的根 H(s)零极点分布图可用pzmap函数画出，调用形式为 pzmap(sys) 表示画出sys所描述系统的零极点图。 * [例1] 用部分分式展开法求X(s)的反变换。 %program6_1 format rat %将结果数据以分数的形式输出 num=[1 2]; den=[1 4 3 0]; [r,p]=residue(num,den) 运行结果为 r = -1/6 ,-1/2 ,2/3 p = -3 ,-1 ,0 故X(s)可展开为 解： * [例2] 用部分分式展开法求X(s)的反变换。 % program6_2 num=[2 3 0 5]; den=conv([1 1],[1 1 2]); %将因子相乘的形式转换成多项式的形式 [r,p,k]=residue(num,den) magr=abs(r) %求r的模 angr=angle(r) %求r的相角 解： * [例2] 用部分分式展开法求X(s)的反变换。 运行结果为 r =-2.0000 + 1.1339i, -2.0000 - 1.1339i, 3.0000 p =-0.5000 + 1.3229i, -0.5000 - 1.3229i, -1.0000 k =2 magr =2.299, 2.2991, 3.0000 angr =2.6258, -2.6258, 0 故X(s)可展开为 * [例3]试画出系统 的零极点分布图，求其单位冲激响应h(t)和频率响应H(j?)，并判断系统是否稳定。 num=[1];den=[1 2 2 1]; sys=tf(num,den); poles=roots(den) figure(1);pzmap(sys); t=0:0.02:10; h=impulse(num,den,t); figure(2);plot(t,h) title(Impulse Respone) [H,w]=freqs(num,den); figure(3);plot(w,abs(H)) xlabel(\omega) title(Magnitude Respone) 解： * 运行结果 * 利用MATLAB进行离散系统的复频域分析 一、部分分式展开的MATLAB实现 二、H(z)的零极点与系统特性的MATLAB计算 * 一、部分分式展开的MATLAB实现 [r,p,k]=residuez(num,den) num,den分别为X(z)分子多项式和分母多项式的系数向量。 r为部分分式的系数；p为极点；k为多项式的系数。若为真分式，则k为空。 * 二、H(z)的零极点与系统特性的MATLAB计算 利用tf2zp函数计算H(z)的零极点，调用形式为 [z,p,k]=tf2zp(b,a) H(z)零极点分布图可用zplane函数画出，调用形式为 zplane(b,a) b和a分别