《matlab作业简要分析》由会员分享，可在线阅读，更多相关《matlab作业简要分析(9页珍藏版)》请在人人文库网上搜索。

1、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:100S=S+m2;if S=1000break;end 每一个环节都对应一个ENDendS=S-m2m=m-1WHILE语句：S=0;m=0;while S0); 求右半平面实轴的极点数P=length(ii) nyquist(num,den);N=0Z=P-Nif Z=0disp(Stable);elsedisp(Unstable。

2、);endz=;p=0 -1 -10;k=10;numo,deno=zp2tf(z,p,k);Gm,Pm=margin(numo,deno)margin函数可以从频率响应数据中计算出幅值裕度，相角裕度以及对应的频率，幅值裕度和相角裕度是针对开环系统而言，指示出系统闭环的相对稳定性，可在当前图形窗口中绘制出带有裕度以及相应频率显示的波特图，其中幅值裕度已分贝为单位re,im.w=nyquist(num,den)mag,phase,w=bode(num,den,w)gm,pm,wcp,wcg=margin(num,den)对应解释：re：奈奎斯特在频率为w时，频率响应的实部；im：奈奎斯特在频率为。

3、w时，频率响应的虚部；w：响应频率(单位：rad/sec)mag：波特图在频率为w时，频率响应的幅值；phase：波特图在频率为w时，频率响应的相位；w：响应频率(单位：rad/sec)gm：幅值裕度；pm：相角裕度；wcp：相角交界频率；wcg：截止频率；对于线性时变系统，如果对取定初始时刻t0的一个非零初始状态x0，存在一个时刻t1,(t1 t0)，和一个无约束的容许控制u(t)，tt0, t1，使状态由 x0 转移到 t1 时 x(t1)=0，则称此x0是在t0时刻为能控的。A=1 0;-1 0;B=1;0;n=length(A); 矩阵A的最大长度Uc=ctrb(A,B); 能控性if。

4、 rank(Uc)=n rank函数是求矩阵的秩,一个向量组的极大线性无关组所含向量的个数称为这个向量组的秩.disp(系统状态能控);elsedisp(系统状态不能控);end对于线性时变系统，如果对取定初始时刻t0的一个非零初始状态x0，存在一个有限时刻t1, (t1 t0)，使对所有的tt0, t1，有 y(t)=0，则称此x0是在t0时刻为不能观测的。A=0 1 0;0 0 1;-2 -4 -3;C=0 0 -1;1 2 1;n=length(A);Vo=obsv(A,C);if rank(Vo)=ndisp(系统状态能观测);elsedisp(系统状态不能观测);endz=;p=0 。

5、-1 -2;k=10;A,b,c=zp2ss(z,p,k); 极点型转状态空间n=length(A);rc=rank(ctrb(A,b); 先判断能控性if rc=ndisp(系统状态能控);P=-2 -1+j -1-j;K=acker(A,b,P) acker函数：A,B为系统的状态空间矩阵，向量P中是期望的闭环极点位置，返回值是增益向量。endt=0:0.05:20; 这边一定要注意上头的引号不要漏掉！(o)要转置成一个列向量！Tf=20;ut=t,10*ones(size(t); ut表示左列是t，右列是全是10的，个数为t阵列长的矩阵t,x,y=sim(Xiti8_6,Tf,ut) s。

6、im(刚才建立的SIM模型名称,要和m文件放在一个文件夹里面,运行时间Tf,选项可以不写,ut);plot(t,y)判断该系统是否为欠阻尼系统，若为欠阻尼系统，计算阶跃响应的过渡过程时间num=100; den=1 4 100;wn=sqrt(num); kesi=den(2)/2/wn;if kesi0); n1=length(ii);jj=find(real(z)0); n2=length(jj); if(n10)disp(System is unstable.);elsedisp(System is stable.);endif(n20)disp(System is a nonminim。

7、al phase one.);elsedisp(Syetem is a minimal phase one.); 分别利用for和while循环语句计算下式： for循环：S=0;for i=0:100S=S+2i;endS (5分) while循环：S=0;i=0;while i1); n=length(ii); (5分)if n0disp(System is unstable.);disp(r(ii);elsedisp(System is stable,); (4分)已知多输入系统的状态方程，采用状态反馈方法将系统闭环极点配置在处。A=0 0 5;1 0 -4;0 1 16;B=2 0;1。

8、 2;0 1; (4分)rc=rank(ctrb(A,B)if rc=length(A)P=-2+2*j -2-2*j -10;K=place(A,B,P) (6分)end随机输入一个自然数，判断能否被16、4、2整除。用if-else结构实现。x=input(please input a natural number );if rem(x,16)=0disp(it can be divided by 16);elseif rem(x,4)=0disp(it can be divided by 4); elseif rem(x,2)=0disp(it can be divided by 2);。

9、elsedisp(it cannot be divided by 16,4 and 2)end已知某一单位负反馈系统的开环传递函数为：绘制出闭环系统的根轨迹；确定并说明使系统产生临界振荡响应的K值，并绘制此时的单位脉冲响应。num=1 6;den=poly(0 -1 -4);rlocus(num,den)k,poles=rlocfind(num,den)nc,dc=cloop(k*num,den);impulse(nc,dc)幅值裕度和相位裕度：The amplitude and phase margin最小相位系统：minimum phase system开环脉冲传递函数：Open loop pulse transfer function临界振荡响应：The critical oscillation response单位负反馈系统：Negative feedback system闭环极点配置：The closed loop poles离散系统：Discrete system欠阻尼：Underdamped 无阻尼：No damping过渡过程时间：Transition time超调量：Overshoot。