《MATLAB数学实验100例题解_1.docx》由会员分享,可在线阅读,更多相关《MATLAB数学实验100例题解_1.docx(25页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、MATLAB数学实验100例题解一元函数微分学实验1一元函数的图形基础实验实验目的通过图形加深对函数及其性质的认识与理解,把握运用函数的图形来观察和分析函数的有关特性与变化趋势的方法,建立数形结合的思想;把握用Matlab作平面曲线图性的方法与技巧.初等函数的图形2作出函数xytan=和xycot=的图形观察其周期性和变化趋势.解:程序代码:x=linspace(0,2*pi,600);t=sin(x)./(cos(x)+eps);plot(x,t);title(tan(x);axis(0,2*pi,-50,50);图象:程序代码:x=linspace(0,2*pi,100);ct=cos(x
2、)./(sin(x)+eps);plot(x,ct);title(cot(x);axis(0,2*pi,-50,50);图象:cot(x)4在区间1,1-画出函数xy1sin=的图形.解:程序代码:x=linspace(-1,1,10000);y=sin(1./x);plot(x,y);axis(-1,1,-2,2)图象:二维参数方程作图6画出参数方程?=tttytttx3cossin)(5coscos)(的图形:解:程序代码:t=linspace(0,2*pi,100);plot(cos(t).*cos(5*t),sin(t).*cos(3*t);图象:极坐标方程作图8作出极坐标方程为10/
3、ter=的对数螺线的图形.解:程序代码:t=0:0.01:2*pi;r=exp(t/10);polar(log(t+eps),log(r+eps);图象:90270分段函数作图10作出符号函数xysgn=的图形.解:程序代码:x=linspace(-100,100,10000);y=sign(x);plot(x,y);axis(-100100-22);函数性质的研究12研究函数)3(log3)(35xexxfx-+=在区间2,2-上图形的特征.解:程序代码:x=linspace(-2,2,10000);y=x.5+3*exp(x)+log(3-x)/log(3);plot(x,y);图象:实验
4、2极限与连续基础实验实验目的通过计算与作图,从直观上揭示极限的本质,加深对极限概念的理解.把握用Matlab画散点图,以及计算极限的方法.深化理解函数连续的概念,熟悉几种间断点的图形特征,理解闭区间上连续函数的几个重要性质.作散点图14分别画出坐标为)10,2,1(),4,(),(3222=+iiiiii的散点图,并画出折线图.解:散点图程序代码:i=1:10;plot(i,i.2,.)或:x=1:10;y=x.2;i=1:10;plot(i.2,4*(i.2)+i.3,-x)数列极限的概念16通过动画观察当n时数列21nan=的变化趋势.解:程序代码:n=1:100;an=(n.2);n=1
5、:100;an=1./(n.2);n=1:100;an=1./(n.2);fori=1:100plot(n(1:i),an(1:i),axis(0,100,0,1)pause(0.1)end图象:函数的极限18在区间4,4-上作出函数xxxxxf-=339)(的图形,并研究)(limxfx和).(lim1xfx解:作出函数xxxxxf-=339)(在区间4,4-上的图形x=-4:0.01:4;y=(x.3-9*x)./(x.3-x+eps);plot(x,y)从图上看,()fx在x1与x时极限为0两个重要极限20计算极限?+xxxxxsin11sinlim)1(0xxex2lim)2(+30s
6、intanlim)3(xxxx-xxx0lim)4(+xxxlncotlnlim)5(0+xxxlnlim)6(20+xxxxxxsincossinlim)7(20-125523lim)8(323+-xxxxxxxxeexxxsin2lim)9(0-xxxxcos110sinlim)10(-?解:1limit(x*sin(1/x)+1/x*sin(x)ans=1(2)limit(x2/exp(x),inf)ans=0(3)limit(tan(x)-sin(8)/x3)ans=NaN(4)limit(xx,x,0,right)ans=1(5)limit(log(cot(x)/log(x),x,0
7、,right)ans=-1(6)limit(x2*log(x),x,0,right)ans=016 (7)limit(sin(x)-x.*cos(x)./(x.2.*sin(x),x,0)ans=1/3(8)limit(3*x.3-2*x.2+5)/(5*x.3+2*+1),x,inf)ans=3/5(9)limit(exp(x)-exp(-x)-2*x)./(x-sin(x)ans=2(10)limit(sin(x)/x).(1/(1-cos(x)ans=exp(-1/3)实验3导数基础实验实验目的深化理解导数与微分的概念,导数的几何意义.把握用Matlab求导数与高阶导数的方法.深化理解和
8、把握求隐函数的导数,以及求由参数方程定义的函数的导数的方法.导数概念与导数的几何意义22作函数71232)(23+-+=xxxxf的图形和在1-=x处的切线.解:作函数71232)(23+-+=xxxxf的图形程序代码:symsx;y=2*x3+3*x2-12*x+7;diff(y)ans=6*x2+6*x-12symsx;y=2*x3+3*x2-12*x+7;f=diff(y)f=6*x2+6*x-12x=-1;f1=6*x2+6*x-12f1=-12f2=2*x3+3*x2-12*x+7f2=20x=linspace(-10,10,1000);y1=2*x.3+3*x.2-12*x+7;y
9、2=-12*(x+1)+20;plot(x,y1,r,x,y2,g)求函数的导数与微分24求函数bxaxxfcossin)(=的一阶导数.并求.1?+baf解:求函数bxaxxfcossin)(=的一阶导数程序代码:symsabxy;y=sin(a*x)*cos(b*x);D1=diff(y,x,1)答案:D1=cos(a*x)*a*cos(b*x)-sin(a*x)*sin(b*x)*b求.1?+baf程序代码:x=1/(a+b);cos(a*x)*a*cos(b*x)-sin(a*x)*sin(b*x)*b答案:ans=cos(a/(a+b)*a*cos(b/(a+b)-sin(a/(a+
10、b)*sin(b/(a+b)*b拉格朗日中值定理26对函数),2)(1()(-=xxxxf观察罗尔定理的几何意义.(1)画出)(xfy=与)(xf的图形,并求出1x与.2x解:程序代码:symsx;f=x*(x-1)*(x-2);f1=diff(f)f1=(x-1)*(x-2)+x*(x-2)+x*(x-1)solve(f1)ans=1+1/3*3(1/2)1-1/3*3(1/2)x=linspace(-10,10,1000);y1=x.*(x-1).*(x-2);y2=(x-1).*(x-2)+x.*(x-2)+x.*(x-1);plot(x,y1,x,y2) (2)画出)(xfy=及其在点
11、)(,(11xfx与)(,(22xfx处的切线.程序代码:symsx;f=x*(x-1)*(x-2);f1=diff(f)f1=(x-1)*(x-2)+x*(x-2)+x*(x-1)solve(f1)ans=1+1/3*3(1/2)1-1/3*3(1/2)x=linspace(-3,3,1000);y1=x.*(x-1).*(x-2);y2=(x-1).*(x-2)+x.*(x-2)+x.*(x-1);plot(x,y1,x,y2)holdonx=1+1/3*3(1/2);yx1=x*(x-1)*(x-2)yx1=-0.3849x=1-1/3*3(1/2);yx2=x*(x-1)*(x-2)y
12、x2=0.3849x=linspace(-3,3,1000);yx1=-0.3849*x.0;yx2=0.3849*x.0;plot(x,yx1,x,yx2)28求下列函数的导数:(1)31+=xey;解:程序代码:symsxy;y=exp(x+1)3);D1=diff(y,1)答案:D1=3*(x+1)2*exp(x+1)3) (2)42lntan(+=xy;解:程序代码:symsx;y=log(tan(x/2+pi/4);D1=diff(y,1)答案:D1=(1/2+1/2*tan(1/2*x+1/4*pi)2)/tan(1/2*x+1/4*pi)(3)xxysinlncot212+=;解
13、:程序代码:symsx;y=1/2*(cot(x)2+log(sin(x);D1=diff(y,1)答案:D1=cot(x)*(-1-cot(x)2)+cos(x)/sin(x)(4)xy2arctan21=.解:程序代码:symsx;y=sqrt(2)*atan(sqrt(2)/x);D1=diff(y,1)答案:D1=-2/x2/(1+2/x2)一元函数积分学与空间图形的画法实验4一元函数积分学(基础实验)实验目的把握用Matlab计算不定积分与定积分的方法.通过作图和观察,深化理解定积分的概念和思想方法.初步了解定积分的近似计算方法.理解变上限积分的概念.提高应用定积分解决各种问题的能力.不定积分计算30求.)1(532?-dxxx解:程序代码:symsxy;y=x2*(1-x3)5;R=int(y,x)答案:R=-1/18*x18+1/3*x15-5/6*x12+10/9*x9-5/6*x6+1/3*x332求.arctan2?xdxx解:程序代码:symsxy;y=x2*atan(x);R=int(y,x)答案:R=1/3*x3*atan(x)-1/6*x2+1/6*log(x2+1)