Двумерная свертка

x = 1:100;
y = complex(cos(x),sin(x));
 
z = zeros(100,100);
s = zeros(11,100);
for i =1 : 10
    z(i,:) = conj(y);
    s(i,:) = (y);
end
 
 
I = conv2(z,s);

lambda = 0.69;
dt = 76.9*10^(-3);
V = 55;
Ti = 30.8 * 10^(-9);
dR = 299792458*Ti/2;
x = -550:V*dt:550;
 
r0 = 3000;%sqrt(h^2 + y^2);
R2 = (r0+(x.^2)./(2*r0));
R = r0:dR:r0+255*dR;
U = 1./(R2.^4);
fi = 0;
A = (U.*exp(1i*(fi + 4*pi*r0/lambda)));
G = 1;
S = exp(1i*((2*pi*x.^2)./(lambda*r0)));
signal2 = A.*G.*S;
 
signal = zeros(13,length(x));
signal1 = zeros(length(R),length(x));
for j = 1:length(R)
   for i=1:length(R2)
       if ( R2(i)>= R(j) && R2(i) < R(j)+dR )
           signal(j,i) = signal2(i);
           signal1(j,i) = signal2(i);
       end
   end
end
 
%моя свертка
I = zeros(length(R),length(x));
F = zeros(length(R)+13,length(x)*2);
F(1:length(R),1:length(x)) = signal1(1:length(R),1:length(x));
 
for i = 1:length(R)
    for j = 1:length(x)      
       J = abs(ifft2(fft2(signal).*fft2(F(i:i+13-1,j:j+length(x)-1))));
       I(i,j) = J(1,1);
    end
    i
end
 
 
%свертка матлаба
I2 = conv2(signal1,signal);
 
mesh(I)
axis([0 300 0 20 ]);
figure,mesh(abs(I2))
axis([0 600 0 40 ]);
figure, mesh(abs(signal1))
axis([0 300 0 20 ]);

lambda = 0.69;  %просто константа
dt = 76.9*10^(-3);         %параметр отвечает за дискрет по оси х (в идеале степень должна быть -6 )
V = 55;                         %параметр отвечает за дискрет по оси х (x = vt)
Ti = 30.8 * 10^(-9);       %просто константа
dR = 299792458*Ti/2;      %параметр отвечает за дискрет по оси y
x = -550:V*dt:550;         % один из параметров по оси х
 
r0 = 3000;%sqrt(h^2 + y^2);      %начало с которого идут отсчеты по оси y
R2 = (r0+(x.^2)./(2*r0));            %значения расстояния (для оси у)
R = r0:dR:r0+255*dR;                 % отсчеты по оси у (расстояние) 256 отсчетов
U = 1./(R2.^4);                         %амплитуда
fi = 0;                                      %константа
A = (U.*exp(1i*(fi + 4*pi*r0/lambda)));           %одна из составляющих синала 
G = 1;                                                        %одна из составляющих синала 
S = exp(1i*((2*pi*x.^2)./(lambda*r0)));           %одна из составляющих синала 
signal2 = A.*G.*S;                                        %сам сигнал, точнее его часть отличная от 0
 
%получился сигнал, которыя по сути зависит от одного параметра х
% но строить его можно по двум параметра х и R
% В итоге получилась парабола в пространстве
% Но по оси y слишком часто взяты отсчеты, поэтому дискреты параметра R2
% нужно уложить в дискреты R
% и создать матрицу сигнал 
 
signal = zeros(13,length(x));      %опорная матрица, заполненная нулями
signal1 = zeros(length(R),length(x));          %матрица сигнала, заполненная нулями
for j = 1:length(R)                                   %проход по дискретам R
   for i=1:length(R2)                                 %цикл по дискретам R2
       if ( R2(i)>= R(j) && R2(i) < R(j)+dR )    %проверка на то, в какойц дискрет попадает точка полученного сигнала
           signal(j,i) = signal2(i);                   %присваивание значения сигнала определенной ячейке матрицы
           signal1(j,i) = signal2(i);                  %присваивание значения сигнала определенной ячейке матрицы
       end
   end
end
 
%моя свертка
I = zeros(length(R),length(x));                  %матрица для результатов 
F = zeros(length(R)+13,length(x)*2);         %матрица по которой буду проходить шаблоном, соответственно по 
                                                            %размеру увеличена на размеры шаблона
F(1:length(R),1:length(x)) = signal1(1:length(R),1:length(x));       %присваиваю ей значения сигнала
 
for i = 1:length(R)       %проход шаблоном по оси у
    for j = 1:length(x)       %проход шаблоном по оси х
       J = abs(ifft2(fft2(signal).*fft2(F(i:i+13-1,j:j+length(x)-1))));       %свертка
       I(i,j) = J(1,1);    %в свертке получается один верный результат, значение которого и берутся для результата
    end
    i
end
 
 
%свертка матлаба
I2 = conv2(signal1,signal);
 
 
%вывод графиков
mesh(I)
axis([0 300 0 20 ]);
figure,mesh(abs(I2))
axis([0 600 0 40 ]);
figure, mesh(abs(signal1))
axis([0 300 0 20 ]);