第3章 数值计算与分析.pptVIP

4. 求标准差和求方差 函数std用于求向量或矩阵中元素的标准差，函数var用于求向量或矩阵中元素的方差。函数std的调用格式如下： · std(X) 或std(X,0) 若输入参数X为向量，返回向量元素的标准差；若为矩阵，则返回一行向量，包括矩阵中每列元素的标准差；其中标准差用n-1正态化，n是数据点数； · std(X,1) 若输入参数X为向量，返回向量元素的标准差；若为矩阵，则返回一行向量，包括矩阵中每列元素的标准差；其中标准差用n正态化； · std(X,Flag,Dim) 返回维数为Dim的标准差，其中，当Flag=0时，标准差用n-1正态化；当Flag不等于0时，标准差用n正态化。 【例30】 a_std=std(a) a_std = 4.0415 2.0817 2.6458 3.2146 a_std2=std(a,1) a_std2 = 3.2998 1.6997 2.1602 2.6247 a_std3=std(a,0,1) a_std3 = 4.0415 2.0817 2.6458 3.2146 a_std4=std(a,0,2) a_std4 = 2.2174 4.5735 2.2174 a_var=var(a) % 求矩阵中每列元素的方差 a_var = 16.3333 4.3333 7.0000 10.3333 3.6.2 相关关系分析函数 在MATLAB中，提供了计算相关系数的函数corrcoef和计算协方差矩阵的函数cov，下面分别说明它们的使用。 1. 相关系数计算 函数corrcoef用于计算相关系数，它的常用调用格式如下： · corrcoef(X) 返回从矩阵X计算所得的相关系数矩阵，其大小与矩阵X一样。矩阵X的每行被看成是一次观测的量，每列对应一个变量。 · corrcoef(X,Y) 输入参数X和Y是列向量，该格式与corrcoef([X Y])的作用一样。 【例31】求相关系数示例。 x = randn(30,4); % 不相关数据 x(:,4) = sum(x,2); % 引入相关关系 r = corrcoef(x) % 计算相关系数 r = 1.0000 -0.0352 0.2673 0.6901 -0.0352 1.0000 -0.5101 0.2617 0.2673 -0.5101 1.0000 0.3504 0.6901 0.2617 0.3504 1.0000 2. 协方差矩阵计算 函数cov用于计算协方差矩阵，它的常用调用格式如下： · cov(X)或cov(X,0) 若输入参数X为向量，返回向量的方差；若为矩阵，其每行被看成一次观测的量，每列作为一个变量，计算返回协方差矩阵，其大小与矩阵X一样；协方差用n-1正态化。 · cov(X,1) 除了协方差用n正态化外，其他同上。 · cov(X,Y)或cov(X,Y,0) X和Y是具有相同元素数的矩阵，等同于cov([X(:),Y(:)])。其中，协方差用n-1正态化。 · cov(X,Y,1) 除了协方差用n正态化外，其他同上。 bb_cov= 0.3678 -0.3908 0.3115 -0.3908 1.0171 0.1093 0.3115 0.1093 0.7376 【例32】向量和矩阵的协方差计算示例。 aa=[2 4 6 8 9]; aa_cov=cov(aa) % 向量的协方差计算 aa_cov = 8.2000 bb=randn(7,3); bb_cov=cov(bb) % 矩阵的协方差计算 3.7 小 结 MATLAB具有强大的科学计算和数值分析功能。科学计算与数值分析是通过数学表达式来实现的，根据所使用的运算符是算术运算符、关系运算符还是逻辑运算符，表达式分为算术表达式、关系表达式和逻辑表达式，同时应注意这些运算符的优先级。 基本矩阵运算包括矩阵的四则运算、矩阵的幂运算、指数运算、对数运算、开方运算以及矩阵求逆、行列式运算等，矩阵运算必须要服从矩阵运算规则。MATLAB也提供了一些的矩阵基本运算函数。 基本数组运算与矩阵运算不同，它是对数组的