侯绍堂matlab实验1报告

1、 数学与软件科学学院 实验报告 学期： 2014 至2015 第 1 学期 2014年10 月 23 日 课程名称:MATLAB程序设计及应用 专业: 信息与计算科学 2012 级 5、6班实验编号：1 实验项目MATLAB2010的基本功能熟悉与应用 指导教师_郭远华_ 姓名：侯绍堂 学号：2012060504 实验成绩：一、实验目的及要求1. 熟悉matlab的基本界面及其快捷键的操作；2. 熟悉常见函数和运算符操作；3. 熟悉数据类型及其运算；4.熟悉数组的各种操作.二、实验内容1. 分别以0、1、45、60、90、pi/3、pi/2、pi为弧度、角度，求余弦。以0，1/2，-1/2，1

2、，-1为余弦值，求反余弦，结果以弧度、角度表示。2. 求正切值1和-1的反正切。以（1,1），（1,-1），(-1,1)， (-1,-1)为参数，用atan2()输出相应的弧度值，并将弧度值转换为角度。比较atan（）和atan2（）的区别。3. 分别求100的自然对数，以2为底，以10为底的对数。4. 输入复数1+2i，求它的共轭值，虚部，实部，相角。5. 求1.2，1.5， 1.8，-1.2，-1.5，-1.8的fix()，floor()，ceil()，round()值。6. 以(17,3) ，（17,-3）为输入参数，调用mod()和rem()求余，并通过help或doc帮助找到它们的数

3、学原理。7. 在历史命令区双击某命令，观察结果。在历史命令区copy命令，在命令输入行修改。通过上下方向键选择输入过的命令。按Esc取消一行的输入。输入一行，按home、end键后观察光标的位置。按住Ctrl+home，Ctrl+end观察命令区的变化。8. 求200的全部素数因子。求15,25的最小公倍数，最大公因数。9. 输入magic(2),magic(3)观察输出结果的行或列的和。生成33，元素为0,10上的整数的两个随机矩阵。分别执行这两个矩阵的+，.*，.。10. 分别输入以下命令clearii=1iclear ii并说明结果产生的原因。输入以下命令观察结果type humpsx=

4、1who xwhos xwhich humpslookfor humps11. 依次输入以下命令，解释结果int8(61)/int8(2)int8(10)+10class(ans)int8(200)int8(-200)int8(1.8)int8(1.2)12. 输入以下命令观察结果3+4i3+4*ib=23+bi3+b*icomplex(6,7)A = complex(8)class(A)13. 手工输入13的逻辑向量。用两个13的向量通过关系运算生成的13逻辑向量。用一个13的向量和标量通过关系运算生成的13逻辑向量。将以上两个逻辑向量相加，解释结果。举例说明&和&的用法和区别。14. 输入

