MATLAB的数值计算 (2).ppt
《MATLAB的数值计算 (2).ppt》由会员分享,可在线阅读,更多相关《MATLAB的数值计算 (2).ppt(74页珍藏版)》请在三一办公上搜索。
1、第二讲 MATLAB的数值计算,matlab 具有出色的数值计算能力,占据世界上数值计算软件的主导地位,数值运算的功能,创建矩阵矩阵运算线性方程组多项式函数线性插值数值分析和傅立叶变换函数优化微分方程的数值解,创建矩阵方法,直接输入法规则:矩阵元素必须用 括住 矩阵元素必须用逗号或空格分隔 在 内矩阵的行与行之间必须 用分号分隔,a=1;b=2;c=3;x=5 b c;a*b a+c c/b x=5.000 2.000 3.000 2.000 4.000 1.500y=2,4,5;3 6 8 y=2 4 5 3 6 8,例:,矩阵元素可以是任何matlab表达式,可以是实数,也可以是复数,复数
2、可用特殊函数i,j 输入。大的矩阵可以用分行输入,回车键代表分号。a=1 2 3;4 5 6 x=2 pi/2;sqrt(3)3+5i,注1 矩阵元素,注2 符号的作用,1)逗号和分号的作用 逗号和分号可作为指令间的分隔符,matlab允许多条语句在同一行出现。分号如果出现在指令后,屏幕上将不显示结果。,2)冒号的作用 用于生成等间隔的向量,默认间隔为1。用于选出矩阵指定行、列及元素。循环语句,注3 只要是赋过值的变量,不管是否在屏幕上显示过,都存储在工作空间中,以后可随时显示或调用。变量名尽可能不要重复,否则会覆盖。注4 当一个指令或矩阵太长时,可用续行,2.用matlab函数创建矩阵,空阵
3、 matlab允许输入空阵,当一项操作无结果时,返回空阵。rand 随机矩阵(01)Randn正态分布的随机矩阵eye 单位矩阵zeros 全部元素都为0的矩阵ones 全部元素都为1的矩阵diag 产生对角矩阵magic魔方阵,创建矩阵方法,例:eye(2,3)zeros(2,3)ans=ans=1 0 0 0 0 0 0 1 0 0 0 0 ones(2,3)ans=1 1 1 1 1 1 V=5 7 2;A=diag(V)A=5 0 0 0 7 0 0 0 2,eye(2)ans=1 0 0 1 zeros(2)ans=0 0 0 0 ones(2)ans=1 1 1 1,在区间20,5
4、0内均匀分布的5阶随机矩阵。命令如下:x=20+(50-20)*rand(5),reshape(A,m,n),它在矩阵总元素保持不变的前提下,将矩阵A重新排成mn的二维矩阵。linspace函数产生行向量。其调用格式为:linspace(a,b,n)其中a和b是生成向量的第一个和最后一个元素,n是元素总数。例 a=linspace(1,10,10)a=1 2 3 4 5 6 7 8 9 10注意:matlab严格区分大小写字母,因此a与A是两个不同的变量。matlab函数名必须小写。,3.用m文件创建矩阵,对于比较大且比较复杂的矩阵,可以为它专门建立一个M文件。下面通过一个简单例子来说明如何利
5、用M文件创建矩阵。,创建矩阵方法,(1)启动有关编辑程序或Matlab文本编辑器,并输入待建矩阵。(2)把输入的内容以纯文本方式存盘(设文件名为mymatrix.m)。(3)在Matlab命令窗口中输入mymatrix,即运行该M文件,就会自动建立一个名为MYMAT的矩阵,可供以后使用。,例 利用M文件建立MYMAT矩阵。,4.用冒号表达式创建矩阵,利用冒号表达式可以线性等间距地建立一个向量来创建矩阵。一般格式是:e1:e2:e3其中e1为初始值,e2为步长,e3为终止值。或者为:(start:step:end)例 a=1:2:10 a=1 3 5 7 9,创建矩阵方法,1.元素的修改(1)直
6、接修改 可用键找到所要修改的矩阵,用键移动到要修改的矩阵元素上即可修改(2)指令修改 可以用A(,)=来修改。,创建矩阵修改,例如 a=1 2 0;3 0 5;7 8 9 a=1 2 0 3 0 5 7 8 9 a(3,3)=0 a=1 2 0 3 0 5 7 8 0,2.矩阵结构的改变,变维 a=1:12;b=reshape(a,3,4)c=zeros(3,4);c(:)=a(:)变向 rot90:旋转;fliplr:左右;flipud:上下;transpose:转置抽取 diag:抽取主对角线;tril:抽取下三角;triu:抽取上三角,合并,如:A=magic(4)A=16 2 3 13
7、 5 11 10 8 9 7 6 12 4 14 15 1 A(2,:)=A=16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1,行列的删除 删除矩阵的某一行或者某一列,只要把该行或者该列赋予一个空矩阵即可。,A(2,3)A(1:k,4)A(:,4)A(:,end)A=1:10B=A(1:3:10)A=1 2 3 4 5 6 7 8 9 10B=1 4 7 10,创建矩阵元素的访问,创建矩阵信息的获取,尺寸信息、元素的数据类型、数据结构,把Matlab工作空间中一些有用的数据长久保存下来的方法是生成mat数据文件。save 将工作空间中所有的变量存到matlab.ma
8、t文件中。,创建矩阵数据的保存与获取,默认文件名,save data将工作空间中所有的变量存到data.mat文件中。save data a b 将工作空间中a和b变量存到data.mat文件中。下次运行Matlab时即可用load指令调用已生成的mat文件。,load load data load data a b mat文件是标准的二进制文件,还可以ASCII码形式保存。,即可恢复保存过的所有变量,加、减(,)相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。,矩阵运算算术运算,2.乘()规则:A矩阵的列数必须等于
9、B矩阵的行数标量可与任何矩阵相乘。a=1 2 3;4 5 6;7 8 0;b=1;2;3;c=a*bc=14 32 23,d=-1;0;2;f=pi*df=-3.1416 0 6.2832 矩阵除的运算在线性代数中没有,有矩阵逆的运算,在matlab中有两种矩阵除运算。,3.两种除法,/和,分别表示左除和右除。如果A矩阵是非奇异方阵,则AB和A/B运算可以实现。AB等效于A的逆右乘B矩阵,而A/B等效于A矩阵的逆左乘B矩阵。对于矩阵来说,左除和右除表示两种不同的除数矩阵和被除数矩阵的关系。对于矩阵运算,一般ABA/B。,A B A自乘B次幂,方阵,1的整数,4.乘方,对于B的其它值,计算将涉及
10、特征值和特征向量,如果B是矩阵,A是标量AB使用特征值和特征向量自乘到B次幂;如A,B都是矩阵,AB则无意义。,a=1,2,3;4,5,6;7,8,9;a2 ans=30 36 42 66 81 96 102 126 150,当一个方阵有复数特征值或负实特征值时,非整数幂是复数阵。,a0.5 ans=0.4498+0.7623i 0.5526+0.2068i 0.6555-0.3487i 1.0185+0.0842i 1.2515+0.0228i 1.4844-0.0385i 1.5873-0.5940i 1.9503-0.1611i 2.3134+0.2717i,5.范数,norm(V)或n
11、orm(V,2):计算矩阵V的 2范数。norm(V,1):计算矩阵V的1范数。norm(V,inf):计算矩阵V的范数。normest:估计矩阵的2范数。,inv 求逆矩阵 trace迹det 行列式的值 rank秩eig 特征值 null化零矩阵diag 对角矩阵 rref约化行阶梯 转置 sqrtm 开方 expm指数 funm一般非线性运算logm对数,6.其它运算,数组运算指元素对元素的算术运算,与通常意义上的由符号表示的线性代数矩阵运算不同。加减(.+,.-)a.+b a.-b,矩阵运算数组运算,对应元素相加减(与矩阵加减等效),2.乘除(,./,.)ab a,b两数组必须有相同的
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB的数值计算 2 MATLAB 数值 计算
链接地址:https://www.31ppt.com/p-6511997.html