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 muestreoL=length(E)/N; % Longitudo del eco sacado del eco compuestoSep=(S*Fs)/v; % Separación entre los ecosAncho=round(2*L+round(Sep*N)); % Ancho y  Alto del ecogramaAlto=L;I=zeros(Alto,Ancho); % EcografiaSE=zeros(N,L); % Subecosfor 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;endsubplot (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    endendsubplot (2,1,2);imshow (I,[-1 1])end