3140504022MATLAB平时作业要点

1、MATLAB程序设计基础专业班级电科14013140504022谢言琴MATLA解序设计基础平时作业说明：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、不合理”2）用if语句实现题目要求MATLAB程序代码如下：x=input（please import results:）;% 输入分数if x100disp（Error:wrong data）% 显示数据错误elseif x=90bank=A %输出等级为Aelseif x=80bank=B %输出等级为Belseif x70bank=C %输出等级为Celseif x60bank=D %输出等级为Delse bank=E %输出等级为E end运行结果如下： Untitledplease import results:123Error:wrong data Untitledplease i

3、mport results:96 bank =A Untitledplease import results:85 bank =B Untitledplease import results:72 bank =C Untitledplease import results:63bank =C Untitledplease import results:6 bank =E(3)用switch语句实现题目要求MATLAB程序代码如下：c=input(please import results:);% 输入分数switch floor(c/10) % 用switch语句实现对分数的求商case9ba

4、nk=A;case8bank=B;case 7 bank=C; case 6 bank=D;case num2cell(0:5)%表示取整后是 0至U 5 bank=E;otherwisebank=输入出错;enddisp(bank);运行结果如下：please import results:123输入出错 Untitledplease import results:95A Untitledplease import results:86B Untitledplease import results:75C Untitledplease import results:65D Untitledp

5、lease import results:3 E Untitledplease import results:-1输入出错2.设计程序，完成10,99之间两位整数的加、减、乘、除四则运算。要求:1)随机产生区间10,99之内的两位随机整数,做相应的运算a=x+y=68 的形式。1);2)在产生符合要求的随机数后，输入一个运算符号3)运算结果显示在命令窗口，其显示类似于“MATLAB程序代码如下：y=10+(99-10)*rand % 产生10到99之间的随机数 a=round(rand*90)+10;disp(a)%产生10到99之间的随机数ab=round(rand*90+10);disp(

6、b)%产生10到99之间的随机数bc=input(请输入一个四则运算符号，用单引号括起来 switch ccase +d=a+b;%做加法运算case -d=a-b;%做减法运算case *d=a*b; %做乘法运算case /;d=a/b; %做除法运算 endd运行结果如下： Untitled y =91.29056719请输入一个四则运算符号，用单引号括起来 +d =86 Untitledy =34.78635996请输入一个四则运算符号，用单引号括起来 -d =-37 Untitledy =95.87512497请输入一个四则运算符号，用单引号括起来*d =2328 Untitledy

7、 =95.18795482请输入一个四则运算符号，用单引号括起来/d =0.65853.求下列分段函数的值并绘制函数曲线。x2 x 一6,x ： 0且 x 吏-3y=x2-5x+6, 0Wx=5y=x.A2-x-1;elseif x10100,并求出n!。MATLAB程序代码如下：s=1;for n=1:100if s10100;n=n-1;break;else s=s*n;endendfprintf(n=%dn ,n,s)运行结果如下： Untitledn=8n=40320MATLAB程序代码如下：s=1;n=1;while s Untitledn=8n=403205.已知f1=1, n=1

8、f2=0, n=2f3=1, n=3fn=f n 1 2fn-2+ fn-3, n3求f1f 50中：X,1)最大值和最小值及它们的位置、各数之和(最大值给变量MAX ,其位置给变量最小值给变量 MIN ,其位置给变量y,各数之和给变量 SUM)。2)正数、零、负数的个数(依次赋值给变量N_positive , N_zero , N_negative)。MATLAB程序代码如下：f(1)=1；f(2)=0；f(3)=1;for n=4:50f(n)=f(n-1)-2*f(n-2)+f(n-3)endMAX,X=max(f)% 将最大值赋值给 MAX,位置给xMIN,Y=min(f)% 将最大值

