雷达系统仿真matlab代码

1、= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = %该程序完成16个脉冲信号的【脉压、动目标显示/动目标检测(MTI/MTD)】% =% 程序中根据每个学生学号的末尾三位(依次为XYZ)来决定仿真参数，034% 目标距离 为3000 8025 9000+(Y*10+Z)*200 8025，4个目标% 目标速度 为50 0 (Y*10+X+Z)*6 100%=% closeall;%关闭所有图形clearall;%清除所有变量clc;%雷达参数%= %C=3.0e8;%光速(m/s)RF=3.140e9/2;%雷达射频1.

2、57GHzLambda二C/RF; %雷达工作波长PulseNumber=16;%回波脉冲数Ban dWidth=2.0e6;% 发射信号带宽带宽B=1/ T,T是脉冲宽度TimeWidth=42.0e-6;%发射信号时宽PRT=240e-6;%雷达发射脉冲重复周期(s),240us对应1/2*240*300=36000米最大无模糊距离PRF=1/PRT;Fs=2.0e6;%采样频率NoisePower=-12;%(dB);%噪声功率(目标为0dB )%-%SampleNumber二fix(Fs*PRT);%计算一个脉冲周期的采样点数 480 ；TotalNumber二SampleNumber

3、*PulseNumber;%总的采样点数 480*16=；BlindNumber二fix(Fs*TimeWidth);%计算一个脉冲周期的盲区-遮挡样点数；%=%目标参数%TargetNumber=4;%目标个数SigPower (1:TargetNumber)=1 1 1 0.25;%目标功率，无量纲TargetDista nee(1:TargetNumber)=3000 8025 15800 8025;% 目标距离,单位 m距离参数为3000 8025 9000+(Y*10+Z)*200 8025DelayNumber (1:TargetNumber)=fix(Fs*2*TargetDis

4、ta nce(1:TargetNumber)/C); %把目标距离换算成 采样点(距离门)fix函数向0靠拢取整Targetvelocity (1:TargetNumber)=50 0 204 100;% 目标径向速度 单位 m/s速度参数为50 0 (Y*10+X+Z)*6 100TargetFd (1:TargetNumber)=2*TargetVelocity(1:TargetNumber)/Lambda;n/.%计算目标多卜勒频移2v/入%;0/%号%-%产生线性调频信%number二fix(Fs*TimeWidth);%回波的采样点数二脉压系数长度二暂态点数目+1if rem(num

5、ber,2)=0%rem 求余nu mber 二nu mber+1;end%把number变为偶数for i=-fix (nu mber/2):fix (nu mber/2)-1Chirp(i+fix( nu mber/2)+1)=exp(j*(pi*(Ba ndWidth/TimeWidth)*(i/Fs)A2);%exp(j*fi)* ，产生复数矩阵 Chirpendcoeff=conj(fliplr(Chirp);%把Chirp矩阵翻转并把复数共轭，产生脉压系数figure(1);%脉压系数的实部plot(real(Chirp);axis(0 90 -1.5 1.5);title(脉压系

6、数实部);% 产生目标回波串产生前3个目标的回波串SignalAII二zeros(1,TotalNumber);%所有脉冲的信号，先填 0for k=1:TargetNumber-1%依次产生各个目标SignalTemp二zeros(1,SampleNumber);% 一个 PRTSig nalTemp(DelayNumber(k)+1:DelayNumber(k)+nu mber)二sqrt(SigPower(k)*Chirp;%一个脉冲的1个目标(未加多普勒速度)(DelayNumber(k)+1):(DelayNumber(k)+number)Sig nal二zeros(1,TotalN

7、umber);for i=1:PulseNumber % 16 个回波脉冲Sig nal(i-1)*SampleNumbe r+1:i*SampleNumber)=Sig nalTemp;%每个目标把 16 个Sig nalTemp 排在一起endFreqMove二exp(j*2*pi*TargetFd(k)*(0:TotalNumber-1)/Fs);%目标的多普勒速度*时间二目标的多普勒相移Sig nal二Sig nal.*FreqMove;Sig nalAII二Sig nalAII+Sig nal; end% %fi=pi/3;%加上多普勒速度后的16个脉冲1个目标%加上多普勒速度后的1

8、6个脉冲4个目标产生第4个目标的回波串%Sig nalTemp二zeros(1,SampleNumber);% 一个脉冲Sig nalTemp(DelayNumber(4)+1:DelayNumber(4)+nu mber)二sqrt(SigP ower(4)*exp(j*fi)*Chirp;%一个脉冲的1个目标(未加多普勒速度)Sig nal二zeros(1,TotalNumber);for i=1:PulseNumberSig nal(i-1)*SampleNumber+1:i*SampleNumber)二Sig nalTemp; endFreqMove=exp(j*2*pi*Target

9、Fd(4)*(0:TotalNumber-1)/Fs);%目标的多普勒速度*时间二目标的多普勒相移Sig nal二Sig nal.*FreqMove;Sig nalAII二Sig nalAII+Sig nal;figure(2);subpIot(2,1,1);pIot(reaI(SignalAII),r- );title( 目标信号的实部 );grid on;zoom on;subpIot(2,1,2);pIot(imag(SignalAII);title(目标信号的虚部);grid on ;zoom on ; %=%产生系统噪声信号%SystemNoise=normr nd(0,10A(No

10、isePower/10),1,TotalNumber)+j* normrnd(0,10A(NoisePower/10),1,TotalNumber);%均值为0,标准差为10A(NoisePower/10)的噪声%=%总的回波信号%Echo二Sig nalAII+SystemNoise;% +SeaClutte r+TerraClutter，加噪声之后的回波for i=1:PulseNumber%在接收机闭锁期，接收的回波为0Echo(i-1)*SampleNumbe r+1:(i-1)*SampleNumber+number)=0;%发射时接收为 0endfigure(3);%加噪声之后的总

11、回波信号subplot(2,1,1);plot(real(Echo),r- );title( 总回波信号的实部，闭锁期为0);subplot(2,1,2);plot(imag(Echo);title(总回波信号的虚部，闭锁期为0);%=时域脉压=%pc_time0=conv(Echo,coeff);%pc_time0 为 Echo 和coeff 的卷积pc_time1=pc_time0 (nu mber:TotalNumber+nu mber-1);%去掉暂态点nu mber-1 个r- );title( 时域脉压结果的时域脉压结果的幅度，无频域脉压figure(4);%时域脉压结果的幅度su

12、bplot(2,1,1);plot(abs(pc_time0),幅度，有暂态点)；%pc_time0的模的曲线 subplot(2,1,2);plot(abs(pc_time1);title( 暂态点);%pc_time1 的模的曲线% = =%Echo_fft=fft(Echo,8192);%理应进行 TotalNumber+number-1点 FFT，但为了提高运算速度，进行了 8192点的FFTcoeff_fft=fft(coeff,8192);pc_fft=Echo_fft.*coeff_fft;pc_freq0=ifft(pc_fft);figure subplot(2,1,1);p

13、lot(abs(pc_freq0(1:TotalNumber+nu mber-1);title(频域脉压结果的幅度，有前暂态点)；subplot(2,1,2);plot(abs(pc_time0(1:TotalNumbe r+nu mber-1)- pc_freq0(1:TotalNumber+number-1),r );title(时域和频域脉压的差别);%去掉暂态点pc_freq1=pc_freq 0(nu mber:TotalNumber+nu mber-1);number-1 个，后填充点若干(8192-number+1-TotalNumber)% = 按照脉冲号、距离门号重排数据=%

14、for i=1:PulseNumberpc(i,1:SampleNumber)=pc_freq1(i-1)*SampleNumber+1:i*SampleNumber);%每个 PRT为一行，每行 480 个采样 点的数据endfigureplot(abs(pc(1,:);title(频域脉压结果的幅度，没有暂态点)；% =MTI (动目标显示)，对消静止目标和低速目标-可抑制 杂波=%for i=1:PulseNumber-1%滑动对消，少了一个脉冲mti(i,:)=pc(i+1,:)-pc(i,:);endfigure(7);mesh(abs(mti);title(MTI result )

15、;% =MTD(动目标检测)，区分不同速度的目标，有测速作用=% mtd=zeros(PulseNumber,SampleNumber);for i=1:SampleNumberbuff(1:PulseNumber)=pc(1:PulseNumber,i);buff_fft=fft(buff);mtd(1:PulseNumber,i)=buff_fft(1:PulseNumber);endfigure(8);mesh(abs(mtd);title(MTD result );%=虚实矩阵转换=% coeff_fft_c=zeros(1,2*8192);for 社1192coeff_fft_c(2*i-1)=real(coeff_fft(i);coeff_fft_c(2*i)=imag(coeff_fft(i);endecho_c=zeros(1,2*TotalNumber);for i=1:TotalNumberecho_c(2*i-1)=real(Echo(i);echo_c(2*i)=imag(Echo(i);end%=以下是为DSP程序提供回波数据、脉压系数=% fo=fope n(F:mystudyVisual_DSP_