数学建模(华中农业大学课件).pptVIP

* * 数 学 实 验 华中农业大学数学建模基地系列课件 微积分实验 实验目的： 掌握MATLAB的极限、导数、积分、微分方程运算； limit(f,x,a)：计算当变量x趋近于常数a时， f(x)函数的极限值； (2) limit(f,x,a,‘right’)：‘right’表示变量x从右边 趋近于a； (3) limit(f,x,a,‘left’)： ‘left’表示变量x从左边 趋近于a + ，则可以用+inf 2.1 极限 例1 求极限 syms x; %定义变量x limit(1/x^2-cot(x)^2,x,0) % ans=2/3 例2 求极限 clear syms a b x; limit((sin(a/x^2)+cos(b/x))^(x^2),x,inf) ans=exp(a-1/2*b^2) 例3 求极限 clear syms x; limit(x^x,x,0,right) ans=1 练习题 （1）求极限 （2）求极限 roots(p)：这是求多项式根的命令，其中p是多项式的系数，按降幂方式排列; (2)fzero(f,x0): 在x0的附近寻找函数f的近似零点 (3) solve(f,x)：求解表达式f的代数方程， 求解变量为x； (4) solve(f1,f2,…,fn,x1,x2,…,xn)：求解符号表 达式f1,f2,…,fn组成的代数方程组，求解变量分别 是x1,x2,…xn。 2.2 代数方程(组)的解 (P34) 例1 求方程ex ? x ? 3 = 0在区间(1, 2)内的一个实根 . syms x; f=exp(x)-x-3; %定义函数f fzero(f,[1,2]) syms x; f=exp(x)-x-3; %定义函数f fzero(f,[1,2]) ans= 1.5.52 例2 求3次方程x3 ? x ? 3 = 0的所有根. roots([1,0,-1,-3]) roots([1,0,-1,-3]) ans = 1.6717 -0.8358 + 1.0469i -0.8358 - 1.0469i clear; dsj=solve(3*x-y=1,x+y=3,x,y) x=simplify(dsj.x) %简化函数 x = 1 y=simplify(dsj.y) %简化函数 y = 2 例3 求解代数方程组3x ? y = 1, x + y = 3. clear; [a,b]=solve(3*x-y=1,x+y=3,x,y) diff(f,x)：以x为自变量，对符号表达式f求一阶导数； 2.3 导数 help diff (P35) diff(f,x,n)：以x为自变量，对符号表达式s求n阶导数。 例1 设y = xe3x, 求y , y (5). syms x; y=x*exp(3*x); y1=diff(y,x); % 1阶导数 y5=diff(y,x,5); % 5阶导数 y1,y5 y1=exp(3*x)+3*x*exp(3*x) y5=405*exp(3*x)+243*x*exp(3*x) 例2 设 求y , y. sym x y=log(x+sqrt(1+x^2)); y1=diff(y,x); % 1阶导数 y2=diff(y,x,2); % 2阶导数 simple(y1),simple(y2)%结果化简 ans =1/(1+x^2)^(1/2) ans =-x/(1+x^2)^(3/2) 例3 设 求yx sym x y z=log(y+sqrt(1+x^2)); zx=diff(z,x); % 1阶导数 simple(zx) %结果化简 (1) fminbnd(f,x1,x2)： 求函数f在区间[x1,x2]上的极小值； (2) fminsearch(‘f’,x0) 求多元函数f在x0附近的极小值 2.4 函数极值 (P38) 例1 观察函数 f (x) = x cos x在[?8, 8]上的图形，研究其在x = 2附近的一个极小值点. syms x; f=x*cos(x); %定义函数f fplot(f,[-8,8]) grid; %显示网格 *