Fractal Estocástico
Publicado por Roberto Júnior 30/10/2006
[ Hits: 9.106 ]
Homepage: http://sites.google.com/site/dibetao
Este programa gera um fractal dito ser aleatório, pois é gerado por um processo estocástico. Uma matriz é usada para armazenar os valores. As dimensões da matriz poderão ser alteradas no fonte do programa.
PS: É necessário o GNUPLOT para gerar o gráfico.
/**
FRACTAL 1
Autor: Roberto Lima
Data: 06/08/2006
Email: dibetao@gmail.com
- Este programa gera um fractal simples.
- Ajuste o numero de linha e colunas da
matriz em LIN e COL respectivamente.
- Ajuste as dimensoes de A[][], para valores
maiores que LIN e COL.
- O programa finaliza assim que ocorre a
percolacao.
- Um arquivo 'pontos.txt' eh gerado, e o grafico
poderá ser visualizado com o GNUPLOT, com o comando
--> plot 'pontos.txt'.
*/
#include <stdio.h>
#include <stdlib.h>
#define LIN 211
#define COL 211
int A[251][251],pp[5],i,j;
//---------------------------------------------------------------------
int vv(int ii,int jj) //verifica vizinhanca
{
int val;
val=0;
if (A[ii+1][jj]==1) val=1;
if (A[ii-1][jj]==1) val=1;
if (A[ii][jj+1]==1) val=1;
if (A[ii][jj-1]==1) val=1;
return(val);
}
//---------------------------------------------------------------------
np()//nova particula
{
int c1=0,pos;
do
{
pos=rand() % 4;
switch (pos)
{
case 0: {
if(pp[2]==1) continue;
i=1; j= (rand() % COL) +1;
break;
}
case 1: {
if (pp[4]==1) continue;
j=1; i= (rand() % LIN) +1;
break;
}
case 2: {
if (pp[1]==1) continue;
i=LIN; j= (rand() % COL) +1;
break;
}
case 3: {
if (pp[3]==1) continue;
j=COL; i= (rand() % LIN) +1;
break;
}
}
if (A[i][j]==0) {
A[i][j]=1;
c1=1;
}
}while (c1!=1);
}
//---------------------------------------------------------------------
mov() //prob movimento
{
int p,ii,jj;
A[i][j]=0;
p=(rand() % 4)+1;
switch (p)
{
case 1: j=j+1; break;
case 2: j=j-1; break;
case 3: i=i+1; break;
case 4: i=i-1; break;
}
if (i>LIN || i<1 || j>COL || j<1) np();
else{
A[i][j]=1;
if (vv(i,j)==1){
if (i==LIN || i==1 || j==COL || j==1)
{
if(i==LIN) pp[1]=1;
if(i==1) pp[2]=1;
if(j==COL) pp[3]=1;
if(j==1) pp[4]=1;
i=0; j=0;
}else np();
}
}
}
//---------------------------------------------------------------------
salva()
{
int ii,jj;
FILE *arq;
arq=fopen("pontos.txt","w+");
for (ii=1; ii<=LIN; ii++)
for (jj=1; jj<=COL; jj++)
if (A[ii][jj]!=0)
fprintf(arq,"%d %d %d\n",ii,jj,A[ii][jj]);
fclose(arq);
}
//---------------------------------------------------------------------
int main()
{
int ii,jj,parar;
float step=0;
srand(time(NULL));
step=0;
parar=0;
for(ii=1;ii<=4;ii++) pp[ii]=0;
for (ii=1; ii<=LIN; ii++) for (jj=1; jj<=COL; jj++) A[ii][jj]=0;
A[LIN/2+1][COL/2+1]=1; //particula no meio da matriz
np();
while (parar!=1)
{
//if(step%10000==0)printf("%f\n",step);//Se necessario verificar loop infinito
if(pp[1]==1 && pp[2]==1 && pp[3]==1 && pp[4]==1)//Quando percolar
{
parar=1;
printf("ACABOU - TOTAL= %f steps\n",step);
salva();
}
else
{
step=step+0.001;
mov();
}
}
}
//FIM
Testar o melhor método de organização C (inserção, bolha e shell-sort)
Gerando uma "number list" de 6 digitos
Cinnamon seria a aposta acertada frente às outras interfaces gráficas mais populares?
KDE Plasma - porque pode ser a melhor opção de interface gráfica
Gentoo: detectando impressoras de rede e como fixar uma impressora por IP
Como o GNOME conseguiu o feito de ser preterido por outras interfaces gráficas
Por que sua empresa precisa de uma PKI (e como automatizar EMISSÕES de certificados via Web API)
Instalando NoMachine no Gentoo com Systemd (acesso Remoto em LAN)
Vou destruir sua infância:) (8)
Interface gráfica com problema (2)
Instalar Linux em notebook Sony Vaio VPCEG13EB (13)









