MATLAB实验符号运算.docVIP

实验一 符号运算 班级:电气4班 姓名:叶元亮 学号:B2012052409 一、实验目的 1、了解符号、数值、字符等数据类型的差别 2、了解符号运算的特点、优缺点 3、掌握符号变量的创建和运算，以及其运算的基本应用 4、掌握基本的符号绘图指令 二、实验内容 1、指出下面的 M1，M2，M3 分别是什么，并上机验证。 取a=1、b=2、c=3、d=4，M1=[a,b;c,d]，M2=[a,b;c,d]，M3=sym([a,b;c,d]); a=1,b=2,c=3,d=4 a = 1 b = 2 c = 3 d = 4 M1=[a,b;c,d] M1 = 1 2 3 4 M2=[a,b;c,d] M2 = [a,b;c,d] M3=sym([a,b;c,d]) M3 = [ a, b] [ c, d] 结论:M1是矩阵，2是字符串，M3是字符变量。 2、下面2种取值情况下，计算并赋给相应情况下的c1、c2，问c1、c2相等吗，为什么？上机验证。 （1） a1=1010; b1=10-10; （2） 将a1、a2作为符号变量赋给a2、b2; a1=1e10; b1=1e-10; c1=(a1+b1-a1)/b1 c1 = 0 a2=sym(a1); b2=sym(b1); c2=(a2+b2-a2)/b2 c2 = 1 结果:c1~=c2，因为c1=0,c2=1,a1、b1是具体的数值，a2、b2是符号变量。 3、符号表达式中自由变量的确定生成符号变量a、b、x、X、Y、k=3、z=，表达式为。 （1）找出f中的全部自由符号变量 （2）在f中确定最优先的自由符号变量 （3）在f中确定2个和3个自由变量时的执行情况 （4）试通过对各符号变量与x的ASCII值做绝对差值，分析自由变量优秀顺序，能得出什么结论？ syms a b x X Y k=sym(3); z=sym(c*sqrt(w)+y*sin(a)); f=a*z*X+(b*x^2+k)*Y; findsym(f) ans = X, Y, a, b, c, w, x, y findsym(f,1) ans = x findsym(f,2) ans = x,y findsym(f,3) ans = x,y,w double(a) ans = 97 double(b) ans = 98 findsym(f,8) ans = x,y,w,c,b,a,X,Y A=abs(97-120) A = 23 B=abs(98-120) B = 22 MATLAB中，x是默认的首选符号变量，其后的次序排列规则是：与x的ASCⅡ码值之差的绝对值小的字母优先，当差绝对值相同时， ASCⅡ码值大的优先。 4、已知，。做符号运算，化简并计算： （1）；； （2）；；； syms x; f=-x+2*x^2+3*x^3+4*x^4; g=5*x^5-6*x^2+9; A=f*g; A=simple(A) A = (-x+2*x^2+3*x^3+4*x^4)*(5*x^5-6*x^2+9) B=f/g; B=simple(B) B = (-x+2*x^2+3*x^3+4*x^4)/(5*x^5-6*x^2+9) C=f\g; C=simple(C) C = (5*x^5-6*x^2+9)/(-x+2*x^2+3*x^3+4*x^4) D=int(f) D = -1/2*x^2+2/3*x^3+3/4*x^4+4/5*x^5 E=int(g,x,1,2) E = 95/2 F=diff(f,x,2) F = 4+18*x+48*x^2 G=limit(B,x,0) G = 0 5、用ezplot命令绘出在[-4,4]之间的图像。 syms x; y=x^2*sin(x); ezplot(y,[-4,4]); 6、用ezsurf命令绘出单位球的曲面图像 syms a b; x=sin(a)*cos(b); y=sin(a)*sin(b); z=cos(a); ezsurf(x,y,z,[0,pi,0,2*pi]); 三：实验总结 通过本次的实验，掌握了符号计算的运用，了解了符号计算的来源，其次在MATLAB中的符号表达式一定要用 单引号括号起来matlab才能识别。 函数ezpl