matlab练习一二答案

1、练习一答案1、一个三位整数各位数字的立方和等于该数本身，则称该数为水仙花数。试编写程序找到,并输出全部水仙花数。for n=100:999a=floor(n/100);b=floor(rem(n,100)/10);c=rem(n,10);if n=a3+b3+c3 s=nendend2、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6 是完数。求1 500之间的全部完数。for i=2:500 s=1; for j=2:(i-1) if (rem(i,j)=0) s=s+j; end end if s=i s end end 3、写一个函数pifun 计算下列级数： F

2、(n)=4*(1-1/3+1/5-1/7+.) 并使用tic 和toc 指令来测量pifun(1000000)的计算时间。n=input(请输入n值：); tic; f=0; for i=1:n f=f+4*(-1)(i+1)/(2*i-1) ; end ftoc4、输入一个百分制成绩，要求输出成绩等级为AE，即90100 为A，8089 为B，7079为C，6069 为D，60 以下为E。m=input(请输入一个百分制成绩：);m=floor(m/10);switch m case 10,9 disp(A); case 8 disp(B); case 7 disp(C); case 6 d

3、isp(D); disp(E); otherwisedisp(E);end5、输入x,y 的值，并将它们的值互换后输出。x = input(x=);y = input(y=);disp (x,y);t=x;x=y;y=t;disp (x,y);6、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price 来表示)：price200 没有折扣；200price500 3%折扣；500price1000 5%折扣；1000price2500 8%折扣；2500price5000 10%折扣；5000price 14%折扣输入所售商品的价格，求其实际销售价格。x=input(price:

4、)if x200 f=x;else if x500 f=x*0.97; else if x1000 f=x*0.95; else if x2500 f=x*0.92; else if x0.5 Ai=i; Aj=j; Aij=i,j end endend8、写一程序求1！+2！+10！sum=0;for m=1:10 i=1; for n=1:m i=i*n; end sum=sum+i;endsum9、编写一个函数，找出矩阵A 中最大元素和最小元素所在的行列号。A=rand(3,5);disp(A);Max=max(max(A);Min=min(min(A);disp(Max);r1,c1=

5、find(A=Max)disp(Min);r2,c2=find(A=Min)10、从键盘输入的10 个整数中，找出第一个能被7 整除的数。若找到，打印此数后退出循环；若未找到，打印“not exist”。a=zeros(1,10);for i=1:10 a(i)=input(please enter number);endfor i=1:10 if (rem(a(i),7)=0) a(i) break endend if (a(i)=0&rem(a(i),7)=0) disp(not exist)end练习二答案1 绘制下列曲线。x=-2:0.1:2;y=exp(-x.2/2)/(2*pi);

6、plot(x,y)1.2t=0:0.1:2*pi;plot(8*cos(t),4*sqrt(2)*sin(t)2分别用plot和fplot函数绘制函数的曲线，分析两曲线的差别(1)x=-2:0.01:2;y=sin(1./x);plot(x,y)(2)fplot(sin(1/x),-2,2)3设y=1/1+e-t,-pi=t=pi,在同一窗口中采用子图的形式绘制不同图形：条形图，阶梯图，杆图和对数坐标图t=-pi:0.1:pi;y=1./(1+exp(-t);plot(t,y);subplot(2,2,1);stairs(t,y);subplot(2,2,2);bar(t,y);subplot

7、(2,2,3);stem(t,y);subplot(2,2,4);loglog(t,y);4绘制三维图X=e-t/20cos(t)Y=e-t/20sin(t) 0t2piZ=t t=0:0.1:2.*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t); x,y=meshgrid(x,y);z=(log(x.2+y.2).*(-10);mesh(x,y,z);5 在同一图形窗口中绘制函数y1=1-sin2(x) (蓝色圆圈), y2=2x+1（绿色点划线），加分格线，x的范围都是0，10，给图形加上标题“y1和y2”，在x 轴上加注“x轴”， 在y 轴上加注“y轴”,在图右侧添加图例，并把“x=5”字符串放置到图形中鼠标所指定的位置上。x=0:0.1:10;y1=1-(sin（x)）.2;y2=2*x+1;plot(x,y1,bo);hold on;plot(x,y2,g-.); xlabel(x)ylabel(y)gridlegend(y1,y2,1)gtext(x=5) 6x=0:0.04:2*pi;y=exp(x./2).*sin(x);y1=exp(x./2).*si