资源描述：
《数学及应用数学数学软件实验》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、《数学软件》实验指导书《数学软件》实验指导书实验一曲线画图实验目的：熟悉matlab 环境，学习和掌握用matlab软件进行曲线绘图。实验内容：画出下列曲线的图形：1、disp('画曲线1')symsxy='x^4-5*x^3+7*x^2-3*x';fplot(y,[-0.5,1.5])2、disp('画曲线2')symsxyf='x^4+y^3=1';ezplot(f,[-1.5,1.5])3、disp('画曲线3')symstx='t^2-t+1';y='t^2+t+1';ezplot(x,y,[-2,2])4、disp('画曲
2、线4')a=2;theta=0:0.01:3*pi;rhol=a*theta;polar(theta,rhol)5、disp('画曲线5')x=-10:0.01:10;y=(x<1).*(2+x)+(x>=1).*(x.*sin(x));plot(x,y)《数学软件》实验指导书第7页《数学软件》实验指导书实验二极限导数积分实验目的：学习和掌握用matlab工具求解极限、导数、积分问题实验内容：1、求函数的极限（1）；symsxdisp('求极限（1）')limit((x^2-6*x+8)/(x^2-5*x+4),x,4)（2）；sy
3、msxdisp('求极限(2)')limit((x+1/x)^x,x,inf)2、求导数：（1）；symsxdisp('求导数(1)')y='3*x^2-5*x+1';yxx=diff(y,x,2)（2）symsatdisp('求导数(2)')x='a*(t-sin(t))';y='a*(1-cos(t))';yx=diff(y,t)/diff(x,t)yxx=diff(yx,t)/diff(x,t)（3）symsabxydisp('求导数(3)')f='b^2*x^2+a^2*y^2-a^2*b^2';yx=-diff(f,x)/
4、diff(f,y)yxx=diff(yx,x)+diff(yx,y)*yx3、求积分：（1）；第7页《数学软件》实验指导书symsaxdisp('求积分(1)')int(x^2/sqrt(x^2+a^2),x)（2）symsxdisp('求积分(2),符号解')int((x*sin(x))^2,x,0,pi)x=0:0.01:pi;disp('求积分(2),数值解')y=(x.*sin(x)).^2;trapz(x,y)《数学软件》实验指导书实验三级数曲面绘图实验目的：学习和掌握用matlab工具求解方程、极值、级数、曲面绘图问题实
5、验内容：1、问题：（1）求方程在限制条件下的根；symsxdisp('1、方程求根(1)')f='x*log(x+sqrt(x^2-1))=sqrt(x^2-1)+0.5*x'x=solve(f,x)（2）求函数的极值：symsxdisp('1、求极值(2)')f='x*sin(x^2-x-2)';fplot(f,[-2,2])disp('极小值点和极小值')[x1,f1]=fminsearch(f,-1.5)[x2,f2]=fminsearch(f,1.5)disp('极大值点和极大值')[x3,f3]=fminsearch('-
6、x*sin(x^2-x-2)',0.5);x3f3=-f32、讨论下列级数的敛散性：；symsndisp('求级数的敛散性')第7页《数学软件》实验指导书a=symsum(log(n)/n^3,n,1,inf)eval(a)disp('收敛')3、画出空间曲面：（1）范围内的图形；[x,y]=meshgrid(-30:0.1:30);disp('画空间曲面(1)')z=10*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2);meshc(x,y,z)（2）双曲抛物面disp('画空间曲面(2)双曲抛物面'
7、)[x,y]=meshgrid(-30:0.1:30);z=(x.^2-y.^2)/3;meshc(x,y,z)《数学软件》实验指导书实验四多元微积分等实验目的：学习和掌握用matlab工具求解多元微积分、微分方程问题实验内容：1、求函数时的全微分.symsxydxdydisp('求函数的全微分')z='log(1+x+y)';dz=diff(z,x)*dx+diff(z,y)*dyx=1y=2dz=subs(dz)2、计算,积分区域D是由抛物线所围成的区域；disp('二重积分')symsxyf1='y^2-x';f2='y-x+
8、2';ezplot(f1,[-2,4.5])holdonezplot(f2,[-2,4.5])[x,y]=solve(f1,f2)symsxy第7页《数学软件》实验指导书s=int(int('x*y','x','y^2','y+2')
