《多元统计分析作业.doc》由会员分享，可在线阅读，更多相关《多元统计分析作业.doc（41页珍藏版）》请在金锄头文库上搜索。
1、重 庆 交 通 大 学学 生 实 验 报 告实验课程名称 MATLAB应用统计实验 开课实验室 理学院实验室 学 院 2011 年级 数学与应用数学 专业班 1班 学 生 姓 名 xxx 学 号 xxxxxx 开 课 时 间 2013 至 2014 学年 第 2 学期假设合理优良中差建模求解全面优良中差结果分析完善优良中差文档清晰优良中差综合成绩教师姓名董志清实验一 Matlab软件的简单熟悉1实验目的：1实验内容：1实验要求：2实验过程：2实验体会：5实验二 统计变量的描述性分析6实验目的：6实验内容：6实验过程：7实验体会：12实验三 聚类分析13实验目的14实验内容15实验过程15实验结
2、果17实验结果分析18实验体会19附录19实验四 判别分析20实验目的20实验内容22实验过程23实验结果23实验结果分析24实验体会24实验五 主成分分析25实验目的25实验内容25实验过程26实验结果27实验结果分析27实验体会28实验六 因子分析29实验目的29实验内容29实验过程30实验结果35实验结果分析36实验体会36实验七 典型相关分析37实验目的37实验内容38实验过程39实验结果39实验结果分析40实验体会 41实验一 MATLAB软件的简单熟悉一、 实验目的熟悉在 MATLAB 中数据的录入、保存和调用方法；熟悉 matlab 中关于矩阵运算和函数运算的各种命令。二、 实验
3、内容1、 矩阵和数组的输入直接输入：MATLAB 中不用描述矩阵的类型和维数，它们由输入的格式和内容决定。小矩阵可以用排列各个元素的方法输入，同一行元素用逗号或空格分开，不同行的元素用分号或回车分开。如 A=1,2,3;4,5,6 （表示键入， 后接回车，然后屏幕上显示输出） A=1 2 3;4 5 6 A=1 2 34 5 6都输入了一个 2*3 的矩阵 A，屏幕上显示结果相同。数组的输入除了向矩阵一样输入以外，还常采用“：”输入方式，如 t=1:5 （从 1 到 5 公差为 1 的等差数组）t =1 2 3 4 5 c=1:2:7（从 1 到 7 公差为 2 的等差数组）c =1 3 5
4、7 linspace(0,1,9) （从 0 到 1 共 9 个等差数组）ans = 0 0.1250 0.2500 0.3750 0.5000 0.6250 0.7500 0.8750 1.0000由语句或函数生成；2 、矩阵的运算MATLAB 中提供了以下矩阵运算符：+ 加法；- 减法； 转置；* 乘法； 乘幂； 左除；/ 右除它们要符合矩阵的运算规律，如果矩阵的行列数不符合运算符的要求将产生错误信息。 a=1 11 3 ;4 12 6 A+a （MATLAB 区分大小写）ans =2 13 68 17 12不指定输出变量，MATLAB 将回应 ans（answer 缩写） B=A-aB
5、=0 -9 00 -7 0A-a 将结果保存在 B 变量中，变量可以随时调用，除非被清除或替代。 左除；/ 右除用法如下：设 A 时可逆矩阵，AX=B 的解是 A 左除 B，即X=AB(当 B 是列向量时， 得到方程出的解)； XA=B 的解是 A 右除 B， 即 X=B/A. 如 A=2 5 -3;3 6 -2;2 4 -3 b=3 1 4 X=AbX =1-1-23 、矩阵的裁剪与拼接从一个矩阵中取出若干行（列）构成新的矩阵称为裁剪，MATLAB 中“：”是非常重要的裁剪工具，如 A(3,:) （A 的第三行） A(:,2) （A 的第二列）将几个矩阵接在一起，称为拼接，左右拼接行数要相同
6、，上下拼接列数要相同，如 E=A,bE =2 5 -3 33 6 -2 12 4 -3 4 F=b;X4、 变量与函数MATLAB 的变量由字母、数字和下划线组成，最多 31 个字符，区分大小写字母，第一个字符必须是字母。当输入一个新变量时，MATLAB 自动建立变量并为其分配内存空间。几个特殊的量：pi 圆周率；Inf 正无穷大指 1/0；NaN 不定数，指 0/0。常用的函数： 三角函数 sin cos tan cot sec csc;其它基本函数 sqrt; exp; log; log10;abs(绝对值)；round( 四舍五入)。这些函数本质上是作用于标量的，当它们作用于矩阵时，使作
7、用于矩阵的每一个元素。例如 x=linspace(0,2*pi,30); （如过语句一分号“；”结束，Matlab 只进行计算，不显示计算的结果。 y=sin(x); plot(x,y)plot 是最常用的画 2 维图形的命令向量函数：max; min; sum; mean(平均值); sort（从小到大排序）这些函数只有作用于向量时才有意义 a=5 2 1 4 3 b=min(a),c=sum(a),e=sort(a)矩阵函数：构造函数 zeros(0 阵)；ones(1 阵);eye(单位阵);rand(随机阵) w=zeros(2,3) u=ones(3) v=eye(3,4) x=ra
8、nd(1,3)矩阵计算函数：求秩 rank(A)；求逆 inv(A)；求迹：trace(A)；求行列式：det(A)；求特征值 eig(A) a=1 11 3 ;4 12 6 ;7 2 9 d=det(a),r=rank(a),t=trace(a),e=eig(a)5、 命令和环境窗口在线帮助可以用命令 help 主题名例如： help sum显示在当前工作区中的所有变量名： whos清除当前工作区中的所有变量： clear把变量储存在文件中： save 文件名调出文件中的变量： load 文件名 save magik clear load magik三、实验要求了解matlab的基本操作。四
9、、实验过程 a=1 2 3;4 5 6;7 8 9a = 1 2 3 4 5 6 7 8 9 a=1,2,3;4,5,6;7,8,9;a = 1 2 3 4 5 6 7 8 9 p=1:0.5:4p = 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000 A=1 3 5;7 9 2;4 6 8A = 1 3 5 7 9 2 4 6 8 a+Aans = 2 5 8 11 14 8 11 14 17 B=a-AB = 0 -1 -2 -3 -4 4 3 2 1 C=2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5; a1=de
10、t(C)a1 = -85 a2=eig(C)a2 = -1.1228 2.5266 5.2981 + 1.3755i 5.2981 - 1.3755i a3=inv(C)a3 = -0.0471 0.5882 -0.2706 -0.9412 0.3882 -0.3529 0.4824 0.7647 -0.2235 0.2941 -0.0353 -0.4706 -0.0353 -0.0588 0.0471 0.2941 a4=rank(C)a4 = 4 a4=rank(C)a4 = 4 a6=C.3a6 = 8 1 -27 -1 27 1 0 343 -1 8 64 -8 1 0 -1 125
11、zeros(3)ans = 0 0 0 0 0 0 0 0 0 eye(3)ans = 1 0 0 0 1 0 0 0 1 ones(3)ans = 1 1 1 1 1 1 1 1 1 rand(3)ans = 0.8147 0.9134 0.2785 0.9058 0.6324 0.5469 0.1270 0.0975 0.9575 randn(3)ans = 2.7694 0.7254 -0.2050 -1.3499 -0.0631 -0.12413.0349 0.7147 1.4897 sin(pi/2)ans = 1 cos(0)ans = 1 tan(pi/4)ans =1.0000 clear x=-pi:0.01:pi; y=sin(x); plot(x,y)图 1 clear syms x f=x*sin(1/x); ezplot(f)图 2实验二 统计数据的描述性分析一、实验目的熟悉在matlab中实现数据的统计描述方法，掌握基本统计命令：样本均值、样本中位数、样本标准差、样本方差、概率密度函数pdf、概率分布函数df、随机数
