matlab 上机实验的主要内容.docVIP

MATLAB 上机实验 上机实验一 主要内容：系统频率响应；系统函数零、极点图；系统的单位脉冲响应。 10-1 ，求系统频率响应；系统函数零、极点图；系统的单位脉冲响应。 b=[0.2 0.3 1];%分子系数 a=[1 0.4 1];%分母系数 figure(1),zplane(b,a);%零极点图 %%求频率响应 b=[0.2 0.3 1];%分子系数 a=[1 0.4 1];%分母系数 figure(2),freqz(b,a);%数字系统频响 %%求系统的单位脉冲响应 b=[0.2 0.3 1];%分子系数 a=[1 0.4 1];%分母系数 figure(3),impz(b,a);%系统的单位脉冲响应 10-2 已知 分析其频率特性，并作、图。 w=[0:1:500]*2*pi/500;%[0，2pi]区域分为501点 X1=1; X2=1-0.9.*exp(-1*j*w); X=X1./X2; magX=abs(X); angX=angle(X).*180./pi; subplot(211);plot(w/pi,magX); xlabel(以pi为单位的频率); title(幅度部分); ylabel(幅度); subplot(212);plot(w/pi,angX); xlabel(以pi为单位的频率); title(相位部分); ylabel(相位); 试验图形： 10-3 已知，画出其幅频特性。 w=[0:1:500]*2*pi/500;%[0,pi]区域分为501点 X1=1-exp(-6*j*w);%分子多项式 magX=abs(X1); angX=angle(X1).*180/pi; subplot(211);plot(w/pi,magX); title(幅度部分); ylabel(振幅); xlabel(单位：pi); subplot(212); plot(w/pi,angX); xlabel(单位：pi); title(相位部分); ylabel(相位); 10-4 分析矩形序列的频率特性，并作出幅频特性图。 N=6 w=[0:1:500]*pi/500;%[0，pi]区域分为501点 X1=1-exp(-1*j*w*N); X2=1-exp(-1*j*w); X22=X2+(X2==0)*eps;%逻辑数组参加运算，使“0”被“机器零“代替 X=X1./X22; magX=abs(X); angX=angle(X).*180./pi; figure(1),subplot(211);plot(w/pi,magX); title(幅度部分);ylabel(幅度); subplot(212);plot(w/pi,angX); xlabel(以pi为单位的频率); title(相位部分);ylabel(相位); %% N=9 N=9; w=[0:1:500]*pi/500;%[0，pi]区域分为501点 X1=1-exp(-1*j*w*N); X2=1-exp(-1*j*w); X22=X2+(X2==0)*eps;%逻辑数组参加运算，使“0”被“机器零“代替 X=X1./X22; magX=abs(X); angX=angle(X).*180./pi; figure(2),subplot(211);plot(w/pi,magX); title(幅度部分);ylabel(幅度); subplot(212);plot(w/pi,angX); xlabel(以pi为单位的频率); title(相位部分);ylabel(相位); 10-5 求横向结构网络的频响图。 w=[0:1:500]*.2*pi/500;%[0,2pi]区域分为501点 X1=1-0.9^8.*exp(-8*j*w); X2=1-0.9.*exp(-1*j*w); X=X1./X2; magX=abs(X); angX=angle(X).*180./pi; subplot(211);plot(w/pi,magX); xlabel(以pi为单位的频率); title(幅度部分); ylabel(幅度); subplot(212);plot(w/pi,angX); xlabel(以pi为单位的频率); title(相位部分); ylabel(相位); 10-6 已知系统函数；求零、极点图；单位脉冲响应；系统频响。 b=[1 0 0 0 0 0 1]; a=[1 0 0 0 0 0 0]; w=[0:1:500]*pi/500; X1=1+exp(-6*j*w); X2=1%分母多项式 X=X1./X2; magX=abs(X); angX=angle(X).*180./pi; figure(1),subplot(221);impz(b,a,1