1、1MATLAB 程序设计基础平时作业说明：1-17 题为必做题，须分别提交纸质文档及与纸质文档对应的电子文档.doc 版。1. 对于每输入的一个百分制成绩，要求输出对应的成绩等级 A、B、C、D、E。其中 A的分值在区间90,100)内；B 的分值在区间80,90)内；C 的分值在区间70,80)内；D 的分值在区间60,70)内；E 的分值在区间0,60)内；要求：1）要有对输入成绩是否合理的判断：输入成绩如果含有小数，只能是“x.5”的形式，其实 x为合理的整数。成绩必须在百分制范围内。如果输入的成绩不合理，应向命令窗口输出报错信息“输入的成绩不合理”2）用 if语句实现题目要求x=inp
2、ut(输入一个数值 :); %输入if mod(x,0.5)=0 %输入if mod(x,0.5)=0 %对 s 赋初值for n=1:100 %默认步长为 1s=s*n; %s=n!if(s10100) break; %如果 s10100,结束运行endenddisp(n=,num2str(n) %输出 ndisp(s=,num2str(s) %输出 s运行结果：n=8s=40320while 循环语句如下：s=1; %对 s 赋初值n=0; %循环次数赋初值while(s3求 f1 f50中：1） 最大值和最小值及它们的位置、各数之和（最大值给变量 MAX，其位置给变量 x，最小值给变量
3、MIN，其位置给变量 y，各数之和给变量 SUM） 。2） 正数、零、负数的个数（依次赋值给变量 N_positive，N_zero，N_negative） 。f(1)=1; f(2)=0;f(3)=1; %输入n=3; %赋值while (n0) %numel 表示返还数组中的元素个数，N_zero=numel(find(f=0) %返还 f 中等于 0 的位置N_negative=numel(find(f0)运行结果：MAX =406631x = 49MIN =-403795y =50SUM =76861N_positive =25N_zero =2N_negative =236. Fib
4、onacci数组的元素满足 Fibonacci规则：ak+2=ak+ak+1,（k=1，2，3，） ；且 a1= a2=1。请设计一段程序，求出该数组中第一个大于 10000的元素对应的 k值以及该元素值。a(1)=1; a(2)=1;k=1;while a(k)=10000 a(k+2)=a(k)+a(k+1); k=k+1;enddisp(k=,num2str(k)disp(a(k)=,num2str(a(k)运行结果：6k=21a(k)=109467. 定义自变量 x的取值区间为5,5，编写一个函数文件，实现绘制函数y=ax3+bx2+cx+d的曲线，其中给定各系数值分别为a=1，b=2
5、， c=5，d=3;a=3， b=7，c=1， d=4;a=0， b=4， c=0， d=5要另外编写调用该函数文件的 m文件，实现调用。function my1(a,b,c,d) %定义主函数，函数文件名为 my1.mx=-5:1:5; %输入 xy=a*x.3+b*x.2+c*x+d;plot(x,y); %绘制曲线xlabel(x 轴); %标志 x 轴ylabel(y 轴); %标志 y 轴title(函数图像 ); %设置图名称a=input(please input a=); %子函数b=input(please input b=);c=input(please input c=)
6、;d=input(please input d=); %利用输出参数 a,b,c,d 向 my1 传递数据my1(a,b,c,d); %调用的是啊 ,a,b,c,d 四个参数图 7（1） 图 7（2）图 7（3） 8. 设计程序，在一个图形窗体中画出函数 y=ex和 y=logx的特性曲线。其中 y=ex中 x取值7区间为(2,2)，y=logx 中 x取值区间为(0.1,5)。要求：1）要有图例、标题（二维图） 、坐标轴标签（x 轴数据和y 轴数据） ；2）每条曲线的属性自己设置，但不能全用默认值。x1=-2:0.1:2;y1=exp(x1);x2=0.1:0.1:5;y2=log(x2);
7、plot(x1,y1,r*-,x2,y2,ko-)legend(函数 y=exp(x)的特性曲线 ,函数 y=logx 的特性曲线)title(二维图 )xlabel(x 轴数据)ylabel(y 轴数据)图 89. 采用 2行 2列子图的方式进行二维图形绘制，要求：1）自变量 x的取值区间为0,22）三个函数分别为：y1=sin(x)，y2=cos(x)，y3=sin(x)+cos(x)3）在第 1幅子图绘制 y1的曲线；在第 2幅子图绘制 y2的曲线；在第 3幅子图绘制y3的曲线；在第 4幅子图中绘制所有 3条曲线4）对每一幅子图，都要有图例、坐标轴标签；曲线的属性自己设置，不能全用默认值
8、。x=0:0.1:2*pi;y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);subplot(2,2,1);plot(x,y1,ro-);legend(sin 曲线 );title(y=sin(x);xlabel(横坐标 1);ylabel(纵坐标 1);8subplot(222);plot(x,y2,g*-);legend(cos 曲线);title(y=cos(x);xlabel(横坐标 2);ylabel(纵坐标 2);subplot(223);plot(x,y3,b-p);legend(sin(x)+cos(x)曲线);title(y=sin(x)+cos(x)
9、;xlabel(横坐标 3);ylabel(纵坐标 3);subplot(224);plot(x,y1,ro-,x,y2,g*-,x,y3,b-p); %同时绘制三个函数的图像legend(sin 曲线 ,cos 曲线,sin(x)+cos(x)曲线);title(三个函数的曲线);xlabel(横坐标);ylabel(纵坐标);图（9）10. 采用 1行 2列子图的方式进行二维图形绘制，要求： 1）在第 1幅子图绘制三维曲线：x0,20y=sin(x)z=cos(x)2）在第 2幅子图绘制三维曲面：x20,20y10,10 4y9xz23） 对每一幅子图，都要有标题、坐标轴标签；曲线的属性自
10、己设置，不能全用默认9值。subplot(121); %在第一个子图中绘图x1=0:0.1:20;y1=sin(x1);z1=cos(x1);plot(x1,y1,ro-,x1,z1,b*-);legend(sin 曲线 ,cos 曲线); %设置图例title(图 1);xlabel(x1);ylabel(y1);subplot(122); %在第二个子图中绘图x2=-20:1:20;y2=0.5*x2;X,Y=meshgrid(x2,y2); %生成 X 和 Y 网络坐标矩阵z=(X.2)/9-(Y.2)/4; %网格化后直接输入 zmesh(X,Y,z); %曲面绘图legend(双曲面
11、);title(图 2);xlabel(x2);ylabel(y2);zlabel(z)grid on;图（10）11. 已知函数 ，其中 x的取值区间为2,2，试编程绘制该二维曲0yx)1(322线，要求：要有标题、坐标轴标签，曲线颜色用红色。p=ezplot(x.2+y.2-1).3-x.2.*y.3,-2 2); %在区间-2 2绘制(x2+y2-1)3-x2y2=0 的图像set(p,color,r); %指定曲线颜色为红色title(曲线图 );xlabel(x 轴);10ylabel(y 轴);图（11）12. 已知原始数据 ， ，插值点为 ，试pi*4:10x)5/xep(*).
12、siny pi*4:3.0xi用 4种不同的插值方法进行插值，并在一幅图的 4个子图中分别绘出 4种插值曲线，并各自注明所使用的插值方法。x=0:1:4*pi;y=sin(x).*exp(-x/5);xi=0:0.3:4*pi; %设置新的插点subplot(221)yi=interp1(x,y,xi,nearset); %使用最近点插值plot(xi,yi)title(最近点插值)subplot(222)yi=interp1(x,y,xi); %默认为线性插值plot(xi,yi)title(线性插值 )subplot(223)yi=interp1(x,y,xi,cubic); %使用三次立方插值plot(xi,yi)title(三次立方插值)subplot(224)yi=interp1(x,y,xi,spline); %使用三次样条插值plot(xi,yi)title(三次样条插值)
