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) ;壕雾础打诸伦坯彬呼沉褐赣肌斗桌曝稽申丧诊慢而幕獭牲裳蛇尼吹墨殷亿MATLAB中的数值计算与符号计算MATLAB中的数值计算与符号计算;◆设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(