雷达原理老师留的作业，干了一整天，收获蛮大，拿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 onfigure(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 onfigure(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 ony=filter(butterLow,mix);%用设计好的截止频率24GHz的低通滤波器滤波
figure(8)
plot(t,real(y))%时域图
title("输出信号时域")
xlabel("t/s")
xlim([0,1e-5])
grid onfigure(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 FrequencyN  = 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{sad}

MATLAB LFCM雷达调频法测距相关推荐

【雷达通信】基于matlab距离角度解耦法MIMO-OFDM雷达波束形成【含Matlab源码 2208期】
⛄一.获取代码方式获取代码方式1: 完整代码已上传我的资源:[雷达通信]基于matlab距离角度解耦法MIMO-OFDM雷达波束形成[含Matlab源码 2208期] 点击上面蓝色字体,直接付费下载 ...
matlab序列负轴,MATLAB中雷达信号回波图为什么会有负半轴
公告: 为响应国家净网行动,部分内容已经删除,感谢读者理解. 话题:MATLAB中雷达信号回波图为什么会有负半轴?哪条产生问题详情:re, plot(0:ts:(16*N-1)*ts,s_pc), ...
怎么用matlab显示噪声,怎么用MATLAB产生噪声调频信号
请问:怎么用matlab产生噪声调频信号? 如果诚心问,就不要问这么业余的问题.以下的BPSK波形成形,高频调制,解调程序,以做参考. %BPSK调制解调系统**************** clea ...
三维激光测距之三角法测距
在激光测距中涉及到三角法测距,但三角法测距的原理是什么呢? 红外发射器按照一定的角度发射红外光束,光束遇到被测物体后将会反射回来,反射回来的光线被CCD检测器检测,其获取一个偏移L.由图示以及已知的各 ...
matlab 坐标不用科学计数法,matlab不用科学计数法
『壹』 matlab中怎么才能不是科学计数法表示结果.比如1.0e+003 * 2.7581,怎么使它显示为2758.1谢谢了,很急啊 format long (小数位14) 或 format sho ...
MATLAB绘制雷达图并导出矢量图到Visio编辑(论文用图)
MATLAB绘制雷达图并导出矢量图到Visio编辑(论文用图) 前言: matlab绘制雷达图雷达图(Radar Chart)又被叫做蜘蛛网图(Spider Chart),适用于显示三个或更多的 ...
MATLAB之牛顿下山法
MATLAB之牛顿下山法算法原理 matlab程序算法原理上一篇博客,我介绍了牛顿法迭代法,接下来我就们接着讲解一下什么是牛顿下山法. 一.迭代公式在牛顿迭代过程中,若满足单调性|f(x(k+ ...
matlab曲面拟合的算法,用Matlab 实现移动曲面拟合法生成DEM
用Matlab实现移动曲面拟合法生成DEM 杜玉军 (武汉大学测绘工程0408班 200431610007 武汉 430079) 摘要:移动曲面拟合法是DEM格网点内插常用的一种方法,利用Matl ...
雷达杂波matlab,基于MATLAB的雷达杂波建模与仿真研究
第 47 卷第 3 期( 总第 185 期) 2018 年 9 月火控雷达技术 Fire Control Radar Technology Vol. 47 No. 3( Series 185) S ...

MATLAB LFCM雷达调频法测距