实验Matlab符号计算说明.docVIP

实验7 Matlab符号计算 实验目的: 掌握定义符号对象的方法; 掌握符号表达式的运算法则以及符号矩阵运算? 掌握求符号函数极限及导数的方法? 掌握求符号函数定积分和不定积分的方法? 实验内容: 已知x=6,y=5,利用符号表达式求 分解因式? (1); (2) 化简表达式 (1); (2) 已知 完成下列运算: (1). (2)B的逆矩阵并验证结果? (3)包括B矩阵主对角线元素的下三角阵? 用符号方法求下列极限或导数? ? ? ,求和? 用符号方法求下列积分? (1) (2) (3) 详细实验内容: 1.已知x=6,y=5,利用符号表达式求 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.用符号方法求下