MatLab在工科数学中的应用.doc

学校代码： 10128 MatLab在工科数学中的应用结课作业学生姓名：梁文隆学 号：201020509049学 院：化工学院专 业：无机非金属材料工程班 级：2010级1班指导教师：杜雅娟 讲师1、 矩阵运算方面的题（参照第一章练习题3）生成矩阵A,通过修改矩阵A的第1行第3列元素为2,把矩阵扩充为4行5列，其中第4行第3列元素为7，其余列为0,删除第2列元素的操作生成矩阵B, 按行的逆顺序取A的1,2,4,5列,且第四行为自然数1到4，生成新的矩阵C;计算B*C，B.* C, B+C,2*B, B和C的行列式，BC, B/C,并找出A中大于3的元素，且将其替换为1。上机实验步骤输入A=0,1,0,2,1; 3,4,6,8,4; 9,7,3,2,2;B=A;B(1,3)=2;B(4,3)=7;B(:,2)=运行结果如下：输入C=A(3:-1:1,1,2,4,5)运行结果如下：输入C(4,:)=1:4运行结果如下：输入B*C运行结果如下：输入 B.*C运行结果如下：输入 B+C运行结果如下：输入 2*B运行结果如下：输入 det(C)运行结果如下：输入 B/C运行结果如下：输入 BC运行结果如下：输入 det(B) 运行结果如下：2、 绘图方面的题，包括二维和三维图形（参照第二章练习题）（1）、二维图形的绘制用ezplot绘制函数的图形。（2）、三维图形的绘制在同一平面中的两个窗口分别用polar绘制心形线和马鞍面z=x*x-y*y的图形，并以不同的角度观察马鞍面。上机实验步骤（1）、在MatLab命令窗口输入 ezplot(2*(t-sin(t),2*(1-cos(t),0,6*pi)运行结果如下： （1）、在MatLab命令窗口输入clfsubplot(1,2,1)theta=linspace(0,2*pi,30);rho=3*(1-cos(theta);polar(theta,rho,b)subplot(1,2,2)x=-3:0.1:3;y=x; X,Y=meshgrid(x,y);Z=X.2-Y.2;surf(X,Y,Z)shading flat运行结果如下：三、M文件程序设计（参照第三章练习题1）已知函数，计算f(-0.5),f(0),f(1),并作出该函数的曲线图，说明x=1.5,x=0,x=1处的连续性和可导性。上机实验步骤（1）在MatLab命令窗口中，点击菜单FileNew-M-file，打开M文件编辑窗口输入程序如下：function y=fun1(x)if (x=-1)&(x0) y=x+1;elseif x1 y=sin(x)+1;elseif xSave，在默认路径下存盘，文件名为fun1.m（3）在MatLab命令窗口中，输入fun1(-0.5)，fun1(0)，fun1(1)得f(-0.5)0.5000，f(0)1，f(1)1 (4) 在MatLab命令窗口中，输入plot(0,fun1(0),*,1,fun1(1),*,1.5,fun1(1.5),*) hold onfplot(fun1,