MATLAB常用函数1、常用矩阵函数和命令d=eig(A)——矩阵特征值[v,d]=eig(A)——矩阵特征值和特征向量det(A)——行列式计算inv(A)——求逆poly(A)——特征多项式trace(A)——对角元素之和orth(A)——正交化2、常用特殊矩阵A=[]——表示空矩阵A=eye(n)——N维单位矩阵A=ones(n,m)——全部元素都为1的矩阵A=zeros(n,m)——全部元素都为0的矩阵A=rand(n,m)——元素为0到1之间均匀分布的随机矩阵A=randn(n,m)——元素为零均值单位方差正态分布的随机矩阵3、图形绘制和修饰(1)坐标轴调整axis([xminxmaxyminymax])——图形的x轴范围限定在[xminxmax],y轴范围限定在[yminymax]axis(‘控制字符串’)——根据字符串控制图形set(gca,‘xtick’,标示向量)——按照向量设置X坐标的刻度标示set(gca,’ytick’,标示向量)——按照向量设置Y坐标的刻度标示set(gca,‘xticklabel’,‘字符串|字符串…')——按字符串设置X坐标的刻度标示set(gca,'yticklabel’,‘字符串|字符串…')——按字符串设置Y坐标的刻度标示(2)文字标示title(‘字符串’)——图形标题xlabel(‘字符串’)——x轴标注ylabel(‘字符串’)——y轴标注text(x,y,‘字符串’)——在坐标处标注说明文字gtext(‘字符串’)——用鼠标在特定处标注说明文字(3)网格控制gridon——在所画的图形中添加网格线gridoff——在所画的图形中去掉网格线(4)图例注解legend(‘字符串1’,‘字符串2’,…,参数)——为区分各条曲线以字符串进行注解说明,参数确定注解在图形中的位置。(5)图形的保持holdon——保持当前图形及轴系的所有特性holdoff——解除命令(6)图形窗口的分割,可把一个绘图窗口分割成几个区域,分别绘图。subplot(m,n,p)——将当前绘图窗口分割成m行n列区域,指定第p个编号区域为当前绘图区。区域编号原则:“先上后下,先左后右”。4、线性系统分析设计(1)模型描述状态空间:写为A、B、C、D四个矩阵的形式;也可用SYS=SS(A,B,C,D)传递函数:用分子、分母的多项式表示;可用SYS=TF(NUM,DEN)建立tf模型零极点:用增益k、分子零点向量z、分母极点向量p表示,同样可用SYS=ZPK(Z,P,K)建立zpk模型部分分式:传递函数没有相同极点时,可与部分分式相互转换:[r,p,k]=residue(num,den)传递函数变换为部分分式形式[num,den]=residue(r,p,k)部分分式变换为传递函数形式闭环系统feedback反馈连接SYS=feedback(SYS1,SYS2,sign)[A,B,C,D]=feedback(A1,B1,C1,D1,A2,B2,C2,D2,sign)[num,den]=feedback(num1,den1,num2,den2,sign)sign=+(-)1反馈极性,缺省-1为负反馈,类似地,有parallel系统并联series系统串联(2)模型转换ss―状态空间、tf―传递函数、zp―零极点[num,den]=ss2tf(a,b,c,d,iu)——状态空间到传递函数[z,p,k]=ss2zp(a,b,c,d,iu)——状态空间到零极点[a,b,c,d]=tf2ss(num,den)——传递函数到状态空间[z,p,k]=tf2zp(num,den)——传递函数到零极点[a,b,c,d]=zp2ss(z,p,k)——零极点到状态空间[num,den]=zp2tf(z,p,k)——零极点到传递函数[r,p,k]=residue(num,den)——传递函数到部分分式[num,den]=(r,p,k)——部分分式到传递函数(3)系统分析(以下模型描述也可以用状态空间表达式的A、B、C、D)step(num,den)——直接画出单位阶跃响应曲线,时间是自动给定。[y,x]=step(num,den,T)——求系统在T时间向量下的单位阶跃响应输出,y是列向量,并返回状态x的变化。Impulse,lism用法类似stepbode(num,den,w)直接绘出bode图。[m,p]=bode(num,den,w)求出幅值和相角,可再用半对数坐标画图Nyquist(num,den,{wmin,wmax})直接绘Nyquist曲线,{wmin,wmax}——频率范围[real,imag,ww]=Nyquist(num,den,{wmin,wmax})——得到实部、虚部,可再绘图。rlocus(num,den)——直接绘出根轨迹[r,k]=rlocus(num,den)——返回增益和对应的复极点[r]=rlocus(num,den,k)——返回给定增益对应的复极点(4)状态反馈设计M=ctrb(A,B)——系统的能控矩阵M=[BABA2B…An-1B]N=obsv(A,C)——求取系统的能观矩阵N=[CCACA2…CAn-1]k=acker(A,B,p)——(单输入)对于期望极点p,求出系统的状态反馈增益阵k