矩阵实验室作业.docVIP

《MatLab在工科数学中的应用》 结课报告 学生姓名： 李 学 院： 能源与动力工程学院 班 级： 环科09 指导教师：曹艳 2011年 12 月 3日 一 学习的主要内容和目的 了解MatLab软件的功能,熟悉MatLab软件的各菜单、工具栏及常用命令的使用。掌握MatLab有关矩阵的创建方法、矩阵的基本运算符、矩阵的计算函数。掌握MatLab的符号运算。熟练掌握二维、三维图形的绘制；掌握简单动画的制作；了解分形几何学，绘制Koch雪花曲线和Minkowski“香肠”曲线。熟炼掌握MatLab程序设计的顺序、分支和循环结构;熟炼掌握脚本M文件和自定义函数的设计和使用；复习高等数学中有关函数极限、导数、不定积分、定积分、二重积分、级数、方程近似求解、常微分方程求解的相关知识.通过作图和计算加深对数学概念：极限、导数、积分的理解.学会用MatLab软件进行有关函数极限、导数、不定积分、级数、常微分方程求解的符号运算；了解数值积分理论,学会用MatLab软件进行数值积分;会用级数进行近似计算.复习线性代数中有关行列式、矩阵、矩阵初等变换、向量的线性相关性、线性方程组的求解、相似矩阵及二次型的相关知识.学会用MatLab软件进行行列式的计算、矩阵的基本运算、矩阵初等变换、向量的线性相关性的判别、线性方程组的求解、二次型化标准形的运算. 二 作图应用 1描点作图(plot) y=sin3x在区间[0,2*pi]画出正弦函数y的图形。 （高数上p154.1（1）） x=0:pi/21:2*pi; y=sin(3*x); plot(x,y) 2显函数作图（fplot） 画y=200*sin(x)/x的图形。（高数上P90、例1） fplot(200*sin(x)/x,[-20,20]) 3隐函数作图（ezplot） 做sqrt(x*x+y*y)=3图，（高数上p36、例3）. ezplot(sqrt(x*x+y*y)=3) 4参数方程作图(ezplot) 在[0,2*pi]上,画x=sint,y=cos2*t,(高数上P123 7题（1))图 输入： ezplot( sin(t), cos(2*t),[0,2*pi]) 5极坐标作图 作r=6a图 输入; theta=linspace(0,2*pi,30); rho=6*theta; polar(theta,rho,b) 6 空间曲面作图 作 z=x2+y2的图 （高数下 P180 3题） 输入： x=-3:0.1:3; y=x; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; surf(X,Y,Z) 三 高等数学应用 1 极限问题 求sin（x)/x的极限（高数 P110 3题） syms x limit(sin(x)/x,x,0) 运行结果 ans =1 2 求导数问题,求y=x*sin(1/x)的导数（高数 P87 5题） syms x; diff(x*sin(1/x),x,1);%求一阶导数 ans diff(x*sin(1/x),x,2);%求二阶导数 ans 运行结果 ans =sin(1/x)-1/x*cos(1/x) ans =-1/x^3*sin(1/x) 3 求不定积分问题，求sqrt(x^2+1)的不定积分(高数上 P205 例1) syms x; int(sqrt(x^2+1),x) 运行结果 ans =1/2*x*(x^2+1)^(1/2)+1/2*asinh(x) 4 求定积分问题求sqrt(x^2+1),x,1,2)二重积分(高数上 P249 例1) syms x; int(sqrt(x^2+1),x,1,2);%二重 ans 运行结果 ans =5^(1/2)-1/2*log(-2+5^(1/2))-1/2*2^(1/2)-1/2*log(2^(1/2)+1) 5 求偏导数问题，F=x/sqrt(x^2+y^2)，求Fx,Fy,Fx2,Fy2，(高数 P80 3) syms x; int(sqrt(x^2+1),x,1,2);%二重 syms F x y Fx Fy Fx2 Fy2 F=x/sqrt(x^2+y^2); Fx=diff(F,x);%求一阶偏微分 Fy=diff(F,y); Fx2=diff(F,x,2);%求二阶偏微分 Fy2=diff(F,y,2); Fx,Fy,Fx2,Fy2 运行结果 Fx =1/(x^2+y^2)^(1/2)-x^2/(x^2+y^2)^(3/2) Fy =-x/(x^2+y^2)^(3/2)*y Fx2 =-3/(x^2+y^2)^(3/2)*x