实验2 数值数组及运算
一、实验目的
（1）掌握一维数组、二维数组、多项式和高维数组的创建
（2）掌握对数组操作的常用函数
二、实验内容
1．使用冒号运算符创建下面各表达式对应的向量。
a) x=linspace(1,10,5); x=1:2.25:10
b) x=linspace(-5,5); x=-5:0.1:5
c) x=logspace(1,3,3); x=10.^[1；3]
d) x=logspace(1,3,5);
2．利用MATLAB 内置函数计算下面的值。
a) cosh(5) cosh(5)
b) sinh(-2) sinh(-2)
c) (e 5+e -5)/2 (exp(5)+exp(-5))/2
d) erf(1.2) ，其中，ηηd e x erf x
⎰-=02
)(（提示：用help 查一下erf 的调用格式）
erf(1.2)
3．用linspace 函数创建下面表达式的对应向量
a) x=0:10; x=linspace(0,10,11)
b) x=0:0.2:10; x=linspace(0,10,51)
c) x=-12:12; x=linspace(-12,12,25)
d) x=10:-1:1; x=linspace(10,1,10)
4．给定行向量[]78910=x 和列向量⎥⎥⎥⎥
⎦
⎤⎢⎢⎢⎢⎣⎡=4321y ，至少使用两种不同的方法求行向量z 。其中i i i y x z -=
x=[10 9 8 7]; y=[1 2 3 4]’; z=x-y ’
5．在MATLAB 中输入下列矩阵，再求出他们的乘积矩阵C ，并将C 矩阵的右下角的2×3子矩阵赋给D 矩阵。
⎥⎥⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎢⎢⎣⎡=49819323753175323321A ，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡++=3498143576255332763441i i B
A=[1 2 3 3;2 3 5 7;1 3 5 7;3 2 3 9;1 8 9 4];
B=[1+4i 4 3 6 7;2 3 3 5 5;2 6+7i 5 3 4;1 8 9 4 3];
C=A*B
D=C(4:5,3:5)
6．根据下面两个矩阵，执行下列的矩阵运算指令。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=136782078451220124A ，⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=087654321B A=[4 12 20;12 45 78;20 78 136];
B=[1 2 3;4 5 6;7 8 0];
I=[1 0 0;0 1 0;0 0 1];
（1） A+5*B 和A-B+I 分别是多少（I 是单位矩阵）
A+5*B A-B+I
（2） A.*B 和A*B 将分别给出什么结果，他们是否相同，为什么？
A.*B A*B
不带点式矩阵相乘，带点式数组相乘
（3） 得出A.^B ，A\B 和A/B 的结果
A.^B A\B A/B
（4） 找出A 矩阵中大于50到100之间的元素，B 矩阵中大于5的元素 A(find((A>50)&(A<100))) B(find(B>5))
7．若[]NaN Inf Inf NaN A 521-=，运行下面的指令，并解释为什么会出现这样的结果。isnan(A)，isfinite(A)，isinf(A)，any(A)，all(A).
8．求出以下矩阵的逆矩阵。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=627943251A ，⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡---=)6sin()2sin()7sin()9sin()4sin()3sin()2sin()5sin()1sin(B inv(A) inv(B) inv(取逆)
9．编写MATLAB 语句来手工输入矩阵A ，然后从矩阵A 得到矩阵B
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A ，⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=321654987B A=[1 2 3;4 5 6;7 8 9];
B=flipud(A) flipud(第一行与最后一行互换)
10． 矩阵可以看成是行向量和列向量的组合。给定行向量[]321=u 和
[]654=v ，用一条语句求2×3矩阵A 。其中u 是第一行，v 是第二行。
u=[1 2 3];
v=[4 5 6];
A=[u;v]
11． 给定矩阵⎥⎥⎦
⎤⎢⎢⎢⎣⎡=71812511C ，用两行语句实现从C 中引用列向量，得出]18211['=s 和]715['=t 。
C=[11 5;2 1;18 7];
s=C(:,1) t=C(:,2)
12． 使用diag 函数和冒号运算符创建下面矩阵
⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡19000
0130000700001)(,1000020000300004)(,4000030000200001)(c b a diag(1:1:4) diag(4:-1:1) diag(1:6:19)
13． 使用diag 函数创建n ×n 对称三对角阵。 ⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡------=2100121001210012D a=[2 2 2 2];b=[-1 -1 -1]
D=diag(a)+diag(b,-1)+diag(b,1)
14． 创建如下矩阵
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=001010100E e=eye(3);E=fliplr(e) fliplr(第一列与最后一列互换) 15． 用一条语句创建如下矩阵（提示：使用矩阵相加以及内置命令tril （取
下三角），ones （全一），eye (单位阵)和zeros （全零））triu(取下三角)
⎥⎥⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎢⎢⎣⎡----------=1111111111101111001110001A a=2*eye(5)-ones(5) b=tril(a) b(:,5)=1
16． 使用rehape 函数和冒号运算符创建下列矩阵。
⎥⎦⎤⎢⎣⎡-----=⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=642024531135)(,24181262216104201482)(B b A a reshape(2:2:24,3,4) reshape(-5:1:6,2,6) reshape(按列加)
