matlab练习一、二答案

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

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

3、);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 price500 3% 折扣； 500 price1000 5% 折扣； 1000 price2500 8% 折扣； 2500 price5000 10% 折扣； 5000price 14% 折扣输入所售商品的价格，求其实际销售价格。x=input(pri

4、ce:)ifx200f=x;elseifx500f=x*0.97;elseifx1000f=x*0.95;elseifx2500f=x*0.92;elseifx0.5Ai=i;Aj=j;Aij=i,jendendend8、 写一程序求 1！ +2！ + +10 ！sum=0;form=1:10i=1;forn=1:mi=i*n;endsum=sum+i;endsum9、 编写一个函数，找出矩阵A 中最大元素和最小元素所在的行列号。A=rand(3,5);disp(A);Max=max(max(A);Min=min(min(A);disp(Max);r1,c1=find(A=Max)disp(M

5、in);r2,c2=find(A=Min)10、从键盘输入的 10 个整数中，找出第一个能被 7 整除的数。若找到，打印此数后退出循环；若未找到，打印“ not exist ”。a=zeros(1,10);fori=1:10a(i)=input(please enter number);endfori=1:10if(rem(a(i),7)=0)a(i)breakendendif(a(i)=0&rem(a(i),7)=0)disp(not exist)end练习二答案1 绘制下列曲线。x3x2(1)y=x-(2) y1 e 23!2(3)x22y 264 (4)xt sin tyt costx=

6、-2:0.1:2;y=exp(-x.2/2)/(2*pi);plot(x,y)1.2t=0:0.1:2*pi;plot(8*cos(t),4*sqrt(2)*sin(t)2 分别用 plot 和 fplot 函数绘制函数 ysin1的曲线，分析两曲线的差别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./(1+exp(-t);plot(t,y);subplot(2,2,1

7、);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.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