matlab实验题汇总要点.docxVIP

湖南科技大学 技能培训3(习 题) 姓 名: 学 号:专 业: 信息与计算科学 学 院: 数学与计算科学学院 指导教师 : 谭 敏 二〇一三 年十二月三十日 第 2 讲： MATLAB 入门 1、用起泡法对 10个数由小到大排序，即将相邻两个数比较 , 将小的调到前头。 解： 代码如下： Untitled1.m clear all ;clc; a=[7 2 1 0 9 4 5 -3 8 6]; n=length(a); for ii=1:n-1 if a(ii+1)=a(ii) t1=a(ii); a(ii)=a(ii+1); a(ii+1)=t1; end for jj=1:n-1 if a(jj+1)=a(jj) t2=a(jj); a(jj)=a(jj+1); a(jj+1)=t2; end end end a 运行结果显示如下： a = 9 8 7 6 5 4 2 1 0 -3 2、有一个 4*5 矩阵 ,编程求出其最大值及其所处的位置。 解： 代码如下： clear;第1页共24页 clc; a=[1 2 3 4 5 34569 67880 12456] max=-1; flage1=0; flage2=0 for i=1:4 for j=1:5 if (a(i,j)max) t=max; max=a(i ,j); a(i,j)=t; flage1=i; flage2=j ； end end end max flage1 flage2 运行结果显示如下： a = 1 2 3 4 5 3 4 5 6 9 6 7 8 8 0 1 2 4 5 6 flage2 = 0 max = 9 flage1 = 2第2页共24页 flage2 = 5 20 3、编程求 n! 。 n 1 解： 代码如下： clear; clc; sum=0; for i=2:11 sum=sum+gamma(i); end sum 运行结果显示如下： sum = 4037913 第 3 讲： MATLAB 作图 1、在同一平面中的两个窗口分别画出心形线和马鞍面。 要求： 1)在图形上加格栅、图例和标注； 2)定制坐标； 3)以不同角度观察马鞍面。 解： 心型线代码如下： clear; clc; i=-pi:0.1:pi; x=2.*(sin(i)-sin(2*i)./2);第3页共24页 y=2.*(cos(i)-cos(i).^2); plot(x,y); axis([-3 3 -4.2 1]);% 制定图形坐标 title( 心形线 );% 当前图顶端加图例 xlabel(x);% 当前图 x 轴加图例 ylabel(y);% 当前图 y 轴加图例 grid on;% 加格栅到当前图上 心型线运行结果显示如下图： 图 3.1（心型线） 马鞍线代码如下： clear; clc; a=10;第4页共24页 b=10; x=-1.5:0.1:1.5; y=-1.5:0.1:1.5; [X,Y]=meshgrid(x,y);%x 行y 列的矩阵 Z=X.^2/a-Y.^2/b; figure(2);% 新建窗口 subplot(221);% 分割当前窗口，作图显示在该分割块上 surf (X,Y,Z); title( 马鞍面 ); xlabel(x); ylabel(y); subplot(222); surf (X,Y,Z); title( 马鞍面 ); xlabel(x); ylabel(y); view(20,30);% 改变视角到（ 20 ，30) subplot(223); surf (X,Y,Z); title( 马鞍面 ); xlabel(x); ylabel(y); view(50,20);% 改变视角到（ 50,20 ） subplot(224); surf (X,Y,Z); title( 马鞍面 ); xlabel(x); ylabel(y); view(60,10);% 改变视角到（ 60,10 ） 马鞍线运行结果显示如下图：第5页共24页 图 3.2 （马鞍面） 2 、 以 不 同 的 视 角 观 察 球 面 x2+y2+z2=r2 和圆柱面 x2+y2=rx 所围区域 . 解： 代码如下： clear; clc; r=2; [x0,y0,z0]=sphere(50); x=r*x0;y=r*y0;z=r*z0; surf(x,y,z);hold on ; f=@(x,y)x.^2+y.^2-r*x; ezsurf(f,[-1.5,2.5],[-2,1.5]); axis equal ; view(