Mudanças entre as edições de "Simulação de Micélio de Fungo"
De Física Computacional
Linha 45: | Linha 45: | ||
+ | </source> | ||
+ | |||
+ | |||
+ | '''Divisão''' | ||
+ | |||
+ | <source lang="python"> | ||
+ | |||
+ | def divisao (x,y): | ||
+ | |||
+ | theta= np.arctan(y/x) | ||
+ | |||
+ | if (x<0) : | ||
+ | theta= theta+ math.pi | ||
+ | |||
+ | angulodivisao= random.random()*math.pi/2 #angulo para divisão de no máximo 90° | ||
+ | |||
+ | angulo1= theta- angulodivisao/2 | ||
+ | angulo2= theta+ angulodivisao/2 | ||
+ | |||
+ | addx1 = r * math.cos(angulo1) | ||
+ | addy1 = r * math.sin(angulo1) | ||
+ | |||
+ | addx2= r * math.cos(angulo2) | ||
+ | addy2= r * math.sin(angulo2) | ||
+ | |||
+ | Ax = x + addx1 | ||
+ | Ay = y + addy1 | ||
+ | |||
+ | Bx = x + addx2 | ||
+ | By= y +addy2 | ||
+ | |||
+ | return (Ax,Ay,Bx,By) | ||
</source> | </source> |
Edição das 17h00min de 20 de maio de 2021
Usando como base o segundo e primeiro capítulo do livro [A Hybrid Mathematical Model of Fungal Mycelia: Tropisms, Polarised Growth and Application to Colony Competition.](https://core.ac.uk/download/pdf/6117416.pdf)
- Mecanismos gerais dos modelos
- Crescimento
- A distribuição de nutrientes ocorre de maneira discreta, não contínua.
- Transporte de nutriente interno # Não será feito por exigir a computação de cada elemento individualmente toda vez
Fungo 1:
- Branching na ponta
- Computar só os da ponta
- Crescimento ocorre só com o nutriente que o ponto final está
Fungo 2:
- Branching lateral
- Morte (implementar idade da linha)
Fungo 3:
- Junção (computação de cada elemento individualmente toda vez) (anastomosis)
Crescimento
def crecimento (x,y):
theta= np.arctan(y/x)
if (x<0) :
theta= theta+ math.pi
aleatorio_theta = random.random()* math.pi/4 - math.pi/8 # angulo de -22.5 até 22.5 (45°)
theta=theta+aleatorio_theta
addx = r * math.cos(theta)
addy = r * math.sin(theta)
fx = x + addx
fy = y + addy
return (fx,fy)
Divisão
def divisao (x,y):
theta= np.arctan(y/x)
if (x<0) :
theta= theta+ math.pi
angulodivisao= random.random()*math.pi/2 #angulo para divisão de no máximo 90°
angulo1= theta- angulodivisao/2
angulo2= theta+ angulodivisao/2
addx1 = r * math.cos(angulo1)
addy1 = r * math.sin(angulo1)
addx2= r * math.cos(angulo2)
addy2= r * math.sin(angulo2)
Ax = x + addx1
Ay = y + addy1
Bx = x + addx2
By= y +addy2
return (Ax,Ay,Bx,By)