MATLAB中数值计算及符号计算.pptVIP

微积分实验(2) MATLAB中的数值计算与符号计算 ;◆函数复合 命令格式: compose(f,g,x), 即求f(g(x); 例子: syms f, x; f=sin(x);compose(f,f,x) ans = sin(sin(x)) syms f, x; f=sin(x); g=1+x^2; compose(f,g,x) ans = sin(1+x^2) ;◆设f(sin(x/2)) = cos(x) + 1, 求f(x)与f(cos(x/2)). syms x, f; old=cos(x); new=1-2*sin(x/2)^2; subs(f(sin(x/2))=cos(x)+1,old,new) ans = f(sin(x/2))=(1-2*sin(x/2)^2)+1 syms u; subs(ans,sin(x/2),u) ans = f(u) = 2-2*u^2 ;Syms x,f,g; f=x/sqrt(1+x^2);g=x; g=simple(compose(f,f,x)) f =x/(1+x^2)^(1/2) g =x/(2*x^2+1)^(1/2) g=simple(compose(g,f,x)) g =x/(3*x^2+1)^(1/2) g=simple(compose(g,f,x)) g =x/(4*x^2+1)^(1/2) g=simple(compose(g,f,x)) g =x/(5*x^2+1)^(1/2) ;syms x,f; f=sin(x)^2/(1+cot(x))+cos(x)^2/(1+tan(x)); diff(f,x) Ans= 2*sin(x)/(1+cot(x))*cos(x)-sin(x)^2/(1+cot(x))^2*(-1-cot(x)^2)-2*cos(x)/(1+tan(x))*sin(x)-cos(x)^2/(1+tan(x))^2*(1+tan(x)^2) simple(ans) ans =-2*cos(x)^2+1 diff(f,x,2) ans=2*cos(x)^2/(1+cot(x))-4*sin(x)/(1+cot(x))^2*cos(x)*(-1-cot(x)^2)-2*sin(x)^2/(1+cot(x))+2*sin(x)^2/(1+cot(x))^3*(-1-cot(x)^2)^2+2*sin(x)^2/(1+cot(x))^2*cot(x)*(-1-cot(x)^2)+2*sin(x)^2/(1+tan(x))+4*cos(x)/(1+tan(x))^2*sin(x)*(1+tan(x)^2)-2*cos(x)^2/(1+tan(x))+2*cos(x)^2/(1+tan(x))^3*(1+tan(x)^2)^2-2*cos(x)^2/(1+tan(x))^2*tan(x)*(1+tan(x)^2) simple(ans) ans =4*sin(x)*cos(x) ;syms x y z f g; f=exp(z)-z+x*y-3; g=jacobian(f,[x, y, z]) g= [ y, x, exp(z)-1] g(1) ans = y g(2) ans = x g(3) ans = exp(z)-1 g(1)/g(3) ans = y/(exp(z)-1) ;syms x y t a; x=a*cos(t)^3; y=a*sin(t)^3; y=simple(diff(y,t)/diff(x,t)) y = -tan(t) simple(diff(y,t)/diff(x,t)) ans = 1/3/cos(t)^4/a/sin(t);syms x; int((1-sin(x)+cos(x))/(1+sin(x)-cos(x))) ans = 2*log(tan(1/2*x))-2*log(tan(1/2*x)+1)-2*atan(tan(1/2*x));先画出图形;再求出驻点(有2个,分别为-1及3) solve(diff(y,x),x) ans = [ -1] [ 3] 求函数的2阶导数 y=diff(y,x,2) y = 12*x-12 ;计算2阶导数在驻点处的函数值 x=-1; eval(y) ans = -24 x=3; eval(y) ans = 24 由高等数学定理知,在-1处取得极大值, 在3处取得极小值 ;这个积分无法计算出理论解,只能计算数值解 syms x; int(sin(sin(x)),1,10) ans = int(sin(sin(x)),x = 1 .. 10) ;第二种方法,用Simpson(辛普森法)计算,其命令是quad(‘f’,a,b,tol),