5、以下命令观察结果piPi0/01/0iskeyword15. 依次输入以下命令并解释结果A = 25；B = 17 ；bitand(A,B)bitor(A,B)bitxor(A,B)B=1；bitand(A,B)B=2；bitand(A,B)B=3；bitand(A,B)16. 输入以下命令观察ticclock date now cuptimetoc17. 分别用ones() zeros() eye() rand() randn() 创建数组18. 手工输入一个22的数组。生成一个33的随机数组，设为A。对A进行以下操作并观察结果。size(A)length(A)A(1,1)A(2,1)A(4

6、)A(1,end)A(end,1)A(1,:)A(:,1)A(:)A(4,4)=10T = A 0.5A(T)A(logical(1 0 1 0 )A(4,:) = A(:,end)=A(3,3) = A(9) = 19. 生成两个随机33数组，分别纵向横向拼接。20. 生成一个22数组A，依次输入以下命令观察并解释结果。A(1 4)A(1 4) = 10 100clear BB(1 3) = A(1 2)21. 生成34的整数数组A，元素在0,10上取值。依次输入以下命令：sum(A)max(A)min(A)sum(sum(A)sum(A(:)max(max(A)min(min(A)sum(

7、A(:)B = reshape(A,4,3)B(:)A(:)比较最后2条命令元素排列的次序。通过rand()和比较运算符随机生成一个34的逻辑数组T，通过T访问A中元素。22. 比较并解释以下命令的结果。1:2:91:2:10linspace(1,9,5)linspace(1,10,5)23. 生成3维整数数组B(3,3,2)，依次输出以下命令并解释结果。BB(:,:,1)B(:,:,2)B(1,:,:)B(:,1,:)24. 通过帮助文档查阅以下函数的功能：numel()、sub2ind()、ind2sub()，并各给出一个实例。25. 写出教材P63页的3题第（6）小题答案。三、实验步骤(

8、该部分不够填写.请填写附页)1、 x=0,1,45,60,90;y=pi/3,pi/2,pi; value1=cosd(x)value1 = 1.0000 0.9998 0.7071 0.5000 0 value2=cos(y)value2 = 0.5000 0.0000 -1.0000A=0,1/2,-1/2,1,-1; acos_value=acos(A)acos_value = 1.5708 1.0472 2.0944 0 3.1416 acosd_value=acosd(A)acosd_value = 90.0000 60.0000 120.0000 0 180.00002、 atan

9、(1) =0.7854 atan(-1) =-0.7854 atan2(1,1) =0.7854 atan2(1,-1)=2.3562 atan2(-1,1)=-0.7854 atan2(-1,-1)=-2.3562 a=0.7854,2.3562, -0.7854, -2.3562; atand_value=atand(a)atand_value = 38.1461 67.0030 -38.1461 -67.0030区别：atan()输入时是一个值，atan2（）是一个坐标3、 X=log(100),Y=log10(100)X = 4.6052Y = 24、x=1+2i; value1=an

10、gle(x),value2=conj(x),value3=imag(x),value4=real(x)value1 = 1.1071value2 = 1.0000 - 2.0000ivalue3 = 2value4 = 15、 x=1.2,1.5, 1.8,-1.2,-1.5,-1.8; v1=fix(x),v2=floor(x),v3=ceil(x),v4=round(x)v1 = 1 1 1 -1 -1 -1v2 = 1 1 1 -2 -2 -2v3 = 2 2 2 -1 -1 -1v4 = 1 2 2 -1 -2 -26、 x=17 17;y=3 -3;ans = 2 -1ans = 2

11、 27、在历史命令区双击某命令，会在命令区运行该命令。在命令区输入命令时。按Esc会直接取消取消一行的输入。随机输入一行命令，按home键光标会直接跳到命令开始、end键则会跳到命令行结尾。按住Ctrl+home，则会跳到命令区的开始，Ctrl+end则会跳到命令区的结尾。8、 factor(200)ans = 2 2 2 5 5x=gcd(15,25),y=lcm(15,25)x = 5y =759、 magic(3) ans = 8 1 6 3 5 7 4 9 2 magic(2)ans = 1 3 4 2magic(3)的每行每列的和都相等，而magic（2）只是每列相等； A=roun

12、d(10*rand(3)A = 8 0 7 10 8 8 7 9 7 B=round(10*rand(3)B = 4 7 0 7 0 1 2 3 8 A+Bans = 12 7 7 17 8 9 9 12 15 A.*Bans = 32 0 0 70 0 8 14 27 56 A2,B2ans = 113 63 105 216 136 190 195 135 170ans = 65 28 7 30 52 8 45 38 67 A.2,B.2ans = 64 0 49 100 64 64 49 81 49ans = 16 49 0 49 0 1 4 9 6410、11、int8(61)/int8

13、(2)int8(10)+10class(ans)int8(200)int8(-200)int8(1.8)int8(1.2)ans = 31ans = 20ans =int8ans = 127ans = -128ans = 2ans =1Int8是代表8位整数，所以61/2=30，int8的范围是-128-127，所以int8(200)=127，int8(-200)=-128;int代表的是整数，所以int8(1.8)=2,int8(1.2)=1;12、ans = 3.0000 + 4.0000ians = 3.0000 + 4.0000ib = 2ans = 3.0000 + 2.0000ia

14、ns = 3.0000 + 2.0000ians = 6.0000 + 7.0000iA = 8ans =double先是通过直接创建的方式来创建复数，后用complex（a,b)函数来创建复数，当complex（a)只有一个参数时，则代表实部为a,虚部为0，用函数class的返回值应为复数；13、 A=1,1,1,B=1,0,1A = 1 1 1B = 1 0 1 ABans = 0 1 014、 piPi=2.00/01/0iskeywordans = 3.1416Pi = 2ans = NaNans = Infans = break case catch classdef continu

15、e else elseif end for function global if otherwise parfor persistent return spmd switch try while变量命名是会区分大小写，pi和Pi代表不同的两个变量。0/0无意义NaN，1/0则是无穷大Inf，iskeyword可以查看系统有哪些预定义的关键字；15、 A = 25;B = 17 ;bitand(A,B)bitor(A,B)bitxor(A,B)B=1;bitand(A,B)B=2;bitand(A,B)B=3;bitand(A,B)ans = 17ans = 25ans = 8ans = 1an

16、s = 0ans = 125=11001,17=10001,bitand925,17)为与操作,ans=10001=17;bitor(25,17)为或操作，ans=11001=25；bitxor(25,17)为异或操作，ans=01000=8；同理，bitand(25,1)=1;bitand(25,2)=0;bitand(25,3)=1.16、ticclock date now cputimetocans = 1.0e+003 * 2.0140 0.0110 0.0130 0.0110 0.0380 0.0362ans =13-Nov-2014ans = 7.3592e+005ans = 34

17、.3750Elapsed time is 0.029479 seconds.17、 A=ones(2,2),B=zeros(2,2),C=rand(3,4),D=randn(3,4),E=eye(3,4)A = 1 1 1 1B = 0 0 0 0C = 0.6787 0.3922 0.7060 0.0462 0.7577 0.6555 0.0318 0.0971 0.7431 0.1712 0.2769 0.8235D = 0.3252 -1.7115 0.3192 -0.0301 -0.7549 -0.1022 0.3129 -0.1649 1.3703 -0.2414 -0.8649 0

18、.6277E = 1 0 0 0 0 1 0 0 0 0 1 0Zeros(m,n),ones(m,n)分别产生填充值为0或1的m行n列的数组；rand（m,n)和rands（m,n)分别产生填充值为01之间的的随机数组，eye（m,n）产生当行等于列时，取值为1，其余时候为0的数组。18、 B=3,4;5,6B = 3 4 5 6 A=round(10*rand(3,3)A = 7 7 1 8 7 5 3 2 10 size(A)length(A)A(1,1)A(2,1)A(4)A(1,end)A(end,1)A(1,:)A(:,1)A(:)A(4,4)=10T = A 0.5A(T)A(l

19、ogical(1 0 1 0 )A(4,:) = A(:,end)=A(3,3) = A(9) = ans = 4 4ans = 4ans = 7ans = 8ans = 0ans = 0ans = 0ans = 7 7 1 0ans = 7 8 3 0ans = 7 8 3 0 7 7 2 0 1 5 10 0 0 0 0 10A = 7 7 1 0 8 7 5 0 3 2 10 0 0 0 0 10T = 1 1 1 0 1 1 1 0 1 1 1 0 0 0 0 1ans = 7 8 3 7 7 2 1 5 10 10ans = 7 3A = 7 7 1 0 8 7 5 0 3 2 10

20、 0A = 7 7 1 8 7 5 3 2 10? Subscripted assignment dimension mismatch.A(3,3）与A（9）不能运行是因为数组只能进行一行或者一列的删除，不能使某一个元素为空。19、 A=rand(3,3),B=rand(3,3)A = 0.3404 0.7513 0.6991 0.5853 0.2551 0.8909 0.2238 0.5060 0.9593B = 0.5472 0.2575 0.8143 0.1386 0.8407 0.2435 0.1493 0.2543 0.9293 C=A BC = 0.3404 0.7513 0.69

21、91 0.5472 0.2575 0.8143 0.5853 0.2551 0.8909 0.1386 0.8407 0.24350.2238 0.5060 0.9593 0.1493 0.2543 0.929320、A=rand(2,2)A = 0.3500 0.2511 0.1966 0.6160 A(1 4)A(1 4) = 10 100clear BB(1 3) = A(1 2)ans = 0.3500 0.6160A = 10.0000 0.2511 0.1966 100.0000B = 10.0000 0 0.196621、A=round(10*rand(3,4)A = 5 6 3

22、 4 4 5 8 6 8 9 8 1 sum(A)max(A)min(A)sum(sum(A)sum(A(:)max(max(A)min(min(A)sum(A(:)B = reshape(A,4,3)B(:)A(:)ans = 17 20 19 11ans = 8 9 8 6ans = 4 5 3 1ans = 67ans = 67ans = 9ans = 1ans = 67B = 5 5 8 4 9 4 8 3 6 6 8 1ans = 5 4 8 6 5 9 3 8 8 4 6 1ans = 5 4 8 6 5 9 3 8 8 4 6 1 T=rand(3,4)0.5T = 0 1 0

23、0 1 0 0 1 1 1 0 0A(T)ans = 4 8 6 9 622、 1:2:91:2:10linspace(1,9,5)linspace(1,10,5)ans = 1 3 5 7 9ans = 1 3 5 7 9ans = 1 3 5 7 9ans =1.0000 3.2500 5.5000 7.7500 10.0000a:b:c表示取值在ab，以b为公差的等差数列，linspace（a,b,c)表示创建一个以a 为下限，b 为上限，c表示向量中元素的个数。23、BB(:,:,1)B(:,:,2)B(1,:,:)B(:,1,:)B(:,:,1) = 0.5285 0.2630 0.

24、7482 0.1656 0.6541 0.4505 0.6020 0.6892 0.0838B(:,:,2) = 0.2290 0.8258 0.0782 0.9133 0.5383 0.4427 0.1524 0.9961 0.1067ans = 0.5285 0.2630 0.7482 0.1656 0.6541 0.4505 0.6020 0.6892 0.0838ans = 0.2290 0.8258 0.0782 0.9133 0.5383 0.4427 0.1524 0.9961 0.1067ans(:,:,1) = 0.5285 0.2630 0.7482ans(:,:,2) = 0.2290 0.8258 0.0782ans(:,:,1) = 0.5285 0.1656 0.6020ans(:,:,2) = 0.2290 0.9133 0.152424、 A=1,2,3A = 1 2 3 numel(A)ans = 325、vec=4 5 2 8 4 7 2 64 2 5