1、 1设有分块矩阵 ，其中 E,R,O,S 分别为单位阵、随机阵、零232SOREA阵和对角阵，试通过数值计算验证 。2S0REE=eye(3);R=rand(3,2);O=zeros(2,3);S=diag(rand(2,1);A=E,RO,S;a=A2-E,R+R*SO,S22.某零售店有 9 种商品的单件进价（元）、售价（元）及一周的销量如表 1.1，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该 10 种商品的总收入和总利润。表 1.1货号 1 2 3 4 5 6 7 8 9 单件进价 7.15 8.25 3.20 10.30 6.68 12.
2、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,753,580,395,2104,1538,810,694;d=(b-a).*c;max(d);min(d);x1,x2=sort(d
3、,1)3.建立一个命令 M-文件：求所有的“水仙花数” ，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153 是一个水仙花数，因为153=13+53+33。function NUM=daffodil()count=0; for i=1:9for j=0:9for k=0:9if i3+j3+k3=100*i+10*j+kcount=count+1;NUM(count)=100*i+10*j+k;elseendendendend4编写函数 M-文件 sq.m：用迭代法求 的值。求平方根的迭代公式为ax)a(21nnxxfunction f=sq(a)format lon
4、gx=a;f=(x+a/x)/2;while abs(f-x)0.000001x=f;f=(x+a/x)/2;endend5.在同一个坐标下作出 y1=ex,y2=1+x,y3=1+x+(1/2)x2,y4= 1+x+(1/2)x2+(1/6)x3 这四条曲线的图形，要求在图上加各种标注，观察到什么现象？发现有什么规律？同时用 subplot 分别在不同的坐标系下作出这四条曲线，为每幅图形加上标题。x=0:0.01:3;y=0:0.01:25;y1=exp(x);plot(x,y1,-m)hold ony2=1+x;plot(x,y2,.c);hold ony3=1+x+x.2/2；plot(
5、x,y3,*k) ；hold ony4=1+x+x.2/2+x.3/6;plot(x,y4,r);hold offlegend(y=e 的 x 次方,y=x+1,y=1+x+x.2/2,y=1+x+x.2/2+x.3/6);title(不知道该取什么名字 )xlabel(x=0:0.01:3)ylabel(y=0:0.01:25);figure(2)x=0:0.01:3;y=0:0.01:25;y1=exp(x);subplot(2,2,1)plot(x,y1,-m);title(y=exp(x)y2=1+x;subplot(2,2,2)plot(x,y2,.c);title(y=1+x)y3
6、=1+x+x.2/2subplot(2,2,3)plot(x,y3,*k)title(y=1+x+x.2/2)y4=1+x+x.2/2+x.3/6;subplot(2,2,4)plot(x,y4,r);title(y=1+x+x.2/2+x.3/6)figure 10 1 2 30510152025 y=exp(x)0 1 2 31234 y=1+x0 1 2 30246810 y=1+x+x.2/20 1 2 3051015 y=1+x+x.2/2+x.3/60 0.5 1 1.5 2 2.5 30510152025x=0:0.01:3y=0:0.01:25不不不不不不不不不y=e不x不不y
7、=x+1y=1+x+x.2/2y=1+x+x.2/2+x.3/6figure 26.作出下列曲面的 3 维图形， ；)sin(2yxzx=-2*pi:pi/10:2*pi;y=-2*pi:pi/10:2*pi;X,Y=meshgrid(x,y);r=sqrt(X.2+Y.2)*piZ=sin(r)mesh(X,Y,Z)shading flat7. 作出函数 y=x4-4x3+3x+5 （x0,6）的图形，找出其在0,6之间的最小值点。x=0:0.01:6;y=x.4-4*x.3+3*x+5;plot(x,y)xmin,ymin=fminbnd(x)(x.4-4*x.3+3*x+5),2,5);text(xmin,ymin,*);
