05实验五 符号计算实验.docVIP

实验五 符号计算实验 一、实验目的： 1、 掌握定义符号对象的方法 syms x y x=6; y=5; z=(x+1)./(sqrt(x+3)-sqrt(y)) z = 9.16311896062463 2.分解因式。 （1）； （2） （1） clear syms x y factor(x.^4-y.^4) ans = (x-y)*(x+y)*(x^2+y^2) （2） clear syms x factor(125.*x.^6+75.*x.^4+15.*x.^2+1) ans = (5*x^2+1)^3 3.化简表达式 （1）；（2） （1） clear syms bata1 bata2 simplify(sin(bata1).*cos(bata2)-cos(bata1).*sin(bata2)) ans = sin(bata1-bata2) （2） clear syms x simplify((4.*x.^2+8.*x+3)./(2.*x+1)) ans = 2*x+3 4.已知 完成下列运算： （1）. （2）B的逆矩阵并验证结果。 （3）包括B矩阵主对角线元素的下三角阵。 （1） p1=sym([0,1,0;1,0,0;0,0,1]); p2=sym([1,0,0;0,1,0;1,0,1]); A=sym([a,b,c;d,e,f;g,h,i]); B=p1*p2*A B = [ d, e, f] [ a, b, c] [ a+g, b+h, c+i] （2） simplify(inv(B)) ans = [ (i*b-c*h)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), (-e*c-i*e+f*b+f*h)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), -(-e*c+f*b)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b)] [ -(i*a-c*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), -(-d*c-i*d+f*a+f*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), (-d*c+f*a)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b)] [ (a*h-b*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), (-d*b-d*h+e*a+e*g)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b), -(-d*b+e*a)/(i*d*b-d*c*h-i*a*e+a*f*h+g*e*c-g*f*b)] simplify(B*inv(B)) ans = [ 1, 0, 0] [ 0, 1, 0] [ 0, 0, 1] （3） B([4,7,8])=[0,0,0] B = [ d, 0, 0] [ a, b, 0] [ a+g, b+h, c+i] 5.用符号方法求下列极限或导数。 。 。 ，求和。 （1） clear syms x; fun=(x.*(exp(sin(x))+1)-2.*(exp(tan(x))-1))./(sin(x)).^3; limit(fun,x,0) ans = -1/2 （2） clear syms x; fun=(sqrt(x)-sqrt(acos(x)))./sqrt(x+1); limit(fun,x,1) ans = 1/2*2^(1/2) （3） clear syms x; y=(1-cos(2.*x))./x; dx1=diff(y) dx1 = 2*sin(2*x)/x-(1-cos(2*x))/x^2 diff(y,2) ans = 4*cos(2*x)/x-4*sin(2*x)/x^2+2*(1-cos(2*x))/x^3 6.用符号方法求下列积分。 （1） （2） （3） （1） clear syms x fun1=1./(1+x.^4+x.^8); int(fun1) ans = 1/6*3^(1/2)*atan(1/3*(2*x-1)*3^(1/2))+1/6*3^(1/2)*atan(1/3*(2*x+1)*3^(1/2))-1/12*3^(1/2)*log(x^2-3^(1/2)*x+1)+1/12*3^(1/2)*log(x^2+3^(1/2)*x+1) (2)