【通用】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)

./,,,,2,311b a b a b a b a i b i a ⨯-+-=+=计算、设有两个复数

6,7,5.4)

cos()sin(2=-=++y x y x y x ，其中、计算的图形。

下分别绘制）同一页面四个坐标系）同一坐标系下（、在（

x y e y x y x y x ln ,,cos ,sin 213====

>> 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

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 + e x 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

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次实验任务

1. 某大学本科生就业情况如下：458人考入研究生，60人签到外企或合资企业，184人签到国营大企业，87人签到私人企业，13人自主创业，画出饼图表示之；

过程：x=[458 60 184 87 13];

pie3(x)

2. 求极限x

x a

x a x )(

lim +-∞→

过程：syms x a;

limit('((x-a)/(x+a))^x',x,inf) 答案：ans=exp(-2*a)

3. 求极限x

x x ln 10

)(tan lim +

→

过程：>> syms x;

>> limit('tan(x)^(1/log10(x))',x,0,'right') 答案：ans=exp （1）

4.设)(',)11(lim )(2t f x

t f x

t x 求+

=∞→

过程：syms x q t;

q = limit ('(1+1/x)^(2*x)',x,inf); df = diff ('q^t',t) 答案：df = q^t*log(q)

5.展开多项式b a b a y 2)1()1(23++-++=；