3_matlab处理数据.pptVIP

Matlab处理数据 天津理工大学理学院 赵瑞 matlab最大的特色可以说是它的矩阵操作了，它提供的相关命 令可以很方便的解决与矩阵有关的数据处理，我们在学习的时候就 可以体会到它的好处. 我们将重点学习一维数值数组和二维数值数组的创建访问和操作, 并将结合中国大学生数学建模比赛的实例加以练习，避免纸上谈 兵。 3.1 一维数组的创建和寻访 一、一维数组的创建 1.一维数组的创建 （1）逐个元素输入法 向量的元素用中括号[]括起来，如果是行向量，则元素之间用逗号或者是空 格隔开；如果是列向量，则元素之间用分号隔开（分号表示换行）。 例如: x=(1,2,3,4,5,6,7,8)T x=[1;2;3;4;5;6;7;8] ;% 一维 8x1 数组 size(x) %查询x的维数 ans= 8 1 二、一维数组的子数组寻访和赋值 访问（输出）数组中某个位置的元素 访问（输出）数组中某些位置的元素 修改某个元素 删除某个元素 添加新元素 例：x=8:-1:1 x = 8 7 6 5 4 3 2 1 三、一维数组的操作 求数组中满足某些条件的元素所在位置或值 求数组的最大值 最小值 平均值 排序 求和 3.2 二维数组（矩阵）的创建和寻访 一、二维数组（矩阵）的创建 二维数组输入必须有以下三个要素： 1.整个输入必须以方括号“[]”把所有元素括起来 2.二维数组的行与行之间必须用分号“;”分隔或者用“回车”分隔 3.二维数组同行元素必须由逗号“,”或空格分隔 【例】输入矩阵 二、二维数组的寻访和赋值 1.访问格式：A(r,c) %输出A的第r行第c列的元素 例： A(2,1) %输出A的第2行第1列的元素 ans= 4 思考：把A的1、3行组成一个新矩阵B， 命令：B=A([1 3],:) 把A的各行改变顺序，按照3、2、1行得到一个新的矩阵 命令：A=A(end:-1:1,:) 删除A的第1、2行， 命令：A([1 2],:)=[] 注意：对于列的操作是一样的道理 下面看一个实例： 附件2 数据 1300多名病人按照4种疗法服药大约每隔8周测试的CD4浓度。 第1列是病人编号，第2列是4种疗法的代码： 第3列是病人年龄，第4列是测试CD4的时刻（周）， 第5列是测得的CD4，取值log(CD4+1). 题目是通过这些数据得出不同年龄段的人应用不同疗法的疗效。 那么首先应该处理这些数据，现阶段我们主要练习用matlab处理 数据。 1.导入数据 首先我们应该把这些数据构成的一个大的矩阵导入matlab，方法 有二： (1)把文本文档(.txt文件)中的矩阵导入matlab的workspace (注意,文本文档中只能是一个矩阵的形式,不能含有其他的文字符号) 首先把文档所在处设为当前的路径，以便matlab能够找到,最好给这个文本重命名，最好是英文名， 比如说我们把“附件2.txt”，命名为data.txt，设置好路径，在command window中输入 load ‘data.txt’ 回车之后就会发现在workspace中出现了一个变量：data 这就完成了文本文档中矩阵的导入 注意:格式是load+空格+单引号+文档名.txt+单引号 (2)File→Import Data 出现以下窗口，选择好文件“打开”就可以了 选择“Next” 选择“Finish”，workspace中就出现了data这个变量 下面双击“data”,就可以看到它的内容 （2）操作练习 导入数据以后，我们做一些矩阵操作的练习 1）把所有使用第一种疗法的病人的数据找出来，生成的新矩阵命名为 p1，并查看这个矩阵的行列数，把p1里cd4浓度出现了0的所在行删除。 命令： place=find(data(:,2)==1); %找出data中第二列的等于1的元素所在行位置 p1=data(place,:); % 把data中第二列的等于1的元素所在行提取出来组成一个新的矩阵p1 p1_size=size(p1); %查询p1的行列数 zeroplace=find(p1(:,5)==0); %找出p1第五列中等于0的元素的位置 p1(zeroplace,:)=[]; % 把p1第五列中等于0的元素所在行删除 下面双击“p1”,就可以看到它的内容 2) 把p1里的病人数据按年龄从大到小的顺序排列 命令： [m,n]=sort(p1(:,3)); %把p1的第三列的元素按从小到大的顺序排列,n返回的是 %排列之后各元素原来的位置 p1=p1(n,:); %把p1各行的顺序改变了，按照第三列的序改变的 %目前的序是从小到