matlab实验一参考模板

1、开课学院、实验室： 实验时间 ： 年 月 日课程名称数学实验实验项目名 称实验一基础实验实验项目类型验证演示综合设计其他指导教师龚劬成 绩实验目的1 熟悉MATLAB软件的用户环境；2 了解MATLAB软件的一般目的命令；3 掌握MATLAB数组操作与运算函数；4 掌握MATLAB软件的基本绘图命令；5 掌握MATLAB语言的几种循环、条件和开关选择结构。基础实验一、实验内容及实验过程（一般应包括实验原理或问题分析，算法设计、程序、计算、图表等， 实验结果及分析）1MATLAB软件的数组操作及运算练习；2直接使用MATLAB软件进行作图练习；3用MATLAB语言编写命令M-文件和函数M-文件。

2、1设有分块矩阵，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证。程序E=eye(3);R=rand(3,2);O=zeros(2,3);V=rand(1,2);S=diag(V);A=E R;O S;X=A*A;Y=E R+R*S;O S*S;.if(X=Y).disp('p=1');.else disp('p=0');.end输出结果：p=1 得证2某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如表1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该10种商品的总收入和总利润。

3、表1.1货号 1 2 3 4 5 6 7 8 9 2 / 9单件进价7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30单件售价11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50销量568 1205 753 580 395 2104 1538 810 694程序 a=7.15 8.25 3.20 10.30 6.68 12.03 16.85 17.51 9.30;.b=11.10 15.00 6.00 16.25 9.90 18.25 20.80 24.15 15.50;.c=568 1205 75

4、3 580 395 2104 1538 810 694;.d=c.*(b-a);e=b.*c;m,im=min(d),M,iM=max(d),.rd,ir=sort(e),sum1=sum(b.*c),sum2=sum(d),运行结果m = 1.2719e+003 im =5 M = 1.3087e+004 iM = 6 rd =1.0e+004 *0.3911 0.4518 0.6305 0.9425 1.0757 1.8075 1.9562 3.1990 3.8398 ir =5 3 1 4 9 2 8 7 6sum1 = 1.4294e+005 sum2 =4.6052e+004所以，第

5、五种商品利润最小，第六种商品利润最大，商品按收入排序如上所示，总收入为1.4294e+005，总利润为4.6052e+0044用subplot分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题，1）概率曲线 ；2）四叶玫瑰线 r=sin2q；3）叶形线 4）曳物线 。程序：x1=linspace(-3,3,100);y1=exp(-x1.*x1);subplot(2,2,1),plot(x1,y1),title('y=exp(-x*x)');.r=-pi:0.01:pi;subplot(2,2,2),polar(r,sin(2*r),'-r'),title

6、('p=sin(2*theda)');.t=linspace(0,9,900);x2=3*t./(1+t.3);y2=3*t.*t./(1+t.3);subplot(2,2,3),plot(x2,y2),title('x=3*t/(1+t3),y=3*t*t(1+t3)');.y3=linspace(-2*pi,2*pi,900);w=sqrt(1-y3.*y3);x3=log(1+w)./y3)-w;x4=log(1-w)./y3)+w;subplot(2,2,4),plot(x3,y3,x4,y3),.title('y=log(1+-sprt(1-x

7、*x)/x)-+sqrt(1-x*x),');Warning: Imaginary parts of complex X and/or Y arguments ignored.Warning: Imaginary parts of complex X and/or Y arguments ignored.图形如图4.1所示疑问：1）为什么将t=linspace(0,9,900)改为t=linspace(-2,9,900)后，运行结果变成了图4.2所示 2）为什么有两个Warning，它们是什么意思图4.1（上上）图4.2（上） 图5.1.1（下）图5.1.2(下下)5作出下列曲面的3维

8、图形，1）；第一次程序:x=-7.5:0.5:7.5;y=x;X,Y=meshgrid(x,y);Z=sin(pi*sqrt(X.2+Y.2);mesh(X,Y,Z);图如图5.1.1第二次程序x=-2:0.01:2;y=x;X,Y=meshgrid(x,y);Z=sin(pi*sqrt(X.2+Y.2);mesh(X,Y,Z);图如图5.1.2所示第三次程序x=-7.5:0.01:7.5;y=x;X,Y=meshgrid(x,y);Z=sin(pi*sqrt(X.2+Y.2);mesh(X,Y,Z);如图5.1.3所示图5.1.3（上）图5.2（下）2）环面： 。程序：u=linspace(

9、0,2*pi,30);v=u;U,V=meshgrid(u,v);X=(1+cos(U).*cos(V);Y=(1+cos(U).*sin(V);Z=sin(U);mesh(X,Y,Z), 图如图5.2所示3）分别作出单位球面在参数为两种不同取值范围的图形,注意坐标轴的单位长度要相等。提示：附加命令rotate3d可实现3维图形旋转。a) ;b) 程序：u1=linspace(0,1.6*pi,30);v1=linspace(0,pi,30);u2=linspace(0,2*pi,30);v2=linspace(0.5*pi,pi,30);.U1,V1=meshgrid(u1,v1);U2,V

10、2=meshgrid(u2,v2);X1=cos(U1).*sin(V1);Y1=sin(U1).*sin(V1);Z1=cos(V1);.X2=cos(U2).*sin(V2);Y2=sin(U2).*sin(V2);Z2=cos(V2);subplot(1,2,1),mesh(X1,Y1,Z1),title('a');subplot(1,2,2),mesh(X2,Y2,Z2),title('b');图如图5.3所示图5.3（上）图5.4（上）图5.5（下）4）z =y2 绕z轴的旋转面图形程序：x=-5:0.01:5;y=x;x,y=meshgrid(x,y

11、);z=x.*x+y.*y;mesh(x,y,z),图如图5.45) y = -,0<x<5 柱面图形程序：x=eps:0.01:5-eps;y=-x;x,y=meshgrid(x,y);z=sqrt(-y);mesh(x,y,z)图如图5.56建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为153=13+53+33。程序：for a=1:9; for b=0:9; for c=0:9; if a*100+b*10+c=a3+b3+c3; d=a*100+b*10+c, end end en

12、dend运行结果：sxh d =153 d = 370 d =371 d =4077编写函数M-文件sq.m：用迭代法求的值。求平方根的迭代公式为迭代的终止条件为前后两次求出的x的差的绝对值小于10-5。程序: function f=sq(a);x=1;for k=1:10000 x=0.5*(x+a/x);endf=x,运行结果：a=4;f=sq(a);f = 28. 求函数的极限、导数或积分：1）当x时；程序：syms x ,f=(x+3.x).(1./x),limit(f,x,1/0,'left') 运行结果：f =(x+3x)(1/x) Warning: Divide by zero. ans =33）求； 程序：syms x,f=(x*x+2*x-1)/(exp(-x)*sin(x)+1),diff(f,x)运行结果：f =(x2+2*x-1)/(exp(-x)*sin(x)+1)ans=(2*x+2)/(exp(-x)*sin(x)+1)-(x2+2*x-1)/(exp(-x)*sin(x)+1)2*(-exp(-x)*sin(x)+exp(-x)*cos(x)5）已知，求；7）;程序：syms x