数学软件与实验第一次上机作业上机时间：2013-4-10 地点：E204班级:071111 学号:07111014 姓名:曹红兴xdhjtang@学号、姓名、MATLAB、第一次作业1.计算三角形三边分别为a，b，c中c边对应内角的角度>> a = 3; b = 3; c = 3;>> acos((a^2+b^2-c^2)/2/a/b) *180/pians =60.0000>> a = 3; b = 4; c = 5;>> acos((a^2+b^2-c^2)/2/a/b) *180/pians =90>> a = 3; b = 4; c = 20;>> acos((a^2+b^2-c^2)/2/a/b) *180/pians =1.8000e+002 -1.9715e+002i2.试分别生成5 阶的单位阵、8 阶均匀分布的随机矩阵及其下三角矩阵，要求矩阵元素为介于10~99之间整数>> C=eye(5,5)C =1 0 0 0 00 1 0 0 00 0 1 0 00 0 0 1 00 0 0 0 1>> N=randsrc(8,8,[10:99])N =59 21 72 34 19 76 25 5266 28 54 28 22 15 45 6812 23 58 60 24 87 22 1265 27 50 67 65 94 12 8542 13 21 47 61 98 94 6014 67 54 28 14 87 37 8654 35 86 95 93 80 36 4127 58 88 17 75 56 39 50>> Z=tril(N)Z =83 0 0 0 0 0 0 091 96 0 0 0 0 0 021 24 81 0 0 0 0 092 97 96 45 0 0 0 066 96 69 68 72 0 0 018 53 13 25 38 54 0 035 82 86 73 95 50 20 059 22 94 12 13 68 54 723.生产列向量x=[1, 3, 10, 15, 20, 25, 30, 35, 40] >> x=[1;3;10;15;20;25;30;35;40]x =13101520253035404.生成以x 的元素为对角线的矩阵A，并输出A 的行数>> A=diag(x)a=length(A)A =1 0 0 0 0 0 0 0 00 3 0 0 0 0 0 0 00 0 10 0 0 0 0 0 00 0 0 15 0 0 0 0 00 0 0 0 20 0 0 0 00 0 0 0 0 25 0 0 00 0 0 0 0 0 30 0 00 0 0 0 0 0 0 35 00 0 0 0 0 0 0 0 40a =95.生成一个与A 同阶的正态分布的随机矩阵B, 要求B矩阵元素为介于0~99之间整数>> B=round(0+rand(9)*(99-0))B =75 10 79 65 25 83 47 38 4674 82 19 16 50 25 35 56 139 69 48 12 69 81 82 8 3365 31 44 49 88 24 58 5 1617 94 64 95 95 92 54 53 7970 3 70 34 54 35 91 77 313 43 75 58 14 19 28 92 5227 38 27 22 15 25 75 13 165 76 67 74 25 61 75 56 606.输出A 与B 的乘积矩阵C>> C=A*B>> C=A*BC =5 996 10 89 30 50 82 7254 84 156 222 21 135 225 99 210660 410 330 730 240 100 620 290 770495 690 255 840 75 1485 135 1110 4351780 1520 420 360 880 660 160 20 1380300 2025 2250 1475 25 725 1925 125 13752940 300 2010 900 2670 180 **** **** 11701855 630 1610 455 665 1050 1855 2100 2102800 1440 3600 840 360 200 440 2080 30807.生成由A 与B 点乘得到的矩阵DD=A.*BD =5 0 0 0 0 0 0 0 00 84 0 0 0 0 0 0 00 0 330 0 0 0 0 0 00 0 0 840 0 0 0 0 00 0 0 0 880 0 0 0 00 0 0 0 0 725 0 0 00 0 0 0 0 02700 0 00 0 0 0 0 0 0 2100 00 0 0 0 0 0 0 0 30808.生成一个由D 的第8、4、2、5 行和第7、1、6、9、2 列组成的子矩阵 E>> E=D([8 4 2 5],[7 1 6 9 2])E =0 0 0 0 00 0 0 0 00 0 0 0 2460 0 0 0 09.求出矩阵E 的最大元素>> m=max(max(E))m =24610.计算的值>> 1.369^2+sin(7*pi/10)*sqrt(26.48)/2.9ans =3.309711.产生一个5阶魔术方阵，并执行如下操作：a)将矩阵的第2行3列元素赋值给变量cb)将由矩阵第2，3，4行第3，5列构成的子矩阵赋值给变量d>> Q=magic(5)c=Q(2,3)d=[Q(2:4,3),Q(2:4,5)]Q =17 24 1 8 1523 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9c =7d =7 1613 2219 312.给出区间[0，1]上的6个等分点数据。
>> L=0:1/5:1L =0 0.2000 0.4000 0.6000 0.8000 1.000013.建立如下矩阵1)错误！未找到引用源。
2)错误！未找到引用源。
>> R=3.14*eye(6)R =3.1400 0 0 0 0 00 3.1400 0 0 0 00 0 3.1400 0 0 00 0 0 3.1400 0 00 0 0 0 3.1400 00 0 0 0 0 3.1400 >> A1=round(10+rand(6)*(10-10))A1(1,1)=0,A1(2,2)=0,A1(3,3)=0,A1(4,4)=0,A1(5,5)=0,A1(6,6)=0A=A1A1 =10 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 10A1 =0 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 10A1 =0 10 10 10 10 1010 0 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 10 A1 =0 10 10 10 10 1010 0 10 10 10 1010 10 0 10 10 1010 10 10 10 10 1010 10 10 10 10 1010 10 10 10 10 10 A1 =0 10 10 10 10 1010 0 10 10 10 1010 10 0 10 10 1010 10 10 0 10 1010 10 10 10 10 1010 10 10 10 10 10 A1 =0 10 10 10 10 1010 0 10 10 10 1010 10 0 10 10 1010 10 10 0 10 1010 10 10 10 0 1010 10 10 10 10 10 A1 =0 10 10 10 10 1010 0 10 10 10 1010 10 0 10 10 1010 10 10 0 10 1010 10 10 10 0 1010 10 10 10 10 0A =0 10 10 10 10 1010 0 10 10 10 1010 10 0 10 10 1010 10 10 0 10 1010 10 10 10 0 1010 10 10 10 10 014.已知函数错误！未找到引用源。
,计算f(-1), f(0.5), f(1.5)。
function y=f(x)n=length(x);for k=1:nif x(k)>=1y(k)=x(k)^2;elseif x(k)>=0y(k)=1;elsey(k)=x(k)+1;endend>> x=[-1,0.5,1.5] y=f(x)x =-1.0000 0.5000 1.5000y =0 1.0000 2.250015.用for-end循环语句求错误！未找到引用源。
和100!。
16.用while-end循环语句求不超过1000的偶数之和与奇数之和。
偶数之和：>> s=0;>> i=0;>> while i<=1000s=s+i;i=i+2;end>> ss =250500>>奇数之和：>> s=0;>> i=1;>> while i<=1000s=s+i;i=i+2;end>> ss =250000>>17.建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。
例如，153是一个水仙花数，因为错误！未找到引用源。
班级学号姓名数学软件与实验-1。