《经济管理数学》
实 验 指 导 册
学院 : 班级: 姓名: 学号:
1
实验一:Matlab基本操作
一、实验目的与要求
1、熟悉Matlab软件的安装、界面、窗口及一些基本功能; 2、掌握Matlab软件的一些基本操作和常用命令; 3、熟悉Matlab软件的一些数值计算功能。 二、实验内容
从Windows中双击Matlab图标,会出现Matlab命令窗口(Command Window),在一段提示信息后,出现系统提示符“>>”.Matlab是一个交互系统,您可以在提示符后键入各种命令,通过上下箭头可以调出以前打入的命令,用滚动条可以查看以前的命令及其输出信息. 如果对一条命令的用法有疑问的话,可以用Help菜单中的相应选项查询有关信息,也可以用help命令在命令行上查询,您可以试一下help、help help和help eig(求特征值的函数)命令. 1、 常用运算符号: 加法 + 减法 - 乘法 * 点乘 .* 点除 ./ 乘幂 左除 \\ 右除 / ^ 例1 计算1.3692?sin??7???26.48?2.9的值。 ?10?输入:>>a=1.369^2+(sin(7*pi/10)*26.48^(1/2))/2.9 结果:>>a = 3.3097 2、 矩阵的输入
Matlab是以矩阵为基本变量单元的,因此矩阵的输入非常方便。输入时,矩阵的元素用方括号括起来,行内元素用逗号分隔或空格分隔,各行之间用分号分隔或直接回车。
?112???例2 输入矩阵 A???103? ,可以在命令窗口中输入
?4?56???>>A=[1 1 2;-1 0 3;4 -5 6] 3 矩阵的结构操作
输入矩阵后,可以对矩阵进行的主要操作包括矩阵的扩充,矩阵元素的提取,矩阵元素的
2
部分删除等。下面对其作简单的介绍 (1) 矩阵的扩充
例如,用下述命令可以在上述矩阵A下面再加上一个行向量: >> A(4,:)=[1 3 2]
下述命令可以在上述矩阵A下面再加上一个列向量: >>A(:,4)=[-1 0 3 2]
可以用下述命令对矩阵进行翻转和旋转 fliplr(A) 矩阵A左右翻转 flipud(A) 矩阵A上下翻转
rot90(A) 矩阵A整体逆时针旋转90度 (2) 矩阵元素的提取
可以用下述命令提取上述矩阵A的第3行第1列的元素; >> A(3,1)
可以用下述命令提取上述矩阵A的第1列和第3列的元素; >> A(:,[1,3])
可以用可以用下述命令提取矩阵的上三角和下三角部分和对角线元素 triu(A) 提取矩阵A的上三角部分 tril(A) 提取矩阵A的下三角部分 diag(A) 提取矩阵A的对角线元素 (3) 矩阵元素的删除
可以用下述命令删除上述矩阵A的第2行的元素; >>A (2,:)=[ ] 4、 特殊矩阵的生成
某些特殊矩阵可以直接调用相应的函数得到,例如: zeros(m,n) 生成一个m行n 列的零矩阵 ones(m,n) 生成一个m行n 列元素都是1的矩阵 eye(n) 生成一个n 阶的单位矩阵 rand(m,n) 生成一个m行n 列的随机矩阵 vander(V) 生成以向量V为基础向量的范得蒙矩阵 magic(n) 生成一个n阶魔方矩阵
2
hilb(n) 生成一个n阶希尔伯特矩阵 invhilb(n) 生成n阶的希尔伯特矩阵的逆矩阵. 例3 随机生成一个6?7的矩阵。 >> rand(6,7)
例4 生成一个以向量(1,2,3,5)为基础向量的范得蒙矩阵. >> vander([1;2;3;5]) 5、数组运算
(1)标量与数组运算 例如 输入: a=[1 2 3 4]; c=2; a1=a+c; a2=a*c; a3=a./c; a4=a.\\c; a5=a.^c; a6=c.^a; 运行结果:
a1 = 3 4 5 6 a2 =2 4 6 8
a3 =0.5000 1.0000 1.5000 2.0000 a4 =2.0000 1.0000 0.6667 0.5000 a5 = 1 4 9 16 a6 = 2 4 8 16 (2)数组与数组运算 例如编写: a=[1 2 3]; b=[2 2 2]; c1=a+b
3
c2=a.*b c3=a./b c4=a.\\b c5=a.^b 结果为:
c1 = 3 4 5 c2 = 2 4 6
c3 =0.5000 1.0000 1.5000 c4 = 2.0000 1.0000 0.6667 c5 = 1 4 9 三、实验作业
?3?781567???1、 输入矩阵A??058?1011?,并提取矩阵A的第3列和第2行元素。
?5?761829???ans = 8 8 6
>> A([2],:) ans =
0 5 8 -10 11 2、生成一个10?12阶的随机矩阵。 3、用matlab生成以下矩阵
4