文档介绍：:..宁德师范学院物理与电气工程系实验报告(2013—2014学年第二学期)课程名称 MATLAB 实验名称 实验一符号运算专业 电气工程及其自动化年级11电气2班 学号 姓名 指导教师 实验日期 2014年5月13日实验一符号运算一、 实验0的1、 了解符号、数值、字符等数据类型的差别2、 了解符号运算的特点、优缺点3、 掌握符号变量的创建和运算,以及其运算的基本应用4、 掌握基本的符号绘图指令二、 实验内容1、指出下面的Ml,M2,M3分别是什么,并上机验证。取a=l、b=2、c=3、d=4,Ml=[a,b;c,d], M2=,[a,b;c,d]’M3=sym(,[a,b;c,d]’);程序:>>a=l,b=2,c=3,d=4a=1b=2c=3d=>>Ml=[a,b;c,d]1 23 4»M2=’[a,b;c,d]’M2-[a,b;c,d]>>M3=sym(’[a,b;c,d]’)M3=[a,b][c,d]结论:Ml是矩阵,M2是向量,M3是数组。a+b-a2、下面2种取值情况下,计算b"-并赋给相应情况下的cl、c2,问cl、c2相等吗,为什么?上机验证。(1)al=10^10;bl=107-10):(2)将al、a2作为符号变量赋给a2、b2;程序:»al二lelO;bl二le-10;cl=(al+bl-al)/bla2=sym(al):b2=sym(bl);c2=(a2+b2-a2)/b2因力cl=0,c2=1,所以cl不等于c2。al、bl是具体的数值,a2、b2是符号变量3、符号表达式中自巾变量的确定生成符号变量a、b、x、X、Y、k=3、cVvv+ysina表达式力/=az^+(如2+(1)找出f中的全部自由符号变量(2)在f中确定最优先的自巾符号变量(3)在f中确定2个和3个自由变量时的执行情况(4)试通过对各符号变景与x的ASCII值做绝对差值,分析自由变景优秀顺序,能得出什么结论?程序:>>symsabxXYk=sym(’a’);z=sym(,c*sqrt(w)+y*sin(a)?):f=a*z*X+(b*x"2+k)*Y;»findsym(f) %找出f中的全部自由符号变量ans=X,Y,a,b,c,w,x,y»findsym(f,1) %在f中确定最优的自由符号变景>>findsym(f,2)%在f中确定2个自由变量时的执行情况ans=x,y»findsym(f,3) %在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二22MATLAB中,x是默认的首选符号变量,其后的次序排列规则是:与x的ASCII码值之差的绝对值小的字母优先,当差绝对值相同时,ASCII码值大的优先。4、已知/U)=-x+2x2+3x3+4*4,g(x)=5x5-6x2+9。做符号运算,化简并计算:(1)f(x)-g(x);f(x)/g(x);f(x)\g(x)(2)程序:>>symsx;»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)(5*x"5-6*x"2+9)/(-x+2*x"2+3*x3+4^x^4)»D-int(f)D=-l/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=05、用ezplot命令绘出>’=x2sin<x)在[-4,4]之间的图像。程序:x=-4::4;y=x."(x);plot(x,y):图像:6、用czsurf命令绘出单位球的曲面图像(选做)>>symsab;ezsurf(sin(a)*cos(b),sin(a)*sin(b),cos(a),[0,2本pi,0,2*pi])图像:三、实验总结在MATLAB中符号表达式一定要用’’单引号括号起来matlab才能标识;’’的内容可以是符号表达式,也可以使符号方程;函数ezplot它无需数据准备,直接画出函数图形,基本调用格式为czplot(f),其中f是字符串或代表数学函数的符号