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

结果：

题目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