Форум программистов, компьютерный форум, киберфорум
Matlab
Войти
Регистрация
Восстановить пароль
 
710 / 529 / 238
Регистрация: 21.10.2012
Сообщений: 1,720
1

Автофокусировка РСА

06.03.2020, 15:21. Просмотров 358. Ответов 1
Метки нет (Все метки)

Здравствуйте. Пытаюсь применить фазоразностный алгоритм автофокусировки в РСА. Алгоритм написал, но, видимо, что-то я неправильно понимаю или что-то не так написал, но результаты получаются практически идентичные (что с фокусировкой, что без неё). Может ли кто-то подсказать, что не так.
Считается это все долго, так как приходится рассчитывать поправки для каждого периода синтезирования. Чтобы просто быстро построить РСА без автофокусировки, то нужно запустить закомментированный код. На рисунках - то, что получается (хоть с фокусировкой, хоть без неё и то, что в идеале бы хотелось увидеть).
Про сам алгоритм оставил ссылку:
http://www.tehnoinfa.ru/radiovidenie/25.html

Кликните здесь для просмотра всего текста
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
clear;
clc;
 
tau_imp = 0.5e-6; % длительность импульса
omega = 190e6; % девиация частоты
fs_now = 200e6; % частота дискретизации
 
tt = (-tau_imp/2 : 1/fs_now : tau_imp/2 - 1/fs_now)'; % время для генерации ЛЧМ сигнала
 
signalQI = exp(1j * 2 * pi * omega .* tt.^2 ./ (2*tau_imp)); % ЛЧМ
signalQI1 = signalQI(end : -1 : 1);
 
txVel = 20; % скорость носителя
Trep = 1 / 1562; % период излучения
minrange = 450; % минимальная наклонная дальность
c = physconst('LightSpeed');
waveLength = c / 9550e6;
 
L = 2 * minrange * tand(8/2); % расстояние, которое попадает в ДН антенны
x = -L/2 : txVel * Trep : L/2; % расстояние на котором проводим синтез (для опорной функции)
 
size0 = 4096 - numel(tt) + 1;
size1 = 4096/2 - numel(tt) + 1;
size2 = 10000;
 
x = repmat(x, [size1 1]);
 
R = minrange + (0:size0-1)' * c/(2*fs_now); % наклонная дальность
R = R(size0 - size1 + 1 : end); % оставляем только то, что будем использовать
 
ref = exp(1j*2*pi/waveLength*x.^2./R); % опорная функция для сжатия по азимуту
 
Tsynt = abs(x(1)/txVel*2); % время синтеза
fs = 1/Trep;
n = size(x, 2)/2;
F = -fs/2: fs/n : fs/2 - fs/n;
slowTime = -Tsynt/2 : Trep : Tsynt/2;
img = zeros(size1, size2-size(x, 2)+1);
 
load('matlab.mat'); % отсчеты комплексной огибающей; каждый столбец - новый период излучения
% строки - отсчеты принятого сигнала
 
trSigMatchFilt = conv2(Ur_filt1, signalQI1, 'valid'); % убираю ЛЧМ сигнал в каждом периоде (сжимаю по дальности)
 
% автофокусировка + сжатие по азимуту
for j = 1:size(trSigMatchFilt, 2)-numel(slowTime)+1
    temp = trSigMatchFilt(:, j:j+numel(slowTime)-1); % беру кусок сигнала за время синтеза
    
   % автофокусировка
    temp1 = conj(temp(:, 1:numel(slowTime)/2)) .* temp(:, numel(slowTime)/2 + 1:end);
    temp2 = fft(temp1, [], 2);
    temp3 = sum(abs(temp2), 1);
    [~, ww] = max(fftshift(temp3, 2));
    coeff = 2*pi*F(ww)/Tsynt;
    fi = exp(-1j .* 2 .* pi * coeff .* slowTime.^2);
    fi = repmat(fi, [size(temp, 1) 1]); % получил корректирующие значения
   
    img(:, j) = sum(conj(ref) .* temp .* fi, 2); % корректировка + сжатие по азимуту
end
 
Z = abs(img);
Z = Z/max(Z(:));
 
 
% быстрая свертка по азимуту
% ref1 = fft(ref, size2, 2);
% Stauf = fft(trSigMatchFilt, [], 2);
% Stauf = Stauf .* conj(ref1);  % взаимная корреляция а не свертка
% img = ifft(Stauf, [], 2);
% Z = abs(img);
% Z = Z/max(Z(:));
% Z(:, end - (size(x, 2) - 1) + 1 : end) = []; 
 
imagesc(Z);
colormap gray;
ax = gca;
ax.YDir = 'normal';
0
Миниатюры
Автофокусировка РСА   Автофокусировка РСА  
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
06.03.2020, 15:21
Ответы с готовыми решениями:

Провести компрессию (сжатие) информации методом главных компонент (РСА) и на основе системы Э. Оя
Здравствуйте! Есть задание провести компрессию (сжатие) информации методом главных компонент (РСА)...

Интересная автофокусировка
Интереса ради создал рекламную компанию, чисто позырить как работает.. :) Сегодня пришло письмо:...

Автофокусировка между объектами
Здравствуйте! Пишу кроссворд. В основном использую jQuery. Сам кроссворд состоит из таблицы и...

Лазерная автофокусировка в смартфонах сейчас необходимость?
Задался вопросом, а является ли сейчас необходимость в лазерной автофокусировки в современных...

1
325 / 196 / 99
Регистрация: 11.01.2019
Сообщений: 783
06.03.2020, 23:23 2
Была такая книга "Радио? Это очень просто!", а РСА? -Это походу очень сложно!
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
06.03.2020, 23:23

Заказываю контрольные, курсовые, дипломные и любые другие студенческие работы здесь.

CryptoAPI | РСА шифрование
Приветствую всех. Пытаюсь шифровать текст при помощи RSA. Имеются следующие вопросы : 1. Дабы...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.