MATLAB结课作业(2016.12) 下载本文

计算应用软件. 工程力学专业

结果:

题目5:设f(x)=x5-4x4+3x2-2x+6;

(1)取x=[-2,8]之间函数的值(取100个点),画出曲线,看它有几个零

点。(提示:用polyval函数)

(2)用roots函数求此多项式的根。 (1)代码:p=[1,-4,3,-2,6];

x=-2:0.1:8; y=polyval(p,x); plot(x,y,'LineWidth',5)

5 / 11

计算应用软件. 工程力学专业

结果:

(2)代码: r=roots(p) n=length(r); for i=1:n

if r(i)>=-2 && r(i)<=8 if (isreal(r(i))) r(i) end end end 结果: ans =

3 ans =

1.6956

6 / 11

计算应用软件. 工程力学专业

题目6:在x=[-10,10],y=[-10,10]范围内画出函数z?代码:

x=-10:0.01:10; y=-10:0.01:10;

[X,Y]=meshgrid(x,y);

Z=sin(sqrt(X.^2+Y.^2))./sqrt(X.^2+Y.^2); mesh(X,Y,Z) 结果:

sinx2?y2x?y22的三维图形。

题目7:由100*rand(30)作为全班成绩S(i),用饼图分别表示Si?60、 60?Si?70、

70?Si?80、80?Si?90、90?Si所占比例。

代码: S=100*rand(30);

k1=0;k2=0;k3=0;k4=0;k5=0; for i=1:900 if S(i)<60 k1=k1+1;

7 / 11

计算应用软件. 工程力学专业

elseif S(i)>=60&S(i)<70 k2=k2+1; elseif S(i)>=70&S(i)<80 k3=k3+1; elseif S(i)>=80&S(i)<90 k4=k4+1; else

k5=k5+1; end end

x=[k1 k2 k3 k4 k5]; y=[1 1 1 1 1];

pie(x,y,{'Si<60','60<=Si<70','70<=Si<80','80<=Si<90','90<=Si'}); 结果:

8 / 11