matlab与信号系统实验题库.docVIP

1．设，试用MATLAB画出该系统的幅频特性和相频特性。w=0:0.025:15; b=[1];a=[0.08,0.4,1]; H=freqs(b,a,w); subplot(2,1,1); plot(w,abs(H));grid; title(H(jw)的幅频特性); subplot(2,1,2); plot(w,angle (H));grid; title(H(jw)的相频特性); 2. 绘制出y（t)=ε（t-1)- ε（t-4)的幅频特性syms t w Gt=sym(heaviside(t-1)-heaviside(t-4)); Fw=fourier(Gt,t,w); FFw=maple(convert,Fw,piecewise); FFP=abs(FFw); ezplot(FFP,[-10*pi 10*pi]);grid; axis([-10*pi 10*pi 0 3.5]) 3.若某连续系统的输入为e(t)，输出为r(t)，系统的微分方程为：’(0_)=5,y(0_)=0;,求全响应。 a=[1 5 6];b=[0 3 2];%定义a，b [A B C D]=tf2ss(b,a);%定义A,B,C,,D的返回值 t=0:0.01:5;%定义时间范围 f=exp(-2*t); %定义函数 X0=[5,0];%系统状态变量X=[x1,x2,…..xn]在t=0时刻的初值。 y=lsim(A,B,C,D,f,t,X0);%绘制函数波形 plot(t,y);%定义横纵坐标 xlabel(t);%定义横标题 title(全响应y(t));%定义总横标题 grid on;%网格 4.求的拉氏变换式，并用MATLAB绘制拉氏变换在s平面的三维曲面图。 syms t s %定义符号变量 ft=sym(2*exp(-t)*Heaviside(t)+5*exp(-3*t)*Heaviside(t)); %定义时间函数f(t)的表达式 Fs=laplace(ft) %求f(t)的拉氏变换式F(s) FFss=abs(Fs); %求出F(s)的模 ezmesh(FFss); %画出拉氏变换的网格曲面图 ezsurf(FFss); %画出带阴影效果的三维曲面图 colormap(hsv); %设置图形中多条曲线的颜色顺序 5.求的拉氏变换式，并用MATLAB绘制拉氏变换在s平面的三维曲面图。 syms t s %定义符号变量 ft=sym(exp(-3*t)*sin(t)*Heaviside(t)); %定义时间函数f(t)的表达式 Fs=laplace(ft) %求f(t)的拉氏变换式F(s) FFss=abs(Fs); %求出F(s)的模 ezmesh(FFss); %画出拉氏变换的网格曲面图 ezsurf(FFss); %画出带阴影效果的三维曲面图 colormap(hsv); %设置图形中多条曲线的颜色顺序 7.已知系统函数，求其零极点图。 8.已知连续系统的系统函数如下，试用MATLAB绘制系统的零极点图，并根据零极点图判断系统的稳定性 num=[1 1 2]; %分子系数，按降幂顺序排列 den=[3 5 4 -6]; %分母系数，按降幂顺序排列 [z,p]=tf2zp (num,den); %求零点z和极点p zplane(z,p) %作出零极点图 9.已知描述离散系统的差分方程为：，且已知系统输入序列为， 求出系统的单位函数响应h(k)在-3 ~10离散时间范围内响应波形。a=[1,-0.5,0.25 -1]; b=[2,2]; impz(b,a,-3:10); title(单位响应) 10.绘制 图形。 t=-1:0.01:10;%定义时间样本向量 t1=0;%指定信号发生的时刻 f1=stepfun(t,t1);%产生左移位的阶跃信号 f2=2-exp(-2*t); f=f2.*f1; plot(t,f) %输出波形 axis([-1,10,0,2.1]) subplot(3,1,1);plot(t,f1);%绘制函数的波形 title(f1) subplot(3,1,2);plot(t,f