东莞理工matlab实验报告1.doc

东莞理工学院实验报告系（院》、专业班级=2013电子卓越1班生名：蔡少委日期：2013-2014学年第2学期，1-3周课程名称：MATLAB与应用学号：201341301218 成绩:实验室名称:计算机房实验名称：MATLAB的符号计算同组同学：无 指导教师：姜鸣1、实验目的1） 熟悉MATLAB的符号运行环境2） 掌握常用的符号计算算符和函数指令3） 掌握MATLAB符号工具求极限、导数、积分和级数和;4） 掌握MATLAB符号工具进行符号矩阵分析；5） 掌握MATLAB符号工具计算线性方程的符号解2、实验设备安装MATLAB 7.7及以上版本软件的计算机3、实验原理应用MATLAB软件进行符号计算的相关指令操作：1） 使用sym和syms指令创建符号对象；2） 使用vpa指令设定符号数字的右效数字位；3） 使用limit、diff、int和symsum指令求极限、导数、积分和级数;4） 使用det、eig、poly指令对符号矩阵进行分析；5） 使用solve指令求代数方程的符号解4、实验内容和步賴OOOOOOOO1）求1分，并进而用该符号解求Ik=()k=()l3k的准确值k=0OO2）对于x〉0，求Ik=02k + l2々十13） （1）通过符号计算求｝= |sinZ的导数（2）然后根据此结果，求&dt dtt=(V4）计算二重积分+5）求符号矩阵A= u21 6/22 a23的行列式值和逆，所得结果应采用“子表达式置换”洁化。

6)求方程x2+y2=l，m = 2的解5、实验数据记录及(分析》讨论(将MATLAB代码写在下面)1，代码:syms x k; fl=x"k;f2=subs(fl,x,-1/3); f3=subs (fl,x, 1/pi); f4=subs(fl,x, * 3 *);sl=symsum (f2, k, 0, inf); s2=symsum(f3,k,0,inf); s3=symsum(f4,k,0,inf);disp (si); disp (s2); disp(s3);结果：3/4pi/(pi - 1)Inf2,代码：syms x positive; syms k;f5=2/(2*k+l)*((x-1)/(x+1))A(2*k+l); s5=symsum(f5,k,0,inf); disp(s5);结果：piecewise([abs(x - 1) < x + 1，2*atanh((x - l)/(x + 1))])3，代码；syms y t;syms d positive;sl=diff(abs(sin (t)));disp (si);s2=limit (si,t, 0, 1 left1); s3=limit(si,t,pi/2); disp (s2); disp(s3);结果：sign(sin(t))*cos(t)-104,代码：syms x y; fl=xA2+yA2;sl=int(int(fl,x,l,x八2) , y, 1,2); disp (si);结果：X八6/3 + (7*xA2)/3 - 8/35代码：syms all al2 al3 a21 a22 a23 a31 a32 a33;A=[all al2 al3;a21 a23 a23;a31 a32 a33]; sl=det(A);s2=inv(A); s3=subexpr(A);disp(si); disp (s2); disp(s3);结果：al2*a23*a31 - al2*a21*a33 - all*a23*a32 + al3*a21*a32 4- all*a23*a33 - al3*a23*a31[(a23*a32 - a23*a33)/(all*a23*a32 + al2*a21*a33 - al2*a23*a31 - al3*a21*a32 - all*a23*a33 + al3*a23*a31)， (al2*a33 - al3*a32)/(al I*a23*a32 + al2*a21*a33 - al2*a23*a31 -al3*a21*a32 - all*a23*a33 + al3*a23*a31)，-(al2*a23 - al3*a23)/(all*a23*a32 + al2*a21*a33 - al2*a23*a31 - al3*a21*a32 - all*a23*a33 + al3*a23*a31)l[(a21*a33 - a23*a31)/(alI*a23*a32 + al2*a21*a33 - al2*a23*a31 - al3*a21*a32 - al I*a23*a33 + al3*a23*a31)，-(all*a33 - al3*a31)/(all*a23*a32 + al2*a21*a33 - al2*a23*a31 - al3*a21*a32 -al I*a23*a33 + al3*a23*a31)，(al l*a23 - al3*a21)/(al I*a23*a32 + al2*a21*a33 - al2*a23*a31 -al3*a21*a32 - all*a23*a33 + al3*a23*a31)][-(a21*a32 - a23*a31)/(al I*a23*a32 + al2*a21*a33 - al2*a23*a31 - al3*a21*a32 - all*a23*a33 + al3氺a23*a31)，(all*a32 - al2*a31)/(al I*a23*a32 + al2*a21*a33 - al2*a23*a31 - al3*a21*a32 -all*a23*a33 + al3*a23*a31)，(al2*a21 - all*a23)/(all*a23*a32 + al2*a21*a33 - al2*a23*a31 -al3氺a21 氺a32 - al l*a23氺a33 + al3*a23*a31)][all，al2，al3][a21,a23, a23][a31, a32, a33]6；代码：syms x y;sl=solve (H+yK，，，x*y=2，)； disp (si);disp (丨 si. x 丨);disp (si. x); disp(si.y);disp(s1.y);结果：x: [4x1 sym] y: [4x1 sym]sl.x((15A(l/2)*i)/2 + l/2)A(l/2)/2 - ((15A(l/2)*i)/2 + l/2)A(3/2)/2 -((15A(l/2)*i)/2 + l/2)A(l/2)/2 + ((15A(l/2)*i)/2 + l/2)A(3/2)/2(1/2 - (15A(l/2)*i)/2)A(l/2)/2 - (1/2 - (15八(l/2)*i)/2)八(3/2)/2 -(1/2 - (15A(l/2)*i)/2)A(l/2)/2 + (1/2 - (15A(l/2)*i)/2)A(3/2)/2sl.y((15A(l/2)*i)/2+ l/2)A(l/2)-((15A(l/2)*i)/2+ 1/2)八(1/2)(l/2-(15A(l/2)*i)/2)A(l/2)-(l/2-(15A(l/2)*i)/2)A(l/2)。