第九讲_MATLAB在电路中应用.pptVIP

本文观看结束！！！ 例2 试画出系统 的零极点分布图，求其单位冲激响应h[k]和频率响应H(ejW) 。 % program 8_2 b =[1 2 1];a =[1 -0.5 -0.005 0.3]; figure(1);zplane(b,a); num=[0 1 2 1]; den=[1 -0.5 -0.005 0.3]; h=impz(num,den); figure(2);stem(h) xlabel(k) title(Impulse Respone) [H,w]=freqz(num,den); figure(3);plot(w/pi,abs(H)) xlabel(Frequency \omega) title(Magnitude Respone) 例1 试用MATLAB画出图示周期三角波信号的频谱。 解：周期信号的频谱为 画三角波信号频谱的MATLAB程序 N=8; n1= -N:-1; %计算n=-N到-1的Fourier系数 c1= -4*j*sin(n1*pi/2)/pi^2./n1.^2; c0=0; %计算n=0时的Fourier系数 n2=1:N; %计算n=1到N的Fourier系数 c2= -4*j*sin(n2*pi/2)/pi^2./n2.^2; cn=[c1 c0 c2]; n= -N:N; subplot(2,1,1); stem(n,abs(cn));ylabel(Cn的幅度); subplot(2,1,2); stem(n,angle(cn)); ylabel(Cn的相位);xlabel(\omega/\omega0); 程序运行结果 例2 求周期矩形脉冲的Fourier级数表示式。并用MATLAB求出由前N项Fourier级数系数得出的信号近似波形。 取A=1, T=2, t=1, w0=p % Gibbs phenomenon t=-2:0.001:2; N=input(Number of harmonics= ); c0=0.5; cN=c0*ones(1,length(t)); %dc component for n=0:2:N-1 % even harmonics are zero cN=cN+cos(pi*n*t)*sinc(n/2); end plot(t,cN); % Gibbs phenomenon N=5 N=15 N=50 N=500 二、用数值积分分析非周期信号频谱 数值函数积分quad8可用来计算非周期信号频谱 F 是一个字符串，它表示被积函数的文件名。 a,b 分别表示定积分的下限和上限 quad8的返回是用自适应Simpson算法得出的积分值 y = quad8(F,a,b) 例3 试用数值方法近似计算三角波信号的频谱 F(jw)= Sa2(w / 2) 解: 图示三角波可表示为 三角波信号频谱的理论值为 例3 试用数值方法近似计算三角波信号的频谱 function y=sf1(t,w);y=(t=-1 t=1).*(1-abs(t)).*exp(-j*w*t); w=linspace(-6*pi,6*pi,512); N=length(w);F=zeros(1,N); for k=1:N F(k)=quad8(sf1,-1,1,[],[],w(k)); end figure(1); plot(w,real(F));title() xlabel(\omega);ylabel(F(j\omega)); figure(2); plot(w,real(F)-sinc(w/2/pi).^2); xlabel(\omega);title(计算误差); 例3 试用数值方法近似计算三角波信号的频谱 运行结果 利用MATLAB进行系统频域分析 连续系统频响特性的计算 周期信号通过系统的响应 离散系统频响特性的计算 一、连续系统频响特性的计算 b 分子多项式系数 a 分母多项式系数 w 需计算的H(jw)的抽样点 (数组w中少需包含两个w的抽样点)。 计算频响的MATLAB函数 H=freqs(b,a,w) 一、连续系统频响特性的计算 例1 三阶归一化的Butterworth低通滤波器的系统函数为 w=linspace(0,5,200); b=[1];a=[1 2 2 1]; h=freqs(b,a,w); subplot(2,1,1); plot(w,abs(