MATLAB符号运算实验.docVIP

MATLAB符号运算实验 用符号函数法求解方程at2+b*t+c=0。 clc;clear; solve(a*t^2+b*t+c) ans = -(b + (b^2 - 4*a*c)^(1/2))/(2*a) -(b - (b^2 - 4*a*c)^(1/2))/(2*a) 用符号计算验证三角等式：sin((1)cos((2)-cos((1)sin((2) =sin((1-(2)、 clc;clear; syms a b expand(sin(a)*cos(b)-cos(a)*sin(b)-sin(a-b)) ans = 0 因为sin(a)*cos(b)-cos(a)*sin(b)-sin(a-b)=0， 所以sin(a)*cos(b)-cos(a)*sin(b)=sin(a-b)。 用符号函数法求矩阵的行列式值、逆和特征根。 clc;clear; syms a11 a12 a21 a22 M=[a11,a12;a21,a22]; [m,n]=size(M) inv(M) [V,x]=eig(M) m = 2 n = 2 ans = [ a22/(a11*a22 - a12*a21), -a12/(a11*a22 - a12*a21)] [ -a21/(a11*a22 - a12*a21), a11/(a11*a22 - a12*a21)] V = [ -(a22 - a11 + (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2))/(2*a21), (a11 - a22 + (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2))/(2*a21)] [1, 1] x = [ a11/2 + a22/2 - (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2)/2, 0] [0, a11/2 + a22/2 + (a11^2 - 2*a11*a22 + a22^2 + 4*a12*a21)^(1/2)/2] 用符号函数法求因式分解： clc;clear; syms x factor(x^4-5*x^3+5*x^2+5*x-6) ans = (x - 1)*(x - 2)*(x - 3)*(x + 1) 5、 ，用符号微分求df/dx。 clc;clear; syms x a f=[a,x^2,1/x;exp(a*x),log10(x),sin(x)]; diff(f,x) ans = [ 0, 2*x, -1/x^2] [ a*exp(a*x), 1/(x*log(10)), cos(x)] 用符号函数法求代数方程组关于x,y的解。 clc;clear; syms x y f1=a*x^2+b*y+c=0; f2=x+y=0; solve(f1,f2) ans = x: (b + (b^2 - 4*a*c)^(1/2))/(2*a) y: -(b + (b^2 - 4*a*c)^(1/2))/(2*a) 7、 符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t的变化范围为[0,2(]。 clc;clear; syms x y t x=sin(3*t)*cos(t); y=sin(3*t)*sin(t); ezplot(x,y,[0,2*pi]) 8、 符号函数绘图法绘制极坐标下sin(3*t)*cos(t)的图形。 clc;clear; syms t ezpolar(sin(3*t)*cos(t))