plotxiang.m

DVD租赁的matlab源程序

p1=polyfit(z,x,6); %以z为自变量将x拟合成六次多项式
p2=polyfit(z,y,6); %以z为自变量将y拟合成六次多项式
zz=[0:0.4:99]'; %重新定义步长更小的高度坐标向量zz
xx=polyval(p1,zz); %用拟合多项求对应于zz的坐标xx
yy=polyval(p2,zz); %用拟合多项求对应于zz的坐标yy
dx1=polyder(p1); %对x的多项式微分多项式
dy1=polyder(p2); %对y的多项式微分多项式
dx2=polyval(dx1,zz); %求x'(z)
dy2=polyval(dy1,zz); %求y'(z)
R2=dx2.^2+dy2.^2;
r1=sqrt(R2);
r2=r1.*sqrt(1+R2);
t=linspace(0,2*pi,251); %定义角度
s=size(t);
X=xx*ones(s)+r*((-dy2./r1)*cos(t)+(dx2./r2)*sin(t));
Y=yy*ones(s)+r*((dx2./r1)*cos(t)+(dy2./r2)*sin(t));
Z=zz*ones(s)-r*(r1./r2)*sin(t);
mesh(X,Y,Z); %三维网线作图
surface(X,Y,Z); %三维表面作图
axis equal; %按实际比例值作图
axis off; %不画坐标轴,使图形更逼真
hidden off;%消隐,去掉曲面上的网线,使图形更逼真
view([35,88]); %选择观察角
shading interp; %使作图时片与片之平滑过渡
colormap hot; %色彩处理(不一定恰当,请读者自己
caxis([-70,450]); %重设颜色,使颜色接近红色尝试)
camlight(200,180); %指定光源位置,这些值为尝试值,不一定合适
lighting gouraud; %设置照明方式