matlab练习一、二答案

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

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

3、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 没有折扣；200 w price500 3% 折扣；500 price1000 5% 折扣； 1000 w price2500 8% 折扣；2500 w price5000 10% 折扣；5000 w price 14% 折扣 输入所售商品的价格，求其实际销售

4、价格。x=input( price: )if x200f=x;else if x500f=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,jendendend8、写一程序求 1！+2！ +, +10！ sum=0;for m=1:10i=1;for n=1:mi=i*n;endsum=sum+i;endsum9、编写一个函数，找出矩阵 A 中最大元素和最小元素所在的行列号。A=rand(3,5);disp(A);Max=max(max(A);Min=min(min(A);dis

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

6、=64x2e_22 二_Lx =tsin t(4)y =t costx=-2:0.1:2;y=exp(-x.A2/2)/(2*pi);Plot(x,y)1.2t=0:0.1:2*pi;plot(8*cos(t),4*sqrt(2)*sin(t)分别用plot和fplot函数绘制函数1y =si n 1的曲线，分析两曲线的差别x(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./

7、(1+exp(-t); plot(t,y); subplot(2,2,1); stairs(t,y); subplot(2,2,2); bar(t,y); subplot(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=tt=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.A2+y.A2).*( -10);mesh(x,y,z);5 在同一图形窗口中

8、绘制函数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).A2;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)x=0:0.04:2*pi; y=exp(x ./2) .*s in (x); y1=exp(x./2).*