Plot matched filter output vs Range of Radar

Question

Muhammad Kashif Latif 2022-4-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1694570-plot-matched-filter-output-vs-range-of-radar

回答： Abhimenyu 2023-12-1

I am trying to create matched filter output vs Range. I created time spaced vector using linspace command. The sampling frequency is 1.25 GHz. I am trying to find Range using R=ct/2. The code is as follows

u_filt=matched_filter(u_rx,h);

Mag= abs((abs(u_filt)).*2);

c=3e8;

fs =1.25e9;

t = linspace(0,length(u_filt),fs);

rangegates = c.*t;

rangegates = rangegates/2;

figure(8)

plot(rangegates,Mag);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Abhimenyu 2023-12-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1694570-plot-matched-filter-output-vs-range-of-radar#answer_1363999

在 MATLAB Online 中打开

Hi Muhammad,

I understand that you want to plot “matched filter output versus the range of the radar”. The “matched filter” can be implemented in various ways. In the code provided by you, a “matched_filter” function is used. This is not an in-built function of MATLAB. MATLAB has “phased.MatchedFilter” object that can be used to implement a “matched filter”. However, if the impulse response of the filter and the signal are already known as in your case, the “conv” function of MATLAB can be used to directly compute the filtered output.

Please refer to the example code below that can be used to plot “matched filter output versus the range of the radar”:

% Assuming u_rx is the received signal and h is the matched filter impulse response 
% Perform matched filtering using convolution 
u_filt = conv(u_rx, fliplr(h)); %"fliplr" function is used to flip the second signal for convolution 
% Take the absolute value of the matched filter output and scale by 2 
Mag = abs(u_filt) * 2; 
% Speed of light 
c = 3e8; 
% Sampling frequency 
fs = 1.25e9; 
% Time vector should be spaced according to the sampling frequency 
t = linspace(0, length(u_filt)/fs, length(u_filt)); 
% Calculate range gates using the formula R = c * t / 2 
rangegates = c * t / 2; 
% Plot the matched filter output versus range 
figure; 
plot(rangegates, Mag); 
xlabel('Range (m)'); 
ylabel('Matched Filter Output'); 
title('Matched Filter Output vs Range'); 