五邑大学《MATLAB程序的设计》选修课大作业.docVIP

.. .. 五 邑 大 学 《MATLAB程序设计》选修课大作业 学生姓名： 学生学号： 学院年级专业： 学院 级 专业 指导老师： 刘战 （副教授） 成绩： 2010年12月 1、求下列联立方程的解 求系数矩阵的秩； 求出方程组的解。 解： 在MATLAB中输入以下指令： a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10] a = 3 4 -7 -12 5 -7 4 2 1 0 8 -5 -6 5 -2 10 rank(a) ans = 4 b=[4;-3;9;-8] b = 4 -3 9 -8 k=a\b k = -1.4841 -0.6816 0.5337 -1.2429 故系数矩阵的秩为4 方程组的解为x=-1.4841 y=-0.6816 z=0.5337 w=-1.2429 2、设把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线； 解： 在MATLAB中输入以下指令： x=linspace(0,2*pi,101); y=cos(0.5+3.*sin(x)./(1+x.^2)); plot(x,y,g-) //这里可以把‘g-’改成其它的，具体请对表基本线型和颜色(见附录) title(函数曲线) xlabel(自变量 x) ylabel(因变量y) text(2,0.5,y=cos[0.5+3sin(x)/(1+x^2)]) grid //这里可以不用grid 也就是没有栅格 用就有 下图就没用 输入结果为： 3、设 （1）取x=[-2,8]之间函数的值（取100个点），画出曲线，看它有几个零点。 （提示：用polyval 函数） 解：经计算可知，当x=-0.2时，f(x)=-68;当x=8时，f(x)=16566 但为了便于观察零点，f(x)的范围取小点。 在MATLAB中输入以下程序： p=[1 -4 0 3 -2 6]; x=linspace(-2,8,100); y=polyval(p,x); plot(x,y,g-); //这里同上 axis([-2,8, -90,800]);//这里-90，800可以适当修改 看出有交点就行了 grid roots(p) //这root(p)可以不要，因为题目不要求求出解 输入结果如图： ans = 3.7999 -1.2607 1.3479 0.0564 + 0.9623i 0.0564 - 0.9623i 经观察和由运算结果可知，有3个零点，分别为x1=-1.2607,x2=1.3479,x3=3.7999。 4、在[-10，10；-10，10]范围内画出函数的三维图形。 解： 在MATLAB中输入以下程序： [X,Y]=meshgrid(-10 : 0.5 :10); a=sqrt(X.^2+Y.^2); Z=sin(a)./a; surf(X,Y,Z); //这里可以把 surf改成mesh 输出结果如图： 5.将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。（10分） ①（曲线图）； ②（曲面图）。 解：在MATLAB中输入以下程序： subplot(2,2,1) subplot(2,2,2) ezplot((cos(x))^(1/2),[-pi./2,pi./2]) xlabel(自变量x) ylabel(因变量y) subplot(2,2,3) x=-2:1:2; y=-4:2:4; ezsurfc(x^2/2^2+y^2/4^2) text(-13,-12,f(x,y)=x^2/2^2+y^2/4^2) subplot(2,2,4) 输出结果如图： 6有一组学生的考试成绩（见表），根据规定，成绩在100分时为满分，成绩在90~99之间时为优秀，成绩在80~89分之间时为良好，成绩在60~79分之间为及格，成绩在60分以下时为不及格，编制一个根据成绩划分等