Matlab 上机实验 数组运算 求阶乘 源程序代码.doc

Matlab实验1
通信11002班李俊良 学号：201006113

第一题 数组的加减乘除和乘方
ü A=[1234]，B=[3456]，求C=A+B，D=A-B，E=A*B，
ü F=A/B，G=A^B
ü 用stem语句画出A，B，C，D，E，F，G。
解:
%ʵÑéÄÚÈÝÒ»Êý×éµÄ¼Ó¼õ³Ë³ýjjcc
clear;clc;%ɾ³ýwordspaceÖеľØÕó²¢½øÐÐÇåÆÁ
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
figure,subplot(2,2,1),stem(A)
title('AµÄÐòÁÐ')
xlabel('Ax'); ylabel('Ay');
subplot(2,2,2),stem(B)
title('BµÄÐòÁÐ')
xlabel('Bx'); ylabel('By');
subplot(2,2,2),stem(C)
title('CµÄÐòÁÐ')
xlabel('Cx'); ylabel('Cy');
subplot(2,2,3),stem(D)
title('DµÄÐòÁÐ')
xlabel('Dx'); ylabel('Dy');
subplot(2,2,4),stem(G)
title('GµÄÐòÁÐ')
xlabel('Gx'); ylabel('Gy');

图1第一题答案

第二题 用MATLAB实现下列序列：

解：
程序如下：
%ʵÑéÄÚÈÝER»­ÐòÁÐ
clear;clc;%ɾ³ýwordspaceÖеľØÕó²¢½øÐÐÇåÆÁ
n=[0:15];%¶¨Òå×Ô±äÁ¿n´Ó0µ½15
x1=0.8.^n;
x2=exp((0.2+3j)*n);
x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);

figure,
%»­µÚÒ»¸ö×Óͼ
subplot(3,1,1),stem(x1) %»­³öx1µÄÐòÁÐͼ
title('x1=0.8^n')
xlabel('n'); ylabel('x1');
%»­µÚ¶þ¸ö×Óͼ
subplot(3,1,2),stem(x2) %»­³öx2µÄÐòÁÐͼ
title('x2=exp((0.2+3j)*n)')
xlabel('n'); ylabel('x2');
%»­µÚÈý¸ö×Óͼ
subplot(3,1,3),stem(x3) %»­³öx3µÄÐòÁÐͼ
title('x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)')
xlabel('n'); ylabel('x3');

图2 第二题图

第三题 举例生成一个3×3的矩阵，分别用
ü 逐个输入
ü 利用冒号表达式“:”生成向量
ü 利用函数生成向量
解：%ʵÑéÄÚÈÝsanÉú³É3*3µÄ¾ØÕó
clear;clc;%ɾ³ýwordspaceÖеľØÕó²¢½øÐÐÇåÆÁ
A=[1 2 3;4 5 6;7 8 9]%Öð¸öÊäÈë·¨

B=[1:3;4:6;7:9]%:Éú³É·¨

C=reshape(A,3,3)%ÀûÓú¯ÊýÉú³ÉÏòÁ¿·¨
运行结果：

A =

1 2 3
4 5 6
7 8 9

B =

1 2 3
4 5 6
7 8 9

C =

1 2 3
4 5 6
7 8 9

第四题 利用for循环求1！+2！+3！+ µµµ +5！的值
解：%ʵÑéÄÚÈÝsiÇóÇ°ÎåÏî½×³ËÖ®ºÍ
clear;clc;%ɾ³ýwordspaceÖеľØÕó²¢½øÐÐÇåÆÁ
jc=1;sum=0;%¶¨Òå½×³ËºÍÇóºÍµÄ³õʼֵ

for n=1:5;
jc=jc*n;%·Ö±ðÇó³önµÄ½×³Ë
sum=sum+jc;%¶ÔËùÇóµÄ½×³Ë½øÐÐÀÛ¼Ó
end
sum=sum%Êä³ö×îÖÕ½á¹û
运行结果：
sum =

153

第五题 已知矩阵A=[-1 5 -4;0 7 8; 3 61 7],b=[8 3 -1;2 5 3; -3 2 0]
ü 求下列表达式的值
ü A+6B和A*A-B+I (I为单位矩阵)
ü A*B,A.*B,B*A
ü A/B,B\A

%ʵÑéÄÚÈÝ5¾ØÕóABºúÂÒÔËËã
clear;clc;%ɾ³ýwordspaceÖеľØÕó²¢½øÐÐÇåÆÁ
A=[-1 5 -4;0 7 8; 3 61 7];B=[8 3 -1;2 5 3; -3 2 0];
I=eye(length(A));%¶¨Ò嵥λ¾ØÕóI£¬Æä´óС¸úAÒ»ÖÂ
z1=A+6*B
z2=A*A-B+I
z3=A*B
z4=A.*B
z5=B*A
z6=A/B
z7=B\A

运行结果：
z1 =

47 23 -10
12 37 26
-15 73 7

z2 =

-18 -217 17
22 533 109
21 867 526

z3 =

14 14 16
-10 51 21
125 328 180

z4 =

-8 15 4
0 35 24
-9 122 0

z5 =

-11 0 -15
7 228 53
3 -1 28

z6 =

1.2234 -0.9255 2.9787
-0.9468 2.3511 -0.9574
4.6170 3.8723 13.8936
第六题A=[1 -1 2 3; 5 1 -4 2;3 0 5 2;11 15 0 9]
ü 求A的对角阵，上三角，下三角，逆矩阵，旋转90度，秩
解：%ʵÑéÄÚÈÝ5¾ØÕóABºúÂÒÔËËã
clear;clc;%ɾ³ýwordspaceÖеľØÕó²¢½øÐÐÇåÆÁ
A=[1 -1 2 3; 5 1 -4 2;3 0 5 2;11 15 0 9] %¶¨Òå¾ØÕóA
I=ones(length(A));%¶¨Ò嵥λ£¨È«1£©¾ØÕóI£¬Æä´óС¸úAÒ»ÖÂ
z1=diag(A)%AµÄ¶Ô½ÇÕó
z2=triu(A)%AµÄÉÏÈý½ÇÕó
z3=tril(A)%AµÄÏÂÈý½ÇÕó
z4=inv(A)%AµÄÄæ¾ØÕó
z5=rot90(A)%AÐýת90¶È
z6=rank(A)%AµÄÖÈ

运行结果：
A =

1 -1 2 3
5 1 -4 2
3 0 5 2
11 15 0 9

z1 =

1
1
5
9

z2 =

1 -1 2 3
0 1 -4 2
0 0 5 2
0 0 0 9

z3 =

1 0 0 0
5 1 0 0
3 0 5 0
11 15 0 9

z4 =

-0.1758 0.1641 0.2016 -0.0227
-0.1055 -0.1016 -0.0391 0.0664
-0.0508 -0.0859 0.1516 0.0023
0.3906 -0.0313 -0.1813 0.0281

z5 =

3 2 2 9
2 -4 5 0
-1 1 0 15
1 5 3 11

z6 =

4

通信11002班李俊良 学号：201006113