matlab上机作业

第一次 上机作业

1．写出完成下列操作的命令：

1) 将矩阵A 的第2～5行中的第1，3，5列元素赋给矩阵B

B=A(:,1:3:5)

2) 删除矩阵A 的第7号元素

A(7)=[]

3) 将矩阵A 的每个元素值加30

A=A+30

4) 求矩阵A 的大小和维数

Whos size(A)

5) 将含有12个元素的向量x 转换成3×4矩阵

Y=reshape(x,3,4)

2．下列命令执行后，L1、L2、L3、L4的值分别是多少？

A=1:9; B=10－A;

L1=A= =B;

L2=A< =5;

L3=A>3&A<7;

L4=find(A>3&A<7);

L1= 0 0 0 0 1 0 0

L2=1 1 1 1 1 0 0

L3=0 0 0 1 1 1 0

L4=4 5 6

3．已知：

A=⎥⎥⎥⎥⎦

⎤⎢⎢⎢⎢⎣⎡---14.35454.9632053256545410778.01023

完成下列操作：

(1)取出A的前3行构成矩阵B，前两列构成矩阵C，右下角3×2子矩阵构成矩

阵D，B与C的乘积构成矩阵E

B=A(1:3,:) C=A(:,1:2) D=A(2:4,3:4) E=B*C

(2)分别求E=10&A<25)。

1.E

0 1

0 0

0 1

2.E&D

1 1

0 1

1 1

3.E|D

1 1

1 1

1 1

4.～E|～D

0 0

1 0

0 0

5.find(A>=10&A<25)

1

5

4．用结构体矩阵来储存5名学生的基本情况数据，每名学生的数据包括学号、姓名、专业和6门课程的成绩。

例如

a(1).x1=1;a(1).x2=’name’;a(1).x3=’grade’;a(1).x4=’number’;a(1).x5=[XX,XX,X];

5．建立单元矩阵B并回答有关问题。

B{1,1}=1;

B{1,2}=`Brenden`;

B{2,1}=reshape(1:9,3,3);

B{2,2}={12,34,2; 54,21,3; 4,23,67};

(1)size(B)和ndims(B)的值分别是多少？

(2)B(2)和B(4)的值分别是多少？

(3)B(3)=[]和B{3}=[]执行后，B的值分别是多少？

1.size(B)

ans=

4 3

ndlims(B)

ans=

2

2.>>B(2)

ans=

7

B(4)

ans=

1

3.B(3)=[]

B=

1 7 1 11 8 5

2 12 9 6 3

6．在[0,3π]区间，求y=sin(x)的值。要求：

(1)消去负半波，即（π,2π）区间内的函数值置0。

(2) （π/3, 2π/3 ）和（7π/3, 8π/3 ）区间内取值均为sin π/3 。

提示：

先根据自变量向量x产生函数值向量y，然后按要求对y进行处理。处理的

思路有两个：一是从自变量着手进行处理，二是从函数值着手进行处理。

(1) x=0:pi/100:3*pi

y=sin(x)

y1=(x2*pi).*y

(2)q=(x>pi/3&x<2*pi/3)|(x<7*pi/3&x<8*pi/3)

qn=-q

y2=q*sin(pi/3)+qn.*y1

7．写出完成下列操作的命令：

1)建立3阶单位矩阵。

a=eye(3)

a=

1 0 0

0 1 0

0 0 1

2)建立5×6随机矩阵A，其元素为[100，200]范围内的随机整数。

A=100+(200-100)*rank(5,6)

3)产生均值为1，方差为0.2的500个正态分布的随机数。

y=1+sqrt(0.2)*randn(50,10)

4)产生和A同样大小的幺矩阵。

ones(size(A))

5)将矩阵A对角线的元素加30。

B=diag(A)+30

6)从矩阵A提取主对角线元素，并以这些元素构成对角阵B.

B=diag(A)

8．使用函数，实现方阵左旋900或右旋900的功能。例如，原矩阵为A, A左旋后得到B, 右旋后得到C。

⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=129631185210741A ，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=321654987121110B ，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=1011

12789456123C A=[1,3,4,10;2,5,8,11;3,6,9,12] >>B=rot90(A)

B=[10,11,12;7,8,9;4,5,6;1,2,3;] >>C=rot90(A,3)

9．求下面线性方程组的解 ⎪⎩

⎪⎨⎧=+=+-=-+83121023224y x z y x z y x

A=[4,2,-1;3,-1,2;12,3,0] B=[2,10,8]’;

x=inv(A)*B

x=

-6.0000

26.6667

27.3333