Problema de optimización del código Matlab. ¿Puede algún experto optimizar este código? El código fuente se ejecuta demasiado lento.

Tengo algunas dudas. Si realiza un bucle de esta manera, el valor de la distancia se actualizará cada vez que se realice el bucle. ¿No es el último valor registrado solo el valor correspondiente a la última posición de lig y col?

[col, lig] = find(fsm == nm); Determine las coordenadas que deben buclearse

for id = 1: length(col)

for area = 1：narea

dist(area) = ...

trace(cohinv(：,：,area)*shiftdim(T3mean(col(id)， lig( id),:,:)));

end

end

Escribir de esta manera no acelerará mucho en principio. . . .

Se recomienda utilizar la interfaz c/c de matlab y utilizar archivos mex para acelerar. . .