雷达系统仿真matlab代码

====================%%该程序完成个脉冲信号的【脉压、动目标显示/动目标检测16MTI/MTD]%====================%%程中依据每个学生学号的末尾三位依次为来确定仿真参数，5XYZ034%目标距离为个目标[300080259000+Y*10+Z*2008025],4%目标速度为[500Y*10+X+Z*6100]%关闭全部图形%清除全部变量%====================%close all;clear all;clc;%============%%雷达参数%%============%C=

3.0e8;m/s%雷达射频RF=

3.140e9/2;

1.57GHz雷达工作波长Lambda=C/RF;%%回波脉冲数PulseNumber=l6;%放射信号带宽带宽八，是脉冲宽度BandWidth=

2.0e6;B=1T%放射信号时宽TimeWidth=42,Oe-6;%雷达放射脉冲重复周期对应夫夫米最大无模PRT=240e-6;s,240us1/2240300=36000糊距离PRF=1/PRT;%采样频率Fs=

2.0e6;;带麋声功率目标为NoisePower=-12;%dB OdB%%;%计算一个脉冲周期的采样点数；SampleNumber=fix Fs*PRT480总的米样点数；To ta lNumber=SampleNumber*Puls eNumber;%480*16=计算一个脉冲周期的盲区-遮挡样点数；BlindNumber=f ixFs*TimeWidth;%%===================================================================================%%目标参数%%===================================================================================招目标个数㊀Targ tNumber=4;%目标功率,无量纲SigPower1:TargetNumber=[

1110.25];%;%目材:星巨离，单位星TargetDistance1:TargetNumber=

[30008025158008025]m巨离参数为[300080259000+Y*10+Z*2008025]DelayNumber1:TargetNumber=fixFs*2*TargetDistance1:TargetNumber/C;%把目标距离换算成采样点距离门函数向靠拢取整fix0目标径向速度单位速度参Targetvelocity1:TargetNumber=

