《2022年MATLAB基础知识考试复习总结 .pdf》由会员分享,可在线阅读,更多相关《2022年MATLAB基础知识考试复习总结 .pdf(9页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、名师总结精品知识点第一章基础系统仿真是根据被研究的真实系统的数学模型研究系统性能的一门学科,现在尤指利用计算机去研究数学模型行为的方法,即数值仿真。集计算,可视化及编程于一身。其主要产品模块构成:() ( ) toolbox ( 3) MATLAB Compiler ( 4)simulink (5) stateflow(6) Real-Time Workshop 。3MATLAB 语言被称为第四代计算机语言。有以下几个主要特点:编程效率高;使用方便;高效方便的科学计算;先进的可视化工具;开放性、可扩展性强;运行时动态连接外部 C或 FORTRAN应用函数;在独立 C或 FORTRAN程序中调用
2、 MATLAB函数;输入输出各种及其他标准格式的数据文件;创建图文并茂的技术文档;特殊应用工具箱;高效仿真工具Smulink 。4 变量命名规则:变量名、函数名对字母大小写敏感;变量名第一个字母必须是英文字母(只能是英文、数字和下连线) 。5 real(z)求复数 Z实部 imag(Z)求复数 Z 虚部 abs(z) 求复数 Z的模angle(Z)求复数 Z 的相角(单位是弧度)callback回校函数mdata=csvread( engdata.txt)clc清除指令窗 clf清除图形窗cd 设置当前工作目录clear清除工作空间保存的变量edit打开 M文件编辑器exit 、quit关闭、
3、退出 MATLAB 6 c=3+5i c=3+5*i (a=3;b=5;c=a+b*i) 1. 3e-4 2.78e23 A. 转置S.*B S./B B.S A.n 7Command History (历史指令)窗记录着用户在命令指令窗中所输入过的所有指令行,且所有这些被记录的指令行都能被复制,并送到指令窗中再运行。8Workspace Browser ( 工作空间浏览器)也叫内存浏览器,他保存了指令窗所使用过的全部变量(除非有意删除), 可通过该浏览器对内存变量进行操作。10 点击 MATLAB 桌面工具条上的?图标,或选择下拉菜单项【Help】 ,都能提供帮助;还提供现场帮助,用鼠标点亮
4、指令并点击右键,在弹出的菜单中选择【 Help On Selection】 。第二章数据及其运算1 简 单 数 组 生 成 方 法 : 逐 个 元 素 输 入 法 ; 冒 号 生 成 法 x=a:intc:b; x=linspace(a,b,n)=a:(b-a)/(n-1):b;logspace(w1,w2,n) 2 diag(3,3,3)产生对角形数组ans =3 0 0 0 3 0 0 0 3 eye(3) 产生单位数组ans =1 0 0 0 1 0 0 0 1 magic(3) 产生魔方数组ans =8 1 6 精品资料 - - - 欢迎下载 - - - - - - - - - - -
5、 欢迎下载 名师归纳 - - - - - - - - - -第 1 页,共 9 页 - - - - - - - - - - 名师总结精品知识点3 5 7 4 9 2 ones(3) 产生全 1 数组ans =1 1 1 1 1 1 1 1 1 rand(3) 产生 0、1 间随机数组ans =0.8147 0.9134 0.2785 0.9058 0.6324 0.5469 0.1270 0.0975 0.9575 zeros(3) 产生全 0 数组ans =0 0 0 0 0 0 0 0 0 zeros(3,2) ans =0 0 0 0 0 0 3 一位数组寻访 x=2.0000 1.04
6、72 1.7321 3.0000+5.0000i; x1=x(3) x1 =1.7321 x2=x(1 2 4) x2 =2.0000 1.0472 3.0000 + 5.0000i x3=x(2:end) x3 =1.0472 1.7321 3.0000 + 5.0000i x4=x(4:-1:1) x4 =3.0000 + 5.0000i 1.7321 1.0472 2.0000 4 二维数组寻访: A(: ,j) 表示 A矩阵的第 j 列全部元素; A(i , : )表示 A 矩阵第 i 行全部元素; A(1:3 ,2:4 )表示对 A矩阵取第 13 行,第 24列中所有元素。例: A=
7、rand(3,5) A =0.1190 0.3404 0.7513 0.6991 0.5472 0.4984 0.5853 0.2551 0.8909 0.1386 0.9597 0.2238 0.5060 0.9593 0.1493 A1=A(1,:) A1 =0.1190 0.3404 0.7513 0.6991 0.5472 A2=A(1:2,2:5) A2 =0.3404 0.7513 0.6991 0.5472 0.5853 0.2551 0.8909 0.1386 A3=A(1,3,2,5) A3 =0.3404 0.5472 0.2238 0.1493 5 非数( Not a N
8、umber ) 指的是 0/0 、/ 、0 x之类的运算,用 NaN或 nan 表精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页,共 9 页 - - - - - - - - - - 名师总结精品知识点示。非数参与运算所得结果也是非数,具有传递性;非数没有大小,不能比较两个非数的大小。6 在 MATLAB 中, 空 数组除了用 表示外,某维或若干维长度均为0 的数组都是 空 数组。7 逻辑运算符: &与| 或 非第三章数据和数组的可视化1 例 t=(0:pi/50:2*pi); % 生成( 101x
9、1)的时间采样列向量k=0.4:0.1:1; %生成( 1x7)的行向量X=cos(t)*k; %生成( 101x7)的矩阵plot(X); %绘制曲线,横坐标为每列元素对应的下标Plot(t,X) %以 t 为横坐标, X为纵坐标,按 Y的列数绘制曲线线型符号 -:-.- - 含义实线虚线点划线双划线色彩符号b g r cmykw含义蓝绿红青品红黄黑白2 grid on 画出分格线 grid off 不画分格线box on使坐标封闭box off使坐标开启title(s)书写图名text(x,y,s) 在(x,y) 处写字符注释xlable(s) 横坐标轴名ylable(s)纵坐标轴名hol
10、d on,hold off多层叠绘axis (0,1,-1,17)坐标范围alpha(0) 完全透明, alpha(0.5)半透明 , alpha(1) 完全不透明plotyy(X1,Y1,X2,Y2) 双纵坐标shadingflat,interp,faceted 3 view(az,el)通过方位角、俯视角设置视点,例view(-82,58)View(vx,xy,xz)通过直角坐标设置视点3 M 文件包括 M脚本文件和 M函数文件; M脚本文件的执行过程与在指令窗中直接输入指令的效果一样,但效率更高;M函数文件与脚本文件类似之处在于他们都是一个扩展名为“ .m”的文件; M函数文件通常由以下
11、几个部分组成:函数定义行、 H1行、函数帮助文本、函数体、注释。legend(s1,s2,)在图的右上角建立图例3 例:t=(0:0.02:2)*pi; x=sin(t);y=cos(t);z=cos(2*t); plot3(x,y,z,b-,x,y,z,bd) box on legend(链,宝石 ) 4 x=-4:4;y=x; X,Y=meshgrid(x,y); Z=X.2+Y.2; subplot(1,3,1),surf(X,Y,Z) %曲面图subplot(1,3,2),mesh(X,Y,Z) %网线图subplot(1,3,3),plot3(x,y,x.2+y.2);box on
12、% 曲线图5 x,y=meshgrid(-3:0.1:3,-2:0.1:2); z=(x.2+2*x).*exp(-x.2-y.2-x.*y); subplot(1,2,1),mesh(x,y,z) axis(-3,3,-2,2,-0.5,1.0) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页,共 9 页 - - - - - - - - - - 名师总结精品知识点title(toushi) hidden off %透视被叠压图形subplot(1,2,2),mesh(x,y,z) title(x
13、iaoyin) hidden on %消隐被叠压图形 axis(-3,3,-2,2,-0.5,1.0) 第四章 MATLAB 编程1for x=array 说明:循环体执行的次数由数组array 的列数决定(commands)例:for x=0.2:0.2:0.8 end b=sqrt(1-x2); end 2while expression说明:当 expression为真时,执行循环体直commands 到为假end 3if expression if expression commands commands end else end 4 ezplot(f) 其中 f=f(x,y)为用符号
14、函数表示的隐函数例:ezplot( x2+x*y+y2 - 10)5x,y=ginput(n)该指令用鼠标从图形上获取n 个点的坐标 (x,y )6 例:用冒泡法对数组A=1 6 8 9 3 5 6进行升序排序解: r=length(A); for i=1:r-1 for j=i+1:r if A(i)A(j) %A(i) sys=zpk(-0.6,-3 -2 -1,5) Zero/pole/gain: 5 (s+0.6) - (s+3) (s+2) (s+1) z,p,k=zpkdata(sys,v) z =-0.6000 p = -3 -2 -1 k =5 num,den=tfdata(s
15、ys,v) num =0 0 5 3 den =1 6 11 6 3 pzmap(sys)注:该指令将在 S平面上用符号“ O ”表示零点,符号“ X”表示极点。4 例:延时模型。 G(s)=e-0.5t(5s+3/s3+6s2+11s+6) num=5 3 ; sys=tf(num,den,inputdelay,tao) den=1 6 11 6; sys=zpk(z,p,k,inputdelay,tao) sys=tf(num,den,inputdelay,0.5) Transfer function: 5 s + 3 exp(-0.5*s) * - s3 + 6 s2 + 11 s +
16、6 5 串联模型: sys=series(sys1,sys2)=sys1*sys2并联模型: sys=parallel(sys1,sys2)=sys1+sys2 反馈连接: sys=feedback(sys1,sys2,sign) 第六章系统时间响应及其仿真1 对连续时间系统进行数字动态仿真,主要有两种方法:基于数值积分的仿真方法;基于离散相似法的仿真方法。2 数值积分方法所得到的离散数值解只是精确解的近似,其误差来自于两个方面,即舍入误差和局部截断误差。舍入误差与计算步长h 成反比,局部截断误差精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 -
17、 - - - - - - - - -第 5 页,共 9 页 - - - - - - - - - - 名师总结精品知识点与 h 成正比,显然选择一个合适的积分步长可使总误差达到最小。3 数值积分常用方法:欧拉法、梯度法、龙格- 库塔法、 Gear 法4 step(sys) step(sys,Tfinal) step(sys,T) 阶跃响应仿真说明:sys 是系统模型; Tfinal为仿真时间,若省略则由系统默认;T 为仿真时间向量,T=T0:dt:Tfinal,dt为连续系统离散化的采样周期, T0 为仿真开始时间。5 impulse(sys) 脉冲响应仿真6U,T=gensig(Type,Ta
18、u)信号发生器 u,t=gensig(square,5,10)说明:Type 为信号类型:sin 、 square 、 pulse ;Tau为信号周期; U为信号值向量;T为与 U对应的时间向量。7 lsim(sys,U,T) 任意输入相应函数8 initial(sys,x0 )用于计算零输入状态下,由初始状态x0 所引起的响应,只能用于状态空间模型。9T,y =ode45( f ,tspan,yo,options)为一种单步显式,变步长RK-45算法,用于求解非刚性微分方程;ode23 也一样; ode113属于多步法,求解非刚性。10 例:已知某系统运动方程及初始条件为:y1=y2y3;
19、y2=-y1y3; y3=-2y1y2 y1(0)=0; y2(0)=0.5; y3(0)=-0.5 求时间区间 t=0,20微分方程的解。解:(1) 建立描述系统微分方程的m-函数文件 rigit.m function dy=rigit(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-2*y(1)*y(2); (2) 编写调用函数 rigit()的 M文件,并执行T,y=ode45(rigit,0,20,0,0.5.-0.5); Plot(T,y(:,1),r,T,y(:,2),b*,T,y(:3),k-.); le
20、gend(y1,y2,y3) 11 已知二阶微分方程 y-(1-y2)y+y=0 y(0)=0,y(0)=1 求时间区间 t=0,20微分方程的解。解: (1)将微分方程表示成一阶微分方程组: y1=y; y1=y2; y2=(1-y12)y2-y1 (2) 建立描述系统微分方程的m-函数文件 vdp.m function dy=vdp(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=(1-y(1)2)*y(2)-y(1) (3)编写 MATLAB 主程序,并执行T,y=ode45(vdp,0,20,0,1); Plot(T,y(:,1),r-,T,y(:,2),b
21、:); legend(y1,y2) 第七章系统频率响应及其仿真精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页,共 9 页 - - - - - - - - - - 名师总结精品知识点1 频率响应是指系统对谐波输入的稳态响应;频率特性是指系统在正弦信号作用下,稳态输出与输入之比对频率的关系特性。G(jw)=Xo(jw)/Xi(jw)=A(jw)ej(w) 2 绘制 G (s)=11(s+1)/s(s2+15s+4)的幅频特性曲线和相频特性曲线解: num=11,11;den=1 15 4 0; w=0
22、.05:0.01:0.5*pi; Gw=polyval(num,j*w)./polyval(den,j*w); mag=abs(Gw); %(p=0) wcwg 不稳定subplot(2,1,1),plot(w,mag); grid on;title(幅频特性)ylable(|G| subplot(2,1,2),plot(w,theta); grid on;title(相频特性)Xlable(omega(rad/s),ylable(deg) 3 多项式函数 Y=polyval(P,X) P: 系数向量 X :自变量4 系统频率响应函数freqs(b,a,w) b、a 为传递函数分子分母系数向量
23、h=freqs(b,a,w0指定正实角频率,返回响应值;h,w=freqs(b,a)自动确定 200 个频率点,返回响应值和对应角频率freqs(b,a,w)绘制对指定正实角频率向量的幅值和相频特性曲线。5nyquist(sys) 绘制 nyquist图nyquist(sys,w) 指定频率范围 w ,绘制 nyquist图nyquist(sys1,sys2, ,sysn,w)在同一坐标系绘制多个nyquist Re,Im,w=nyquist(sys)返回频率响应的实部和虚部及对应w,不绘图6 bode(sys) 绘制 bode 图bodemag(sys)仅绘制幅频 bode 图7 margi
24、n(sys)绘制 Bode图,标出幅值与相位裕度。kg,y,wg,Wc=margin(sys) 返回幅值裕度 kg,相位裕度 y, 相位穿越频率 Wc和幅值穿越频率 Wg ,不绘制 Bode 图。第八章 控制系统的综合校正1 控制系统的性能指标大体上可分为两类,即频域指标和时域指标。2已知系统传递函数为G(s)=5(s+1)/(s3+4s2+6s+5),求最大超调量 Mp,调整时间 ts 和峰值时间 tp.解: num=5 5;den=1 4 6 5;sys=tf(num,den); finalvalue=polyval(num,0)/polyval(den,0); y,t=step(sys)
25、; %y()=lim sY(s)=limG(s) yp,k=max(y);tp=t(k); G(s)=Y(s)/(1/s)=sY(s) Mp=100*(yp-finalvalue)/finalvalue; len=length(t); while(y(len)0.98*finalvalue)&(y(len)1.02*finalvalue) len=len-1;end; ts=t(len);step(sys) 3 已知单位负反馈系统的开环传递函数为Gk=10/s(s+1)(s+5), 求其单位斜坡输入时,系统稳态误差。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢
26、迎下载 名师归纳 - - - - - - - - - -第 7 页,共 9 页 - - - - - - - - - - 名师总结精品知识点解:G k=zpk(,0 -1 -5,10); Xi=zpk(,0 0,1); sys=1/(1+Gk); Es=sys*Xi; %dcgain(sys)=lim sys(s) ess=dcgain(tf(1 0,1)*Es); %计算稳态偏差 ss=lim sE(s) t=0:0.05:10;xi=t; E(s)=Xi(s)/(1+G(s)H(s) y=lsim(sys*Gk,xi,t); E(s)=H(s)E1(s) plot(t,xi,r-.,t,y,
27、t,xi-y,k:) legend(shuru,shuchu,wucha,0); xlable(t(s),ylable(fuzhi、wucha) 4 频域指标包含: (1)开环频域指标: 主要指幅值穿越频率Wc,幅值裕度频率 Kg,相位裕度 等: (2)闭环频域指标:谐振峰值Mr,谐振频率Wr,闭环频宽Wcc等。5 PID 校正例:设校正前开环传递函数为G(s)=4/(2s+1)(0.5s+1)(0.05s+1),设计滞后控制器使系统相位裕度为60 度,开环增益为 49。解: (1)计算改变增益前后,系统开环对数频率特性Gp=tf(1,2,1)*tf(1,0.5,1)*tf(1,0.05,1)
28、*4; Gp1=Gp*49/4: figure(1),margain(Gp) figure(2),margain(Gp1) (2)计算改变增益后,具有希望相位裕度的系统开环对数幅频穿越频率Wc W=logspace(-1,2,100); 稳态误差: ess=lim e(t)=lim sE1(s) mag,ph=bode(Gp1,W); 稳态偏差: ss=lim (t)=lim sE(s) mag=reshape(mag,100,1);GB(s)=XO(s)/XI(s)=G(s)/(1G (s)H(s) ph=reshape(ph,100,1); wc=interp1(ph,W,-110) (3
29、) 确定校正环节mag110=interp1(ph,mag,-110); Beta=mag110 T1=6/Wc;BT1=Beta*T1; Gc=tf(T1,1,BT1,1) (4) 系统校核sys=Gc*Gp1 figure(3),margin(sys) (5) 时间响应比较figure(4) subplot(2,1,1),step(feedback(Gp1,1,-1) subplot(2,1,2),step(feedback(sys,1,-1), r)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页,共 9 页 - - - - - - - - - - 文档编码:KDHSIBDSUFVBSUDHSIDHSIBF-SDSD587FCDCVDCJUH 欢迎下载 精美文档欢迎下载 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 9 页,共 9 页 - - - - - - - - - -