Cómo implementar este algoritmo usando Matlab

En el futuro, este tipo de preguntas es mucho más apropiado para Stack Overflow.

Pero veo dos problemas potenciales. Pero admitiré que no conozco el algoritmo

Primero es que el bucle for interno no hace nada más que calcular la opción j = 4 y almacenarla. Si desea la suma de esos valores, debe hacerlo correctamente.

En segundo lugar, su operador de comparación para deltaE analiza todos los DeltaEs, no el más reciente. Esto es probablemente lo que realmente está causando su problema. Intente ejecutar el clode a continuación para una explicación simple.

A = unos (5,1);
B = -ones (5,1);
C = [A, B];

si (A <0)
disp (‘Hola, estoy roto’)
fin

si (B <0)
disp (‘Hola, estoy trabajando pero no muy bien definido’)
fin

si (C <0)
disp (‘Hola, estoy ignorando ciertos valores’)
fin

Puedo agregar una cosa a la respuesta de William, ya que sé el algoritmo que estás tratando de implementar:

se supone que tienes una red de valores de giro, algo así como N = 100; giros = 2 * pi * rand (N); Luego tomas un punto al azar usando algo como coordenadas = randi ([1, N], 1,2); Luego, compara esos puntos con sus vecinos, teniendo en cuenta sus condiciones de contorno.