Matlab数学实验报告
w .. . .. 数学实验报告
：
班级：
学号：
第一次实验任务
过程： 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)
./,,,,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====
>> 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 + 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
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次实验任务
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 1
)(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++-++=； 过程：syms a b;
F=(a+b)^3+(b-1)^2+a+2*b; expand(F)
答案： ans = a^3+3*a^2*b+3*a*b^2+b^3+b^2+1+a
6.分解因式123245-++-=x x x x y
过程：syms x;
F=x^5+3*x^4+2*x^2+x-1; factor(F)
答案：ans= x^5+3*x^4+2*x^2+x-1
7.求方程0123=+-x x 的根 过程：syms x;
C=solve(‘x^3-2*x+1=0’) 答案C=
1 1/2*5^(1/2)-1/
2 -1/2-1/2*5^(1/2) 8.
过程：>> syms x;
>> f=sqrt(x+sqrt(x+sqrt(x))); >> df=diff(f) 答
案
：
df=
1/2/(x+(x+x^(1/2))^(1/2))^(1/2)*(1+1/2/(x+x^(1/2))^(1/2)*(1+1/2/x^(1/2)))
.
的导数求函数x x x y +
+
=