Matlab数学实验报告.doc

数学实验报告 姓名： 班级： 学号： 第一次实验任务 过程： a=1+3i; b=2-i; 结果： a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程： x=-4.5*pi/180; y=7.6*pi/180; 结果： sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得：对于matlab中的角度计算应转为弧度。 （1）过程： x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果： （2）过程： subplot(2,2,1) plot(x,y1) subplot(2,2,2) plot(x,y2) subplot(2,2,3) plot(x,y3) subplot(2.2.4) subplot(2,2,4) plot(x,y4) 结果： 心得：在matlab中，用subplot能够实现在同一页面输出多个坐标系的图像，应注意将它与hold on进行区别，后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B，计算（1）AB，（2）对B中每个元素平方后得到的矩阵C，（3）sinB，（4） A的行列式，（5）判断A是否可逆，若可逆，计算 A的逆矩阵，（6）解矩阵方程AX=B，（7）矩阵A中第二行元素加1，其余元素不变，得到矩阵D,计算D。 过程： A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10); 结果： (1) A*B = 22 28 49 64 76 100 (2) C=B.^2 C = 1 4 9 16 25 36 (3) sin(B) ans = 0.8415 0.9093 0.1411 -0.7568 -0.9589 -0.2794 (4) a=det(A) a = 0 (5) 由det(A)存在，且det（A）= 0，故A不可逆 (6) X无解 (7) l=A(2,1),l1=l+1,m=A(2,2),m1=m+1,n=A(2,3),n1=n+1 D=A,D(2,1)=l1,D(2,2)=m1,D(2,3)=n1 D = 1 2 3 5 6 7 7 8 9 心得：熟练掌握matlab中对数组的操作方法尤为重要，其中，除了要牢记一些十分方便的函数外，应注意对数组的每一个元素进行操作时，应在运算符前加“.”号，如.^ .* 等。 6、设y=(x2 + ex cosx+[x])/x,分别计算x=1, 3, 5, 7.4时y的值。其中[x]表示x的取整函数。 过程：x=1,y=(x^2+exp(x)*cos(x)+fix(x))/x 结果：x = 1 y = 3.4687 x = 3 y = -2.6282 x = 5 y = 14.4198 x = 7.4 y = 105.2995 7、已知某地区1-12月份的平均气温为5，8，13，20，25，28，31，33，27，24，19，10，请绘图表示。 过程:： T=[5,8,13,20,25,28,31,33,27,24,19,10];bar(T) 结果： 第2次实验任务 某大学本科生就业情况如下：458人考入研究生，60人签到外企或合资企业，184人签到国营大企业，87人签到私人企业，13人自主创业，画出饼图表示之； 过程：x=[458 60 184 87 13]; pie3(x) 答案： 求极限 过程：syms x a; limit(((x-a)/(x+a))^x,x,inf) 答案：ans=exp(-2*a) 求极限 过程： syms x; limit(tan(x)^(1/log10(x)),x,0,right) 答案：ans=exp（1） 4.设 过程：syms x q t; q = limit ((1+1/x)^(2*x),x,inf); df = diff (q^t,t) 答案：df = q^t*log(q) 5.展开多项式； 过程：syms a b; F=(a+b)^3+(b-1)^2+a+2*b; expand(F) 答案： ans