CHNO3(数值计算).pptVIP

第三章 基本数值计算 第一节 数据类型 （1）变量。分数值变量和字符变量。区分大小写、变量名不 超过31位、不能有中文(包括文件名）。 （2）常量。计算机中不变的量。如i,j,pi,NaN（不确定）, Inf（无穷大）等。常量可任意调用。 （3）字符变量。将字符串作为变量。有三种方法表示字符串。 1）用单引号‘’。 2）用函数sym(‘’)。 3）用命令syms。注意函数与命令的用法区别。(exno9t) 第二节 矩阵构造及运算 二、矩阵的加、减、乘和除 （1）进行矩阵加、减、乘和除的条件： 1）进行矩阵加减时，参与运算的矩阵必须同维。 2）A/B时，A、B列数必须相同。A\B时， A、B行数 必须相同。 （2）注意矩阵的左乘与右乘不同。（A*B 不等于 B*A） （3）注意矩阵的左除与右除不同。（A/B不等于 A\B） （4）标量与矩阵的数运算和数学函数对矩阵的运算等于 对矩阵的每一个元素的运算。 （exno12t） 练习3-1 求解四元线性方程 定义三个矩阵：系数矩阵、未知数矩阵、常数矩阵： 上述方程可变为求解 （exno13t） 三、对矩阵的操作 1.对矩阵元素的操作 （1）提取矩阵A的第n行第m列的元素，a=A(n,m)。 A([1:3],4) A([1,3],4) （2）提取矩阵A的第n行的所有元素， b=A(n,:)。 A(b):提取逻辑阵b中为1的A阵元素 （3）提取矩阵A的第m列的所有元素， c=A(:,m)。 （4）将矩阵A的第n行第m列的元素重新赋值b，A(n,m)=b。 （5）将矩阵A的第n行的所有元素重新赋值b, A(n,:)=b。 （6）将矩阵A的第m列的所有元素重新赋值b A(:,m)=b。 （7）将矩阵A的第n行第m列的元素删除， A(n,m)=[ ]。 （8）将矩阵A的第n行的所有元素删除， A(n,:)=[ ]。 （9）将矩阵A的第m列所有元素删除， A(:,m)=[ ]。(命令窗练习) 元素的行列为与的关系 2.对矩阵A的部分元素操作(exno3r) 3.对矩阵A的部分分析 （1）方阵A的行列式值，可表达为 det(A)。 （2）矩阵A的秩，可表达为 rank(A)。 （3）行列式值不为零的方阵A，求逆矩阵可表达为inv(A)。 （4）矩阵A的转置阵，可表达为A’。 （5）矩阵A的特征向量与特征值： [V,D]=eig(A) 式中，等号左边是函数对矩阵A的作用后返回的变量信息。其中， V和D分别为A的特征向量与特征值。这里，注意MATLAB的函数发生 作用后返回的变量信息这种形式，特别注意中、小括号的使用。 四、矩阵的数组运算(matrix,array) 运算中，往往需要在矩阵之间进行对应元素的运算，这就是矩阵的数组运算。矩阵运算与矩阵的数组运算区别主要反映在乘和除两种运算中。 1.乘法 （1）矩阵运算 A*B 。（2）矩阵的数组运算 A.*B (称为数乘)。 2.除法 （1）矩阵运算 A/B或A\B。 （2）矩阵的数组运算A./B (称为数除)。 矩阵的两种不同运算方式提醒读者特别需要引起高度重视。 练习3-2 已知 A=[1 3 5;2 4 6;7 8 9] B=[1 2 1;2 3 2;3 4 3] 比较A*B和A.*B以及A/B和A./B的区别。(命令窗练习) 练习3-3 已知 x=0:0.5:8 按如下函数关系求y序列 练习3-4 exno7t 五、矩阵元素的关系运算与逻辑运算 1.矩阵元素之间的关系运算 （1）元素之间的关系运算有 《 小于 》 大于 《= 小于等于 》= 大于等于 == 等于 ~= 不等于 （2）运算结果：满足条件时为 1，不满足条件时为 0，特别注意 等于的运算符号。 2.矩阵元素之间的逻辑运算 （1