实验2 数值数组及运算

1、实验2数值数组及运算i=一、实验目的掌握一维数组、二维数组、多项式和高维数组的创建掌握对数组操作的常用函数二、实验内容使用冒号运算符创建下面各表达式对应的向量。x=linspace(1,10,5); x=1:2.25:10 x=linspace(-5,5); x=-5:0.1:5x=logspace(1,3,3); x=10.”1； 3x=logspace(1,3,5);利用MATLAB内置函数计算下面的值。cosh(5)cosh(5)sinh(-2)sinh(-2)(e5+e-5)/2(exp(5)+exp(-5)/2d) erf(1.2)，其中，erf (x) = I x e -门2d门(

2、提示：用 help 查一下 erf 0的调用格式)erf(1.2)用linspace函数创建下面表达式的对应向量x=0:10;x=linspace(0,10,11)x=0:0.2:10;x=linspace(0,10,51)x=-12:12;x=linspace(-12,12,25)x=10:-1:1;x=linspace(10,1,10)4.给定行向量x =110 9 8 7和列向量j =,至少使用两种不同的方法求行向量2。其中z = x. yx=10 9 8 7; y=1 2 3 4 ; z=x-y在MATLAB中输入下列矩阵，再求出他们的乘积矩阵C,并将C矩阵的右下 角的2X3子矩阵赋给

3、D矩阵。12331 + 4.4367一2357233551357,B =26 + 7.534323918941894311A=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*BD=C(4:5,3:5)根据下面两个矩阵，执行下列的矩阵运算指令。-4 12 20 -1 2 3一A =12 45 78，B =4 5 620 78 1367 8 0A=412 20;12 4578;2078 136B=1 2 3;4 5 6;7 8 0;I=1 0 0;0 1 0;0

4、0 1;A+5*B和A-B+I分别是多少(I是单位矩阵)A+5*B A-B+IA.*B和A*B将分别给出什么结果，他们是否相同，为什么？A.*B A*B不带点式矩阵相乘，带点式数组相乘得出A.B AB和A/B的结果A.B AB A/B找出A矩阵中大于50到100之间的元素，B矩阵中大于5的元素 A(find(A50)&(A5)若A = 1 2 NaN Inf - Inf 5 NaN，运行下面的指令，并解释为什么 会出现这样的结果。isnan(A)，isfinite(A)， isinf(A)，any(A)，all(A).求出以下矩阵的逆矩阵。一 1- 5 - 2一sin(1)sin(-5)sin

5、(-2)A =34 - 9，B =sin(3) sin(4)sin(9)-726sin(-7)sin(2)sin(6)inv(A) inv(B)inv(取逆)9.编写MATLAB语句来手工输入矩阵A，然后从矩阵A得到矩阵B1 2 3-7 8 9一A=456，B =4 5 67 8 91 2 3A=12 3;4 56;7 89;B=flipud(A)flipud(第一行与最后一行互换)10.矩阵可以看成是行向量和列向量的组合。给定行向量u = 1 2 3和v =如5 6，用一条语句求2X3矩阵A。其中u是第一行，v是第二行。u=1 2 3;v=4 5 6;A=u;v11. 给定矩阵 C 二s =

6、 11 21121818和 t = 5，用两行语句实现从C中引用列向量，得出1 7。12.1；18 7;t=C(:,2)C=11 5;2s=C（:,1）1000-4000-1000 一0200,(b)0300&）070000300020001300004000100019(a)diag(1:1:4)diag（4:-1:1）13. 使用diag函数创建nXn对称三对角阵。diag(1:6:19)使用diag函数和冒号运算符创建下面矩阵 TOC o 1-5 h z III0 HYPERLINK l bookmark113 o Current Document D = -12- 100T2-1 HY

7、PERLINK l bookmark116 o Current Document 00-12a=2 2 2 2;b=-1 -1 -1D=diag(a)+diag(b,-1)+diag(b,1)创建如下矩阵0 0 1E = 0 10 e=eye(3);E=fliplr(e) fliplr(第一列与最后一列互换)1 0 0用一条语句创建如下矩阵（提示：使用矩阵相加以及内置命令tril （取下三角），ones(A100011100111101-1111111111）,eye（单位阵）和zeros （全零）triu（取下三角）a=2*eye(5)-ones(5) b=tril(a) b(:,5)=13

8、546使用rehape函数和冒号运算符创建下列矩阵。281420(a) A =4101622,(b)B =6121824reshape(2:2:24,3,4)reshape(-5:1:6,2,6) reshape(按列加)对下列每个矩阵(a)找出每个矩阵各列中绝对值最大的元素；(b)找出每个矩阵各行中绝对值最大的元素;(c)找出每个矩阵中绝对值最大的元素;(d)找出每个矩阵中绝对值最小的元素;(9 + 2)C = (3 -1)(-3 + 4)(-2 - 2)-(9 + 22)(-2 - 22)一(1 +1)，D =(3 T)(1 +1)(7 + 4) _(-3+42)(7 + 42)-2-21

9、174_ 9-2-一 2 - 2-注意：C，D和E矩阵由下列矩阵运算得到：A =31和=-11-3744C=A+BD=A+B.*abs(B)E=A./Babs(取绝对值)max(最大)min(最小)先求C：1)c=C(:,1)max(real(c)d=C(:,2)max(real(d)2)e=C(1,:)max(real(e)f=C(2,:)max(real(f)求D：1)c=D(:,1)max(real(c)d=D(:,2)max(real(d)2)e=D(1,:)max(real(e)f=D(2,:)max(real(f)求E：1)c=E(:,1)max(real(c)d=E(:,2)ma

10、x(real(d)2)e=E(1,:)max(real(e)f=E(2,:)max(real(f)思考下面的程序，并解释为什么第一条语句A*B是合法的，而第二句不 合法，这个错误是由那条语句造成的。A=ones(3,2);B=2*ones(2,3);A*B;A(2,3)=2;A*B;?Error using = = *Inner Matrix dimension must agree.解释下列程序中使用.*数组运算符造成的错误。如果不使用循环创建向 量w,怎样可以使得w(i) = u(i)*v(i)。 u=0:3;v=(3:T:0) w=u.*v;?Error using = = *Inner

11、 Matrix dimension must agree.为什么下列程序中A(k,k)和A。,i)表达式的出错信息各不相同? clear all A=ones(3,3); A(k,k)? Undefined function or variable k .A(i,i)? Index exceeds matrix dimensions.定义多项式x3 -12x2 + 25x + 6，并计算x = 2时多项式的值poly2sym(方程组)a=1 -12 25 6; p=poly2sym(a) x(2)求下面矩阵的特征多项式 TOC o 1-5 h z HYPERLINK l bookmark107 o Current Document 2 -100 HYPERLINK l bookmark110 o Current Document D 二-12-10 HYPERLINK l bookmark60 o Current Document 0-12-1 HYPERLINK l bookmark64 o Current Document 00-12D=2 -1 0 0;-1 2 -1 0;0 -1 2 -1;0 0 -1 2;PA=poly(D)PPA=poly2str(PA, s)已知一多项式的零点为2, -3, 1 +