《matlab中期作业_附源代码》.pdfVIP

一、题目如下 1. 在同一个图中绘制基本函数 sin(x)和 cos(x)函数，注意使用不通颜色的曲线区分，并有图 题、纵横坐标等信息；自己举例绘制三位矢量图、三维曲面图及等值线图。 2. 利用 LU 分解法计算线性方程组的解。 4x+5.3y-5.6z-3m-3.4n=100.06 5x-2.1y+3.2z+4m-8n=-75.72 2x-4y-7.2z-5m-2.4n=98.2 5x-3y-8z+2.3m+3n=57.1 4.2x-3y-2n=3.72 3 2 3. 采用不动点迭代法计算非线性方差x +4x -10=0，在区间[1,2]上的一个根。 4. 利用函数的递归调用，求n ！。 5. 用筛选法求某自然数范围内的全部素数。 （素数是大于1，且除了 1 和它本身以外，不能被其他任何整数所整除的整数。用筛选法求 素数的基本思想是：要找出 2~m 之间的全部素数，首先在 2~m 中划去2 的倍数(不包括 2) ， 然后划去 3 的倍数(不包括 3)，由于 4 已被划去，再找 5 的倍数 (不包括 5)，…，直到再划 去不超过的数的倍数，剩下的数都是素数。） 二、1~5 题解答及源程序如下： 1.1 在同一个图中绘制基本函数 sin(x)和 cos(x)函数：（源程序m 文件见 Ti_1_1 ） ①做y sin x 和y cos x 图如下： ②源程序Ti_1_1 ： clc; x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1, g,x,y2,r-.); title(正弦曲线和余弦曲线 ); xlabel(横坐标X); ylabel(纵坐标Y); legend(sin(x),cos(x)); 1.2绘制三维矢量图：（源程序m文件见Ti_1_2 ） ①作z x 2 y 2 的三维矢量图如下： ②源程序Ti_1_2 ： clc;clear; [x,y]=meshgrid(-4:0.2:4); z=x.^2+y.^2; [x1,y1,z1]=surfnorm(x,y,z); quiver3(x,y,z,x1,y1,z1); title( z=x^2+y^2三维矢量图 ); xlabel( X轴 ); ylabel( Y轴 ); zlabel( Z轴 ); 1.3绘制三维曲面图：（源程序m文件见Ti_1_3 ） ①作z sin x cos y 的三维网格曲面图如下： ②源程序Ti_1_3 ： clc;clear; [x,y]=meshgrid(0:0.15:3*pi); z=sin(x)+cos(y); mesh(x,y,z); title(z=sinx+cosy三维网格曲面图形 ); xlabel( X轴 ); ylabel( Y轴 ); zlabel( Z轴 ); 1.4绘制等值线图：（源程序m文件见Ti_1_4 ） 2 2 ①作z e(x y ) 的三维等值图如下： ②源程序Ti_1_4 ： clc;clear; [x,y]=meshgrid(-4:0.2:4); z=exp(-x.^2-y.^2); contour3(x,y,z,20); title( z=exp(-x^2-y^2)三维等值图 ); xlabel( X轴 ); ylabel( Y轴 ); zlabel( Z轴 ); 2.利用LU分解法计算线性方程组的解：（源程序m文件见Ti_2 ） ①求解结果如下： ②源程序Ti_2 ： clear;clc; fprintf(*********求解线性方程组LU分解后U矩阵如下 :*********\n,); fprintf(4x+5.3y-5.6z-3m-3.4n=100.06\n,); fprintf(5x-2.1y+3.2z+4m-8n=-75.72\n,); fprintf(2x-4y-7.2z-5m-2.4n=98.2\n,); fprintf(5x-3y-8z+2.3m+3n=57.1\n,); fprintf(4.2-3y-2n=3.72\n\n,); fprintf(*********要求解线性方程组系数矩阵A如下：*********,); A=[4