Fractal Estocástico
Publicado por Roberto Júnior 30/10/2006
[ Hits: 8.764 ]
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
Crud em C++ orientado a objetos com banco de dados MySQL
Exibi os números primos de um numero recebido pelo usuário, sem estrutura de repetição
Enviar mensagem ao usuário trabalhando com as opções do php.ini
Meu Fork do Plugin de Integração do CVS para o KDevelop
Compartilhando a tela do Computador no Celular via Deskreen
Como Configurar um Túnel SSH Reverso para Acessar Sua Máquina Local a Partir de uma Máquina Remota
Configuração para desligamento automatizado de Computadores em um Ambiente Comercial
Compartilhamento de Rede com samba em modo Público/Anônimo de forma simples, rápido e fácil
Cups: Mapear/listar todas as impressoras de outro Servidor CUPS de forma rápida e fácil
Criando uma VPC na AWS via CLI
Tem como instalar o gerenciador AMD Adrenalin no Ubuntu 24.04? (6)