计算机图形学程序88243.docVIP

实验一 function midpointline(x0,y0,x1,y1) x0=1; y0=2; x1=7; y1=6; x=x0; y=y0; dx=x1-x0; dy=y1-y0; grid on; hold on; plot(x,y,o); d0=dx-2*dy; d1=-2*dy; d2=2*(dx-dy); while(xx1) if(d0 = 0) x=x+1; d0=d0+d1; elseif(d0 0) x=x+1; y=y+1; d0=d0+d2; end hold on; plot(x,y,o); end end %Bresenham画圆法 function BresenhamCircle(x0,y0,R) figure; R=4; x0=8; y0=8; x=0; y=R; p=3-2*R; grid on;hold on; for x=0:y plot(x+x0,y+y0,.); hold on; plot(y+x0,x+y0,.);%关于y=x的对称点 hold on; plot(x+x0,-y+y0,.);%关于x轴的对称点 hold on; plot(y+x0,-x+y0,.); hold on; plot(-x+x0,y+y0,.);%关于y轴(未平移之前)的对称点 hold on; plot(-y+x0,x+y0,.); hold on; plot(-x+x0,-y+y0,.); hold on; plot(-y+x0,-x+y0,.); hold on; if p=0 p=p+4*(x-y)+10; y=y-1; else p=p+4*x+6; end end set(gca,XTick,[x0-R:1:x0+R]); set(gca,YTick,[y0-R:1:y0+R]); axis equal; 实验二 x=meshgrid(-10:0.5:10); y=x; z=sin((x.^2+y.^2).^(1/2)).*(x.^2+y.^2).^(1/2); h=surf(x,y,z); axis([-15,15,-15,15,-15,15]); try while 1 rotate(h,[0 0 1],1); pause(0.05); end end clc; clear; M=moviein(16); for j=1:16 sphere(j); axis equal; axis off; M(j)=getframe; end movie(M,2); x=meshgrid(-10:0.5:10); y=x; z=sin((x.^2+y.^2).^(1/2)).*(x.^2+y.^2).^(1/2); surf(x,y,z); h=gca; axis vis3d; for i=0:inf view(h,i,30); drawnow; end clc; clear; y=rand(1,30); axis off; h1=axes(position,[0 0 1 1]); plot(y); axis off; h=axes(position,[0 0 0.1 0.1]); while 1 for i=1:30 set(h,position,[i/30,y(i),0.1,0.1]); sphere(30); axis square off; pause(0.01); end end p=peaks(50); h=axes(position,[0 0 1 1],visible,off); for i=1:20 p1=p/i; surf(p1); set(h,zlim,[0 10]); axis off; pause(0.1); end clc;clear; M=moviein(16); for j=1:16 rectangle(position,[4,5,15,10],curvature,j/16); M(j)=getframe; axis equal; end movie(M,2); A=imread(text.jpg)