数理统计matlab,MATLAB数理统计方法

MATLAB数理统计方法

MATLAB数理统计方法 一 MATLAB基础 MATLAB的意思是Matrix laboratory，是进行科学计算的重要工具．启动MATLAB后，出现如下图所示的界面，在缺省状态，呈现3个窗口．右边的窗口(Command Window)为命令窗口；左边的两个窗口分别为启动平台(Launch Pad)和命令历史(Command history)窗口，也可以切换到工作空间(Workspace)和当前目录浏览器窗口(Current Directory)． 可以在命令窗口通过键盘输入要执行的命令并按回车键确认，也可以点击菜单“File”—“New”—“M-file”先建立.m文件(也称为m-文件或M-文件)，在该文件中依次输入要执行命令，形成一个程序，然后执行这个程序． 程序必须被储存后方可执行．可以点击菜单“Debug”—“Run”执行程序，也可以按F5键执行程序． M-文件的名字可以由英文字母a-z及A-Z，数字0-9和下划线组成，但必须以英文字母打头！否则将会产生重大错误． 输入矩阵的最简单的方法是把矩阵的元素直接排列在方括号“[]”中，每行内的元素用空格或逗号分开，行与行之间用分号分开，多个空格被视为一个空格．例如输入 a = [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 大矩阵可以分行输入，用回车键代替分号，这样的输入形式更接近我们平时使用的矩阵格式．例如 a = [1 2 3 4 5 6 7 8 9] 可以先将一个矩阵输入到Excel工作表，即先建立数据文件．然后打开“file”—“Import Data”菜单，即启动导入数据导航，按提示打开Excel工作表，将数据导入到MATLAB的工作空间中去． 矩阵也称为数组．只有一行或一列的数组称为一维数组，有多个行和多个列的数组称为二维数组．一维数组x的第i个元素记为x(i)，二维数组a的第i行第j 列的元素记为a(i,j)． MATLAB采用双精度储存变量和数值计算，但能以多种格式输出数据．例如 x = [4/3 1.2345e–6] 则在几种常见的的格式下，输出结果分别为 at short 1.3333 0.0000 at short e 1.3333e+000 1.2345e–006 at short g 1.3333 1.2345e–006 at long 1.33333333333333 0.00000123450000 at long e 1.333333333333333e+000 1.234500000000000e–006 at long g 1.33333333333333 1.2345e–006 at rat 4/3 1/810045 应当先执行格式命令，然后再输出数据．MATLAB的缺省格式为at short． 要显示一个变量的内容，只需在命令窗口或程序中键入该变量的名字．例如在命令窗口依次执行a=1/3 ，at long，a这三个命令，结果如下： >> a=1/3 a = 0.3333 >> at long >> a a = 0.33333333333333 若最大的元素大于1000或小于0.001，则显示short或long格式时会加上一个比例．在命令后加上分号“；”，则屏幕上不会立即显示出结果，这在运算大的数据量时十分有用，如下命令产生100*100的魔方矩阵，但并不在屏幕上显示． A = magic(100); 如果一个命令很长，想另起一行接着输入命令，须要在末尾加上“.”，如: s = 1 – 1/2 + 1/3 – 1/4 + 1/5 – 1/6 + 1/7 . – 1/8 + 1/9 – 1/10 + 1/11 – 1/12; 可用who或whos来察看当前工作空间中有哪些变量．若要从工作空间中删除所有的变量，用 clear 也可以根据需要一次删除若干个变量，例如要从工作空间中删除x1，x2两个变量，用 clear x1 x2 你可以将工作间保存为一个二进制的.mat文件，以备以后调用．命令 save june10 将工作空间保存到文件june10.mat．也可只保存工作间中的部分变量值，例如要保存变量x,y,z到文件june10.mat，使用命令 save june10 x y z 重载时只需输入 load june10 要想将变量a更名为b，可使用命令 b=a; clear a MATLAB支持矩阵的加法、减法、乘法、转置、求逆等各种运算．命令a 和inv(a)分别计算矩阵a的转置矩阵和逆矩阵．MATLAB的算数运算符为： + 加 - 减 * 乘 ^ 幂 / 左除 \ 右除 对于两个标量a, b来说，a/b=b\a=a÷b，而对于两个矩阵a, b来说， a/b=a*inv(b) b\a=inv(b)*a 两个同维数组相加减，等于其对应元素相加减．一个数组与一个标量相加减，则等于数组的各元素分别与这个常数相加减． 用符号“. / ”表示两个数组的除法．若x, y是同维数组，则x ./ y表示x的元素分别除以y的对应元素得到的数组． z=x ./ y 即z(i, j)=x(i, j)/y(i, j)．x ./ y 与 y .\ x 相等，都表示x除以y，但运算是在对应元素间进行的，与矩阵的除法是不同的． 矩阵的乘方用“^”符号表示，a^p的意思是a的p次方． 数组的乘方用 “.^” 符号表示．若x=[x1, x2, …, xn], y=[y1, y2, …, yn]是同维数组，则 z = x .^ y=[x1^y1, x2^y2, …, xn^yn] 若c是一个标量，则 z = x .^ c=[x1^c, x2^c, …, xn^c] 可以使用help命令寻求帮助．例如