第5章MATLAB数值计算(免费阅读).pptVIP

第5章 MATLAB数值计算 第5章 MATLAB数值计算 本章要点 数据的分析与统计 数值插值 曲线拟合 求解常微分方程 教学时数 2学时 在实际工程技术中，人们一般通过数值解来获得问题的解，主要原因： 1、解析解不存在 2、解析解存在但不实用 与符号计算相比，数值计算在科研和工程中的应用更为广泛。MATLAB也正是凭借其卓越的数值计算能力而称雄世界。 随着科研领域、工程实践的数字化进程的深入，具有数字化本质的数值计算就显得愈益重要。 5.1 数据分析 MATLAB：面向矩阵进行运算 矩阵的每列代表不同的被测变量，相应的行代表被测向量的观测值，通过对矩阵元素的访问进行数据统计分析。 一、 数据统计 函数名称 功能 函数名称 功能 max(x) 找x各列的最大元素 min(x) 找x各列的最小元素 mean(x) 求x各列的平均值 sum(x) 求x各列元素之和 median(x) 找x各列的中间值元素 sort(x) 使x的各列元素按递增排序 prod(x) 求x各列元素之积 【例】对矩阵 求各列的最大元素、中值和平均值。 clear A=[4 8 -9;11 -12 4;-8 0 5;0.6 5 10]; maxA=max(A) maxA = 11 8 10 medA=median(A) medA = 2.3000 2.5000 4.5000 meanA=mean(A) meanA = 1.9000 0.2500 2.5000 二、离差和相关 函数名称 功能 函数名称 功能 var(x) x各列的方差 cov(x，y) 两个矩阵x和y的协方差 std(x) x各列的标准差 corrcoef(x) x的自相关阵 range(x) x各列的极差 corrcoef(x，y) 两个矩阵x和y的互相关系数，结果为方阵 cov(x) x的协方差阵 corr2(x，y) 两个矩阵x和y的相关系数 【例】建立一个3×4的随机矩阵，求方差、标准差、极差、协方差和自相关阵 A=rand(3,4) A = 0.9501 0.4860 0.4565 0.4447 0.2311 0.8913 0.0185 0.6154 0.6068 0.7621 0.8214 0.7919 B=var(A) B = 0.1293 0.0429 0.1616 0.0301 C=std(A) C = 0.3596 0.2070 0.4020 0.1736 D=range(A) D = 0.7190 0.4053 0.8029 0.3472 E=cov(A) E = 0.1293 -0.0725 0.0819 -0.0293 -0.0725 0.0429 -0.0301 0.0237 0.0819 -0.0301 0.1616 0.0323 -0.0293 0.0237 0.0323 0.0301 F=corrcoef(A) F = 1.0000 -0.9731 0.5664 -0.4688 -0.9731 1.0000 -0.3614 0.6596 0.5664 -0.3614 1.0000 0.4624 -0.4688 0.6596 0.4624 1.0000 5.2 数值运算 一、 多项式 一个n次多项式有n+1个系数，在MATLAB中用一个长度为n+1的行向量表示一个n次多项式，其中多项式的各元素按降幂顺序排列，缺少的系数要补0. MATLAB表示的系数向量为 如果知道多项式的根为 则 多项式相关函数 名称 函数格式 说明 创建多项式 P=[ a0 a1 a2 …an-1 an] P为多项式（以下各函数中P均为多项式），a0 a1 a2 … an-1 an为按降幂顺序排列的多项式系数 P=poly(A) A为向量。创建以向量A中元素为根的多项式 求根 roots(P) 求该多项式的根，以列向量的形式给出 求值 polyval(P，A) 当A为标量时，求多项式P在自变量x=A时的值；当A为向量时，求x分别等于A中每个元素时