matlab作业简要分析..docx

t=0:0.1:30;y=1-2*exp(-t).*sin(t);plot(t,y);xlabel(Time);ylabel(Amplitude);title(Decaying-oscillating Exponential);FOR语句：S=0;for m=1:100 S=S+m^2; if S=1000 break; end 每一个环节都对应一个ENDendS=S-m^2m=m-1WHILE语句：S=0;m=0;while S1000 m=m+1; S=S+m^2;endS=S-m^2m=m-1K=6;Z=[-3];P=[-1;-2;-5];[num,den]=zp2tf(Z,P,K) ZP零极点模型 TF传递函数模型 SS状态空间模型tf(num,den)[A,B,C,D]=zp2ss(Z,P,K)ss(A,B,C,D)num1=[10];den1=[1 1];num2=[1];den2=[2 0.5];num3=[540];den3=[1];num4=[0.1];den4=[1];[na,da]=series(num1,den1,num2,den2);[nb,db]=feedback(na,da,num4,den4);[nc,dc]=series(num3,den3,nb,db);[num,den]=cloop(nc,dc); cloop:单位反馈系统闭环函数tf(num,den)num1={1 1;0 1}; 分子写一个矩阵den1={[1 1],[1 2];1,[1 0]}; 分母写一个矩阵G1=tf(num1,den1) 传递函数矩阵num2={1 1;1 0};den2={[1 3],[1 1];[1 1],1};G2=tf(num2,den2)Gs=series(G1,G2) 串联Gp=parallel(G1,G2) 并联num=7*[2 3];den=conv([1 0],conv([1 0],conv([3 1],conv([1 2],conv([1 2],[5 0 3 8])))));tf(num,den)numo=[20];deno=[1 8 36 40 0];[numc,denc]=cloop(numo,deno); cloop:单位反馈系统闭环函数t=0:0.1:10; y=step(numc,denc,t); 单位阶跃响应曲线，t为仿真时间plot(t,y); Mp=((max(y)-1)/1)*100; disp([超调量Mp= ,num2str(Mp),%]); disp:屏幕输出函数 num2str:把数值转换成字符串wn=6;zeta=[0.2:0.2:1.0,2.0];figure(1);hold on hold on:在当前图的轴（坐标系）中画了一幅图，再画另一幅图时，原来的图还在for i=zeta num=wn^2;den=[1 2*i*wn wn^2]; step(num,den); gtext(num2str(i)); gtext:图片注释，标题endnumo=[1];deno=conv([1 0],conv([0.1 1],[0.2 3]));[numc,denc]=cloop(numo,deno); cloop:单位反馈系统闭环函数t=0:0.1:30;r=t;lsim(numc,denc,r,t) [y,x]=lsim(num,den,u,t)：u为系统输入信号，t为仿真时间，针对线性时不变模型，给定任意输入，得到输出响应。num=[1 5];den=conv([1 0],conv([1 2],[1 3]));rlocus(num,den); K=rlocfind(num,den)%移动鼠标到根轨迹与实轴上的交点，得到增益K Rlocus格式：[r,k]=rlocus(n,d)[r,k]=rlocus(g)[r,k]=rlocus(n,d,k)[r,k]=rlocus(g,k) 函数可计算出或画出SISO系统的根轨迹，其中g(或n,d)为对象模型，输入变量k为用户自已选择的增益向量，当k缺省时则为系统自动生成增益向量k,返回变量r为根轨迹各个点构成的复数矩阵.如果在函数调用中不返回任何参数，则rlocus函数在当前窗口中画出系统的根轨图。numo=[1 2];deno=conv([1 4 3], [1 4 3]);figure(1);rlocus(numo,deno);K0=rlocfind(numo,deno)%移动鼠标到根轨迹与虚轴上的交点，得到增益K0，使闭环系统稳定K的取值范围为：0KK0。K1=10;[numc1,denc1]=cloop(K1*numo,deno);figure(2); impulse(numc1,denc1); impulse是脉冲信