资源描述:
《MatLab在工科数学中的应用.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、学校代码:10128《MatLab在工科数学中的应用》结课作业学生姓名:梁文隆学号:201020509049学院:化工学院专业:无机非金属材料工程班级:2010级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中
2、大于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
3、)运行结果如下:一、一、绘图方面的题,包括二维和三维图形(参照第二章练习题)(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*p
4、i,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)shadingflat运行结果如下:三、M文件程序设计(参照第三章练习题1)已知函数,计算f(-0.5),f(0),f(1),并作出该函数的曲线图,说明x=1.5,x=0,x=1处的连续性和可导性。上机实验步骤(1)在MatLab命令窗口中,点击菜单File->New->M-file,打
5、开M文件编辑窗口输入程序如下:functiony=fun1(x)if(x>=-1)&(x<0)y=x+1;elseifx<1y=sin(x)+1;elseifx<=2y=x.*x;end(2)在M文件编辑窗口中,点击File->Save,在默认路径下存盘,文件名为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
6、.5,fun1(1.5),'*')holdonfplot('fun1',[-1,2])运行结果如下:由图像观察得各点的连续性和可导性:在x=1.5处,连续可导;在x=0处,连续可导;在x=1处,不连续不可导。