线性调频信号回波处理程序.doc

%程序：线性调频信号回波处理 %LFMEchoProc.m clc; %清楚命令窗口信息 clear; %清除工作区的变量 close all; %清除绘图 %****************参数设置*************% fr=3e9; %射频 f0=10e6; %中频 pw=256e-6; %脉宽 bw=2e6; %带宽 pri=1e-3; %脉冲重复周期 fs=70e6; %采样率 R=100e3; %距离 v=0; %速度 Pt=10*log10(1e6); %功率 lamda=10*log10(0.1);%波长 G=30; %天线增益dB L=6; %损耗dB Fn=3; %噪声系数dB RCS=10*log10(2); %RCS %***************信号参数计算*************% Pr=Pt+2*G+2*lamda+RCS-3*10*log10(4*pi)-4*10*log10(R)-L; %回波功率计算 A=sqrt(10^(Pr/10)); %回波信号幅度 tou=2*R/3e8; %收发时间 fd=2*v*fr/3e8; %回波多普勒频率 k=bw/pw; %调频带宽频率变化率 %************************发射信号产生*******************% t=(0:1/fs:pri).; %时间 sampleTbegin=ceil(0*fs+1); %采样时间起始点 sampleTEnd=ceil(pw*fs+1); %采样时间结束点 sTrans=zeros(length(t),1); %产生发射信号 sTrans(sampleTbegin:sampleTEnd,1)=... sin(2*pi*(f0*t(sampleTbegin:sampleTEnd,1)+k*t(sampleTbegin:sampleTEnd,1).^2/2)); %图像展示 figure; plot(t*1e6,sTrans); %发射信号时域波形 xlabel(时间/us); ylabel(幅度); freTrans=fftshift(fft(sTrans)); %发射信号频域波形 f=linspace(-fs/2,fs/2,length(t)); figure; plot(f/1e6,abs(freTrans)); xlabel(频率/MHz); ylabel(幅度); %***********************回波信号产生*******************% beginsample=ceil(tou*fs); %采样起始点 endsample=ceil((tou+pw)*fs); %采样结束点 sEcho=zeros(length(t),1); %产生回波信号 sEcho(beginsample:endsample,1)=... A*sin(2*pi*((f0+fd)*(t(beginsample:endsample,1)-tou)+k*(t(beginsample:endsample,1)-tou).^2/2)); figure; pl