混沌分形之迭代函数系统(IFS).docx

混沌分形之迭代函数系统( IFS) - 电脑资料
IFS 是分形的重要分支，
算法：
.设定一个起始点(x0, y0)及总的迭代步数。
.以概率P选取仿射变换 W形式为
X1=a*x0+b*y0+e
Y1=c*x0+d*y0+f;m[3][6]=;
}
voidIterateValue(floatx,floaty,floatz,float&outX,float&o utY,float&outZ)const
{
floata,b,c,d,e,f;//' 仿射变幻中的系数
floatR=(float)rand()/RAND_MAX;
if(R<=m[0][6])
{
a=m[0][0];b=m[0][1];c=m[0][2];d=m[0][3];e=m[0][4];f=m[0] [5];
}
elseif(R<=m[0][6]+m[1][6])
{
a=m[1][0];b=m[1][1];c=m[1][2];d=m[1][3];e=m[1][4];f=m[1] [5];
}
elseif(R<=m[0][6]+m[1][6]+m[2][6])
a=m[2][0];b=m[2][1];c=m[2][2];d=m[2][3];e=m[2][4];f=m[2] [5];
}
else
{
a=m[3][0];b=m[3][1];c=m[3][2];d=m[3][3];e=m[3][4];f=m[3] [5];
}
utX=(a*x*cosf(c/180))-(b*y*sinf(d/180))+e*FRACTAL_RADIUS;
utY=(a*x*sinf(c/180))+(b*y*cosf(d/180))+f*FRACTAL_RADIUS; utZ=z;
}
boolIsValidParamA()const{returntrue;}
voidSetParamA(floatv)
{
m_ParamA=v;
m[2][4]=v;
m[2][5]=v;
m[3][4]=1-v;
m[3][5]=v;
}
private:
floatm[4][7];//' 存放 IFS 码
};
复制代码
( 3)芦苇
复制代码
// 芦苇
classIFSBulrush:publicFractalEquation
{
public:
IFSBulrush()
{
m_StartX=;
m_StartY=;
m_StartZ=;
m_ParamA=;
floatk=m_ParamA*;
//'IFS 码赋值
m[0][0]=;m[0][1]=;m[0][2]=0;m[0][3]=0;m[0][4]=0; m[0][5]=0;m[0][6]=;
m[1][0]=;m[1][1]=;m[1][2]=k;m[1][3]=k;m[1][4]=1; m[1][5]=k/1600;m[1][6]=;
m[2][0]=;m[2][1]=;m[2][2]=0;m[2][3]=0;m[2][4]=0. 5f;m[2][5]=0;m[2][6]=;
voidIterateValue(floatx,floaty,floatz,float&outX,float&o utY,float&outZ)const
{
floata,b,c,d,e,f;//' 仿射变幻中的系数
floatR=(float)rand()/RAND_MAX;
if(R<=m[0][6])
{
a=m[0][0];b=m[0][1];c=m[0][2];d=m[0][3];e=m[0][4];f=m[0] [5];
}
elseif(R<=m[0][6]+m[1][6])
{
a=m[1][0];b=m[1][1];c=m[1][2];d=m[1][3];e=m[1][4];f=m[1]
[5];
}
else
{
a=m[2][0];b=m[2][1];c=m[2][2];d=m[2][3];e=m[2][4];f=m[2] [5];
}
utX=(a*x*cosf(c/180))-(b*y*sinf(d/180))+e*FRACTAL_RADIUS;
utY=(a*x*sinf(c/180))+(b*y*cosf(d/180))+f*FRACTAL_RADIUS;
utZ=z;
boolIsValidParamA()const{returntrue;}
voidSetParamA(floatv)
{
m_ParamA=v;