Fractal Estocástico
Publicado por Roberto Júnior 30/10/2006
[ Hits: 8.947 ]
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
Sintetizador sonoro que toca redez vouz 4
Servidor de Backup com Ubuntu Server 24.04 LTS, RAID e Duplicati (Dell PowerEdge T420)
Visualizar câmeras IP ONVIF no Linux sem necessidade de instalar aplicativos
Atualizar Debian Online de uma Versão para outra
Instalar driver Nvidia no Debian 13
Redimensionando, espelhando, convertendo e rotacionando imagens com script
Debian 13 Trixie para Iniciantes
Convertendo pacotes DEB que usam ZSTD (Padrão Novo) para XZ (Padrão Antigo)
Debian 13 "trixie&... lançado! (44)
Rust é o "C da nossa geração"? (3)
Gerenciador de arquivos é finalizado ao abrir pasta específica (2)