数学软件实验报告实验四.doc

-数学软件实验报告实验四MATLAB的符号方程求解与符号绘图一实验目的MATLAB不仅能够进展符号微积分运算和符号积分变换，还能够进展求解符号代数方程组、求解符号微分方程组，还有符号绘图功能、图示化函数计算器以及MAPLE命令的使用和MAPLE语言接口。本次实验的目的是学会使用MATLAB求解符号代数方程和方程组、求解符号解微分方程和微分方程组，学会使用函数计算器的各种组合运算，熟悉MAPLE命令的使用和MAPLE语言接口。二实验要求熟练掌握函数solve和dsolve的使用，并熟悉这两个函数各个参数的变化，学会使用函数maple、mapleinit、mhelp和precread函数的使用过程，并熟悉with和readlib两个maple函数的功能和使用方法，掌握MAPLE库函数的分类。熟练掌握二维符号绘图、极坐标符号绘图、三维符号曲线绘图、三维符号网格绘图、三维符号曲面绘图、符号等高线图的绘制等，了解符号函数图形计算器的使用。三实验容第三章 3.7 符号代数方程求解线性方程组的符号解法>> *,y,z=solve('10*-y=9','-*+10*y-2*z=7','-2*y+10*z=6')* = 473 - 475y = 91 - 95 z = 376 - 475非线性方程符号的解法>> *0=0.5 0.5;>> fsolve('fc',*0)ans = 0.5803 0.18733.8 符号微分方程求解>> dsolve('D*=-a*')ans = _C1 e*p(-a t)3.9 符号函数的二维图符号函数的建议绘图函数ezplotezplot('erf(*)') 绘制函数图函数fplot*=0;.05;1;fplot('tan(*),sin(*),cos(*)',2*pi*-1 1 -1 1);第四章 图形处理功能4.1 二维图形根本绘图命令y=rand(100,1);>> plot(y)*=0:0.01*pi:pi;>> y=sin(*).*cos(*);>> plot(*)*=0:0.01*pi:pi;y=sin(*'),cos(*');>> plot(*',*',y)*=1:0.1*pi:2*piy=sin(*)z=cos(*)plot(*,y,'-k',*,z,'-.rd')特殊的二维图形函数二维特殊图形函数4.2 三维图形根本绘图命令三维罗线图>> *=0:pi/50:10*pi>> y=sin(*)>> z=cos(*)>> plot3(*,y,z)2绘制参数为矩阵的三维图>> *,y=meshgrid(-2:0.1:2,-2:0.1:2)>> z=*.*e*p(-*.2-y.2)>> plot3(*,y,z)3网图函数使用mesh函数绘制三维面图>> *=-8:0.5:8;y=*'；>> a=ones(size(y)*;>> b=y*ones(size(*);>> c=sqrt(a.2+b.2)+eps;>> z=sin(c)./c;>> mesh(z)使用mesh函数绘制Hilbert矩阵三维面图>> z=hilb(10);>> mesh(z)使用meshchanshu 绘制三维面图>> *,Y=meshgrid(-4:0.5:4);>> Z=sqrt(*.2+Y.2);>> meshc(Z)特殊的三维图形函数*=2,4,6,8;>> pie3(*,0,0,1,0)*,Y=meshgrid(-4:0.5:4);Z=sqrt(*.2+Y.2);>> surfc(*,Y,Z)*,Y=meshgrid(-4:0.5:4);>> contour3(peaks(*,Y),25)*=0:pi/20:pi*3;>> r=5+cos(*);>> a,b,c=cylinder(r,30);>> mesh(a,b,c)a,b,c=sphere(40);t=abs(c);surf(a,b,c,t);a*is('equal');a*is('square');colormap('hot')四、实验总结在本次试验中，我做了符号代数方程的求解，符号微分方程求解，符号函数的二维图，二维图形、三维图形的处理功能，用到了ezplot,fplot,sphere,contour3,mesh等语言，绘制了一些二维三维的图形，对图形的绘制有了一定的了解。 这次实验根本没有出现问题，相对于前几个实验，感觉图形绘制更加好做一些。但是在语言运用上还是有一些不熟悉，觉得如果没有书可能就会分不清要用哪个语言函数，需要加深对语言函数的理解和应用，做到能记住那些语言函数的功能。. z.