《2022年MATLAB知识点资料讲解 .pdf》由会员分享，可在线阅读，更多相关《2022年MATLAB知识点资料讲解 .pdf（25页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、资料收集于网络，如有侵权请联系网站删除word 可编辑1，计算矩阵的特征值和特征向量【V,D】=eig（A）其中 A 指的是矩阵V 指的是特征向量组成的矩阵D 指的是特征值所组成的矩阵A=0 1 0;1 0 1;0 1 0; V,D=eig(A) V = 0.5000 -0.7071 0.5000 -0.7071 0.0000 0.7071 0.5000 0.7071 0.5000 D = -1.4142 0 0 0 -0.0000 0 0 0 1.4142 2，求线性方程的解1231231232314354578950 xxxxxxxxx这是一个非齐次方程方程，对于线性代数来讲，很难求解，但
2、是利用MATLAB来求A=1 2 3;3 -5 4;7 8 9; 1 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑B=14 5 50; 2 x=AB 3 B=14;5;50; 4 这 4 个语句的区别，特别是2.3 两句的区别， 2 语句中 B 表示的是 1 行 3 列，但是 4 语句中的表示的3 行 1 列。同时， X=AB很 X=A/B也是很有能区别的，具体的区别如
3、下：X=AB表示的是 A*X=B 的解X=A/B 表示的是 X*A=B 的解，在这里，具体解释如下；A*X=B，在线性代数中， X=A(-1)B , 所以 X=AB。X*A=B的解为 X=BA(-1),就是 X=B/A。3，一元方程求根这个在计算中是个难题，一元多次方程式很难求解的，对于纯粹的手工计算来说，所以有MATLAB比较方便，具体语句如下：231.350.6680.4360.69552fxxxxP=0.69552 0.436 0.668 1.35; X=roots(P) X = 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - -
4、 - 名师精心整理 - - - - - - - 第 2 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑0.2817 + 1.2456i 0.2817 - 1.2456i -1.1902 + 0.0000i 对于这个方程来说，要注意的是，在构成向量时，一定要从告辞往低次排列， 中间缺少次数的， 用 0 代替，求解语句是 roots （P） 。4，图形处理功能A, 可以绘制函数图像，具体的语句如下： x=linspace(0,6); 这个表示 X轴，在 0 到 6 内取 100 个点，这 是 默 认 的 数 值100 。 当 然 还 有
5、 可 以 设 定 数 值 的 方 法linspace(0,6,100)表示的是在 0 到 6 内取等间隔取100 个点。 y1=sin(2.*x); plot(x,y1)这是绘制函数sin （2X） 的图像，在这里要注意的，这个乘号的作用，其中还有一个“点号”。这是和普通的数学中的乘法不同的地方。y1=sin(2.*x);y2=sin(x.2);y3=(sin(x).2 这是其他的图像，分别是 sin2x sin(x2) (sinx)2，但是在这3 个式中，用 MATLAB表达时，总会用到“点号”这个符号，这是要注意的地方。绘图用plot 函数就可以吧。名师资料总结 - - -精品资料欢迎下载
6、 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑B，在同一个窗口绘制多个函数的图像，但是并不是在同一个坐标系中绘制。具体语句如下：例题：用四种方法描述cos(x)*sin(y)图形，分别采用以下的集中函数，分别是 surf，mesh ,meshc,waterfall 进行函数图像的绘制具体语句如下：x=linspace(-10,10,100); y=linspace(-10,10,100); x1,y1=mesh
7、grid(x,y); x 是 n 为向量，y 是 m 维向量，如 x1=1 ;2 ;3;y1=4; 5 则x1,y1 就产生一系列坐标点，(1,4)，(1,5) (2,4)，(2,5) (3,4)，(3，5) 一般在作 3D 图像的时候遇到，这是为了把XY 变成矩阵，如果不这样做，就会产生报错。但是在前面做2D 图像时就不会用到meshgrid函数，这是在做3D名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 25 页 - - - - - - - - - 资料收集于网络，如
8、有侵权请联系网站删除word 可编辑图像时要注意到的地方。z=cos(x1).*sin(y1); subplot(2,2,1);surf(x1,y1,z); subplot(2,2,2);mesh(x1,y1,z); subplot(2,2,3);meshc(x1,y1,z); subplot(2,2,4);waterfall(x1,y1,z)当没有用 surf,mesh ,meshc, waterfall 时，图像时画不出来的。当我们确实 waterfall（x1,y1,z）语句时，绘制的图像具体如下，第四幅图缺失了图像。5,矩阵的加法具体的语句如下， 只要注意不要把矩阵的表示方法搞错了就可
9、以了a=1,2,3;4,5,6;7,8,9; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑 b=1,1,1;2,2,2;3,3,3; a+b ans = 2 3 4 6 7 8 10 11 12 6，矩阵的乘法在这之前，首先要注意的是A*B 和 A.*B 是不同的计算，所以不能够乱用。 A*B 表示的是我们学习的数学中的算法。但是A.*B表示的是两个矩阵中相同的位置相乘
10、，是不同的用法，要注意。具体的算法如下；a=1,2,3;4,5,6;7,8,9; c=1 2 3;4 5 6;7 8 9; a*c ans = 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑30 36 42 66 81 96 102 126 150 a.*c ans = 1 4 9 16 25 36 49 64 81 7，超定系统首先要了解什么是超定系统，超定系统是对于解
11、决多元方程而言的，同时方程的个数要多于自变量的个数，这个就是超定方程x1+2x2=1 2x1+3x2=2 3x1+4x2=3如何去解这个方程，具体的语句如下：a=1 2;2 3;3 4;b=1;2;3; ab 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑ans = 1.0000 0.0000 有超定方程，那么必然就有欠定方程，欠定方程的定义是，方程的个数小于自变量的个数
12、。如下；x1+2x2+3x3=1 2x1+3x2+4x3=2 a=1 2 3;2 3 4; b=1;2;ab ans = 1 0 0 8，矩阵的幂运算a=1 2 3;2 3 4; C=a1.5; 和普通计算的区别不是很大9，矩阵的转置名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑分为共轭转置和非共轭转置两种情况，共轭转置用语句“a ”,非共轭转置用“ a. ”.有个在前面
13、用过的逗号运算符. a=1+2i 2+4i; a(共轭转置，不仅进行了转置，还进行了共轭的处理) ans = 1.0000 - 2.0000i 2.0000 - 4.0000i a. ans =1.0000 + 2.0000i 2.0000 + 4.0000i 10，矩阵的关系和逻辑运算 a=0 -1 2; b=-3 1 2; a a ab ans = 1 0 0 a=b ans = 1 0 1 a=b ans =0 0 1 a=b ans = 1 1 0 也就是比较两个矩阵的关系，当矩阵的关系是真的时候，则输出为 1。当矩阵的关系是假的时候，就会输出为0。当然，还存在或与非 3 种逻辑运算。
14、11，矩阵的函数运算具体的使用函数如下：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑inv矩阵求逆det行列式的值eig矩阵的特征值和特征矢量diag对角矩阵trace矩阵的迹rank矩阵的秩pinv矩阵的伪逆norm矩阵或矢量的范数求解方程的解法，主要是利用rank 函数， rank 函数是求解这个矩阵的秩的。 来判断这个方程是不是有一个解，当系数矩阵和增广矩阵有相
15、同的秩的时候，只有一个解， 但是当他们的秩是不同的时候，要根据法则来判断，具体的事例如下：c = 5 2 -9 -18 -9 -2 2 -7 6 7 3 29 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑 rank(c) ans = 3 12，矩阵的分解A，特征值得分解利用 eig 函数，事例 c,d=eig(A),这个在前面已经解释过了，参看前面的解释。B.奇异值分
16、解利用 svd 函数， u,s,v=svd(A) C,三角分解l,u=lu(A),其中 l 代表下三角， u 代表上三角D，Cholesky（乔里斯基）分解如果 A 为 n 阶对称正定矩阵， 则存在一个非奇异的上三角实矩阵L，使：LLAT。当限定 L的对角元素为正时，这种分解时唯一的，称为Cholesky分解。使用 chol（a）就可以了。E，QR分解（也称为正交分解）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 25 页 - - - - - - - - - 资料收集
17、于网络，如有侵权请联系网站删除word 可编辑实矩阵 A可以写成AQR的形式， 其中 Q 为正交阵， R为上三角阵。规定若R的对角元为正数，则分解唯一。q,r=qr(A) 13,矩阵的特殊操作特殊矩阵：空阵： 创建空阵。全 0 阵：各个元素都为零的矩阵，函数zeros。a=zeros(M,N)，生成 M 行 N 列的矩阵。a=zeros(size(B) ，生成与矩阵B维数相同的矩阵。单位阵：对角线元素为1，其他元素为0 的矩阵，函数eye(m,n)。全 1 阵：各个元素都为1 的矩阵，函数ones。随机阵：矩阵元素由随机数构成的矩阵。函数 rand、randn。rand(M,N)， 生成 M
18、行 N列随机矩阵，矩阵元素值在区间(0,1)之间。randn(M,N)，生成 M 行 N 列随机矩阵，矩阵元素值服从正态分布N(0,1)。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑14，变维操作一般先排列，在排行。a=1:12 a = 1 2 3 4 5 6 7 8 9 10 11 12 reshape(a,3,4) ans = 1 4 7 10 2 5 8 11
19、3 6 9 12 15，矩阵的翻转对矩阵进行左右、上下翻转、旋转等操作。fliplr ：左右翻转flpdim：第 n 维翻转flipud ：上下翻转名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑rot90：逆时针旋转9016，矩阵的抽取函数 diag 实现矩阵对角元素的抽取：c=diag(a,n)，c 为抽取矩阵 a 的第 n 条对角线所创建的元素矢量。a=diag(c
20、,n)，创建对角矩阵a，使矢量 c 成为 a 的第 n 条对角线矢量。 n=0 或不指定 n 时，为主对角线。函数 tril 实现下三角矩阵抽取：c=tril(a,n)，抽取矩阵 a 的第 n 条对角线下面的部分，包括第 n 条对角线。函数 triu 实现上三角矩阵抽取：c=triu(a,n)，抽取矩阵 a 的第 n 条对角线上面的部分，包括第 n 条对角线。17，如何取出矩阵A 中的数组b=A(3:5,1:2)表示的是把A矩阵的 3 到 5 行， 1 到 2 列取出来。18，如何取出数组A中的某个元素名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - -
21、 - - - - - 名师精心整理 - - - - - - - 第 15 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑x=rand(1,5) x = 0.7577 0.7431 0.3922 0.6555 0.1712 x(3) ans = 0.3922 x(1:3)取出数组的元素的前3 个ans = 0.7577 0.7431 0.3922 x(1 3 5)取出数组的元素的第1,3,5 个ans = 0.7577 0.3922 0.171 符号运算功能18，创建符号矩阵利用 sym( ),方括号中是符号。查找函数 findsym
22、(A,N),在矩阵 A 中找到 N 个与 X相近的数字。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑19，求导数Diff(f,x)这句话表示的是F函数中， 对 X进行求导数。 具体的语句如下：Syms x n %定义两个自变量字符f=xn %定义函数 f 的定义式Diff(f,x) %对 f 中 x 求导数ans = n*x(n - 1) diff(f,n) %定义对
23、 f 中的 n 求导数ans = xn*log(x) 20，关于多项式的乘法和除法具体的语句如下：a=1 -5 3 -4 2; b=1 2 -5 3; conv(a,b) %乘法语句 conv() ans = 1 -3 -12 30 -36 33 -22 6 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑 deconv(ans,b) %除法语句 deconv() ans
24、 = 1 -5 3 -4 2 21，因式分解函数是 factor（F），F表示的是函数。具体的例子如下：sym x; f=x9-1; factor(f) ans = x - 1, x2 + x + 1, x6 + x3 + 1 22，算式展开利用 expand 函数，具体的例子如下：syms x y f=(x+1)5; expand(f) ans =x5 + 5*x4 + 10*x3 + 10*x2 + 5*x + 1 f=sin(x+y); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - -
25、- 第 18 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑 expand(f) ans =cos(x)*sin(y) + cos(y)*sin(x) 但是在这里要注意的关于sym 和 syms 函数的区别， sym 适合定义一个符号，但是syms适合定义多个运算符号。23，合并同类项Collect 函数，具体的语句如下：syms x t f=x*x*(x-6)+12*t; expand(f) ans = t*x3 - 6*t*x2 + 12*t*x collect(f) ans = t*x3 - 6*t*x2 + 12*t*x f
26、actor(f) ans = t, x, x2 - 6*x + 12 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑24，分式的通分syms x y f=x/y+y/x; numden(f) ans = x2 + y2 n,d=numden(f) n = x2 + y2 d =x*y 25，讲多项式写成嵌套的形式sym x ans = x f=x3+6x2+11x-6;
27、 f=x3+6x2+11x-6; 错误 : 不应为MATLAB 表达式。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑是不是想输入 : f = x3 + 6*x2 + 11*x - 6; horner(f) ans = x*(x*(x + 6) + 11) - 6 26，替换函数subs( 函数式 ，要替换的字符，带入的字符) 具体的语句如下：Syms a b Subs
28、(a+4*b,a,3) Ans=3+4*b 还有另外一种写法：就是不用syms 函数来定义字符Subs( a+4*b ,a,3) Ans=3+4*b 我们在没有使用syms 函数式，所以a,b 就必须要用这个符名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑号要定义。27，精度的控制可控精度运算：函数vpa vpa(S) ， 显示符号表达式S在当前精度 D 下的值。D 是
29、使用 digits函数设置的数值精度。vpa(S,D) ，显示符号表达式S在精度 D 下的值， D 不是当前的精度值，而是临时使用digits 函数设置为 D 位精度。设定所用数值的精度：函数digits digits(D)，设置数值的精度为D 位。digits，在命令窗口显示当前设定的数值精度。具体的语句如下：a=1/4,exp(1);log(3),7/3 a = 0.2500 2.7183 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 22 页，共 25 页 - - - - -
30、 - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑1.0986 2.3333 vpa(a,10) ans = 0.25, 2.718281828 1.098612289, 2.333333333 a=sym(1/4,exp(1);log(3),7/3) 警告 : Support of strings that are not valid variable names or define 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 23 页，共 25 页 -
31、- - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑a number will be removed in a future release. To create symbolic expressions, first create symbolic variables and then use operations on them. In symconvertExpression (line 1536) In symconvertChar (line 1441) In symtomupad (line 1198) In sym (line 177) a =
32、1/4, exp(1) log(3), 7/3 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 24 页，共 25 页 - - - - - - - - - 资料收集于网络，如有侵权请联系网站删除word 可编辑 vpa(a,10) ans = 0.25, 2.718281828 1.098612289, 2.333333333 在这里，有两种写法，sym 和没有 sym，一个是有精度的，一个是没有精度的。可以从上面的例子可以看出来的。没有sym，系统会自动把按照精度来写。下面来谈谈
33、vpa 和 digits 的区别；首先，这两个函数的表达形式，Vpa(a,10)或者 vpa(a),这是临时定义的。我们会定义一个数据的精度digits（10），这个 digits 所定义的数据精度是对所有的函数有用的，但是，当出现vpa(10)时，精度会临时的变化，会遵循这个临时的vpa()所带的精度。但是在经过 vpa 语句过后，再用vpa（a），这个 a 表示的是算式 ，这个精度是前面digits()定义的精度。28，符号微积分名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 25 页，共 25 页 - - - - - - - - -