大样本数据的处理方.docVIP

第一章 大样本数据的处理方法 实验1.1 矩阵的基本运算 1.1.1实验目的 （1）熟练掌握矩阵的基本运算并能利用索引向量进行矩阵的合成与分解. （2）灵活运用sort命令，解决有关排序的实际问题. （3）熟练掌握矩阵的特征值与特征向量的计算方法. 1.1.2 实验背景知识介绍 在线性代数课程中我们学习了矩阵的加减法、数乘、矩阵的乘法、矩阵的秩以及矩阵的转置与求逆运算.为了运用Matlab 软件进行矩阵的基本运算，我们首先回顾以下基本概念： 1．矩阵的基本运算 设有矩阵，，则有以下的基本运算： (1) 矩阵的加减法 (2) 矩阵的数乘 ， (3) 矩阵的转置 (4) 矩阵的乘法 (5) 矩阵的点乘与点除 ， (6) 矩阵可逆的充分必要条件 阶矩阵可逆 2．矩阵的特征值与特征向量 设有阶矩阵和常数，如果存在非零向量，使得 则称是的一个特征值，非零向量称为矩阵对应于特征值的一个特征向量. 特征值与特征向量具有以下性质： (1) 有相同的特征值； (2) 不同的特征值对应的特征向量线性无关； (3) 阶矩阵的个特征值满足： ， 其中称为矩阵的迹. 1.1.3实验内容 1．矩阵运算的Matlab实现 在Matlab中，任何矩阵（向量）采取按行输入每个元素：同一行中的元素用逗号或者用空格符来分隔，且空格个数不限；换行时用分号进行分隔，所有元素处于方括号[ ]内. 例如，矩阵，在Matlab中的输入方法为： A=[1,5,1,0,1;2,6,0,1,1;3,7,1,0,1;4,8,0,1,1]; % 此处分号的作用在于上述输入的结果并不显示，如果要显示结果，可以去掉分号或者将分号改成逗号. 如果要从矩阵中取出部分元素组成新的矩阵，可以有以下方法： (1) 取出矩阵的第行与第列： A(i,:) , A(:,j) , % 此处冒号的作用表示取出行（列）的所有元素； (2) 取出矩阵的第i,j,k行元素组成的矩阵： 方法1索引向量法： 首先输入索引向量a=[i,j,k];然后输入 B=A(a,:), 方法2 矩阵拼接法： B=[A(i,:);A(j,:);A(k,:)] 为了方便实现矩阵的基本运算，我们将Matlab命令列于下表： 表1.1 矩阵的基本运算 命令 功能 命令 功能 A+B 矩阵对应元素相加 std(A) 计算A的各列标准差 A-B 矩阵对应元素相减 range(A) 计算矩阵A的各列极差 A.*B 矩阵对应元素相乘 rank(A) 计算矩阵A的秩 A./B 矩阵对应元素相除 A.^n 矩阵A的各元素n次方 A*B 矩阵A左乘矩阵B A^n 方阵A的n次方 k*A 常数k乘矩阵A的各元素 flipud(A) 将矩阵A的元素上下翻转 A 矩阵A转置运算 fliplr(A) 将矩阵A的元素左右翻转 inv(A) 矩阵A的求逆运算 exp(A) 将矩阵A的各元素做指数 mean(A) 计算A的各列均值 log(A) 对矩阵A的各元素取对数 var(A) 计算A的各列方差 sqrt(A) 对A的各元素求算术平方根 abs(A) 对A的各元素取绝对值 sum(A) 计算A的各列元素之和 【例1.1】设矩阵 (1) 取出矩阵的第一行与第四列； (2) 从矩阵中去掉一、三两列后剩余元素组成的矩阵； (3) 计算矩阵各行的均值与方差； (4) 计算矩阵各列的极差与标准差； (5) 将矩阵进行左右翻转和上下翻转. 题目分析：从矩阵中去掉一、三两列后剩余元素组成的矩阵也就是从中取出第2，4，5列元素组成的矩阵，由于Matlab中的命令多数是针对列元素进行计算的，所以对矩阵的行元素进行计算只需先将矩阵转置即可. 计算程序： A=[1,5,1,0,1;2,6,0,1,1;3,7,1,0,1;4,8,0,1,1]; A1=A(1,:), A2=A(:,4) 方法1：B=[A(:,2),A(:,4:5)] 方法2：c=[2,4,5];B=A(:,c) mean (A),var(A), range (A), std (A) flipud (A), fliplr (A) 结果说明：在利用拼接法时要注意行拼接与列拼接之间的不同之处 【例1.2】长江流域主要站点水质检测数据如下（2003-06—2003-10），将数据用矩阵表示，并将各站点水质状况用17个矩阵分别表示. 表1.2 长江流域主要站点水质检测数据(2003-06—2003-10） 序号 点位名称 断面情况 主要监测项目(单位:mg/L) pH* DO CODMn NH3-N 1 四川攀枝花 干流 7.6? 6.8? 0.2? 0.1? 2 重庆朱沱? 干流（川-渝省界）? 7.63? 8.41? 2.8? 0.34? 3 湖北宜