[500204100];%m/s数为[500Y*10+X+Z*6100]%计算目标多TargetFd1:TargetNumber=2*Targetvelocity1:TargetNumber/Lambda;卜勒频移入2v/======================================================================%%产生线性调频信号%======================================================================%==============%夫回波的采样点数=脉压系数长度=暂态点数目㊀number=f ixFs TimWidth;%+1求余if remnumber,2~=0%remnumber=number+l;把变为偶数㊀end2numb rfor i=-fixnumber/2:fixnumber/2-1人Chirpi+fixnumber/2+1=expj*pi*BandWidth/TimeWidth*i/Fs2;%ex产生复数矩阵pj*fi*,Chirp end;%把矩阵翻转并把复数共粗，产生脉压系数coef f=conj fliplrChirp110⑴；％脉压系数的实部㊀figur1脉压系数实部1;plot realChirp;axis[090-

1.

51.5];title%----------------------------------------------产生目标回波串----------------------------------%%产生前个目标的回波串%3;%全部脉冲的信号，先填SignalAll=zeros1,TotalNumber0%依次产生各个目标for k=l:TargetNumber-1一个SignalTemp=zeros1,SampleNumber;%PRT SignalTempDelayNumberk+1DelayNumberk+number=sqrtSigPowerk*Ch irp;%一个脉冲的个目标未加多普勒速度1DelayNumberk+1:DelayNumberk+numberSignal=zeros1,TotalNumber;个回波脉冲for i=l:PulseNumber%16㊀%每个目标Signali-1*SampleNumber+1i*SampleNumber=SignalT mp;把个排在一起16SignalTempend目标的多普FreqMove=exp j*2*pi*TargetFd k*0:TotalNumber-1/Fs;%勒速度*时间=目标的多普勒相移__加上多普勒速度后的个脉冲个目标Signal=Signal.*FreqMove;%161加上多普勒速度后的个脉冲个目标SignalAll=SignalAl1+Signal;%164end%%产生第个目标的回波串%4fi=pi/3;一个脉冲SignalTemp=zeros1,SampleNumber;%SignalTempDelayNumber4+1DelayNumber4+number=sqrtSigPower4*ex pj一个脉冲的个目标未加多普勒速度*fi*Chirp;%1Signal=zeros1,TotalNumber;for i=l:PulseNumberSignali-1*SamplcNumbcr+l:i*SampleNumber-SignalTemp;end目标的多普勒FreqMove=exp j*2*pi*TargetFd4*0:TotalNumber-1/Fs;%速度大时间=目标的多普勒相移Signal=Signal.*FreqMove;SignalAll=SignalAll+Signal;figure2;目标信号的实部1;11subplot2,1,1;plot realSignalAll,r-;title*gridon;zoom on;目标信号的虚部f1subplot2,1,2;plot imagSignalAll;title;grid on;zoomon;======================================================================%==============%%产生系统噪声信号%%====================================================================================%A ASystemNoise=normrnd0,10NoisePower/10,1,TotalNumber+j*normrnd0,10NoisePower/10,1,TotalNumber;%均值为标准差为人的噪0,10NoisePower/10声%====================================================================================%%总的回波信号%======================================================================%==============%力噪声之后的回波Echo=SignalAll+SystemNoise;%+SeaClutter+TerraClutter,H for%在接收机闭锁期，接收的回波为i=l:PulseNumber0%放射时接收㊀Echoi-1*Sampl Number+1:i-1*SampleNumber+number=0;为加噪声之后的总回波信号1end figure3;%subplot2,1,1;plot realEchoz「总回波信号的实部，闭锁期为fr-;title01;subplot2,1,2;plot imagEcho;「总回波信号的虚部,闭锁期为寸域脉压1title0;===============================4%为㊀=================================pc_time0=conv Echo,coef f;%pc_tim0%和的卷积;%去掉暂Echo coef f pc_timel=pc_time0number:TotalNumber+number-1态点个;%同•域脉压结果的幅度㊀numb r-l figure4subplot2,1,1;plot abs「时域脉压结果的幅度，有暂态点的模的曲线1f1pc_time0,r-;title;%pc_time0时域脉压结果的幅度，无暂态点1subplot2,1,2;plot abspc_timel;title的模的曲线%===二==========================页域脉压1㊀;%pc_tim19=================================%;%理应进行点但为了提高运㊀Echo_fft=f ftEcho,8192TotalNumb r+number-l FFT,算速度，进行了点的8192FFT coeff_fft=fftcoeff,8192;pc_fft=Echo_fft,*coeff_fft;pc_freq0=ifftpc_fft;figure5;1㊀㊀subplot2,1,1;plot abspc_f reqO1:TotalNumb r+numb r-l;title频域脉压结果的幅度，有前暂态点D;subplot2,1,2;plotabspc_timeO1:TotalNumber+number-1-时域和频域脉压的差另卜;111pc_f reqO1:TotalNumber+number-1,r;title去掉暂态点个,pc_f reql=pc_f reqOnumber:TotalNumber+number-1;%number-1后填充点若干8192-number+l-TotalNumber%================依据脉冲号、距离门号重排数据=================================%for i=l:PulseNumberpci,1:SampleNumber=pc_freqli-每个为——行，每行、采样点的数据1*SampleNumber+l:i*SampleNumber;%PRT4804endfigure6;频域脉压结果的幅度,没有暂态点1;1plot abspc1,;title动目标显示，对消静止目标和低速目标--可抑制杂波%================MTI=================================%%滑动对消，少了一个脉冲for i=l:PulseNumber-1mti i,:=pc i+1,:-pc i,:;endfigure7;11meshabsmti;titleMTI result;动目标检测，区分不同速度的目标，有测速作用%================MTD=================================招mtd=zerosPulseNumber,SampleNumber;for i=l:SampleNumberbuff1:PulseNumber=pc1:PulseNumber,i;buff_fft=fftbuff;mtd1:PulseNumber,i=buff_fft1:PulseNumber;end1figure8;meshabsmtd;titleMTD result*;%=======================================虚实矩阵转换=======================================受coeff_fft_c=zeros1,2*8192;for i=l:8192coeff_fft_c2*i-l=realcoeff_ffti;coeff_fft_c2*i=imagcoeff_ffti;endecho_c=zeros1,2*TotalNumber;for i=l:TotalNumberecho_c2*i-l=realEchoi;echo_c2*i=imagEchoi;end===========================以下是为程序供应回波数据、脉压系数%DSP===============================%频域1f1%fo=fopen*F:\my study\Visual_DSP_test\test_l\coeff_fft_c.dat,wt;%脉压系数一—一—一%for i=l:2*81921%fprintffo,%f,\r\n,coeff_fft_ci;z%end%fclosefo;%次回波1%fo=fopen*F:\my study\Visual_DSP_test\test_l\echo_c.dat*,,wt;%16的一一一一%fori=l:2*TotalNumber1%fprintffo,*%f,\r\n echo_ci;%endz%fclosefo;。