Matlab上机操作2——符号计算.docVIP

1、%创建数值常量和符号常量 a1=2*sqrt(5)+pi %创建数值常量 a2=sym(2*sqrt(5)+pi) %按符号表达式创建符号常量 a3=sym(2*sqrt(5)+pi) %按最接近的两个正整数表示符号常量 a4=sym(2*sqrt(5)+pi,d) %按最接近的十进制浮点数表示符号常量 a31=a3-a1 %数值常量和符号常量的计算 a5=2*sqrt(5)+pi %字符串常量 2、%创建符号变量和符号表达式 f1=sym(a*x^2+b*x+c) %创建符号表达式 syms a b c x %创建符号变量 f2=a*x^2+b*x+c %创建符号表达式 syms(a,b,c,x) %注意这里是syms，不是sym sym(a,b,c,x) 错误使用 sym 输入参数太多。 f3=a*x^2+b*x+c; %创建符号表达式 syms x y real %创建实数符号变量 z=x+i*y; %创建z为复数符号变量 real(z) %复数z的实部是实数x sym(x,unreal); %清除符号变量的实数特性 real(z) %复数z的实部 3、%比较符号矩阵与字符串矩阵 A=sym([a,b;c,d]) %创建符号矩阵 B=[a,b;c,d] %创建字符串矩阵 C=[a,b;c,d] %创建数值矩阵 C=sym(B) %转换为符号矩阵 whos 4、%计算符号矩阵的行列式值、非共轭转置和特征值 syms a11 a12 a21 a22 A=[a11 a12;a21 a22] %创建符号矩阵 det(A) %计算行列式 A. %计算非共轭转置 eig(A)%计算特征值 5、%符号表达式的代数运算 f=sym(2*x^2+3*x+4) g=sym(5*x+6) f+g %符号表达式相加 f*g %符号表达式相乘 6、%对符号表达式进行任意精度控制并用三种运算方式表示同一符号常量 a=sym(2*sqrt(5)+pi) digits %显示默认的有效位数 vpa(a) %用默认的位数计算并显示 vpa(a,20) %按指定的精度计算并显示 digits(15) %改变默认的有效位数 vpa(a) %按digits指定的精度计算并显示 a1 =2/3 %数值型 a2 = sym(2/3) %有理数型 digits a3 =vpa(2/3,32) %VPA型 format long a1 7、%符号变量与数值变量进行转换 a1=sym(2*sqrt(5)+pi) b1=double(a1) %转换为数值变量 a2=vpa(sym(2*sqrt(5)+pi),32) b2=numeric(a2) %转换为数值变量 b2=numeric(a2) 未定义函数或变量 numeric。 b3=eval(a1) %将括号内的字符串视为语句并运行 whos 8、%符号表达式中的自由符号变量和符号变量 f=sym(a*x^2+b*x+c) findsym(f) %得出所有的符号变量 g=sym(sin(z)+cos(v)) findsym(g,1) %得出第一个符号变量 9、%多项式符号表达式的化简 f=sym(x^3-6*x^2+11*x-6) %多项式形式 g= sym((x-1)*(x-2)*(x-3)) %因式形式 h= sym( x*(x*(x-6)+11)-6) %嵌套形式 pretty(f) collect(g) f1=sym(x^3+2*x^2*y+4*x*y+6) collect(f1,y) %按y来合并同类项 expand(g) horner(f) 秦久韶 factor(f) y=sym(cos(x)^2-sin(x)^2) simplify(y) simple(y) f=sym(x^3-6*x^2+11*x-6) f = x^3 - 6*x^2 + 11*x - 6 g=sym((x-1)*(x-2)*(x-3)) g = (x - 1)*(x - 2)*(x - 3) h=sym(x*(x*(x-6)+11)-6) h = x*(x*(x - 6) + 11) - 6 pretty(f) 3 2 x - 6 x + 11 x - 6 collect(g) ans = x^3 - 6*x^2 + 11*x - 6 f1=sym(x^3+2*x^2*y+4*x*y+6) f1 = x^3 + 2*y*x^2 + 4*y*x + 6 collect(f1,y) ans = (2