佛山科学技术学院Matlab实训报告工科.docVIP

Matlab 实训报告 班级： 姓名： 学号： 一、、…………………… 1） 和的图像，区间取[0,5]。 2)用Matlab命令求函数的极限：=?，=？，=？（注：开方函数用sqrt） 3）用surf，mesh绘制曲面，。 4) 求函数的一阶和二阶导数，然后在同一坐标系下作出一阶和二阶导函数的图形。 5）在同一坐标系中，取自变量区间[15，200]，用不同的颜色或线型分别画出下面四个函数的图形： 观测当增大时图形的走向。(注：可选用’r’,’g’,’b’颜色) 6）将函数（）的图形分别向左和向右平移2个单位；然后将函数图形向上和向下平移2个单位。 7）求矩阵的行列式、逆矩阵、特征值和特征向量。 8)将函数（）的图形以原点为中心，逆时针旋转60度角。 9）用莱布尼茨级数公式计算的近似值。 选做题： 10）对莱布尼茨级数公式利用while循环编写计算(的近似值的Matlab程序，要求pi与公式计算出近似值的误差的绝对值小于0.0001。 三、 和的图像，区间取[0,5]。 解题： x=0:0.01:5;y1=sin(5*x);y2=cos(3*x);y3=sin(5*x).*cos(3*x);y4=sin(5*x)+cos(3*x);subplot(2,2,1)plot(x,y1)subplot(2,2,2)plot(x,y2)subplot(2,2,3)plot(x,y3)subplot(2,2,4)plot(x,y4) 第二题2)用Matlab命令求函数的极限：（1）=?，（2）=？，（3）=？ 解题： 实验程序为： （1）syms x; limit(((1+x)/x)^2*x,x,inf) （2）syms x; limit((1-sqrt(1-x^2)/exp(x)-cos(x)),x,0) （3）syms x; limit(((3*x^2+5)/5*x+3)*sin(2/x),x,inf) 结果为： （1） （2） （3） 第三题，用surf，mesh绘制曲面，。 解题： 实验程序： x=-5:5;y=-5:5; [x,y]=meshgrid(x,y);z=2*x.^2+y.^2; surf(x,y,z) xlabel(横坐标x);ylabel(纵坐标y);zlabel(函数值z) 结果如下： 第四题，求函数的一阶和二阶导数，然后在同一坐标系下作出一阶和二阶导函数的图形。 解题：程序如下， 一阶：syms x; diff(2*x.^3-12*x+6,x,1) 输出结果： 二阶：syms x; diff(2*x.^3-12*x+6,x,2) 输出结果： 作图： x=-4:0.1:4; y1=6*x.^2 - 12; y2=12*x; plot(x,y1,x,y2,:) 输出结果： 第五题，在同一坐标系中，取自变量区间[15，200]，用不同的颜色或线型分别画出下面四个函数的图形： 观测当增大时图形的走向。(注：可选用’r’,’g’,’b’颜色)。 解题：程序如下； x=15:0.1:200; y1=((x+1)/x).^x; plot(x,y1,-r); hold on; y2=((x+1)/x).^(x+1); plot(x,y2,-g); hold on; y3=((x+2)/(x+1)).^x; plot(x,y3,-b); hold on; y4=2.7183; plot(x,y4,+k); legend(y1=((x+1)/x).^x,y2=((x+1)/x).^(x+1),y3=((x+2)/(x+1)).^x,y4=2.7183); hold off 输出结果： 第六题，将函数（）的图形分别向左和向右平移2个单位；然后将函数图形向上和向下平移2个单位。 解题：程序如下， x=-4:0.1:4;y=exp(3-x.^2); x1=x+2; x2=x-2; y1=y+2; y2=y-2; plot(x,y,x1,y,:,x2,y,-.,x,y,x,y1,:,x,y2,-.) xlabel(x);ylabel(y) 输出结果： 第七题，求矩阵的行列式、逆矩阵、特征值和特征向量。 解题：程序如下， a=[2 7 3;3 9 4;1 5 3]; det(a)——计算行列式 输出结果： inv(a)——计算逆矩阵 输出结果： eig(a)——求特征值 输出结果： [x,d]=eig(a)——求特征向量 输出结果： 其中，X是矩阵，它的列是相应的特征向量，d的对角线元素是特征值。 第八题，将函数（）的图形以原点