Estudo_Sinais_Sonar_Radar

N = 49;
n = [0:N];
Tx = 1/1000;
Ty = [0:0.00005:0.001];
f0 = 2*%pi*1500;
x = sin(f0*n*Tx);
h = sin(f0*Tx*(N - 1 - n));
a0 = 0.8;
//w = zeros(length(Ty), length());
w = zeros(length(Ty), 2*N + 1);
xv = zeros(length(Ty), length(n));
Axx = zeros(length(Ty), 2*N + 1);
phi = zeros(length(Ty), 2*N + 1);

for i=1:length(Ty)
   for j=1:length(n)
      xv(i,j) = sin(f0*j*Ty(i));
   end
end

for i=1:length(Ty)
    w(i,:) = convol(h,xv(i,:))
end

for i =1:length(n)
    phi(:, i) = w(:, n(i) + N - 1);
end

for i=1:length(Ty)
    Axx(i, :) = phi(i, :)*sqrt(Tx*Ty(i))
end

plot(Axx)


//for n = 0:N
//    for k = 0:100
//        Ea = (Axx[n*T][a0^k])^2
//end
//end
//f = ((Axx)^2)/Ea;
//for n = 0:999
//    for k = 0:100
//        nmd = n*f[n][k];
//    end
//end
//
//for n = 0:999
//    for k = 0:100
//        (enx)^2 = ((nmd - n)^2)*f[n][k];
//    end
//end
//
//for n = 0:999
//    for k = 0:100
//        kmd = k*f[n][k];
//    end
//end
//for n = 0:999
//    for k = 0:100
//        (ekx)^2 = ((kmd - k)^2)*f[n][k];
//    end
//end