MATLAB刘卫国版本

海 南 大 学

MATLAB 实验 (刘卫国版本)

题 目： MATLAB 实验 学 号： 2009050******

姓 名： **** 年 级： 09电气 学 院： 机电工程学院 系 别： 电气工程系 专 业： 电气工程及其自动化

1

实验一 MATLAB运算基础

二（1）二、实验内容

1.先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 ⑴z2sin85?1?1?e2；

程序为>> z1=2*sin(85*pi/180)/(1+exp(2)) z1 =

0.2375 ⑵z12ln(x?1?x2)，其中x???21?2i?2???0.455??； >> x=[2 1+2i;-0.45 5];

>> z2=1/2*log(x+sqrt(1+x^2)) z2 =

0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i

⑶ze0.3a?e?0.3a3?2sin(a?0.3)?ln0.3?a2,a??3.0,?2.9,?2.8,?,2.8,2.9,3.0程序为>> a=(-3.0:0.1:3.0);

>> z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) z3 =

Columns 1 through 3

0.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i

Columns 4 through 6

0.7913 + 3.1416i 0.7822 + 3.1416i 0.7602 + 3.1416i

2

Columns 7 through 9

0.7254 + 3.1416i 0.6784 + 3.1416i 0.6196 + 3.1416i

Columns 10 through 12

0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780 + 3.1416i

Columns 13 through 15

0.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i

Columns 16 through 18

-0.0771 + 3.1416i -0.2124 + 3.1416i -0.3566 + 3.1416i

Columns 19 through 21

-0.5104 + 3.1416i -0.6752 + 3.1416i -0.8536 + 3.1416i

Columns 22 through 24

-1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416i

Columns 25 through 27

-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i

Columns 28 through 30

-37.0245 -3.0017 -2.3085

Columns 31 through 33

-1.8971 -1.5978 -1.3575

Columns 34 through 36

-1.1531 -0.9723 -0.8083

Columns 37 through 39

-0.6567 -0.5151 -0.3819

3

Columns 40 through 42

-0.2561 -0.1374 -0.0255

Columns 43 through 45

0.0792 0.1766 0.2663

Columns 46 through 48

0.3478 0.4206 0.4841

Columns 49 through 51

0.5379 0.5815 0.6145

Columns 52 through 54

0.6366 0.6474 0.6470

Columns 55 through 57

0.6351 0.6119 0.5777

Columns 58 through 60

0.5327 0.4774 0.4126

Column 61

0.3388

?t2 ⑷z?0?t?1?t24??1 1?t?2，其中t?0:0.5:2.5

??t2?2t?1,2?t?3程序为>> t=(0:0.5:2.5);

>> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(2<=t&t<3).*(t.^2-2.*t+1) z4 =

0 0.2500 0 1.2500 1.0000 2.2500

4

?1234?4??13?1?????，求下列表达式的值：

32.已知：A?34787,B?20???????3657???3?27?? ⑴A?6*B和A?B?I(其中I为单位矩阵)；⑵A*B和A.*B; ⑶A^3和A.^3;⑷A/B和A\\B;⑸[A,B]和[A([1,3]),:);B^2]. 程序为>> A=[12 34 -4;34 7 87;3 65 7];

>> B=[1 3 -1;2 0 3;3 -2 7]; >> I=[1 0 0;0 1 0;0 0 1]; >> A+6*B

ans =

18 52 -10 46 7 105 21 53 49

>> A-B+I

ans =

12 31 -3 32 8 84 0 67 1

>> A^3

ans =

37226 233824 48604 247370 149188 600766 78688 454142 118820

>> A.^3

ans =

1728 39304 -64 39304 343 658503 27 274625 343

>> A/B

5

ans =

16.4000 -13.6000 7.6000 35.8000 -76.2000 50.2000 67.0000 -134.0000 68.0000

>> B\\A

ans =

109.4000 -131.2000 322.8000 -53.0000 85.0000 -171.0000 -61.6000 89.8000 -186.2000

>> [A,B]

ans =

12 34 -4 1 3 -1 34 7 87 2 0 3 3 65 7 3 -2 7

>> [A([1,3],:);B^2]

ans =

12 34 -4 3 65 7 4 5 1 11 0 19 20 -5 40

?12345??30?678910??17?6???3.设有矩阵A和B：A??1112131415?,B??023???1617181920???97???2122232425???413⑴求它们的乘积C；⑵将矩阵C的右下角3?2子矩阵赋给D；

16?9???4? ?0?11??⑶查看MATLAB工作空间的使用情况。

程序为>> A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25]; >> B=[3 0 16;17 -6 9 ;0 23 -4;9 7 0;4 13 11]; >> C=A*B

C = 93 150 77

6

258 335 237 423 520 397 588 705 557 753 890 717 >> D=C([3,4,5],[2,3])

D = 520 397 705 557 890 717 >> whos

Name Size Bytes Class Attributes A 5x5 200 double B 5x3 120 double C 5x3 120 double D 3x2 48 double 4.完成下列操作：

