eco.m

code for show a ecography using beam forming

function I=eco(E,N,S,v,Fs)
% I=eco(E,N,S,v,FS)
% calcula el eco en base a la señal de los trasductores introducida en un
% solo vector:
% E Señal del eco compuesto
% N numero de ecos en la señal
% S Separación de los transductores en mm
% v velocidad del sonido del medio
% Fs Frecuencia de muestreo

L=length(E)/N; % Longitudo del eco sacado del eco compuesto
Sep=(S*Fs)/v; % Separación entre los ecos
Ancho=round(2*L+round(Sep*N)); % Ancho y  Alto del ecograma
Alto=L;
I=zeros(Alto,Ancho); % Ecografia
SE=zeros(N,L); % Subecos

for M=0:N-1
    SE1(M+1,:)=E(1+M*L:L+M*L); % SubEco
    [m,p]=max(SE1(M+1,:)); % Maximo en el Subeco, donde llega la onda de presión principal
    SE(M+1,1:L-p+1)=SE1(M+1,p:L)/m; %Subeco recortado para empezar cuando llega el maximo
    %SE=SE1;
end

subplot (2,1,1);
plot (SE')

for Y = 1:Ancho % Recorre el ancho de la matriz
    for X = 1:Alto % Recorre el alto de la matriz
        MEDIA=0;
        for M=0:N-1
            P=round(sqrt((Y-L-(M*Sep))^2+X^2)); % distancia desde el punto de la matriz hasta la posición del transductor
            if P<L
                MEDIA=MEDIA+SE(M+1,P)*1/N;
            end
        end
            I(X,Y)=MEDIA;
        end
    end
end

subplot (2,1,2);
imshow (I,[-1 1])

end