绘图答案

求方程的正根

解答：

syms x;

f=x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x;

g=solve(f);

for i=1:length(g)

if(‘g(i)’>0)

g(i)

end

end

2. 在同一坐标轴中绘制下列两条曲线并标注两曲线的交叉点。（1）y=2x - 0.5

（2）。

解答：

x1=-pi/2:0.01:pi/2;

y1=2*x1-0.5;

t=0:0.01:pi;

x2=sin(3*t).*cos(t);

y2=sin(3*t).*sin(t);

x=[];

y=[];

for i=1:length(x1)

for j=1:length(x2)

if abs(y1(i)-y2(j))<1e-2&abs(x1(i)-x2(j))<1e-2

x=[x,x1(i)];

y=[y,y1(i)];

end

end

end

plot(x1,y1,x2,y2,x,y,'rd')

3. 分别用plot和fplot函数绘制函数的曲线，分析两曲线的差别。

4. 设，在同一图形窗口采用子图的形式绘制不同图形：条形图、阶梯

图、杆图和对数坐标图。

5. 在[0 2π]范围内绘制以Y轴为对数的二维曲线图: y=|1000sin(4x)|+1

6. 绘制下列极坐标图。

7. 绘制下列三维图线。

解答：

(1) t=-0:0.1:2*pi; plot3(exp(-t/20).*cos(t),exp(-t/20).*sin(t),t);

(2) t=-0:0.1:1; plot3(t,t.^2,t.^3);

8. 绘制下列曲面图。

解答：

(1) x=-5:0.1:5; [X,Y]=meshgrid(x); Z=5*ones(size(X)); surf(X,Y,Z);

(2) x=-3:0.1:3; [X,Y]=meshgrid(x); Z=-5./(1+X.^2+Y.^2); surf(X,Y,Z);

9. 播放一个直径不断变化的球体。

解答：

[x,y,z]=sphere(20);

m=moviein(100);

for i=1:100

xi=x*(1+i/100);

yi=y*(1+i/100);

zi=z*(1+i/100);

surf(xi,yi,zi);

m(:,i)=getframe;

end

movie(m)

10. 用柱坐标系绘制一个圆锥曲面。

11. 绘制曲面。

解答：[x,y]=meshgrid(-6:0.03:6);

z=sin(sqrt(x.^2+y.^2));

i=y

z(i)=NaN;

surf(z)

shading interp

12. 绘制由函数形成的三维立体图，并通过改变视点获得该图形在各

个坐标平面上的平面投影。

解答：[x,y]=meshgrid(-1:0.03:1);

z=x.^2+y.^2;

i=x.^2+y.^2>1;

z(i)=NaN;

subplot(2,2,1)

surf(z);

shading interp

subplot(2,2,2)

surf(z);

shading interp

view(0,0);

title('yoz投影');

subplot(2,2,3)

surf(z)

shading interp

view(90,0);

title('xoz投影');

subplot(2,2,4)

surf(z)

shading interp

view(0,90);

title('xoy投影');

13. 绘制一个单位球体在第一卦限的部分。

14. 在同一图形窗口中绘制函数y1=1-sin2(x) (蓝色圆圈), y2=2x+1（绿色点划线），加分格线，x的范围都是[0，10]，给图形加上标题“y1和

y2”，在x 轴上加注“x轴”，在y 轴上加注“y轴”,在图右侧添加图例，并把“x=5”字符串放置到图形中鼠标所指定的位置上。

解答：x=0:0.1:10;

y1=1-sin(x).^2;

y2=2*x+1;

plot(x,y1,'blueo');

hold on;

plot(x,y2,'g-.');

xlabel('x')

ylabel('y')

grid

legend('y1','y2',4)

gtext('x=5')

15. 在同一图形窗口中绘制椭圆x^2+4y^2=4及其内切圆和外切圆

x^2+y^2=1, x^2+y^2=4。

解答： t=0:pi/99:2*pi;

x=cos(t);

y=sin(t);

plot(x,y)

axis square

x=2*cos(t);

y=2*sin(t);

hold on

plot(x,y)

x=2*cos(t);

y=sin(t);

plot(x,y)