日志
二次B样条法绘制曲线函数
|
void byangtiao(int *x,int *y,int n,unsigned int k)
{
unsigned int i,j;
float t,t1,t2,a,b,c,tx,ty;
*x=*(x+1);*y=*(y+1);
*(x+n+1)=*(x+n);*(y+n+1)=*(y+n);
t=1.0/k;
setcolor(13);
moveto((*x+(*(x+1)))/2.0,(*y+(*(y+1)))/2.0);
for(i=0;i<n;i++)
{
for(j=1;j<k;j++)
{
t1=j*t;
t2=t1*t1;
a=(t2-2*t1+1)/2.0;
b=t1-t2+1/2.0;
c=t2/2.0;
tx=a*(*(x+i))+b*(*(x+i+1))+c*(*(x+i+2));
ty=a*(*(y+i))+b*(*(y+i+1))+c*(*(y+i+2));
lineto(tx,ty);
}
}
}
{
unsigned int i,j;
float t,t1,t2,a,b,c,tx,ty;
*x=*(x+1);*y=*(y+1);
*(x+n+1)=*(x+n);*(y+n+1)=*(y+n);
t=1.0/k;
setcolor(13);
moveto((*x+(*(x+1)))/2.0,(*y+(*(y+1)))/2.0);
for(i=0;i<n;i++)
{
for(j=1;j<k;j++)
{
t1=j*t;
t2=t1*t1;
a=(t2-2*t1+1)/2.0;
b=t1-t2+1/2.0;
c=t2/2.0;
tx=a*(*(x+i))+b*(*(x+i+1))+c*(*(x+i+2));
ty=a*(*(y+i))+b*(*(y+i+1))+c*(*(y+i+2));
lineto(tx,ty);
}
}
}
