Matlab数学实验报告
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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 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++-++=;
过程: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