第3章数值运算.pptVIP

第1页，共31页，星期日，2025年，2月5日3.1矩阵的构造通过直接输入矩阵的元素构造矩阵：用中括号[]把所有矩阵元素括起来同一行的不同数据元素之间用空格或逗号间隔用分号（；）指定一行结束可分成几行进行输入，用回车符代替分号数据元素可以是表达式，系统将自动计算结果*第2页，共31页，星期日，2025年，2月5日例：输入矩阵A、B的值A=[1234;5678;9101112;13141516]B=[1,sqrt(25),9,13;2,6,10,7*23+sin(pi),7,11,15;4,abs(-8),12,16]*第3页，共31页，星期日，2025年，2月5日3.2矩阵下标与子矩阵提取A(m,n) 提取第m行，第n列元素A(:,n) 提取第n列元素A(m,:) 提取第m行元素A(m1:m2,n1:n2) 提取第m1行到第m2行和第n1列到第n2列的所有元素A(m:end,n) 提取从第m行到最末行和第n列的子块A(:) 得到一个长列矢量，该矢量的元素按矩阵的列进行排列*第4页，共31页，星期日，2025年，2月5日例使用下标得到矩阵元素和子块X=[1230;5608;901112;0141516]X(2,3)％提取矩阵第二行第三列的元素X(2,:)％提取矩阵第二行所有元素X(2:3,1:3)％提取第二行到第三行和第一列到第三列的子块X(2:end,1)％提取从第二行到末行和第一列的子块X(:)％得到一个长列矢量，该矢量的元素按X的列排列*第5页，共31页，星期日，2025年，2月5日例：修改矩阵A中元素的数值A=[1234;5678;9101112;13141516]A(1,1)=0A(2,2)=A(1,2)+A(2,1)A(4,4)=cos(0)*第6页，共31页，星期日，2025年，2月5日3.3矩阵的算术运算矩阵的加减运算：＋(加)、－(减)运算规则：相加、减的两矩阵必须有相同的行和列；两矩阵对应元素相加减。允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。矩阵乘法：*(乘)运算规则：A矩阵的列数必须等于B矩阵的行数。标量可与任何矩阵相乘。*第7页，共31页，星期日，2025年，2月5日矩阵除法矩阵除法：/(右除)、\(左除) X=A\B是方程A*X=B的解；X=B/A是方程X*A=B的解。*第8页，共31页，星期日，2025年，2月5日a^p——a自乘p次*方阵1的整数如a,p都是矩阵，a^p则无意义。矩阵的乘方：^(乘方)a=[1,2,3;4,5,6;7,8,9];a^2ans=303642668196102126150※当一个方阵有复数特征值或负实特征值时，非整数幂是复数阵。第9页，共31页，星期日，2025年，2月5日对于p的其它值,计算涉及特征值和特征向量。例：b=a^0.5b=0.4498+0.7623i0.5526+0.2068i0.6555-0.3487i1.0185+0.0842i1.2515+0.0228i1.4844-0.0385i1.5873-0.5940i1.9503-0.1611i2.3134+0.2717i*第10页，共31页，星期日，2025年，2月5日矩阵转置运算规则：如果矩阵A是实数阵，则矩阵A的转置是将第i行第j列的元素与第j行第i列的元素互换；如果矩阵A是复数阵，则A转置后的元素由A对应元素的共轭复数构成。X=[123;560;9011];XY=[2309;5601];YZ=[3+7i9-2i;5+2i1];Z*第11页，共31页，星期日，2025年，2月5日3.4矩阵的关系运算关系运算符：(小于)、=(小于或等于)、(大于)=(大于或等于)、==(等于)、～=(不等于)关系运算符的运算法则：矩阵的维数一致。将两个矩阵的对应元素进行比较关系运算结果。若比较关系的结果为真，置为1，否则置0*第12页，共31页，星期日，2025年，2月5日3.5矩阵的逻辑运算必须是两个同维矩阵或其中一个矩阵为标量才能进行MATLA