MATLAB程序设计基础 作业1

MATLAB程序设计基础 作业1

1. 求下列表达式的值，观察MATLAB工作空间使用情况并保存全部变量

2sin85?(1) z1?

1?e2?21?2i?12?(2) z2?ln(x?1?x) 其中x???0.45? 25??e0.3a?e?0.3a0.3?asin(a?0.3)?ln(3) z3?，a??3,?2.9,?2.8,?,2.8,2.9,3

22（提示：用冒号表达式生成a向量，求各点函数值时用点乘运算）

?1234?4??13?1?????3?，求下列表达式的值： 2. 已知A??34787?，B??20?3?27??3657?????(1) A?6*B 和A?B?I（其中I为单位阵） (2) A*B 和A.*B (3) A^3 和A.^3 (4) A/B 和B\\A

(5) [A,B] 和[A([1,3],:);B^2] 3. 设有矩阵A和B

?1234??6789A??11121314??16171819?21222324?5?016??3???10?17?69??15?，B??023?4? ???20?70??9?41311?25???? (1) 求它们的乘积C。

(2) 将矩阵C的右下角3×2子矩阵赋给D。 (3) 查看MATLAB工作空间的使用情况。 4. 完成下列操作：

(1) [100,999]之间能被21整除的数的个数（提示：先利用冒号表达式，再利用find和length函数，通过Help查询用法）。

(2) 建立一个字符串向量，删除其中的大写字母（提示：利用find函数和空矩阵）。

?I3?3?A?5. 设有分块矩阵?0?2?3R3?2??，其中I,R,0,S分别为单位矩阵，随机矩阵，S2?2???IR?RS??零矩阵和对角矩阵，试通过数值计算验证A??2?0?。 S??21．(1)

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

z1 = 0.2375

----------------------------------------------------------------------------------------------------------- (2)

>> x=[2 1+2i;0.45 5] x =

2.0000 1.0000 + 2.0000i 0.4500 5.0000 >> z2=log10(x+sqrt(1+x*x))/2 z2 =

0.3191 + 0.0096i 0.3895 + 0.1589i 0.1979 0.5031 + 0.0019

------------------------------------------------------------------------------------------------------------ (3) a=-0.3:0.1:3

z3=(exp(0.3*a)-exp(-0.3*a).*sin(a+0.3))/2+log10(0.3+a)/2 2．(1) >> a+6*b ans =

18 52 -10 46 7 105 21 53 49 >> a-b+eye(3) ans =

12 31 -3 32 8 84 0 67 1

----------------------------------------------------------------------------------------------------------

(2) >> a*b ans =

68 44 62 309 -72 596 154 -5 241

>> a.*b ans =

12 102 4 68 0 261 9 -130 49

(3) >> a^3 ans =

37226 233824 48604 247370 149188 600766 78688 454142 118820

>> a.^3 ans =

1728 39304 -64 39304 343 658503 27 274625 343

------------------------------------------------------------------------------------------------------------

(4) >> a/b 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

------------------------------------------------------------------------------------------------------------

(5) >> [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

================================================================

3. (1) >> c=a*b c =

93 150 77 258 335 237 423 520 397 588 705 557 753 890 717

----------------------------------------------------------------------------------------------- (2) >> d=c(3:5,2:3) d =

520 397 705 557

890 717 (3)略

================================================================

4. (1) >> n=100:999

length(find(rem(n,21)==0)) ans = 43

(2) >> str='dengZWingM#*%VC509'

>> index=find(str>='A'&str<='Z') >> str(index)=0 str =

deng ing #*% 509

================================================================

5. >> i=eye(3)

>> r=rand(3,2) >> o=zeros(2,3) >> s=diag([1 2]) >> a=[i r;o s] a =

1.0000 0 0 0.8147 0.9134 0 1.0000 0 0.9058 0.6324 0 0 1.0000 0.1270 0.0975 0 0 0 1.0000 0 0 0 0 0 2.0000 >> a*a ans =

1.0000 0 0 1.6294 2.7401 0 1.0000 0 1.8116 1.8971 0 0 1.0000 0.2540 0.2926 0 0 0 1.0000 0

0 0 0 0 4.0000 >> [i,r+r*s;o,s*s] ans =

1.0000 0 0 1.6294 2.7401 0 1.0000 0 1.8116 1.8971 0 0 1.0000 0.2540 0.2926 0 0 0 0

0 1.0000 0 0 0 4.0000

本文来源：https://www.bwwdw.com/article/om8f.html