⑴求[100,999]之间能被21整除的数的个数； ⑵建立一个字符串向量，删除其中的大写字母。 程序为>> A=100:1:999; K=find(rem(A,21)==0); length(K) ans = 43

>> ch='aHDJLK143663CFHI'; >> H=find(ch>='A'&ch<'Z'); >> ch(H)=[ ] ch = a143663

实验二 MATLAB矩阵分析与处理二（1）E=eye(3);

R=rand(3,2); O=zeros(2,3); S=diag([2,3]); A=[E,R;O,S]; A^2

B=[E,(R+R*S);O,S^2] ans =

1.0000 0 0 2.8504 1.9439

7

0 1.0000 0 0.6934 3.5652 0 0 1.0000 1.8205 3.0484 0 0 0 4.0000 0 0 0 0 0 9.0000 B =

1.0000 0 0 2.8504 1.9439 0 1.0000 0 0.6934 3.5652 0 0 1.0000 1.8205 3.0484 0 0 0 4.0000 0 0 0 0 0 9.0000

（2）>> H=hilb(5) P=pascal(5) Hh=det(H) Hp=det(P) Th=cond(H) Tp=cond(P) H =

1.0000 0.5000 0.3333 0.2500 0.2000 0.5000 0.3333 0.2500 0.2000 0.1667 0.3333 0.2500 0.2000 0.1667 0.1429 0.2500 0.2000 0.1667 0.1429 0.1250 0.2000 0.1667 0.1429 0.1250 0.1111 P =

1 1 1 1 1 1 2 3 4 5 1 3 6 10 15 1 4 10 20 35 1 5 15 35 70 Hh =

3.7493e-012

8

Hp = 1 Th =

4.7661e+005 Tp =

8.5175e+003

（3）>> A=fix(10*rand(5)) H=det(A) Trace=trace(A) Rank=rank(A) Norm=norm(A) A =

7 8 0 7 6 4 5 6 4 3 9 2 3 3 5 4 6 8 1 1 4 8 5 1 6 H =

-6498

Trace = 22 Rank = 5

9

Norm =

23.8478

（4）>> A=[-29,6,18;20,5,12;-8,8,5] [V,D]=eig(A) %数学意义略 A =

-29 6 18 20 5 12 -8 8 5 V =

0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050 D =

-25.3169 0 0 0 -10.5182 0 0 0 16.8351

（5）A=hilb(4) A(:,1)=[] A(4,:)=[]

B=[0.95,0.67,0.52]'; X=inv(A)*B

B1=[0.95,0.67,0.53]'; X1=inv(A)*B1 N=cond(B) N1=cond(B1)

Na=cond(A) %矩阵A为病态矩阵 A =

1.0000 0.5000 0.3333 0.2500

10

0.5000 0.3333 0.2500 0.2000 0.3333 0.2500 0.2000 0.1667 0.2500 0.2000 0.1667 0.1429 A =

0.5000 0.3333 0.2500 0.3333 0.2500 0.2000 0.2500 0.2000 0.1667 0.2000 0.1667 0.1429 A =

0.5000 0.3333 0.2500 0.3333 0.2500 0.2000 0.2500 0.2000 0.1667 X =

1.2000 0.6000 0.6000 X1 =

3.0000 -6.6000 6.6000 N = 1 N1 = 1

11

Na =

1.3533e+003

（6）A=[1,4,9;16,25,36;49,64,81] B=sqrtm(A)

C=sqrt(A) %sqrtm函数是以矩阵为单位进行计算，sqrt函数是以矩阵中的元素进行计算 A =

1 4 9 16 25 36 49 64 81 B =

0.6344 + 1.3620i 0.3688 + 0.7235i 0.7983 - 0.4388i 1.4489 + 1.1717i 2.7697 + 0.6224i 3.2141 - 0.3775i 4.3578 - 1.6237i 5.7110 - 0.8625i 7.7767 + 0.5231i C =

1 2 3 4 5 6 7 8 9

实验三 选择程序结构设计

二、（1） x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0]; y=[]; %建立存放所有y值的矩阵 for x0=x

if x0<0&x0~=-3 y=[y,x0*x0+x0-6];

elseif x0>=0&x0<5&x0~=2&x0~=3 y=[y,x0*x0-5*x0+6]; else

y=[y,x0*x0-x0-1]; end end

x %输出所有x y

12

x =

-5.0000 -3.0000 1.0000 2.0000 2.5000 3.0000 5.0000 y =

14.0000 11.0000 2.0000 1.0000 -0.2500 5.0000 19.0000

（2）x=input('请输入一个百分制成绩：'); if x>100|x<0

disp('您输入的成绩不是百分制成绩，请重新输入。'); else

if x<=100&x>=90 disp('A');

elseif x<=89&x>=80 disp('B');

elseif x<=79&x>=70 disp('C');

elseif x<=69&x>60 disp('D'); else disp('E'); end end

请输入一个百分制成绩：89 B

（3）n=input('请输入员工工号：'); h=input('该员工工作时数是：'); if h>120

x=(h-120)*84*(1+0.15)+120*84; elseif h