雷达系统建模与仿真报告Word下载.docx
《雷达系统建模与仿真报告Word下载.docx》由会员分享,可在线阅读,更多相关《雷达系统建模与仿真报告Word下载.docx(32页珍藏版)》请在冰豆网上搜索。
其概率密度函数如下:
lamade=4*201+1;
a=6;
b=10;
x
(1)=11;
M=2^36;
end;
%%%%%%%%%%%%%%%%%%%%%
i=2:
y(i)=(b-a)*x(i)+a;
n=5:
0.1:
11;
hist(y,n),axis([a-1b+10max(hist(y,n))+20]);
mean(y)
var(y)
上面程序中取a=6,b=10.即(6,10)区间上的均匀分布。
均值=8.0070方差=1.3311
2.2高斯分布:
高斯分布的概率密度函数如下;
其产生方法是在均匀分布随机数的基础上通过函数变换法来产生。
产生步骤是①产生均匀分布的随机数。
②产生服从标准正态分布的随机数。
③由标准正态分布产生一般正态分布。
标准正态分布
其部分程序如下:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
i=1:
u(i)=sqrt(-2*log(x(i))).*cos(2.*pi.*y(i));
v(i)=sqrt(-2*log(x(i))).*sin(2.*pi.*y(i));
n1=-5:
0.2:
5;
n2=-5:
subplot(1,2,1);
hist(u,n1);
subplot(1,2,2);
hist(v,n2);
mean(u)
var(u)
mean(v)
var(v)
均值=-0.0182方差=0.9910
一般正态分布
a=2;
b=2;
i=1:
v=b*u+a;
n=-10:
10;
hist(v,n);
均值=2.0301方差=4.0482
2.3指数分布:
服从正态分布的信号通过线性检波器后其包络强度(功率)服从指数分布。
其概率密度函数为:
其产生方法亦有:
①在均匀分布随机数的基础上产生指数分布。
②在正态分布随机数的基础上产生该分布。
产生程序分别如下:
程序1(部分)
lamade1=1
y(i)=-log(x(i))./lamade1;
n=0:
hist(y,n);
运行结果:
均值=1.0140方差=1.0292
程序2(部分)
s(i)=(u(i).*u(i)+v(i).*v(i));
0.3:
25;
hist(s,n);
mean(s)
ar(s)
2.4瑞利分布:
在雷达系统中载带信号的包络服从瑞利分布。
正态随机过程在其杂波载频
上可以表示为:
其中
、
是服从
的相互独立的随机过程,检波器的包络幅度(电压):
服从瑞利分布
。
瑞利分布的概率密度函数为:
①在均匀分布随机数的基础上产生瑞利分布。
程序1(部分):
segma=2;
y(i)=segma*sqrt(-2*log(x(i)));
var(y)
均值=2.5239方差=1.7417
程序2(部分):
s(i)=sqrt(u(i).*u(i)+v(i).*v(i));
var(s)
均值=1.2537方差=0.4317
2.5广义指数分布
概率密度函数为:
式中:
s-信噪比
部分程序如下:
s=8;
h(i)=u(i)+sqrt(2*s);
z(i)=h(i).*h(i)+y(i).^2;
60;
hist(z,n);
均值=17.1432方差=69.0430
2.6广义瑞利分布
-信噪比
a=1;
10240;
s(i)=sqrt((u(i)+a).^2+v(i).^2);
n=-1:
15;
均值=1.5539方差=0.6022
2.7韦布尔分布
韦布尔分布模型的性能介于瑞利分布模型与对数一正态分布模型之间.海浪杂波和地面杂波都可以用它来表示,并且在一个相当宽的条件围它能精确地表示实际的杂波分布。
韦布尔分布的概率密度函数为:
a-形状参数;
b-比例参数;
x0-位置参数;
该分布是在服从瑞利分布随机数的基础上用变换法产生的,其产生源程序(部分)及直方图如下:
a=3;
m=5;
y(i)=m+b*(-log(x(i)).^(1/a));
y=m+b*((-log(x)).^(1/a));
hist(y,60);
均值=6.7896方差=0.4212
2.8对数-正态分布
对数一正态分布模型可以用来表示高分辨率雷达在观察角小于
时,观察到的海浪杂波,在低观察角时观察到的地面杂波也可用对数一正态分布模型,这类杂波通常是形状不规则的大反射体,例如远洋舰船,较大的空间飞行器,或者SAR雷达观察到的城市等等。
其概率密度函数是:
均值
,方差
其产生源程序及直方图如下:
%%%%%%%%%%%%%%%%%
a=0.4;
b=0.4;
v=sqrt(b)*u+a;
%%%%%%%%%%%%%%%%%%
L=exp(v);
hist(L,100);
mean(L)
var(L)
均值=1.8499方差=2.2399
2.9Swerling分布
雷达系统中两次回波幅度之差服从Swerling1型。
可表示雷达反射回波功率或截面积或信噪比。
产生源程序(部分)如下:
rp=10;
lamade1=1/rp;
hist(y,100);
其直方图如下:
均值=10.1319方差=102.9242
2.10.t分布随机数
调用MATLAB函数trnd()产生服从t分布的随机数
均值=0.0506方差=7.7854
程序如下:
%产生满足t分布的随机数
clear;
clc;
y=trnd(2,1,2048);
%自由度为2
hist(y,50);
mean=mean(y);
var=var(y);
设计报告二产生两种相关随机变量
一、根据课本123页,产生相关随机序列的步骤:
(1)对功率谱进行采样,得到序列
(2)产生独立的
区间均匀分布的随机相位矢量序列
,其总体均匀功率谱等于1,即
=1。
(3)然后,给每个随机相位矢量乘以比例系数,得
(4)最后取逆离散傅立叶变换得到相关随机序列
1
二、两种相关随机序列
(1)功率谱密度为均匀分布随机变量:
clear
fs=512;
%设频率为512
x=rand(1,fs);
b=7;
z=(b-a)*x+a;
Sz=fft(z,fs);
N=length(Sz);
Pz=abs(Sz.^2)/N;
fs1=512;
x1=rand(1,fs);
c=2*pi;
d=0;
z=(c-d)*x1+d;
xn=sqrt(Pz).*z;
xk=ifft(xn,fs);
e=mean(xk);
%均值
d=std(xk)^2;
%方差
Pz1=abs(xk.^2)/N;
0.05:
100;
subplot(2,1,1);
plot((0:
511),xk);
axis([0,100,0,2]);
title('
均匀分布功率谱采样产生的相关随机变量'
);
subplot(2,1,2);
511),10*log10(Pz1));
相关随机变量的功率谱'
(2)功率谱密度为正态分布随机变量:
y=rand(1,fs);
z=sqrt(-2*log10(x)).*cos(2*pi*y);
x1=rand(1,fs1);
%产生0-2pai的均匀随机变量
xk=ifft(xn,fs1);
n=--256:
255;
plot((-256:
255),xk);
axis([-256,255,-0.5,0.5]);
正态分布功率谱采样产生的相关随机变量'
)
设计报告三雷达系统仿真
一.概述:
从一般意义上讲,系统仿真可以理解为在对一个已经存在或尚不存在但正在开发的系统进行研究的过程中,为了了解系统的在特性,设计构造即能反映系统特征又能符合系统实验要求的系统模型,并在该系统的模型上进行实验,以达到了解或设计系统的目的。
系统仿真本质上由三个要素构成,即系统,系统模型和实验。
系统是问题的本源,是分析问题的目的,实验是解决问题达到目的的手段,而系统模型则是连接系统和实验之间的桥梁。
二.任务要求及其各部分模型及仿真
仿真面向整个雷达系统,包含:
信号的混频、相干检波、脉冲压缩、MTI、相干积累、恒虚警等。
仿真雷达系统的原理框图如下:
此作业中从中放开始仿真,经过相干检波,然后在进行1/4抽取,之后在进行脉冲压缩(频域),取模后进行固定对消,然后进行相参积累,最后进行恒虚警处理。
(1)参数要求:
发射信号:
线性调频信号(17个周期):
脉宽100us带宽5MHz
周期1ms
回波:
目标回波:
杂波回波噪声
一个距离单元Rt=20Kmfd=250Hz
杂波回波:
四个距离单元Rc=25-25.12Km
噪声:
高斯白噪声
中放:
中频频率:
30MHz带宽:
10MHz
采样频率:
40MHz
相干检波:
参考信号频率30MHzLPF6MHz
1/4抽取
脉冲压缩:
频域实现
固定对消:
一次MTI
相参积累:
16点FFT
CFAR:
选大单元平均前16后16
(2)脉冲压缩信号的选取
线性调频信号是研究最多和应用最广的脉冲压缩信号,它容易产生和处理,多普勒频移容限大,可得到的脉压比围大。
发射的线性调频信号表示为:
(3)产生回波信号:
回波信号=目标回波+杂波回波+噪声
1产生目标回波(程序如下):
fd=250;
%多普勒频移
Rt=20*10^3;
%目标距离(一个距离单元)
Delay_s=fix(2*Rt*fs/c);
%目标回波延时点数
fori=1:
n
S_send(i,:
)=a.*exp(j*2*pi*((0+K*t/2).*t));
temp1(1,:
)=a.*exp(j*2*pi*((f01+K*t/2).*t+fd*i*T));
hbxh_s(i,Delay_s+1:
Num_T)=temp1(1,1:
Num_T-Delay_s);
hbxh_z(i,:
)=hbxh_z(1,:
end;
2产生杂波回波(程序如下):
Rc=25*10^3;
%杂波距离(4个距离单元)
Delay_c=fix(2*Rc*fs/c);
%杂波延时
Delay_c1=0.2*10^(-6)*fs;
temp2(1,:
)=a.*exp(j*2*pi*((f01+K*t/2).*t));
temp3(1,Delay_c+1+(0)*Delay_c1:
Num_T)=temp2(1,1:
Num_T-Delay_c-(0)*Delay_c1);
temp3(2,Delay_c+1+
(1)*Delay_c1:
Num_T-Delay_c-
(1)*Delay_c1);
temp3(3,Delay_c+1+
(2)*Delay_c1:
Num_T-Delay_c-
(2)*Delay_c1);
temp3(4,Delay_c+1+(3)*Delay_c1:
Num_T-Delay_c-(3)*Delay_c1);
hbxh_z(1,:
)=temp3(1,:
)=hbxh_z(1,:
)+temp3(2,:
)+temp3(3,:
)+temp3(4,:
3产生高斯白噪声(程序如下):
hbxh_n=0.2*randn(n,Num_T);
(4)中放:
中放模型为:
(相当于通过一个有增益的带通滤波器)
程序如下:
Mhz=1*10^6;
fsx=40*Mhz;
rfgain=5;
%增溢
fcuts=[3*Mhz,5*Mhz,15*Mhz,17*Mhz];
%带通滤波器f0=30MHzB=10MHz
mags=[0,1,0];
devs=[0.0080.050.008];
[m,Wn,beta,ftype]=kaiserord(fcuts,mags,devs,fsx);
hh=fir1(m,Wn,ftype,kaiser(m+1,beta),'
noscale'
hh=rfgain*hh;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
f_hh=fft(hh,5*length(hh));
figure(3),plot([0:
length(f_hh)-1]/length(f_hh),abs(f_hh));
title('
带通滤波器的频域图'
fork=1:
y_zf_out(k,:
)=conv(hbxh(k,:
),hh);
y_zf_out=y_zf_out(:
25:
40024);
figure(4);
subplot(211);
plot(real(y_zf_out(1,:
)));
通过带通滤波器(中放)后的信号'
subplot(212);
plot(abs(fft(y_zf_out(1,:
))));
通过带通滤波器(中放)后的信号(频域)'
(5)相干检波:
输入信号与一个本地振荡信号进行乘法运算,然后通过低通滤波器,再进行1/4抽取。
其模型为:
源程序如下:
(部分)
nn=[0:
length(y_zf_out)-1];
forl=1:
i1(l,:
)=y_zf_out(l,:
).*cos(2.*pi.*f0.*nn*1/fs);
q1(l,:
).*sin(2.*pi.*f0.*nn*1/fs);
end
f_i1=fft(i1(1,:
));
f_q1=fft(q1(1,:
figure(5),plot([0:
length(f_i1)-1/length(f_i1)],abs(f_i1)),title('
低通滤波器之前的i信号(频域)'
figure(6),plot([0:
length(f_i1)-1/length(f_i1)],abs(f_q1)),title('
低通滤波器之前的q信号(频域)'
aa=[1100];
desired_faa=[06-16+120]*10^6;
%MHz
faa=desired_faa/fs*2;
%频率值
low_filter=remez(m,faa,aa);
%n为阶数
figure(7),plot(low_filter),title('
低通滤波器的时域图'
f_filterl=fft(low_filter,512);
figure(8),plot([0:
511]/512,abs(f_filterl));
低通滤波器的频域图'
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
forkk=1:
i2(kk,:
)=conv(i1(kk,:
),low_filter);
q2(kk,:
)=conv(q1(kk,:
i2=i2(:
q2=q2(:
figure(9),subplot(2,1,1),plot(real(i2(1,:
))),title('
相干检波后输出的i路信号'
subplot(2,1,2),plot(real(q2(1,:
相干检波后输出的q路信号'
%时域图
f_i2=fft(i2(1,:
f_q2=fft(q2(1,:
figure(10),subplot(2,1,1),
plot(abs(f_i2)),title('
i路信号(频域)'
subplot(2,1,2),plot(abs(f_q2));
q路信号(频域)'
仿真结果:
低通滤波前的I路信号:
低通滤波前的Q路信号:
(6)脉冲压缩:
线性调频是最早和发展最充分的脉冲压缩技术。
它通过对雷达的载波频率进行调制以增加雷达的发射带宽并在接收时实现脉冲压缩。
可以在时域也可以在频域进行,当雷达接收机中的脉冲压缩滤波器与发射波形相匹配时,此时的相关性最强将输出幅度的峰值,可以明显改善信噪比。
此次仿真在频域进行:
其模型为:
Y(n)
加窗(汉明窗)
部分源程序如下:
y_my=i_cq+j*q_cq;
%将i和q成一路复信号
pp_filter1=fliplr((S_send));
%匹配滤波器(时域)
forcount=1:
f_y_my(count,:
)=fft(y_my(count,:
%将信号变到频域
pp_filter=pp_filter1(1:
4:
end);
%对匹配滤波器进行抽取
f_pp_filter=fft(pp_filter);
my_out(count,:
)=f_y_my(count,:
).*f_pp_filter.*hamming(10000).'
;
my_out1(count,:
)=ifft(my_out(count,:
figure(13);
plot(abs(my_out1(1,:
脉压后的信号(时域)'
(7)固定对消:
对于动目标显示雷达来讲,它应将杂波抑制,而让各种速度的目标回波信号通过。
实际能做到的杂波抑制滤波器,只能使滤波特性的凹口宽度基本和杂波的梳状谱的宽度相当。
常用的相消器有一次相消器、二次相消器、多次相消梳状滤波器和抑制运动杂波滤波器等。
本仿真采用的是一次相消器,一次对消器的原理及频域特性如下图示:
foraa=1:
n-1(对消后减少一个周期)
mti_out1(aa,:
)=my_out1(aa+1,:
)-my_out1(aa,:
end
figure(14);
mti_out2=abs(mti_out1(:
:
%对消后取模
plot(mti_out2(1,:
一