Matlab练习及答案

Matlab 练习一
1、一个三位整数各位数字的立方和等于该数本身，则称该数为水仙花数。

试编写程序找到并输出全部水仙花数。

2、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6 是完数。

求[1 500]之间的全部完数。

3、写一个函数pifun 计算下列级数：
F(n)=4*(1-1/3+1/5-1/7+...)
并使用tic 和toc 指令来测量pifun(1000000)的计算时间。

4、输入一个百分制成绩，要求输出成绩等级为A~E，即90~100 为A，80~89 为B，70~79 为C，60~69 为D，60 以下为E。

5、输入x,y 的值，并将它们的值互换后输出。

6、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price 来表示)：
price<200 没有折扣；200≤price<500 3%折扣；500≤price<1000 5%折扣；
1000≤price<2500 8%折扣；2500≤price<5000 10%折扣；5000≤price 14%折扣
输入所售商品的价格，求其实际销售价格。

7、由指令A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5 的元素的位置。

8、写一程序求1！+2！+ (10)
9、编写一个函数，找出矩阵A 中最大元素和最小元素所在的行列号。

10、从键盘输入的10 个整数中，找出第一个能被7 整除的数。

若找到，打印此数后退出循环；若未找到，打印“not exist”。

function shuixian(n) %n为用户输入的一个整数n=3;
W=zeros(1,n);
ss=0;
for m=10^(n-1):10^n-1
for i=1:n
W(i)=fix(rem(m,10^(n+1-i))/10^(n-i));
end
W=W.^n;
for j=1:n
ss=ss+W(j);
end
if m==ss
disp(m);
end
ss=0;
W=zeros(1,n);
end
2.
for m=1:500
sum=0;
for n=1:m/2
if rem(m,n)==0
sum=sum+n;
end
end
if sum==m
m
end
end
3.
n=input('请输入n值：');
tic;
f=0;
for i=1:n
f=f+4*((-1)^(i+1))/(2*i-1) ;
end
f
toc
4.1
m=input('请输入一个百分制成绩：');
m=ceil(m/10);
switch m
case {10,9}
disp('A');
case 8
disp('B');
case 7
disp('C');
case 6
disp('D');
case {5,4,3,2,1,0}
disp('E');
otherwise
disp('输入有误！')
end
4.2
a=input('请输入一个百分制成绩：');
if a>=90 & a<=100
disp('A');
else if a>=80 & a<90
disp('B'); else if a>=70 & a<80 disp('C'); else if a>=60 & a<70 disp('D'); else if a>=0 & a<60 disp('E'); else
disp('输入数据有误，请核对......');
end
end
end
end
end
5．
x = input('x=');
y = input('y=');
disp ([x,y]);
t=x;
x=y;
y=t;
disp ([x,y]);
x=input('price:')
if x<200
f=x;
else if x<500
f=x*0.97;
else if x<1000
f=x*0.95;
else if x<2500
f=x*0.92;
else if x<5000 f=x*0.90;
else
f=x*0.86;
end
end
end
end
end
f
7.
A=rand(3,5)
for i=1:3
for j=1:5
if A(i,j)>0.5
Ai=i;
Aj=j;
Aij=[i,j]
end
end
end
sum=0;
for m=1:10
i=1;
for n=1:m
i=i*n;
end
sum=sum+i;
end
sum
10.
a=zeros(1,10);
for i=1:10
a(i)=input('please enter number'); end
for i=1:10
if (rem(a(i),7)==0)
a(i)
break
end
end
if (i==0&&rem(a(i),7)~=0)
disp('not exist')
end。