第七讲MATLAB在电路中的应用解说.ppt

二、周期信号通过系统的响应 例2 周期方波通过RC系统的响应。 二、周期信号通过系统的响应 例2 周期方波通过RC系统的响应。 %p5_2 Periodic signal pass LTI system file123.m T=4;w0=2*pi/T;RC=0.1; t= -6:0.01:6;N=51; c0=0.5;xN=c0*ones(1,length(t)); %dc for n=1:2:N % even harmonics are zero H=abs(1/(1+j*RC*w0*n)); phi=angle(1/(1+j*RC*w0*n)); xN=xN+H*cos(w0*n*t+phi)*sinc(n*0.5); end plot(t,xN); xlabel([time RC=,num2str(RC)]);grid; set(gca,xtick,[-5 -3 -1 0 1 3 5]); 二、周期信号通过系统的响应 例2 周期方波通过RC系统的响应。 三、离散系统频率响应的计算 计算频率响应的MATLAB函数 b 分子的系数 a 分母系数 w 抽样的频率点(至少2点)， w在0~2p之间 幅频特性: abs, 相频特性: angle h = freqz(b,a,w) b=[1]; w=linspace(0,2*pi,512); h2=freqz(b,a2,w); plot(w/pi,abs(h1),w/pi,abs(h2),:); 三、离散系统频响特性的计算 a1=[1 -0.9]; a2=[1 0.9]; h1=freqz(b,a1,w); legend(\alpha=0.9,\alpha=-0.9); %file124.m 三、离散系统频响特性的计算 利用MATLAB进行连续系统的复频域分析 部分分式展开的MATLAB实现 H(s)的零极点与系统特性的MATLAB计算 一、部分分式展开的MATLAB实现 [r,p,k]=residue(num,den) num,den分别为F(s)分子多项式和分母多项式的系数向量。 r为部分分式的系数，p为极点，k为多项式 二、H(s)的零极点与系统特性的MATLAB计算 计算多项式根roots的函数可用于计算H(s)的零极点。 r=roots(N) %计算多项式N的根 H(s)零极点分布图可用pzmap函数画出，调用形式为 pzmap(sys) 表示画出sys所描述系统的零极点图。 例1 试画出系统 的零极点分布图，求其单位冲激响应h(t)和频率响应H(j?)，并判断系统是否稳定。 %file125.m 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(Magnitu) 利用MATLAB进行离散系统的Z域分析 部分分式展开的MATLAB实现 H(z)的零极点与系统特性的MATLAB计算 一、部分分式展开的MATLAB实现 [r,p,k]=residuez(num,den) num,den分别为F(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分别为H(z)分子多项式和分母多项式的系数向量。 返回值z为零点、p为极点、 k为增益常数。 例1 试画出系统 的零极点分布图，求其单位冲激响应h[k]和频率响应H(ejW) 。 % file126.m b =[1 2 1];a =[1 -0.5 -0.005 0.3]; figure(1);zplane(b,a); n