没有合适的资源?快使用搜索试试~ 我知道了~
脉冲雷达仿真
3星 · 超过75%的资源 需积分: 48 87 下载量 105 浏览量
2018-04-12
16:29:35
上传
评论 12
收藏 10KB TXT 举报
温馨提示
试读
19页
仿真脉冲多卜勒雷达的信号处理。设脉冲宽度为各学生学号末两位数,单位为μs,重复周期为200μs,雷达载频为10GHz,输入噪声为高斯白噪声。目标模拟分单目标和双目标两种情况,目标回波输入信噪比可变(-35dB~10dB),目标速度可变(0~1000m/s),目标幅度可变(1~100),目标距离可变(0~10000m),相干积累总时宽不大于10ms。单目标时,给出回波视频表达式;脉压和FFT 后的表达式;仿真给出脉压和FFT 后的输出图形;通过仿真说明各级处理的增益,与各级时宽和带宽的关系;仿真说明脉压时多卜勒敏感现象和多卜勒容限及其性能损失(脉压主旁比与多卜勒的曲线)。双目标时,仿真出大目标旁瓣盖掩盖小目标的情况;仿真出距离分辨和速度分辨的情况。 注意: 1、白噪声加入采用randn函数;2、整体产生回波,再整体加时延和多卜勒信息;3、通过数据计算出输出信噪比;
资源推荐
资源详情
资源评论
clear all;
close all;
clc;
Trp=200e-6; %脉冲重复周期us
tm=36e-6; % 脉冲宽度6us
N=1024; % 相干积累时间小于10ms,Nmax=10ms/200us=50
f0=10e9;%载波频率
fs=10e5;%采样频率
ts=1/fs;%一个采样点时间
c=3e8; % 光速
namna=c/f0;%波长
R0=5000; %距离
v=100; %目标径向速度
A=2; %回波幅度
NTrp=fs*Trp;%一个脉冲时间内采样点个数
Ntm=fs*tm;%一个脉宽内采样点个数
SampleNumber=NTrp;%一个脉冲时间内采样点个数
PulseNumber=N;%总脉冲数
TotalNumber=N*NTrp;%总采样点个数
t=0:ts:N*Trp-ts;% time
i=1:N;
ti=2*(R0-(i-1).*v*Trp)/c;%delay of different Trp
%plot(real(Uo))
for i=1:N
tempu(i)=A*exp(j*2*pi*(i-1)*Trp*(2*v)/namna);
close all;
clc;
Trp=200e-6; %脉冲重复周期us
tm=36e-6; % 脉冲宽度6us
N=1024; % 相干积累时间小于10ms,Nmax=10ms/200us=50
f0=10e9;%载波频率
fs=10e5;%采样频率
ts=1/fs;%一个采样点时间
c=3e8; % 光速
namna=c/f0;%波长
R0=5000; %距离
v=100; %目标径向速度
A=2; %回波幅度
NTrp=fs*Trp;%一个脉冲时间内采样点个数
Ntm=fs*tm;%一个脉宽内采样点个数
SampleNumber=NTrp;%一个脉冲时间内采样点个数
PulseNumber=N;%总脉冲数
TotalNumber=N*NTrp;%总采样点个数
t=0:ts:N*Trp-ts;% time
i=1:N;
ti=2*(R0-(i-1).*v*Trp)/c;%delay of different Trp
%plot(real(Uo))
for i=1:N
tempu(i)=A*exp(j*2*pi*(i-1)*Trp*(2*v)/namna);
for k=1:NTrp
if(k>ti(i)*fs && k<(ti(i)+tm)*fs)
u((i-1)*NTrp+k)=tempu(i);% 回波的产生
else
u((i-1)*NTrp+k)=0;
end
end
end
[N0,wn]=buttord(1/36,1/20,0.5,30);
[b0,a0]=butter(N0,wn,'low');%巴特沃兹低通滤波器
%s12=filter(b0,a0,s1);
noise=randn(1,length(u(:)));%初始噪声,带宽为1MHZ
fnoise=filter(b0,a0,noise);%滤波噪声,带宽=信号带宽=27.8khz
xishu=sum(fnoise.^2)/length(fnoise(:))
SNR1=3;%信噪比
%sigPower1 = sum(abs(u(:)).^2)/length(u(:));
noisePower1=A*A/(10^(SNR1/10))/xishu;%限制噪声功率,使其满足snr
s12=u+sqrt(noisePower1)*fnoise;%加噪声的回波信号
%s1=awgn(u,SNR1,'measured');
%{
figure(8)
plot(real(s1));
figure(9)
plot(abs(fft(u)));
%}
if(k>ti(i)*fs && k<(ti(i)+tm)*fs)
u((i-1)*NTrp+k)=tempu(i);% 回波的产生
else
u((i-1)*NTrp+k)=0;
end
end
end
[N0,wn]=buttord(1/36,1/20,0.5,30);
[b0,a0]=butter(N0,wn,'low');%巴特沃兹低通滤波器
%s12=filter(b0,a0,s1);
noise=randn(1,length(u(:)));%初始噪声,带宽为1MHZ
fnoise=filter(b0,a0,noise);%滤波噪声,带宽=信号带宽=27.8khz
xishu=sum(fnoise.^2)/length(fnoise(:))
SNR1=3;%信噪比
%sigPower1 = sum(abs(u(:)).^2)/length(u(:));
noisePower1=A*A/(10^(SNR1/10))/xishu;%限制噪声功率,使其满足snr
s12=u+sqrt(noisePower1)*fnoise;%加噪声的回波信号
%s1=awgn(u,SNR1,'measured');
%{
figure(8)
plot(real(s1));
figure(9)
plot(abs(fft(u)));
%}
剩余18页未读,继续阅读
资源评论
- 王林森2019-06-24资料还行。可以看下。
- moodybat2021-12-09没什么参考意义
jinqu——✌
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功