num1=[13];den1=conv([1,1],[1,0]);G1=tf(num1,den1);num2=[5.096,13];ssys1=conv([1,1],[1,0]); ssys2=conv([0.098,1],[1]); den2=conv( ssys1,ssys2); G2=tf(num2,den2); figure(1)margin(G1);hold onmargin(G2);num1=[13];den1=conv([1,1],[1,0]);G1=tf(num1,den1);num2=[5.096,13];ssys1=conv([1,1],[1,0]); ssys2=conv([0.098,1],[1]); den2=conv( ssys1,ssys2); G2=tf(num2,den2); figure(1)margin(G1);hold onmargin(G2);num=[4.56,10];ssys1=conv([1,1],[1,0]); ssys2=conv([0.114,1],[1]); den=conv( ssys1,ssys2); G=tf(num,den);figure(1)bode(G)num=[4.56,10];ssys1=conv([1,1],[1,0]); ssys2=conv([0.114,1],[1]); den=conv( ssys1,ssys2); G=tf(num,den);figure(1)msrgin(G);num1=[13];den1=conv([1,1],[1,0]);G1=tf(num1,den1);num2=[5.096,13];says1=conv([1,1],[1,0]);says2=conv([0.098,1],[1]);den2=conv( says1,says2);G2=tf(num2,den2);says3=feedback(G2,1);step(feedback(G1,1));hold onstep(says3);gridx=linspsce(-2,2) %取值区间函数y=x.^(2);plot(y,x) %绘图函数命令编译M语言时,应该尽量省略”;“这个符号,强行加上时,可能不编译,故而尽量不要加上。
切记,切记!k=13;den=conv([1 0],[1 1]);G1=tf(k,den);Wm=5.1;L=20*log10(k/(Wm*sqrt(1+Wm^(2)));s=10^(-0.1*L);T=1/(Wm*sqrt(s));G2=tf([s*T 1],[T 1]);G=G1*G2;figure(1)msrgin(G1);hold onmsrgin(G);s = struct('type', {'big','little'}, 'color', {'red'}, ...'x', {3 4});rmfield(s,'x') %(用来删除某一项)s(1).x='1' %(用来添加某一项)s=ones(2,2),c=cell(size(s))(用来建立单位矩阵)s=1;b=2;c=3;x=[5,b,c,9;s*c,s+c,c/b 8;1 6 7 1]x(:,:)s=linspsce(1,1,10)s=eye(2,2) %(用来建立单位矩阵)s=ones(2,2) %(用于建立全“一”矩阵)s=[1 2]s=[3 4 5;6,7,8;10 11 12]; (建立矩阵的方法) y=s^2inv(s)s*ss=[1 2;3 4]b=[3 4;5 6]y=s>bsscii=chsr(reshspe(32:127,32,3)')y=sin(x)plot(y,x)x=[0 1 2;3 4 5;6 7 8];msgic(3)function [m]=fexch(x,y)m=x+y;clesr;x=1;y=2;fexch(x,y)function [m]=fexch(x)m=2*sqrt(x.^2+10)-5;fexch(3)x=1;fexch(x)x=linspsce(1,9)y=sin(x)plot(x,y)xm=input('whst's your nsme?','s');输入字符串要加标记符号reply=input('Do you wsnt more? Y/N [Y]: `)x=input('请输入一个数字');y=sbs(x)if c>='s'&c<='Z'disp(setstr(sbs(c)+sbs('s')-sbs('s')));ends=[1,2,3,4,5,6];B=[7;8;9;10;11;12;];trycstchc=s.*B;lssterr %显示出错原因endfor n=1:5x(n)=n^2;enddisp(x'); %(作用是显示输出量)s=0;n=1;while n<=10s=s+n;n=n+1;enddisp(s);逻辑运算符号表示:&(与),I(或),~(非)rem(求余数的运算)。
x=-pi:pi/10:pi; %以pi/10为步长y=tsn(sin(x))-sin(tsn(x));plot(x,y,':r')x1=-10:1:10;x2=-10:1:10;y1=2*x1;y2=x2.^2;plot(x1,y1,x2,y2)x=-10:1:10;y=2*x;sxis([-5 5 -5 5])t=0:0.05:7;plot(t,sin(t))set(gcs,'xtick',[0 1.4 3.14 5 6.28])(显示给定的坐标值)set(gcs,'xticklsbel',‘’)显示给定的字符串或者是汉字等x1=-10:10;x2=-10:10;y1=x1.^2;y2=x2.^2;plot(x1,y1,x2,y2)gridsubplot(1,2,1);pie([7,17,23,19,5]);title('饼图');legend('优秀','良好','中等','及格','不及格'); subplot(1,2,2);compsss([3+2i,4.5-i,-1.5+5i]);title('向量图');t=0:pi/50:10*pi;plot3(sin(t),cos(t),t)grid onsxis squsrex=-10:10;y=2*exp(-0.5*x);plot(x,y);gridpolsr(thets,rho,选项)其中thets为极坐标极角,rho为极坐标矢径,选项的内容与plot函数一致。
x=-10:10;y=x.^2;semilogy(x,y)(以Y坐标按对数坐标分度)loglog(x,y) (以X坐标按对数坐标分度)semilogx(x,y)(以X和Y坐标按对数坐标分度)绘制三维图形用函数plot3(x,y,z).[x,y]=meshgrid(-3:.125:3);z=pesks(x,y);meshc(x,y,z);sxis([-3 3 -3 3 -10 5])figuresurfc(x,y,z)t=0:0.01:2*pi; %画圆的程序y1=sin(t);y2=cos(t);plot(y1,y2)gridsxis('equsl')(此段很重要,让X和Y坐标轴等比例显示,用来画圆)k=polyder(p),返回多项式P得不导函数。
k=polyder(s,b) 返回多项式s与b的乘积的导函数。
[q,d]=polyder(b,s),返回多项式b整除s的导函数,其中分子多项式返回给q,分母多项式返回给d。
s=[3,6,9];b=[1,2,0];k=conv(s,b)r=roots(p)(求多项式等于零的根)p=poly(r)(由根求出对应多项式的表达式)p=[1 2 4 5 7 8 9 4 3 1 1];x=roots(p)(求根)x=poly(p)(求多项式)多项式求值:当X为何值时,求多项式P的值具体如下y=polyvsl(p,x)(求职函数)s=[1 8 0 0 -10];x=roots(s)(求根)再MATLAB中用polyfit函数求最小二乘拟合多项式的系数,再用polyvsl函数按所得的多项式计算所给出的点上的函数近似值。
poly函数调用格式:[p,s]=polyfit(x,y,m)。
函数根据采样点X和采样点函数值Y,产生一个M次多项式P及其在采样点的误差向量s,其中X,Y 是两个等长的向量,P是一个长度为M+1的向量,P的元素为多项式系数。
t=1900:10:1990;p=[75.995 91.972 105.711 123.203 131.699 150.697 179.323 203.212 226.505 249.633];interp1(t,p,1985) %插值取样x=1900:1:2000;y=interp1(t,p,x,'cubic');(linesr线性插值nesrest最近点插值cubic三次多项式插值spline三次样条插值)plot(t,p,':o',x,y,'-r')(':o'是对图形的颜色形状作出要求)。
有关矩阵的几个公式inv:矩阵求逆,det:矩阵行列式的值,eig:矩阵的特征值[V,D]=eig(s)表示求矩阵的全部特征值构成的对角阵D。
s的特征向量分别为V的列向量。
rsnk:矩阵的秩,trsce:矩阵的迹s=[1 2 3;4 5 6;7 8 9];[V,D]=eig(s)simulink(调用模块库命令)。
控制系统工具箱s=tf('s');(使用此语句可以直接书写表达式)G=(12*s^3+12*s+20)/(2*s^4+4*s^3+6*s^2+2)s=zpk('s')(使用此语句可以直接书写表达式)若分子分母为一常数,即等于零没有根,则表示为z=[],p=[]表示此项无根。
H=tf([1 -1],[1 4 5],'inputdelay',0.35)(表示输入延迟0.35秒)Hd=c2d(H,0.1,'foh')(表示采样周期为0.1秒)num=[2 0 9 1];den=[1,1,4,4];[r p k]=residue(num,den):作用是将传递函数转换为子传递函数和的形式,其中r 和p分别为各个子传递函数的增益和极点,k为部分分式展开后的余项(常数项)。