matlab简明教程 第二章【DOC精选】.doc

第二章 MATLAB在微积分问题求解中的应用 2.1 微分问题的MATLAB求解 1. 函数作图 MATLAB函数画图可通过ezplot或fplot等函数实现。 1）ezplot ezplot函数的调用格式如下 ezplot(f,[a,b]) 功能：表示在区间[a,b]绘制y=f(x)的函数图，当区间缺省时默认区间[-2*pi,2*pi]。 ezplot(x,y,[tmin,tmax]) 功能：在区间tmin t tmax上绘制参数方程x = x(t)，y = y(t)的图形当区间缺省时默认区间[-2*pi,2*pi]。 例1 ezplot(sin(x)) 图2.1.1 例2 ezplot(t*cos(t),t*sin(t),[0,4*pi]) 图2.1.2 2) fplot fplot 函数的调用格式如下 fplot(fun,lims) 功能：绘制函数fun在区间lims上的图形。 例３fplot(tan(x),[-pi/4 pi/4]) 图2.1.3 2 极限的符号运算 极限是高等数学中基本概念之一，在微积分中，很多概念是用极限定义的，例如导数和定积分。因此，掌握极限的运算对学好高等数学是极为重要的。在MATLAB中，极限的求解可由limit函数来实现，limit函数的格式及功能见表2.2.1。 表2.2.1 1limit函数的格式及功能 MATLAB 函数名 函数功能 双 侧 极 限 limit(F,x,a) 计算 limit(F,x,inf) 计算 limit(F) 计算 limit(F,a) 计算 （x为默认自变量） 单 侧 极 限 limit(F,x,a, right) 计算 limit(F,x,a, left) 计算 limit(F(x),x,inf, left) 计算 limit(F(-x),x,inf, left) 计算 因为数列实际上就是定义在正整数集合上的函数，因此数列的极限可看成时的特殊函数的极限；多元函数的极限可化为累次极限实现。 例1 求下列数列的极限 1） 2） 3） 4） 5） 6） 7） 8） 9） 解：syms n a r1=limit(sqrt(n^2+a^2)/n,n,inf,left) 输出 r1 =1 r2=limit(sqrt(n^2+3)-sqrt(n^2-3),n,inf,left) 输出r2 =0 r3=limit(3^n*sin(pi/3^n),n,inf,left) 输出r3 =pi r4=limit((2^n-3^(n+1))/(3^n-2^(n+1)),n,inf,left) 输出r4 =-3 r5=limit(sin(pi*sqrt(n^2+1)),n,inf,left) 输出r5 =1 .. 1 r6=limit(((n-1)/(n+1))^n,n,inf,left) 输出r6 =exp(-2) r7=limit((n-1)^2/(n+1),n,inf,left) 输出r7 =Inf r8=limit((-1)^n,n,inf,left) 输出r8 =-1 .. 1 r9=limit((-2)^n,n,inf,left) 输出r9 =NaN 例2 求下列函数的极限 1） 2） 3） 4） 5） 6） 7） 8） 9） 解： syms x h t f1=limit((sin(x+h)-sin(x))/h,h,0) 输出f1 =cos(x) f2=limit(1/(1-x)-3/(1-x^3),x,1) 输出 f2 =-1 f3=limit(x*sin(1/x)) 输出 f3 =0 f4=limit(t/(x-2),3) 输出f4 =t f5=limit(abs(x)/x,x,0,left) 输出f5 =-1 f6=limit((1+t/(-3*x))^(-x),x,inf,left) 输出f6 =exp(1/3*t) f7=limit(((2*x+3)/(2*x+1))^(x+1),x,inf) 输出f7 =exp(1) f8=limit(x*sin(1/(x-1)),x,1) 输出f8 =-1 .. 1 f9=limit(x*sin(x),x,inf) 输出f9 =NaN 例3 求下列函数的极限 1） 2） 解： syms x y; p1=limit(limit((2-sqrt(x*y+4))/(x*y),x,0),y,0) 输出p1 =-1/4 p2=limit(limit(log(x+exp(y))/sqrt(x^2+y^2),x,1),y,0) 输出p2 =log(2) 3. 一阶微商的计算 由导数的定义可知，一切导数的问题，都可以用极限的方法求得