MATLAB入门与矩阵运算讲义.ppt

四、矩阵的分解 对矩阵X进行QR分解，就是把X分解为一个正交矩阵Q和 一个上三角矩阵R的乘积形式。QR分解只能对方阵进行。 调用格式： [Q,R]=qr(X)：产生一个一个正交矩阵Q和一个上三角矩阵R，使之满足X=QR。? [Q,R,E]=qr(X)：产生一个一个正交矩阵Q、一个上三角矩阵R以及一个置换矩阵E，使之满足XE=QR。 满秩矩阵的QR分解 如果矩阵X是对称正定的，则Cholesky分解将矩阵X分解成一个下三角矩阵和上三角矩阵的乘积。设上三角矩阵为R，则下三角矩阵为其转置，即X=RR。 调用格式： R=chol(X)：产生一个上三角阵R，使RR=X。若X为非对称正定，则输出一个出错信息。? [R,p]=chol(X)：这个命令格式将不输出出错信息。当X为对称正定的，则p=0，R与上述格式得到的结果相同；否则p为一个正整数。 对称正定矩阵的Cholesky分解 五、解方程组 第一步:定义变量syms x y z ...; 第二步:求[x,y,z,...]=solve(eqn1,...,eqnN,var1,...varN); 第三步:求出n位有效数字的数值x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);... syms?x?y;? [x,y]=solve(x^2+3*y+1=0,y^2+4*x+1=0);? x=vpa(x,4); y=vpa(y,4); 例：解二（多）元二（高）次方程组：? x^2+3*y+1=0? y^2+4*x+1=0 解多元多次方程组 2、线性方程组的线性求解 考虑给出的线性方程组Ax=B，k=rank(A),r=rank(A B) （1）当k=r=n，方程组有唯一解,用matlab语言可以立即得出该方程组的解X=inv(A)*B或X=A^(-1)*B。 (2) 当k=rn,方程组有无穷多组解，此时可以在matlab语言中由null(A)先求出化零空间Z,其调用格式为Z=null(sym(A))，通过MATLAB语言xo=pinv(A)*B求出方程的一个特解。即 Z= null(sym(A)) x0=sym(pinv(A))*B syms a1 a2 ；x=Z*[a1;a2]+x0 MATLAB入门与矩阵计算 * 低级语言包括机器语言和汇编语言。机器语言就是计算机指令的集合，它与计算机同时诞生，是第一代的计算机语言；汇编语言是用符号来表示计算机指令，被称为第二代语言 * 低级语言包括机器语言和汇编语言。机器语言就是计算机指令的集合，它与计算机同时诞生，是第一代的计算机语言；汇编语言是用符号来表示计算机指令，被称为第二代语言 * 低级语言包括机器语言和汇编语言。机器语言就是计算机指令的集合，它与计算机同时诞生，是第一代的计算机语言；汇编语言是用符号来表示计算机指令，被称为第二代语言 * 低级语言包括机器语言和汇编语言。机器语言就是计算机指令的集合，它与计算机同时诞生，是第一代的计算机语言；汇编语言是用符号来表示计算机指令，被称为第二代语言 * 低级语言包括机器语言和汇编语言。机器语言就是计算机指令的集合，它与计算机同时诞生，是第一代的计算机语言；汇编语言是用符号来表示计算机指令，被称为第二代语言 MATLAB入门与矩阵计算 制作人：慕黎明、耿寒 MATLAB 基本内容 Matlab 的特点与功能 Matlab 具有很强的数值计算功能 Matlab 以矩阵作为数据操作的基本单位， 但无需预先指定矩阵维数（动态定维） 按照 IEEE 的数值计算标准进行计算 提供十分丰富的数值计算函数，方便计算，提高效率 Matlab 命令与数学中的符号、公式非常接近， 可读性强，容易掌握 Matlab 是一个交互式软件系统 输入一条命令，立即就可以得出该命令的结果 Matlab 的特点与功能 Matlab 符号计算功能 Matlab 和著名的符号计算语言 Maple 相结合 Matlab 的编程功能 Matlab具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。通过 Matlab 进行编程完成特定的任务 Matlab 的绘图功能 Matlab提供丰富的绘图命令，很方便实现数据的可视化 Desktop操作桌面简介 命令窗口 工作区窗口 命令历史窗口 菜单栏 工具栏 文件目录浏览器 文件概括 新建 打开 剪切 粘贴 恢复 Help：打开 MATLAB 帮助； Current Directory：设置当前目录 MATLAB 的工具栏 复制 撤消 变量命名原则 Matla