Matlab上机题及答案.doc

1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。
for m=100:999
m1=fix(m/100); %求m的百位数字
m2=rem(fix(m/10),10); %求m的十位数字
m3=rem(m,10); %求m的个位数字
if m==m1*m1*m1+m2*m2*m2+m3*m3*m3
disp(m)
end
end

2. 从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。
sum=0;
n=0;
val=input('Enter a number (end in 0):');
while (val~=0)
sum=sum+val;
n=n+1;
val=input('Enter a number (end in 0):');
end
if (n > 0)
sum
mean=sum/n
end

3. 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求[1,500]之间的全部完数。
for m=1:500
s=0;
for k=1:m/2
if rem(m,k)==0
s=s+k;
end
end
if m==s
disp(m);
end
end

4. 从键盘上输入数字星期，在屏幕上显示对应英文星期的单词。
function week
n=input('input the number:');
if isempty(n)
errror('please input !!')
end
if n>7|n<1
error('n between 1 and 7')
end

switch n
case 1
disp('Monday')
case 2
disp('Tuesday')
case 3
disp('Wednesday')
case 4
disp('Thursday')
case 5
disp('Friday')
case 6
disp('Saturday')
case 7
disp('Sunday')
end

5. 某公司销售电脑打印机的价格方案如下：
（）如果顾客只买一台打印机，则一台的基本价格为$150。
（）如果顾客购买两台以上打印机，则第二台价格为$120。
（）第三台以后，每台$110。
写一段程序分别计算出购买1--10台打印机所需的钱数。打印机台数可以在程序开始处指定，或通过input命令读入。运行程序，计算出购买10台打印机的总价格。
写出程序，生成分别购买1--10台打印机所需价格的图表（使用fprintf命令输出图表，不允许手算）。
x=input('请输入购买的打印机台数：');
for m=1:x
if m<=1
y(m)=150*m;
elseif m<=2
y(m)=150+120*(m-1);
else y(m)=150+120+110*(m-2);
y(1,m)=y(m);
end
end
y(x)
plot(1:m,y,'r*--')

6. 有一组学生的考试成绩（见表），根据规定，成绩在100分时为满分，成绩在90~99之间时为优秀，成绩在80~89分之间时为良好，成绩在60~79分之间为及格，成绩在60分以下时为不及格，编制一个根据成绩划分等级的程序。
学生姓名
王
张
刘
李
陈
杨
于
黄
郭
赵
成 绩
72
83
56
94
100
88
96
68
54
65

作业4 绘图练习
1、有一组测量数据满足，t的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线。图中添加标题（title）和图例框（图注legend）。

2.20 在5.1题结果图中添加标题，并用箭头线标识出各曲线a的取值。

2.21 在5.1题结果图中添加标题和图例框。

2、绘制三维曲线：x=sin(t); y=cos(t); z=cos(2*t);并设计标题和图注。
>>t=(0:0.02:2)*pi;
>>x=sin(t);y=cos(t);z=cos(2*t);
>>plot3(x,y,z,'b-',x,y,z,'bd')
>>view([-82,58]),box on
>>xlabel('x'),ylabel('y'),zlabel('z')
>>legend('链','宝石','Location','best')

3、在同一窗口中画出一个正六变形、半径为1的圆以及球。并标明各自的标题

4、，当x和y的取值范围均为-2到2时，用建立子窗口的方法在同一个图形窗口中绘制出三维线图（plot3）、网线图(mesh)、表面图(surf)和带渲染效果(shading interp)的表面图,标明图形的标题，使图像完美。

>> [x,y]=meshgrid([-2:.2:2]);
>> z=x.*exp(-x.^2-y.^2);
>> mesh(x,y,z)
>> subplot(2,2,1),plot3(x,y,z)
>> title('plot3(x,y,z)')
>> subplot(2,2,2),mesh(x,y,z)
>> title('mesh(x,y,z)')
>> subplot(2,2,3),surf(x,y,z)
>> title('surf(x,y,z)')
>> subplot(2,2,4),surf(x,y,z),shading interp
>> title('surf(x,y,z),shading interp')