当前位置:文档之家› 浅析Matlab数学实验报告

浅析Matlab数学实验报告

数学实验报告姓名:班级:学号:第一次实验任务过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000ia-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 2849 6476 100(2) C=B.^2C =1 49 1625 36(3) sin(B)ans =0.8415 0.90930.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+1D=A,D(2,1)=l1,D(2,2)=m1,D(2,3)=n1D =1 2 35 6 77 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 = 1y = 3.4687x = 3y = -2.6282x = 5y = 14.4198x = 7.4y = 105.29957、已知某地区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 ax 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 xt 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-17.求方程0123=+-x x 的根过程:syms x;C=solve(…x^3-2*x+1=0?)答案C=11/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)))9. 求不定积分过程: syms x;F=1/(sqrt(2*x+3)+sqrt(2*x-1));Nf=int(F,x,)答案:Nf = -1/12*(2*x-1)^(3/2)+1/12*(2*x+3)^(3/2)过程: syms x;F=1/(1+cos(x));Nf=int(F)答案:Nf = tan(1/2*x).的导数求函数x x x y ++=dx x x ?-++12321.cos 11?+dx x10 求定积分过程:syms x;F=1/(x*sqrt(log(x)*(1-log(x)));Nf=int(F,x,exp(1/2),exp(3/4))答案:ans=asin(2*log()-96*log(2)-1)-asin(2*log(7)-100*lo g(2)-1)11. 解方程组 ???=-=+1382y x y x 过程:syms x y;eq1=sym(…2*x+y=8?);eq2=sym(…x-3*y=1?);[x,y]=solve(eq1,eq2)答案:x =25/7y = 6/712. 求和∑=20121k k 过程:syms k;symsum(1/k^2,k,1,20)答案:ans = 41/2013. .156x 2cos 次泰勒多项式处的在求π=x 过程:syms x; taylor(cos(2*x),x,15,pi/6)答案:ans=1/2-3^(1/2)*(x-1/6*pi)-(x-1/6*pi)^2+2/3*3^(1/2)*(x-1/6*pi)^3+1/3*(x-1/6*pi)^4-2/15*3^(1/2)*(x-1/6*pi)^5-2/45*(x-1/6*pi)^6+4/315*3^(1/2)*(x-1/6*pi)^7+1/315*(x-1/6*pi)^8-2/2835*3^(1/2)*(x-1/6*pi)^9-2/14175*(x-1/6*pi)^10+4/155925*3^(1/2)*(x-1/6*pi)^11+2/467775*(x-1/6*pi)^12-4/6081075*3^(1/2)*(x-1/6*pi)^13-4/42567525*(x-1/6*pi)^1414. 编写猜数游戏程序:首先由计算机随机产生一个 [1,100] 之间的一个整数,然后由用户猜测所产生的这个数。

根据用户猜测的情况给出不同的提示,如果猜测的数大于产生的数,则显示“High” ,小于则显示“ Low ” ,等于则显示“You won !”,同时退出游戏。

用户最多有 7 次机会。

过程:.)ln 1(ln 43?-e e x x x dxa=fix(rand(1)*100);for i=1:1:7x=input('please input the number you guess');if x>afprintf('high');continue;end。

相关主题