9、赋值给 MIN,位置给ySUM=sum(f)% 和赋值给 SUMN_positive=length(find(f0)% 找到正数的个数和位置，赋值给 N_positiveN_zero=length(find(f=0)%找到零的个数和位置赋值给N_zeroN_negative=length(find(f0)%找到负数的个数和位置赋值给N_negative运行结果如下：MAX =406631X =49MIN = -403795Y =50SUM =76861N_positive =25N_zero =2N_negative =236. Fibonacci数组的元素满足 Fibonacci规则：ak+

10、2=ak+ak+1, (k=1 , 2, 3,)；且 a= a2=1 请设计一段程序，求出该数组中第一个大于10000的元素对应的k值以及该元素值。MATLAB程序代码如下：a1=1；a2=1；i=3；a3=a1+a2；while a3 Untitled i=21 a(i)=10946a(21)=109467 .定义自变量x的取值区间为-5,5,编写一个函数文件，实现绘制函数y=ax3+bx2+cx+d的曲线，其中给定各系数值分别为a=T, b=2, c= - 5, d=3;a=3, b= - 7, c=1,d=4;a=0,b=4,c=0, d= - 5要另外编写调用函数文件的m文件，实现调用

11、。MATLAB程序代码如下：function mydraw（a,b,c,d） % 建立函数 mydrawx=-5:0.1:5;y=a*x.A3+b*x.A2+c*x+d;plot(x,y);mydraw(-1,2,-5,3) mydraw(3,-7,1,4) mydraw（0,4,0,-5）结果截图如下：图2系数取第组的曲图4系数取第组的曲线8 .设计程序，在一个图形窗体中画出函数y=ex和y=logx的特性曲线。其中y=ex中x取值区间为（-2,2）, y=logx中x取值区间为（0.1,5）。要求：1）要有图例、标题（二维图）、坐标轴标签（x轴数据和y轴数据）;2）每条曲线的属性自己设置，

12、但不能全用默认值。MATLAB程序代码如下：x1=-2:0.1:2;y1=exp(x1);x2=0.1:0.1:5;y2=log(x2);plot(x1,y1,r-,x2,y2,k-d);title(二维图);xlabel(x 轴数据);ylabel(y 轴数据);legend(eAx曲线,log(x)曲线);%设置图例结果截图如下：图5函数y=ex和y=logx的特性曲线9 .采用2行2列子图的方式进行二维图形绘制，要求：1)自变量x的取值区间为0,2可2)三个函数分别为：y1=sin(x) , y2=cos(x) , y3=sin(x)+cos(x)3)在第1幅子图绘制y1的曲线；在第2幅

13、子图绘制y2的曲线；在第3幅子图绘制y3的曲线；在第4幅子图中绘制所有3条曲线4)对每一幅子图，都要有图例、坐标轴标签；曲线的属性自己设置，不能全用默认值。MATLAB程序代码如下：x=0:0.1:2.*pi;y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);subplot(221);plot(x,y1,r-p);legend(sin(x);xlabel(横坐标 1);ylabel(纵坐标 1);subplot(222);plot(x,y2,go-);legend(cos(x);xlabel(横坐标 2);ylabel(纵坐标 2);subplot(223);plot(

14、x,y3,k-d);legend(sin+cos);xlabel（横坐标 3）;ylabel（纵坐标 3）;subplot（224）;plot（x,y1,x,y2,x,y3,ro-）;legend（所有）;xlabel（横坐标 4）;ylabel（纵坐标 4）;结果截图如下：i-Ok、匕* Vinr IwimvE Tuahi tugp Windkm Hiip*& l . ? 0 n d例学标TMW4 e 惜学ha24*54圣制 Untitled2 ra =3 rb=3%rank(A)=rank(B)=未知数个数=3,方程组有唯一解R =2-101 00010001 x=inv(A)*b x =2-101 4x1 -6x3 3x4 = 52) 3x1 -x2 -3x3 4x4 = 4x1 5x2 -9x3 -8x4 = 015A=4 0 -6 3;3 -1 -3 4;1 5 -9 -8;b=5;4;0;B=A b;ra=rank(A)rb=rank(B) format ratR=rref(B) 运行结果如下： Untit