实验一熟悉MATLAB 环境
[实验目的]
1．熟悉MATLAB 主界面，并学会简单的菜单操作。
2．学会简单的矩阵输入与信号输入。
3．掌握部分绘图函数。
[实验原理]
MATLAB 是以复杂矩阵作为基本编程单元的一种程序设计语言。它提供了各种矩阵的运算与操作，并有较强的绘图功能。
用户第一次使用MATLAB 时，建议首先在屏幕上键入demo 命令，它将启动MATLAB 的演试程序，用户可在此演示程序中领略MATLAB 所提供的强大的运算与绘图功能。也可以键入help 进行进一步了解。MATLAB 启动界面如图所示：
操作界面主要的介绍如下：
指令窗（Command Window），在该窗可键入各种送给MATLAB 运作的指令、函数、表达
式，并显示除图形外的所以运算结果。
历史指令窗（Command History），该窗记录已经运行过的指令、函数、表达式；允许
用户对它们进行选择复制、重运行，以及产生M 文件。
工作空间浏览器（Workspace Browser），该窗口罗列出MATLAB 工作空间中所有的变量
名、大小、字节数；并且在该窗中，可对变量进行观察、编辑、提取和保存。
其它还有当前目录浏览器（ Current Directory Browser ）、M 文件编辑/ 调试器
（Editor/Debugger）以及帮助导航/浏览器（Help Navigator/Browser）等，但通常不随操
作界面的出现而启动。
利用File 菜单可方便对文件或窗口进行管理。其中File | New 的各子菜单，M-file（M 文件）、Figure（图形窗口）、或Model（SIMULINK 编辑界面）分别可创建对应文件或模块。Edit 菜单允许用户和Windows 的剪切板交互信息。
MATLAB 语言最基本的赋值语句结构为：变量名列表=表达式。表达式由操作符或其它字符，函数和变量名组成，表达式的结果为一个矩阵，显示在屏幕上，同时输送到一个变量中并存放于工作空间中以备调用。如果变量名和“=”省略，则ans 变量将自动建立，例如键入：1900/81，得到输出结果：ans =23.4568。MATLAB 中变量命名的原则要求必须以英文字母开头，文件夹名字中可以还包括、下划线和数字，不要使用其他的字符，更不要单纯使用数字或者中文名命名，有时在运行MATLAB 中一些莫名的错误可能就是不规范的命名引起的。这种规则包括将来为自己编写的脚本文件、函数文件命名以及为使用的变量命名也应遵循这个规则。
[实验内容]
一．用户工作目录和当前目录的建立和设置
1．为管理方便，每个用户在使用MATLAB 前，尽量为自己建立一个专门的工作目录，即“用户目录”，用来存放自己创建应用文件。例如首先打开资源管理器，在E 驱动器下可以根据自己喜好建立一个新文件夹，但应注意：该文件夹必须以英文开头，文件夹名字中可以包括字母、下划线和数字，不要使用其他的字符，更不要单纯使用数字或者中文名命名，有时在运行MATLAB 中一些莫名的错误可能就是不规范的命名引起的，这与MATLAB 中为变量和文件名命名原则是一致的。尽管MATLAB\work允许用户存放用户文件，但最好把它仅作为临时工作目录来使用。2．为了用户运行方便，在开始工作前可把用户定义的目录设置为当前目录，方法是可直接在MATLAB 命令窗口桌面上方通过目录设置栏来实现。
3．把自己的工作目录加到MATLAB 搜索路径上，这样可以在任何情况下方便地调用自己的编写文件。MATLAB 工作时，基本搜索过程为：首先在工作空间，即MATLAB 内存中进行检查，看输入的指令是不是变量；如不是，则检查输入指令是不是内建函数（比如sin 函数等）；如不是，则在当前目录上，检查是否有与输入指令相同的M文件存在；如还没有，则在MATLAB 定义的搜索路径其他目录中，检查是否有该M文件存在。设置方法是通过File 菜单下设置路径对话框进行。
二．熟悉简单的矩阵输
1．从屏幕上输入矩阵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
2．试用回车代替分号，观察输出结果
A=[1,2,3
4,5,6
7,8,9]
A =
1 2 3
4 5 6
7 8 9
3．输入size(A)，观察结果。
size(A)
ans =
3 3
4．输入矩阵B=[9，8，7；6，5，4；3，2，1]；回车
输入矩阵C=[4，5，6；7，8，9；1，2，3]；回车，然后分别键入A B C 回车观察结果。
B=[9,8,7;6,5,4;3,2,1]
B =
9 8 7
6 5 4
3 2 1
C=[4,5,6;7,8,9;1,2,3]
C =
4 5 6
7 8 9
1 2 3
5．选择File|new 菜单中的M-file，输入B=[9 ，8，7；6，5，4；3，2，1]；保存为B.M 文件，退出编辑环境。此时在工作环境中使用B 命令就可调出B 矩阵。
6．再试着输入一些矩阵，矩阵中的元素可为任意数值表达式，但注意:
矩阵中各行各列的元素个数需分别相等，否则会给出出错信息。
D=[1,2,3,4;4,5,6,7;7,8,9,4;5,2,3,1]
D =
1 2 3 4
4 5 6 7
7 8 9 4
5 2 3 1
7．输入who 和whos 观察结果，了解其作用。
who
Your variables are:
A B C D ans xqxqxq1234 xqxqxq1235
whos
Name Size Bytes Class
A 3x3 72 double array
B 3x3 72 double array
C 3x3 72 double array
D 4x4 128 double array
ans 1x2 16 double array
xqxqxq1234 1x1 1 logical array
xqxqxq1235 1x3 6 char array
Grand total is 49 elements using 367 bytes
二．常用基本命令的使用
数据文件的存取训练：利用save 和load 指令可以方便地把当前工作空间的数据变量保
存为数据文件。接题一，输入：save a 回车，表示把当前工作空间的所有变量保存为a.mat 文件。输入save a A C 回车，表示把当前内存中的变量A、C 保存为a.mat 文件。类似的还有load 命令，实现保存数据变量的装载，执行clear 命令，然后运行load a 回车，观察工作空间的变化。save、load、clear 命令还有一些其他的应用格式，具体可参看帮助文件。
save a
save a A C
三．基本序列运算
1．数组的加减乘除和乘方运算
A=[1,2,3;4,5,6;7,8,9]
A =
1 2 3
4 5 6
7 8 9
B=[7,8,9;4,5,6;1,2,3]
B =
7 8 9
4 5 6
1 2 3
A*B
ans =
18 24 30
54 69 84
90 114 138
A+B
ans =
8 10 12
8 10 12
8 10 12
A-B
ans =
-6 -6 -6
0 0 0
6 6 6
A\B
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 2.203039e-018.
ans =
-0.8004 0.7984 -0.3999
-8.3992 -12.5967 -11.2003
8.1996 10.7984 10.6001
2．输入A=[1 2 3]，B=[4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B,并用stem 画出A，B，C，D，E，F，G。再输入一些数组，进行类似运算。
A=[1 2 3]
A =
1 2 3
B=[4 5 6]
C=A+B
C =
5 7 9
B =
4 5 6
D=A-B
D =
-3 -3 -3
E=A.*B
E =
4 10 18
F=A./B
F =
0.2500 0.4000 0.5000 G=A.^B
G =
1 3
2 729
stem(A)
stem(B)
stem(C)
stem(D)
stem(E)
stem(F)
stem(G)
3． 在命令窗口用plot 指令粗略描绘下列各函数的波形（其中对于连续信号可取时间间隔为
0.001，可参看下面的a 的实例来实现，数学函数MATLAB 实现参见附录） a.
t e t f --=3)( 30< 实现方法，在命令窗口执行一下指令，可简单描绘出函数曲线： t=0:.001:3; y=3-exp(-t); plot(t,y) b. t t e e t f 235)(--+= 30< y=5*exp(-t)+3*exp(-2*t) plot(t,y) c. t e t f t π2sin )(-= 30< y=exp(-t).*sin(2*pi*t); plot(t,y) d. at at t f /sin )(= ? 2π < t < 2π （π在MATLAB 中用pi 来实现） clear t=-2*pi:0.001:2*pi; a=2; y=sin(a*t)./(a*t); plot(t,y) f. k e k f =)( 50<≤k t=0:0.001:5; k=2; y=exp(k); plot(t,y) 三．M 脚本文件编写初步 用M 文件编辑器编写脚本文件重新描绘上面函数波形.clear a b c 实验二数值数组创建、应用及可视化 [实验目的] 1．掌握二维数组的创建、寻访，区分数组运算与矩阵运算的区别。 2．掌握标准数组生成函数和数组构造技法。 3．进一步熟悉M 脚本文件编写的方法和技巧。 [实验原理] 数值数组和数组运算MATALB 核心内容，是MATLAB 最重要的一种内建数据类型。通常，数组是由一组实数或复数排成的长方阵列(Array ),它可以是一维的“行”或“列”，可以 是二维的“矩形”，也可是三维的若同维矩形的堆叠，甚至是更高的任意维。而数组运算是指无论在数组上施加什么运算(加减乘除或函数),总认定那种运算对被运算数组中的每个元 素(Element)平等地时实施同样的操作。这使得计算程序简单、易读，使程序指令更接近教 科书上的数学计算公式，并提高了程序的向量化程度，提高计算效率，节省计算机开销。 一维数组的创建可采用逐个元素输入法，这是最简单，但又最通用的构造方法，如： x=[2 pi/2 sqrt(3) 3+5i];另外有规律的产生数组可以采用冒号生成法，通用格式是x = a : inc : b，a 是数组的第一个元素，inc 是采样点之间的步长。若（b-a）是inc 的整数倍，则生成数组的最后一个元素是b,否则小于b。或者采用定数线性采样法，该法是在设定“总点数”下，均匀采样生成一维“行”数组。格式为x = linspace ( a , b , n )，a、 b 分别是生成数组的第一个各最后一个元素,n 是采样总点数，该指令生成（1×n）数组。 二维数组是由实数或复数排列成矩形而构成的。从数据结构上看，矩形和二维数组没 有什么区别。当二维数组带有线性变换含义时，该二维数组就是矩阵。二维数组的创建也可采用直接输入法，或者利用构造M 文件创建和保存数组。 除此以外，还可以采用MATLAB 提供的标准函数生成我们需要的数组，诸如zeros、ones、rand、eye、diag、magic 等，详请见附录。 一维数组元素的寻访和标识采用X(index)方法，只是要注意MATLAB 中第一个元素下 标index 是1，而不是C 语言中的0。二维数组元素的标识和寻访可分为“全下标”标识和 “单下标”标识，“全下标”标识，即指出是“第几行，第几列”的元素，如A(3,5)表示二 维数组A 的第三行第五列元素。该标识法的优点是几何概念清楚，引述简单，在MATLAB 的 寻址和赋值中最为常用。“单下标”标识，顾名思义就是只用一个下标来指明元素在数组中 的位置，首先对二维数组的所有元素进行“一维编号”。“一维编号”是指：先设想把二维数组的所有列，按先左后右的次序、首尾相接排成“一维长列”，然后自上往下对元素未知进行编号，其优点是简洁、方便，特别是如果碰到对二维数组进行诸如for 循环操作时可以减少循环次数，提高编程效率。“全下标”标识和“单下标”标识可以通过sub2ind 和ind2sub 指令进行转换，详见帮助。另外，不论二维数组还是一维数组还可以采用“逻辑1”标识， 这种方法常用于寻找数组中所有大于或小于某值的元素的问题中。比如X(abs(X)>3)可以找 出数组X 中所有绝对值大于3 的元素。另外还可以借助ones, zeros, rand, randn 和cat, repmat, reshape 等函数直接或间接构作高维数组，详见附录或MATLAB 帮助。 Matlab 中的许多函数可以直接对任意维的数组直接运算，相当于对数组中的每个元素 分别进行运算。比如Y=sin(X)可以直接得到与数组X 中每一个元素相对应的的正弦值，这 大大简化了编程。可以执行数组运算的常用函数已经列表于附录中，以供参考。但值得注意 的是，虽然从外观形状和数据结构上看，二维数组和（数学中的）矩阵）没有区别，但矩阵 作为一种变换或映射算子的体现，矩阵运算有着明确而严格的数学规则。数组运算是MATLAB 软件所定义的规则，其目的是为了数据管理方便、操作简单、指令形式自然和执行计算的有 效。为了区别数组和矩阵运算，在易混淆的地方，数组运算在运算符前加一小黑点“.”以 示区别，比如Y=A.*B，代表的示数组A 和数组B 对应元素相乘，而Y=A*B，则表示内维相同 的矩阵A 和B 的乘积。由此也可看出，在执行数组与数组的运算时，参与运算的数组必须同 维，运算所得结果也总与原数组同维。 本节只涉及数组可视化方法的简单实现。通常，对于离散数据可采用stem 命令或者使 用plot 绘点的方法，而对于连续函数可直接采用plot 指令来实现。 [实验内容] 一．数组的创建和寻访 1．一维数组 在命令窗口执行下面指令，观察输出结果，体味数组创建和寻访方法，%号后面的为注 释，不用输入。 1．一维数组 在命令窗口执行下面指令，观察输出结果，体味数组创建和寻访方法，%号后面的为注释，不用输入。 rand('state',0) %把均匀分布伪随机发生器置为0 状态 x=rand(1,5) %产生(1×5) 的均布随机数组 x(3) %寻访数组x 的第三个元素。 x([1 2 5]) %寻访数组x 的第一、二、五个元素组成的子数组。 x(1:3) %寻访前三个元素组成的子数组 x(3:end) %寻访除前2 个元素外的全部其他元素。end 是最后一个元素的下标。 x(3:-1:1) %由前三个元素倒排构成的子数组 x(find(x>0.5)) %由大于0.5 的元素构成的子数组 x([1 2 3 4 4 3 2 1]) %对元素可以重复寻访，使所得数组长度允许大于原数组。 x(3) = 0 %把上例中的第三个元素重新赋值为0 x([1 4])=[1 1] %把当前x 数组的第一、四个元素都赋值为1 x[3]=[]; %空数组的赋值操作 rand('state',0) x=rand(1,5) x = 0.9501 0.2311 0.6068 0.4860 0.8913 x(3) ans = 0.6068 x([1 2 5]) ans = 0.9501 0.2311 0.8913 x(3:end) ans = 0.6068 0.4860 0.8913 x(3:-1:1) ans = 0.6068 0.2311 0.9501 x(find(x>0.5)) ans = 0.9501 0.6068 0.8913 x(find(x>0.5)) ans = 0.9501 0.6068 0.8913 x([1 2 3 4 4 3 2 1]) ans = 0.9501 0.2311 0.6068 0.4860 0.4860 0.6068 0.2311 0.9501 x(3)=0 x = 0.9501 0.2311 0 0.4860 0.8913 x([1 4])=[1 1] x = 1.0000 0.2311 0 1.0000 0.8913 x([3])=[] x = 1.0000 0.2311 1.0000 0.8913 x([3])=[] x = 1.0000 0.2311 0.8913 2．二维数组 （1)在命令窗口执行下面指令，观察输出结果 a=2.7358; b=33/79; C=[1,2*a+i*b,b*sqrt(a);sin(pi/4),a+5*b,3.5+i] C = 1.0000 5.4716 + 0.4177i 0.6909 0.7071 4.8244 3.5000 + 1.0000i M_r=[1,2,3;4,5,6],M_i=[11,12,13;14,15,16] M_r = 1 2 3 4 5 6 M_i = 11 12 13 14 15 16 CN=M_r+i*M_i CN = 1.0000 +11.0000i 2.0000 +12.0000i 3.0000 +13.0000i 4.0000 +14.0000i 5.0000 +15.0000i 6.0000 +16.0000i (2)利用M 文件创建和保存下面的数组（分号后面是一个空格后加3 个小数点，表示 换行续写） A=(-4:1:5) A = -4 -3 -2 -1 0 1 2 3 4 5 B=reshape(A,2,5) B = -4 -2 0 2 4 -3 -1 1 3 5 B(find(abs(B)>3)) ans = -4 4 5 4)在命令窗口执行下面指令，体会二维数组的子数组寻访和赋值 A=zeros(2,4) 数学实验报告 班级: 学号： 姓名： 实验序号：1 日期：年 月 日 实验名称：特殊函数与图形 ◆ 问题背景描述：绘图是数学中的一种重要手段，借助图形，可以使抽象的对象得到 明白直观的体现，如函数的性质等。同时，借助直观的图形，使初学者更容易接受新知识，激发学习兴趣。 ◆ 实验目的：本实验通过绘制一些特殊函数的图形，一方面展示这些函数的特点属性， 另一方面，就 Matlab 强大的作图功能作一个简单介绍。 实验原理与数学模型： 1、 球2222x y z R ++= ，x=Rsin φcos θ, y= Rsin φsin θ, z= cos φ, 0≤θ≤2π , 0≤φ≤π 环面 222222222()4(),(cos )cos ,x y z a r a x y x a r φθ+++-=+=- (cos )sin ,sin ,02,02y a r z r φθφφπθπ=-=≤≤≤≤ 2、 平面摆线：2 22 31150,(sin ),(1cos ),0233 x y x a t t y a t t π+-==-=-≤≤ 3、 空间螺线：（圆柱螺线）x=acost , y=asint , z=bt ;（圆锥螺线）22 cos ,sin ,x t t y t t z t === 4、 椭球面sin cos ,sin sin ,cos ,02,0x a y b z c φθφθφθπφπ===≤<≤≤ 双叶双曲面3 tan cos ,tan sin ,sec ,02,22 x a y b z c π φθφθφθπφπ===≤<- << 双曲抛物面2 sec ,tan 2 u x au y bu z θθ=== 实验所用软件及版本：mathematica(3.0) 主要内容（要点）： 1、 作出下列三维图形（球、环面） 2、 作出下列的墨西哥帽子 3、 作出球面、椭球面、双叶双曲面，单叶双曲面的图形 4、 试画出田螺上的一根螺线 5、 作出如图的马鞍面 实验1 熟悉Matlb环境及基本操作 实验目的： 1．熟悉Matlab环境，掌握Matlab的主要窗口及功能； 2．学会Matlab的帮助使用； 3．掌握向量、矩阵的定义、生成方法和基本运算； 4．掌握Matlab的基本符号运算； 5．掌握Matlab中的二维图形的绘制和控制。 实验内容： 1．启动Matlab，说明主窗口、命令窗口、当前目录窗口、工作空间窗口、历史窗口、图形窗口、M文件编辑器窗口的功能。 2．实例操作Matlab的帮助使用。 3．实例操作向量、矩阵的定义、生成方法和基本运算。 4．实例操作Matlab的基本符号运算。 5．实例操作Matlab中的二维图形绘制和控制。 实验仪器与软件： 1．CPU主频在2GHz以上，内存在512Mb以上的PC； 2．Matlab 7及以上版本。 实验讲评： 实验成绩： 评阅教师： 年月日 实验1 熟悉Matlab环境及基本操作 一、Matlab环境及主要窗口的功能 运行Matlab安装目录下的matlab.exe文件可启动Matlab环境，其默认布局如下图： 其中， 1．主窗口的功能是：主窗口不能进行任何计算任务操作，只用来进行一些整体的环境参数设置，它主要对6个下拉菜单的各项和10个按钮逐一解脱。 2．命令窗口的功能是：对MATLAB搜索路径中的每一个M文件的注释区的第一行进行扫描，一旦发现此行中含有所查询的字符串，则将该函数名及第一行注释全部显示在屏幕上。 3. 历史窗口的功能是：历史窗口显示命令窗口中的所有执行过的命令，一方面可以查看曾经执行过的命令，另一方面也可以重复利用原来输入的命令行，可以从命令窗口中直接通过双击某个命令行来执行该命令， Matlab 与多元统计分析 胡云峰 安庆师范学院 第三章习题 3.1对某地区的6名2周岁男婴的身高、胸围、上半臂进行测量。得样本数据如表3.1所示。 假设男婴的测量数据X （a ）（a=1，…，6）来自正态总体N 3(μ,∑) 的随机样本。根据以往的资料，该地区城市2周岁男婴的这三项的均值向量μ0=（90,58,16）’，试检验该地区农村男婴与城市男婴是否有相同的均值向量。 表3.1 某地区农村2周岁男婴的体格测量数据 1．预备知识 ∑未知时均值向量的检验： H 0：μ=μ0 H 1：μ≠μ0 H 0成立时 122)(0,)(1)(1,) ()'((1)))()'()(,1)(1)1(,) (1)P P X N n S W n n X n S X n X S X T p n n p T F P n p n p μμμμμ---∑--∑??∴----=-----+∴-- 当 2 (,)(1) n p T F p n p p n α-≥--或者22T T α≥拒绝0H 当 2 (,)(1) n p T F p n p p n α-<--或者22T T α<接受0H 这里2 (1) (, )p n T F p n p n p αα-= -- 2．根据预备知识用matlab 实现本例题 算样本协方差和均值 程序x=[78 60.6 16.5;76 58.1 12.5;92 63.2 14.5;81 59.0 14.0;81 60.8 15.5;84 59.5 14.0]; [n,p]=size(x); i=1:1:n; xjunzhi=(1/n)*sum(x(i,:)); 《数字信号处理》实验报告 学院信息学院专业电子信息工程班级电子08 姓名学号时间2010.10.25 实验一熟悉MATLAB环境 一、实验目的 1、熟悉MATLAB（2008a版本）的主要操作命令。 2、学会离散信号及系统的表达和运算。 3、掌握绘图及查看结果的命令。 4、用MATLAB编程并学会创建函数。 二、实验内容 1、用M语句实现如下向量数组（序列），并由stem语句绘离散序列图（即杆图）。（并标注横轴、纵轴和标题） a. x(n)=0.8n ；n取0－15。 图a 指数序列 b. x(n)=e（0.2＋3j）n ；n取0－15。 图b 复指数序列：模相角 c. x(n)=3cos（0.125πn＋0.2π）＋2sin（0.25πn＋0.1π） n取0－15。 图C 复合正弦序列 d. 把c.的x(n)周期化，周期16点，绘4个周期。 图d－16点周期化 e. 把c.的x(n)周期化，周期10点，绘4个周期。 图e－10点周期化 2、序列的运算： 已知x(n)= [1，－1，3，5 ]，计算y(n)及p(n)并绘杆图。（y(n)为有限长7点,P(n)为9点）。 a. y(n)= 2x(n＋2）－ x(n－1)－2 x(n)； 图a y(n)序列 b. P(n)= ∑nx(n－k），求和k＝1，2，～5 图b P(n)序列 3、编一个用户函数名stepshift（n0，n1，n2）实现阶跃序列 U[n－n0]。其中n0在起点n1和终点n2之间位置可选。并自选3个入口参数绘图。 M文件程序如下所列。 function u = stepsh(N0,N1,N2); u=[zeros(1,N1),zeros(1,N0-N1-1),ones(1,N2-N0+1)]; stem(u); 实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1 例1-2 例1-3 例1-4 3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。 5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？ 多元统计分析中的应用研究 , 摘要：许多实际问题往往需要对数据进行统计分析，建立合适的统计模型，过去一般采用SAS 、SPSS软件分析，本文给出 Matlab软件在多元统计分析上的应用, 主要介绍Matlab 在聚类分析、判别分析、主成份分析上的应用,文中均给以实例, 结果令人满意。 关键词：Matlab软件；聚类分析；主成份分析 Research for application of Multivariate Statistical Analysis Abstract:Many practice question sometimes need Statistical Analysis to data.,and establish appropriate Statistical model SAS and SPSS software were commonly used in foretime ,this paper give the application of Matlab software in Multivariate Statistical Analysis,mostly introduce the application of Matlab software in priciple component analysis and cluster analysis and differentiate analysis.The example are given in writing and the result are satisfaction. Key words: Matlab software; cluster analysis; priciple component analysis 0 引言 许多实际问题往往需要对数据进行多元统计分析, 建立合适的模型, 在多元统计分析方面, 常用的软件有SAS 、SPSS 、S-PLUS等。我们在这里给出Matlab在多元统计分析上的应用, 在较早的版本中, 统计功能不那么强大, 而在Matlab6.x版本中, 仅在统计工具中的功能函数就达200多个, 功能已足以赶超任何其他专用的统计软件,在应用上Matlab具有其他软件不可比拟的操作简单,接口方便, 扩充能力强等优势, 再加上Matlab的应用范围广泛, 因此可以预见其在统计应用上越来越占有极其重要的地位，下面用实例给出Matlab 在聚类分析、主成份分析上的应用。 1 聚类分析 聚类分析法是一门多元统计分类法，其目的是把分类对象按一定规则分成若干类，所分成的类是根据数据本身的特征确定的。聚类分析法根据变量（或样品或指标）的属性或特征的相似性，用数学方法把他们逐步地划类，最后得到一个能反映样品之间或指标之间亲疏关系的客观分类系统图，称为谱系聚类图。 聚类分析的步骤有：数据变换，计算ｎ个样品的两两间的距离，先分为一类，在剩下的ｎ－１个样品计算距离，按照不同距离最小的原则，增加分类的个数，减少所需要分类的样品的个数，循环进行下去，直到类的总个数为１时止。根 《数字信号处理》实验报告 学院诚毅学院专业电子信息工程班级电子109 姓名学号时间2012.10. 实验一熟悉MATLAB环境 一、实验目的 1、熟悉MATLAB的主要操作命令。 2、学会简单的矩阵输入和数据读写。 3、掌握简单绘图命令。 4、用MATLAB编程并学会创建函数。 5、观察离散系统的频率响应 二、实验内容 1、数组的运算：已知两个数组A＝[1 2 3 4] ；B＝[3 4 5 6]；求以下数组： C＝A+B; D＝A－B；E＝A .*B; F＝A ./B; G＝A .^B; 用stem语句画出其中6个数组（向量）的离散序列图（即杆图）。 图－1 A数组图－2 B数组图－3 C数组 图－4 E数组图－5 F数组图－6 G数组 2、序列的产生：（绘杆图，标注横轴、纵轴和标题） a. x(n)=0.8n ；n取0－15。 b. x(n)=e n（0.2＋3j）；n取0－15。 图a 指数序列图b 复指数序列：模相角 c. x(n)=3cos（0.125πn＋0.2π）＋2sin（0.25πn＋0.1π） n取0－15。 图C 复合正弦序列 d. 把c.的x(n)周期化，周期16点，绘4个周期。 e. 把c.的x(n)周期化，周期10点，绘4个周期。 图d－16点周期化图e－10点周期化 3、序列的运算： x(n)= [1，－1，3，5 ] 计算y(n)及p(n)并绘杆图（y(n)为有限长7点,P(n)为9点）。 a. y(n)= 2x(n＋2）－ x(n－1)－2 x(n)； b. p(n)= ∑nx(n－k），求和k＝1，2，～5 图a y(n)序列图b P(n)序列 4、绘时间函数的图形：（在x轴和y轴以及图形上方应加上适当的标注）。 a. x(t)=sin(2πt)；t取0－10秒 b. y(t)=cos(100πt) sin(πt) t 取0－4秒 图a x(t) 正弦信号图b y(t) 调幅信号 实验一小球做自由落体运动内容：一小球竖直方向做自由落体，并无损做往返运动。程序： theta=0:0.01:2*pi x=cos(theta) y=sin(theta) l=1 v=1 while l<10 for t=1:10 y=y+(-1)^l*v*t plot(x,y,[-1,1],[-56,2],'.') axis equal pause(0.1) end l=l+1 end 结果： -50 -40 -30 -20 -10 收获：通过运用小球自由落体规律，及（-1）^n 来实现无损往 返运动！ 实验二 旋转五角星 内容：一个五角星在圆内匀速旋转 程序：x=[2 2 2 2 2 2] y=[0 4/5*pi 8/5*pi 2/5*pi 6/5*pi 0] y1=2*sin(y) x1=2*cos(y) theta=0:4/5*pi:4*pi x2=2*cos(theta) y2=2*sin(theta) plot(x,y,x1,y1,x2,y2) axis equal theta1=theta+pi/10 x2=2*cos(theta1) y2=2*sin(theta1) plot(x2,y2) axis equal theta=0:4/5*pi:4*pi for rot=pi/10:pi/10:2*pi x=2*cos(theta+rot) y=2*sin(theta+rot) plot(x,y) pause(0.1) end 结果： -2 -1.5-1-0.500.51 1.52 -2-1.5-1-0.500.511.5 2 收获：通过theta1=theta+pi/10，我们可以实现五角星在圆内匀速 旋转！ 实验三 转动的自行车 内容：一辆自行车在圆内匀速转动 程序：x=-4:0.08:4; y=sqrt(16-x.^2); theta1=-pi/2:0.01*pi:3*pi/2; x3=0.5*cos(theta1); y3=0.5*sin(theta1); theta=-pi/2+0.02*pi for k=1:100 实验一 Matlab基础知识 一、实验目的： 1.熟悉启动和退出Matlab的方法。 2.熟悉Matlab命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握Matlab各种表达式的书写规则以及常用函数的使 用。 二、实验内容： 1.求[100，999]之间能被21整除的数的个数。（rem） 2.建立一个字符串向量，删除其中的大写字母。（find） 3.输入矩阵，并找出其中大于或等于5的元素。（find） 4.不采用循环的形式求出和式 63 1 2i i= ∑ 的数值解。(sum) 三、实验步骤： ●求[100，199]之间能被21整除的数的个数。（rem） 1.开始→程序→Matlab 2.输入命令： ?m=100:999; ?p=rem(m,21); ?q=sum(p==0) ans=43 ●建立一个字符串向量，删除其中的大写字母。（find） 1.输入命令： ?k=input('’,’s’)； Eie48458DHUEI4778 ?f=find(k>=’A’&k<=’Z’)； f=9 10 11 12 13 ?k(f)=[ ] K=eie484584778 ●输入矩阵，并找出其中大于或等于5的元素。（find） 1.输入命令： ?h=[4 8 10;3 6 9; 5 7 3]； ?[i,j]=find(h>=5) i=3 j=1 1 2 2 2 3 2 1 3 2 3 ●不采用循环的形式求出和式的数值解。（sum） 1.输入命令： ?w=1:63； ?q=sum(2.^w) q=1.8447e+019 实验二 Matlab 基本程序 一、 实验目的： 1. 熟悉Matlab 的环境与工作空间。 2. 熟悉M 文件与M 函数的编写与应用。 3. 熟悉Matlab 的控制语句。 4. 掌握if,switch,for 等语句的使用。 二、 实验内容： 1. 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 2. 编程完成，对输入的函数的百分制成绩进行等绩转换，90~100为优，80~89为良，70~79为中，60~69为及格。 3. 编写M 函数文件表示函数 ，并分别求x=12和56时的函数值。 4. 编程求分段函数 2226;03 56;0532 1;x x x x y x x x x x x x +-<≠=-+≤<≠≠-+且且及其它，并求输入x=[-5.0,-3.0,1.0,2.0,2.5,3.0,3.5]时的输出y 。 三、 实验步骤： 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 1. 打开Matlab ，新建M 文件 2. 输入命令： 51022-+x Matlab多元统计分析程序 1. 主成分分析M程序 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 主成分分析 % % 设对变量x1,x2,...,xp进行n次观测，得到n×p数据矩阵x=x(i,j)， % 本程序对初始数据进行主成分分析，要求先请将观测矩阵输入到变 % 量x,再运行本程序。 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % 确定观测矩阵x 的尺寸，以便数据标准化. % [n,p]=size(x); % % 数据处理方式设置，即是否先将数据标准化. % fprintf('\n 1---使用原始数据直接计算距离') fprintf('\n 2---使用标准化后的数据计算距离') k=input('请输入你的选择（1～2）'); % % 数据标准化 % switch k case 1 xs=x; case 2 mx=mean(x); xs=(x-repmat(mx,n,1))./repmat(stdr,n,1); end % % 主成分分析，返回各主成分pc，所谓的z－得分score，x的协方差 % 矩阵的特征值latent和每个数据点的Hotelling统计量tsquare. % [pc score latent tsquare]=princomp(xs) 2. 典型相关分析M程序 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % 在运行本程序之前，请先把数据输入/导入到MATLAB 的 % 内存空间，并存放在变量x 中，每行存放一个样本。 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % 确定观测矩阵x 的尺寸 % [n,h]=size(x); % % 输入基本参数 % p=input('\n第一组变量的个数p = ? '); fprintf('\n1--使用样本协方差矩阵计算典型相关变量') fprintf('\n2--使用样本相关矩阵计算典型相关变量') ctl=input('\n请输入你的选择'); % % 默认的显著性水平为alpha=0.05，可以改变下面语句中的alpha值。 % alpha=0.05; % % 按要求计算样本协方差矩阵或样本相关矩阵 % switch ctl case 1 st=cov(x); case 2 实验一 MA TLAB 环境熟悉及基本操作 1. 命令窗口运行入门 1-1 求2 ))39(212(2 -*+的算术运算结果 1-2 简单矩阵A 的生成 1-3 绘制图形比较函数曲线5432x x x x 、、、，要求（一）四个曲线在一个图形文件窗口（二） 子图 2. 命令窗口常用的基本操作 1-4 若用户想计算51) 2sin(21+=πy ，51) 2cos(22+=πy 的值 1-5 变量的保存：save 和load 命令，将已生成的变量保存，名为×××，存入目录C ：\my_dir 3. 系统的在线帮助 1-6 help sqrt/help syntax/help (回车) 1-7 lookfor line(查找与直线、线性问题有关的函数) 1-8 help 菜单 实验二 数组(矩阵)及其运算 2-1 用户不能清除的固定变量ans 、eps 、pi 、i 、j 、NaN 、Inf 等 如>>eps Ans=2.2204e-016 2-2 矩阵的创建 （1）>>a=[1 2 3;4 5 6;7 8 9] >>a=[1:2;4:6;7:9]; a1=1:6 （2）利用MA TLAB 内部函数产生矩阵 >>b=eye(3) >>c=ones(2,5) >>d=zeros(3,2) >>e=linspace(-3,6,10) %注意三个参数的含义 >>f=logspace(1,6,5) %注意三个参数的含义 >>R=rand(1,6) （3） 访问元素 >>aa=a(1,2)+a(2,3) >>bb=a(4)+a(8) >>a(end-2:end,2) %在下标中可以直接用end 表示这一维的最后一个元素 （4） 矩阵操作 >>A=[1,2,3;4,5,6;7,8,9] >>B=diag(A) >>C=diag(B) >>D=rot90(A) >>E=reshape(D,1,9) %注意这三个参数的含义 2-3 矩阵运算 （1）算术运算 + - * / \ ^ ‘ 作业一 1.2 分析2016年经济发展情况 排名省gdp 占比累计占比 1 广东79512.05 10.30 10.30 2 江苏76086.2 9.86 20.17 3 山东67008.2 8.68 28.85 4 浙江4648 5 6.02 34.87 5 河南40160.01 5.20 40.08 6 四川32680.5 4.24 44.31 7 湖北32297.9 4.19 48.50 8 河北31827.9 4.12 52.62 9 湖南31244.7 4.05 56.67 10 福建28519.2 3.70 60.37 11 上海27466.2 3.56 63.93 12 北京24899.3 3.23 67.16 13 安徽24117.9 3.13 70.28 14 辽宁22037.88 2.86 73.14 15 陕西19165.39 2.48 75.62 16 内蒙古18632.6 2.41 78.04 17 江西18364.4 2.38 80.42 18 广西18245.07 2.36 82.78 19 天津17885.4 2.32 85.10 20 重庆17558.8 2.28 87.37 21 黑龙江15386.09 1.99 89.37 22 吉林14886.23 1.93 91.30 23 云南14869.95 1.93 93.22 24 山西12928.3 1.68 94.90 25 贵州11734.43 1.52 96.42 26 新疆9550 1.24 97.66 27 甘肃7152.04 0.93 98.59 28 海南4044.51 0.52 99.11 29 宁夏3150.06 0.41 99.52 30 青海2572.49 0.33 99.85 31 西藏1150.07 0.15 100.00 将2016各省的GDP进行排名，可以发现，经济发达的的地区主要集中在东部地区。西部gdp的占比较小。作出2016各省的gdp直方图如下： Matlab实验报告 实验二图像处理 一、实验目的 （1）通过应用MA TLAB语言编程实现对图像的处理，进一步熟悉MATLAB软件的编程及应用； （2）通过实验进一步掌握图像处理的基本技术和方法。 二、实验内容及代码 ㈠.应用MA TLAB语言编写显示一幅灰度图像、二值图像、索引图像及彩色图像的程序，并进行相互之间的转换 首先，在matlab页面中的current directory下打开存放图像的文件夹。 1.显示各种图像 ⑴显示彩色图像： ①代码：>> mousetif=imread('tif.TIF'); >> image(mousetif) 显示截图： ②代码：>> mousetif=imread('tif.TIF'); >> imshow(mousetif) 显示截图： ③代码：mousetif=imread('tif.TIF'); subimage(mousetif) 显示截图： 显示截图： ⑵显示二值图像 ①代码：>> I=imread('单色bmp.bmp'); >> imagesc(I,[0 2]) 显示截图： ②代码：>> I=imread('单色bmp.bmp'); >> imshow(I,2) 显示截图： ③代码：>> I=imread('单色bmp.bmp'); >> subimage(I) 显示截图： ⑶显示灰度图像 ①代码：>> I1=imread('256bmp.bmp'); >> imagesc(I1,[0,256]) 显示截图： 代码：>> I1=imread('256bmp.bmp'); >> colormap(gray); >> subplot(1,2,1); >> imagesc(I1,[0,256]); >> title('灰度级为[0 256]的mouse.bmp图'); >> subplot(1,2,2); >> imagesc(I1,[0,64]); >> colormap(gray); >> title('灰度级为[0 64]的mouse.bmp图'); 显示截图： 实验四浸出制剂的制备 实验一 MATLAB集成环境使用实验目的: 1.熟悉MATLAB的基本环境 2.掌握MATLAB的基本命令 3.掌握简单线性方程组的求解 4.掌握MATLAB的行编辑和页面编辑 实验内容与要求 1.观察MATLAB 的基本环境,熟悉菜单和工具栏 2.MATLAB基本命令的使用 a.运行demos, 观看演示程序 b.工作空间管理命令: who , whos , clear , which , quit ,exit workspace c.路径编辑命令: path, cd ,addpath rmpath d.显示格式命令: type wall.m , more on , type coneball.m , more off ,type coneball.m e.窗口清理命令: coneball2 clc cla clf close f. 操作系统命令: copyfile pwd dir delete mkdir g. 联机帮助命令:help help . helpwin lookfor zoom which rguizoom help rguizoom what ver signal 3.简单线性方程组的求解 以学号 208040401 为例 A=[ 2 0 8 0 4 0 4 0 1] b=[ 208 40 401] x==A\b 4.行编辑与页面编辑 a.键入magic(4) 使用,调出该命令,并改变为magic(6)执行. b.执行coneball命令 edit命令, 打开coneball.m文件 改变文件中的半径和高, 存盘后,重新执行该命令,比较前后的不同. 页脚内容1 MATLAB的统计工具箱中的多元统计分析中提供了聚类分析的两种方法： 1.层次聚类hierarchical clustering 2.k-means聚类 这里用最简单的实例说明以下层次聚类原理和应用发法。 层次聚类是基于距离的聚类方法，MATLAB中通过pdist、linkage、dendrogram、cluster等函数来完成。层次聚类的过程可以分这么几步： (1) 确定对象（实际上就是数据集中的每个数据点）之间的相似性，实际上就是定义一个表征对象之间差异的距离，例如最简单的平面上点的聚类中，最经常使用的就是欧几里得距离。 这在MATLAB中可以通过Y=pdist（X）实现，例如 >> X=randn(6,2) X = -0.4326 1.1892 -1.6656 -0.0376 0.1253 0.3273 0.2877 0.1746 -1.1465 -0.1867 1.1909 0.7258 >> plot(X(:,1),X(:,2),'bo') %给个图，将来对照聚类结果把 >> Y=pdist(X) Y = Columns 1 through 14 1.7394 1.0267 1.2442 1.5501 1.6883 1.8277 1.9648 0.5401 2.9568 0.2228 1.3717 1.1377 1.4790 1.0581 Column 15 2.5092 例子中X数据集可以看作包含6个平面数据点，pdist之后的Y是一个行向量，15个元素分别代表X 的第1点与2-6点、第2点与3-6点,......这样的距离。那么对于M个点的数据集X，pdist之后的Y 将是具有M*(M-1)/2个元素的行向量。Y这样的显示虽然节省了内存空间，但对用户来说不是很易懂，如果需要对这些距离进行特定操作的话，也不太好索引。MATLAB中可以用squareform把Y转换成方阵形式，方阵中 实验一 matlab 的开发环境与基本运算 一、实验学时 4学时（课内2学时，课外2学时） 二、实验类型 验证性实验 三、实验目的 1、熟悉matlab 的开发环境 2、掌握矩阵、变量、表达式的各种运算 四、所需设备及软件 1、安装了windows xp/win7/win8/win10的计算机 2、matlab 开发工具 五、实验基本原理 1、熟悉MATLAB 环境： MATLAB 桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器。 2、掌握MATLAB 常用命令 表1 MATLAB 常用命令 3、MATLAB 变量与运算符 1）变量命名规则 2）MATLAB 的各种常用运算符 表2 MATLAB 算术运算符 表3 MATLAB 关系运算符 表4 MATLAB 逻辑运算符 表5 MATLAB 特殊运算 4、MATLAB 的一维、二维数组的寻访（约定大写字母为矩阵） 表6 子数组访问与赋值常用的相关指令格式 指令 含义 指令 含义 A=a 把标量a 赋值给A 的每个元素 a+A 标量a 分别与A 元素之和 a-A 标量a 分别与A 元素之差 a.*A 标量a 分别与A 元素之积 a *A 标量a 分别与A 每个元素之积 A.^n A 的每个元素自乘n 次 A ^n A 为方阵，自乘n 次 A.*B 对应元素相乘 A*B 内维相同，矩阵相乘 指令 含义 指令 含义 Exp(A) 以自然数e 为底，分别以A 的元素为指数，求幂 log(A) 对A 的各元素求对数 Sqrt(A) 对A 的各元素求平方根 Diag(A ) A 矩阵的对角线元素 Rank 矩阵的秩 Det 方阵的行列式值 Ones 产生全1的数组或矩阵 Zeros 产生全0的数组或矩阵 Rand 产生均匀分布的随机数组或矩阵 Randn 产生正态分布的随机数组 Size 返回指定矩阵的行数和列数 Rot90 矩阵逆时针旋转90度 Factor 因式分解 Expand 展开表达式 六、实验内容 1、建立一个以学号命名的文件夹；将该文件夹添加到matlab 的路径管理器中 2、启动matlab ，熟悉 命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器 3、学习使用help 命令 4、使用窗口命令：clc 、close 、close all 、clear ，了解其功能和作用，观察各个窗口[command window 、command history 、workspace]的变化结果。要求用文字说明。 5、练习“五”中各运算符、函数的应用。注意观察运算结果。在实验报告中必须完成以下内容。 （1）矩阵运算 12344347873657A -????=??????，131203327B -?? ??=?? ??-?? 求下列表达式的值： MATLAB实验报告 一、实验名称 实验4图形绘制（1） 二、实验目的： 熟悉和掌握MA TLAB基本的二维图形绘制函数。 三、实验内容： 1．绘制简单的二维图形 2．一个坐标系绘制多幅图形 3.图形标识和坐标控制 4.交互式图形指令 四、回答问题： （本次实验未预留问题） 五、遇到的问题及解决： 遇到了求y=lnx时，输入“y=ln(x)”不被软件识别的问题，查看常用数学函数表后改为y=log(x)成功解决。 在求10x时不知道用什么函数，函数表里也查不到，在老师的点拨下用“y=10.^x”解决。 在绘图时发现默认线型不够明显，查表后使用尖三角、叉号代替默认线型。 六、体会： 本次实验我学会了利用MATLAB绘制图形的基本方法，以及相应的备注方法。 难点是了解各种函数的具体作用并熟练掌握。 体会是：多学多练，孰能生巧，日积月累，必有提高。 思考题： 1.在同一坐标系绘制t3,-t2,t2sint在[0,2π]内的曲线图。 x=0:pi/50:2*pi; y1=t.*t.*t; y2=-t.*t; y3=t.*t.*sin(t); plot(t,y1,'^k',t,y2,'.k',t,y3,'xk'); legend('\ity=t^3','\ity=-t^2','\itt^2*sint'); 2.在一幅图中画出4幅子图，分别绘制sin2x，tanx，lnx，10x的图形，并加上适当的图形注释。注意：把函数变成MATLAB对应的形式。 x=0:pi/50:2*pi; y1=sin(2*t); y2=tan(x); y3=log(x); y4=10.^x; subplot(2,2,1) plot(x,y1); legend('y=sin2x'); subplot(2,2,2) plot(x,y2) legend('y=tanx'); subplot(2,2,3) plot(x,y3) 实验一 MATLAB 的实验环境及基本命令 一实验目的： 1. 学习了解MATLAB 勺实验环境 2. 在MATLA 系统命令窗口练习有关 MATLAB^令的使用。 二实验步骤 1. 学习了解MATLAB 勺实验环境： 在 Windows 桌面上，用 mouse 双击MATLAB!标，即可进入 MATLA 啲默认界面: 图1-1 MATLAB 系统命令窗口 MATLAB 勺默认界面由三部分组成： ① Comma nd Win dows 命令窗口 ② Comma nd History 历史窗口 ③ Workspace 工作空间 2. 练习MATLA 療统命令的使用。 ①表达式 MATLAB 勺表达式由变量、数值、函数及操作符构成。实验前应掌握有关变量、数值、函数及操作 符的有关内容及使用方法。 三实验内容 练习1-1 :计算下列表达式： 要求计算完毕后，键入相应的变量名，查看并记录变量的值。 a=(1+{10)/2 b=|3+5| c=sin(e±) d = 250 /sin ( ’) 实验过程及结果：a=(1+sqrt(10))/2 a =2.0811 b= abs(3+5*i) b =5.8310 c= si n(exp(-2.3)) c=0.1001 d=250/s in (pi/6) d=500.0000 ②.向量运算：n 维向量是由n 个成员组成的行或列数组。在 MATLA 中，由分号分隔的方括号中的元 素产生一个列向量；由逗号或空号分隔的方括号中的元素产生一个列向量；同维的向量可进行加减运 算，乘法须遵守特殊的原则。 练习1-2 已知：X=[2 ; -4 ; 8] 求：Y=R ; P=5*R; E=X* Y ; S=X * Y 实验过程及结果：R=[tan(pi/4),sqrt(9),-5] R=1.0000 3.0000 - Y=R' Y=1.0000 P=5*R P =5.0000 15.0000 -25.0000 R=1.0000 3.0000 -5.0000 m a t l a b与多元统计分 析 Company Document number：WTUT-WT88Y-W8BBGB-BWYTT-19998 Matlab 与多元统计分析 胡云峰 安庆师范学院 第三章习题 对某地区的6名2周岁男婴的身高、胸围、上半臂进行测量。得样本数据如表所示。假设男婴的测量数据X （a ）（a=1，…，6）来自正态总体N 3(,∑) 的 随机样本。根据以往的资料，该地区城市2周岁男婴的这三项的均值向量0= （90,58,16）’，试检验该地区农村男婴与城市男婴是否有相同的均值向量。 表 某地区农村2周岁男婴的体格测量数据 解 1．预备知识 ∑未知时均值向量的检验： H 0：=0 H 1：≠0 H 0成立时 122)(0,)(1)(1,) ()'((1)))()'()(,1)(1)1(,) (1)P P X N n S W n n X n S X n X S X T p n n p T F P n p n p μμμμμ---∑--∑??∴----=-----+∴-- 当 2 (,)(1) n p T F p n p p n α-≥--或者22T T α≥拒绝0H 当 2 (,)(1) n p T F p n p p n α-<--或者22T T α<接受0H 这里2(1) (, )p n T F p n p n p αα-= -- 2．根据预备知识用matlab 实现本例题 算样本协方差和均值 程序x=[78 ;76 ;92 ;81 ;81 ;84 ]; [n,p]=size(x); i=1:1:n; xjunzhi=(1/n)*sum(x(i,:)); y=rand(p,n); for j=1:1:n y(:,j)= x(j,:)'-xjunzhi'; y=y; end A=zeros(p,p); for k=1:1:n; A=A+(y(:,k)*y(:,k)'); end xjunzhi=xjunzhi' S=((n-1)^(-1))*A 输出结果xjunzhi = S = 然后u=[90;58;16]; t2=n*(xjunzhi-u)'*(S^(-1))*(xjunzhi-u) f=((n-p)/(p*(n-1)))*t2 输出结果t2 = f = 所以21()'()T n X S X μμ-=--= 实验一熟悉MATLAB开发环境 1、实验目的 1）熟悉M A TLAB 的开发环境，掌握常用菜单的使用方法； 2）熟悉M A TLAB 工作界面的多个常用窗口包括命令窗口、历史命令窗口、当前工作目录窗口、工作空间览器窗口、变量编辑窗口和M文件编辑／调试窗口等； 3）了解M A TLAB 的命令格式。 2、实验内容 2.1 MATLAB的启动 方式一：Windows的“开始”菜单方式； 方式二：桌面图标方式。 2.2 MATLAB 的窗口布局 图1 MATLAB 7.x 工作环境 2.3 使用命令窗口（Command Window） 命令窗口如图 S1-2 所示，在命令窗口中输入： a=[1 2;3 4] a = 1 2 3 4 b=1/3 b = 0.3333 c=a*b c = 0.3333 0.6667 1.0000 1.3333 2.4 使用标点符号来修改命令行 (1) %添加注释 c=a+b % c为矩阵 a+b 的和 (2) ;不显示计算结果 a=[1 2;3 4]; (3)…续行 d=123456 + … 111111 (4) ,用作数组元素的分隔 a=[1,2;3,4] a = 1 2 3 4 2.5 使用操作键 ↑：向前调回已输入过的命令行； ↓：向后调回已输入过的命令行； Esc：清除当前行的全部内容。 2.6 使用 Format 设置数值的显示格式 format long：显示 15 位长格式 format long c c = 0.33333333333333 0.66666666666667 1.00000000000000 1.33333333333333 format rat：显示近似有理式格式 format rat c c = 1/3 2/3 1 4/3 2.7 设置命令窗口的外观 在 MATLAB 的界面选择菜单“File”---“Preferences” 2.8 使用控制命令 clc %清空命令窗口的显示内容 clear beep on %发出 beep 的声音 beep 2.9 历史命令窗口（Command History） 历史命令窗口在 MATLAB界面的左下侧，在历史命令窗口中可以看到本次启动MATLAB 的时间和已经输入的命令。 2.10 工作空间窗口（Workspace） 工作空间窗口在 MATLAB界面的左上 保存变量 Save filename 读回： Load filename 2.11 当前目录浏览器窗口（Current Directory Browser） （略）matlab实验报告
实验1熟悉matlab环境和基本操作
matlab与多元统计分析
实验一 熟悉MATLAB环境(实验报告)2010
MATLAB实验报告50059
matlab与应用多元统计分析
熟悉MATLAB环境(实验报告)
matlab实验报告
参考答案Matlab实验报告
Matlab多元统计分析程序
MATLAB环境熟悉及基本操作
多元统计分析
matlab实验报告
实验1 MATLAB集成环境使用
MATLAB的统计工具箱中的多元统计分析中提供了聚类分析的两种方法
实验一 matlab的开发环境与基本运算
MATLAB程序设计实验报告
MATLAB的实验环境及基本命令
matlab与多元统计分析
实验一 熟悉 MATLAB开发环境
相关主题
文本预览