MATLAB基础及其在数学建模中的应用讲述.ppt

2 积分 int(f) 函数f 对符号变量x 或接近字母x 的符号变量求不定积分 int(f,t) 函数f 对符号变量t 求不定积分 int(f,a,b) 函数f 对符号变量x 或接近字母x 的符号变量求从a 到b 的定积分（a,b控制上下线） int(f,t,a,b) 函数f 对符号变量t 求从a 到b 的定积分 例 syms a x f=sin(a*x) g=int(f) gg=int(f,a) f =sin(a*x) ff =sin(x^3) g =-1/a*cos(a*x) gg =-1/x*cos(a*x) 例 syms a x f=sin(a*x) g=int(f,0,pi) f = sin(a*x) g =-(cos(pi*a)-1)/a 3 极限 limit(f) 当符号变量x 或最接近字母x 的符号变量趋于0 时函数f 的极限 limit(f,t,a) 当符号变量t 趋于a 时函数f 的极限 左右极限的求法 例：sym x limit(1/x) limit(1/x,x,0,‘left’) left左极限 limit(1/x,x,0,‘right’) 右 ans =x ans =NaN ans =-inf ans =inf 例 syms x t a f=sin(x)/x g=limit(f) limit((cos(x+a)-cos(x))/a,a,0) limit((1+x/t)^t,t,inf) f =sin(x)/x g =1 ans =-sin(x) ans =exp(x) 4 求级数和 symsum(s,t,a,b) 表示s 中的符号变量t 从a 到b 的级数和t 缺省时设定为x 或最接近 x 的字母 例 syms x k symsum(1/x,1,3) ans =11/6 s1=symsum(1/x^2,1,inf) s2=symsum(x^k,k,0,inf) s1 =1/6*pi^2 s2 =-1/(x-1) 5 泰勒 Taylor 多项式（n阶泰勒展式） taylor(f,n,a) 函数f 对符号变量x 或最接近字母x 的符号变量在a点的n-1 阶泰勒多项式，n 缺省时值为6 ，a 缺省值为0 例 taylor(sin(x)) ans =x-1/6*x^3+1/120*x^5 f=log(x) s=taylor(f,4,2) f =log(x) s =log(2)+1/2*x-1-1/8*(x-2)^2+1/24*(x-2)^3 三 解方程 1 代数方程 solve(f,t) 对f 中的符号变量t 解方程f=0 ，t 缺省值为x 或最接近x 的字母 例 syms a b x c f=a*x^2+b*x+c s=solve(f) ss=solve(f,b) f =a*x^2+b*x+c s = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))] ss = -(a*x^2+c)/x 注 求解形如f(x)=q(x)形式的方程则需要用单引号把方程括起来 例 s=solve(cos(2*x)+sin(x)=1) s = [ 0] [ pi] [ 1/6*pi] [ 5/6*pi] 求解方程组 例 [x,y]=solve(x^2+x*y+y=3,x^2-4*x+3=0) x = [ 1] [ 3] y = [ 1] [ -3/2] 2 微分方程 dsolve(S,s1,s2,…,x) 其中S 为方程，s1,s1,s3,…为初始条件，x 为自变量，方程S 中用D 表示求导数，D2,D3,…表示二阶三阶等高阶导数，初始条件缺省时给出带任意常数C1,C2,..的通解，自变量缺省值为t ，也可求解微分方程组 例 dsolve(Dy=1+y^2) ans =tan(t+C1) y=dsolve(Dy=1+y^2,y(0)=1,x) y =tan(x+1/4*pi) x=dsolve(D2x+2*D1x+2*x=exp(t),x(0)=1,Dx(0)=0) x =1/5*exp(t)+3/5*exp(-t)*sin(t)+4/5*exp(-t)*cos(t) S=dsolve(Df=3