雷達原理老師留的作業，干了一整天，識訓蠻大，拿word寫的報告，夜已經深了（21點），懶得轉成LaTeX了，就直接截圖放上來了，
一、問題提出
LFCM（線性調頻連續波）雷達，發射信號帶寬1GHz,時寬10us. 3個目標徑向距離分別為100m,105m,105.1m，分析中頻信號帶寬，仿真產生回波信號及去斜處理程序，對中頻輸出做傅里葉變換，觀察距離維回波特性并進行討論分析，在此基礎上，討論該類雷達性能指標和信號引數之間的關系，
二、問題解決
1.線性調頻信號（LFM）原理：
2.線性調頻信號（LFM）產生：
3.雷達回波產生：
4.下混頻原理：
5.下混頻實作：
三、心得感悟
雖說截圖寫博客確實不好，但是確實太困了，如果有需要word版的可以評論或者私信找我，這些地方我也學的不是很懂，要是有發現可以改進的地方，歡迎評論區指正！
附全部代碼：
fs=20e10;
T=10e-6;%時寬
B=1e9;%帶寬
mu=B/T;%調頻率
n=round(T*fs);%采樣點個數
t=linspace(0,T,n);
f0=24e9;%載頻
s=exp(2j*pi*(f0*t+0.5*mu*t.^2));
plot(t,real(s))%時域圖
title("發射信號時域")
xlabel("t/s")
xlim([0,1e-9])
figure(2)
S=fftshift(fft(real(s))./n);
f=linspace(-fs/2,fs/2-1,n);%頻域橫坐標，注意奈奎斯特采樣定理，最大原信號最大頻率不超過采樣頻率的一半
plot(f,abs(S))%頻域圖
title("發射信號頻域")
xlim([-3e10,3e10])
xlabel("f/Hz")
grid on
figure(3)
w=f0+mu.*t;
plot(t,w)%線性調頻圖
grid on
title("發射信號調制圖")
xlabel("t/s")
ylabel("f/Hz")
pos=[100,105,105.1];%三個目標的位置/m
timedelay=2*pos/(3e8);%雷達回波時延
r1=exp(2j*pi*(f0*(t-timedelay(1))+0.5*mu*(t-timedelay(1)).^2));
r2=exp(2j*pi*(f0*(t-timedelay(2))+0.5*mu*(t-timedelay(2)).^2));
r3=exp(2j*pi*(f0*(t-timedelay(3))+0.5*mu*(t-timedelay(3)).^2));
figure(4)
subplot(411)
plot(t,real(r1))%時域圖
title("100m回波信號時域")
xlim([0,1e-5])
xlabel("t/s")
subplot(412)
plot(t,real(r2))%時域圖
xlim([0,1e-5])
title("105m回波信號時域")
xlabel("t/s")
subplot(413)
plot(t,real(r3))%時域圖
xlim([0,1e-5])
title("105.1m回波信號時域")
xlabel("t/s")
r=r1+r2+r3;%回波信號加在一起
subplot(414)
plot(t,real(r))%時域圖
xlim([0,1e-5])
title("全部回波信號時域")
xlabel("t/s")
figure(5)
R=fftshift(fft(real(r))./n);
f=linspace(-fs/2,fs/2-1,n);%頻域橫坐標，注意奈奎斯特采樣定理，最大原信號最大頻率不超過采樣頻率的一半
plot(f,abs(R))%頻域圖
title("回波信號頻域")
xlim([-3e10,3e10])
xlabel("f/Hz")
grid on
%回波信號乘以本振信號做下混頻
mix=real(s).*real(r);
figure(6)
plot(t,real(mix))%時域圖
title("混頻信號時域")
xlabel("t/s")
xlim([0,1e-5])
grid on
figure(7)
MIX=fftshift(fft(mix)./n);
f=linspace(-fs/2,fs/2-1,n);%頻域橫坐標，注意奈奎斯特采樣定理，最大原信號最大頻率不超過采樣頻率的一半
plot(f,abs(MIX))%頻域圖
title("混頻信號頻域")
xlim([-6e10,6e10])
ylim([0,0.3])
xlabel("f/Hz")
grid on
y=filter(butterLow,mix);%用設計好的截止頻率24GHz的低通濾波器濾波
figure(8)
plot(t,real(y))%時域圖
title("輸出信號時域")
xlabel("t/s")
xlim([0,1e-5])
grid on
figure(9)
Y=fftshift(fft(y)./n);
f=linspace(-fs/2,fs/2-1,n);%頻域橫坐標，注意奈奎斯特采樣定理，最大原信號最大頻率不超過采樣頻率的一半
plot(f,abs(Y))%頻域圖
title("輸出信號頻域")
xlim([-10^8,10^8])
ylim([0,0.3])
xlabel("f/Hz")
grid on
其中butterLow濾波器代碼：
function Hd = butterLow
%BUTTERLOW Returns a discrete-time filter object.
% MATLAB Code
% Generated by MATLAB(R) 9.5 and Signal Processing Toolbox 8.1.
% Generated on: 21-Apr-2021 11:52:59
% Butterworth Lowpass filter designed using FDESIGN.LOWPASS.
% All frequency values are in Hz.
Fs = 2e+11; % Sampling Frequency
N = 7; % Order
Fc = 24000000000; % Cutoff Frequency
% Construct an FDESIGN object and call its BUTTER method.
h = fdesign.lowpass('N,F3dB', N, Fc, Fs);
Hd = design(h, 'butter');
% [EOF]
留個贊再走啊，不然都沒有信心寫下一篇博客了\emoji{said}
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/278961.html
標籤：其他