Koch

最新推荐文章于 2021-03-16 18:53:48 发布

_NULL_Coder

最新推荐文章于 2021-03-16 18:53:48 发布

阅读量366

点赞数

文章标签：递归 Koch

本文链接：https://blog.csdn.net/weixin_40446473/article/details/79542436

版权

科赫曲线点的求法。

思路是运用递归。

易知n阶曲线有四段n-1阶曲线构成。

如b是由四段a构成的，c是由四段b构成的。//长度每次缩短三倍

可知，循环体如下：

pl[j].x=X+abs(d*cos(k*Pi));//d为线段的长度，k控制方向
pl[j].y=Y+d*sin(k*Pi);

X=pl[j].x,Y=pl[j].y;

代码如下：

# include<iostream>
# include<cmath>
# include<cstring>
using namespace std;
long double Pi=acos(-1);
int n,j=1;
double X=0.0,Y=0.0;
double d=100.0,k;
struct KOCH
{
	double x,y;
}pl[2000];


int koch(int i)
{
	if(i==0)
	{
		pl[j].x=X+abs(d*cos(k*Pi));
		pl[j].y=Y+d*sin(k*Pi);
		X=pl[j].x,Y=pl[j].y;
		j++;
		return 1;
	}
	d/=3;
	for(int num=0;num<3;num++)
	{
		switch(num)
		{
			case 0:
				k=0.0;
				break;
			case 1:
			    k=1.0/3;
				break;
			case 2:
			    k=1.0+1.0/3;
				break;
			case 3:
			    k=1.0/3;				
		}
		koch(i-1);
	}	
}


int main()
{
	memset(pl,0,sizeof(pl));
	scanf("%d",&n);
	koch(n);
	printf("0.00000000 0.00000000\n");
	for(int i=1;;i++)
	{
		if(pl[i].x==0&&pl[i].y==0)
		   break;
		else
		  printf("%.8lf %.8lf\n",pl[i].x,pl[i].y);   
	}
	printf("100.00000000 0.00000000\n");
}