答案详解实验指导书Matlab软件应用与开发分析.docVIP

《MATLAB软件应用与开发》实验指导书 实验1 MATLAB绘图 (一)、实验类型：验证型 (二)、实验类别：基础实验 (三)、每组人数：1 (四)、实验要求：选修 (五)、实验学时：3个学时 (六)、实验目的： 学会使用MATLAB进行二维图形以及三维图形绘制。 (七)、预备知识： 本实验所用MATLAB命令 利用冒号表达式生成行向量a=初始值:步长:终值 二维显函数图形绘制：直角坐标系plot(x,y)，极坐标系polar(theta,rho) 二维隐函数图形绘制：ezplot(隐函数表达式) 三维曲线绘制：plot3(x,y,z) 三维曲面绘制：三维表面图[x,y]=meshigrid(v1,v2); %生成网格数据 surf(x,y,z) 三维网格图[x,y]=meshigrid(v1,v2);mesh(x,y,z) (八)、内容与要求： 1.1 已知迭代模型，试写出求解该模型的M函数。如果取迭代初值为x(0)=0，y(0)=0，那么请进行30000次迭代求出一组x和y向量，然后在所有的x(k)和y(k)坐标处点亮一个点（使用红色点作图），并使用蓝色*号和黑色圆圈分别标注出起点和终点，最后绘制出所需的图形。提示：所绘制出的图形又称作Henon引力线图，它将迭代出来的随机点吸引到一起，最后得出貌似连贯的引力线图。 1.2 按要求绘制出下列图形 （1）y=sinx，其中x的定义域为(-2pi,2pi)，步长0.1，蓝色点线，尺寸为6个标准单位，标注出x，y轴，图形名称； （2）y=cos(2x)+5x^2+ln(x)，其中x的定义域为(-pi,pi)，红色虚线，尺寸为4个标准单位，标注出x，y轴，图形名称； （3）将y1=sinx，y2=cos(2x)+5x^2+ln(x)绘制在同一直角坐标系内，其中y1为黄色*，尺寸为4个标准单位，y2位黑色实线，尺寸为6个标准单位。 1.3 分别选取合适的θ范围，将下列极坐标方程以2行2列的形式绘制在同一极坐标系内： （1），（2）， （3），（4） 1.4 分别绘制出下列隐函数的图形： （1）x^2+y^2=3xy^2，（2）x^3-x^2=y^2-y， （3）e^(-(x+y)^2)+π/2sin(5x+2y)=0， （4）(x^2-y^2+xy)e^(-x^2-y^2-xy)=0， （5）分别将（1）和（2）以及（3）和（4）画在同一图形内。 1.5 请分别使用surf函数和mesh函数绘制出xy和sinxy的三维表面图和三维网格图。 1.6 在图形绘制语句中，若函数为不定式NaN，则相应的部分不绘制出来。试利用该规律绘制z=cos(xy)+3xy^2的三维表面图，并剪切下x^2+y^2≤0.25的部分。 (九)、实验操作 习题1.1 x(1,1)=0;y(1,1)=0; for k=1:29999 x(1,k+1)=1+y(1,k)-1.4*x(1,k)^2; y(1,k+1)=0.3*x(1,k); end plot(x(1,1),y(1,1),b*,...... x(1,:),y(1,:),r.,...... x(1,end),y(1,end),ko); xlabel(x轴);ylabel(y轴);title(Henon引力线图) 习题1.2 （1） x1=-2*pi:0.1:2*pi; y1=sin(x1); plot(x1,y1,b.,LineWidth,6); xlabel(x轴);ylabel(y轴); title(y=sin(x)); （2） x2=-pi:0.1:pi; y2=cos(2*x2)+5*x2.^2+log(x2); plot(x2,y2,r-,LineWidth,4); xlabel(x轴);ylabel(y轴); title(y=cos(2x)+5x^2+ln(x)); Warning: Imaginary parts of complex X and/or Y arguments ignored （3）plot(x1,y1,y*,LineWidth,4); hold on; plot(x2,y2,k,LineWidth,6); xlabel(x轴);ylabel(y轴); Warning: Imaginary parts of complex X and/or Y arguments ignored 习题1.3 theta1=-2:0.01:2; theta2=-4*pi:0.01:4*pi; theta3=-2*pi:0.01:2*pi; theta4=-10:0.01:10; rho1=