MATLAB数学实验100例题解_3.docx
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)./(sin(x)+eps);plot(x,ct);title('cot(x)');axis(0,2*pi,-50,50);图象：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/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);图象：实验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: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(+30sintanlim)3(xxxx-xxx0lim)4(+xxxlncotlnlim)5(0+xxxlnlim)6(20+xxxxxxsincossinlim)7(20-125523lim)8(323+-xxxxxxxxeexxxsin2lim)9(0-xxxxcos110sinlim)10(-?解：1>>limit(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,'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求导数与高阶导数的方法.深化理解和把握求隐函数的导数,以及求由参数方程定义的函数的导数的方法.导数概念与导数的几何意义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-12>>symsx;y=2*x3+3*x2-12*x+7;>>f=diff(y)f=6*x2+6*x-12>>x=-1;f1=6*x2+6*x-12f1=-12>>f2=2*x3+3*x2-12*x+7f2=20>>x=linspace(-10,10,1000);y1=2*x.3+3*x.2-12*x+7;y2=-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+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及其在点)(,(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)>>holdon>>x=1+1/3*3(1/2);>>yx1=x*(x-1)*(x-2)yx1=-0.3849>>x=1-1/3*3(1/2);>>yx2=x*(x-1)*(x-2)yx2=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+=;解：程序代码：>>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)
