《电子科大毕设》由会员分享，可在线阅读，更多相关《电子科大毕设（8页珍藏版）》请在金锄头文库上搜索。
1、第6 章 软件测试6.1 接收机处理模块6.1.1 程序名：IQchannal函数 IQchannal 用于完成对目标回波进行下变频 IQ 两路正交分解工作。输 入参数为目标回波数据；输出参数为IQ两路信号数据。The Receive Wave after Down-Conversin 图 6.3 载波频率为 25MHz ，采样频率为 20MHz 的 13 位巴克码信号下变频波形 程序编码 IQchannal图The Receive Wave after Down apn-E.w形function Srt=IQchannal(Srt,Trec,f0) global FsN=length(Srt
2、);% % 混频% t1=linspace(Trec(1),Trec(2),N); t1=Trec(1):(Trec(2)-Trec(1)/N:(Trec(2)-(Trec(2)-Trec(1)/N);for i=1:NXbi(i)=Srt(i)*sin(2*pi*f0*t1(i);Xbq(i)=Srt(i)*(cos(2*pi*f0*t1(i);end% IQ相抽取if Fs20e6D=3;Rbi=zeros(1,D),Xbi;Rbq=zeros(1,D),Xbq;Dlength=floor(length(Xbi)/D);clear Xbi Xbqfor i=1:Dfor j=1:Dleng
3、thXbi(i,j)=Rbi(D*(j-1)+D+2-i); Xbq(i,j)=Rbq(D*(j-1)+D+2-i);endend%低通滤波器load ham_filter.mat ham_length=length(ham_h);for i=1:Dfor j=1:floor(ham_length/D)H(i,j)=ham_h(D*(j-1)+i);endendham_length=length(H(i,:);ham_n=floor(ham_length/2);ham_zeros=zeros(1,ham_n);for i=1:DWbi(i,:)=filter(H(i,:),1,0,Xbi(i,
4、:),ham_zeros);Zbi(i,:)=Wbi(i,ham_n+1:end);Wbq(i,:)=filter(H(i,:),1,0,Xbq(i,:),ham_zeros);Zbq(i,:)=Wbq(i,ham_n+1:end);endfor j=1:DlengthYbi(j)=sum(Zbi(:,j);Ybq(j)=sum(Zbq(:,j);endelse%低通滤波器load ham_filter.matham_length=length(ham_h);ham_n=floor(ham_length/2);ham_zeros=zeros(1,ham_n);Wbi=filter(ham_h,
5、1,0,Xbi,ham_zeros);Ybi=Wbi(ham_n+1:end);Wbq=filter(ham_h,1,0,Xbq,ham_zeros);Ybq=Wbq(ham_n+1:end);end%生成IQ相的复信号。Srt=complex(Ybi,Ybq);6.1.2 程序名：match_functionmatch_function 包括 match_lfm、match_dc 函数，用于完成对线性调频信号 和巴克编码信号的脉冲压缩(匹配滤波)工作。数据参数为目标回波下变频信号 数据；输出参数为目标回波脉压后数据10。(qp)llp匚一 100-10-2Q50-60,6.4The-.Wav
6、e aftePulse Compression带-0- (qp) Bpnl_ubmw(qp)-gn-FIAIThe Wave after Pul.se- Gcimpres.ion-40波器单位冲激响应15-2Time2.63W5t=linspace(O,T,NO);K=B/T;f1=f0+B/2;h=cos(2*pi*fl*t-pi*K*t42);h=IQchannal(h,0,T,f0);Sot=conv(Srt,h);程序编码 match_dcfunction Sot=match_dc(Srt)global B f0 T N0 %匹配滤波器单位冲激响应 h=send_dc(B,T,l,N0
7、,0,f0); h=IQchannal(h,0 T,f0);h=fliplr(conj(h);Sot=conv(Srt,h);6.1.3 程序名： Measure_informationMeasure_information 函数用于测量目标参数(距离，角度等)。输入参数为 回波下变频波形数据和回波脉压波形数据；输出参数为目标的距离，方位角，俯 仰角等。EjlllclLI 5一 PQwnE石戎eu苍芯Lnfunctiondist,Velocity,Azimuth_m,Elevation_m,track_flog=Measure_information(Trec,Add IQ,CF_in,Sot
8、_Add_IQ,Azimuthchannal,Elevationchannal,track_flog,T0,f0)global T N0 Ts K wavelength R Velocity_function N=length(CF_in);% % 接受波形实部t1=linspace(Trec(1),Trec(2),N);if track_flog=1% 距离检测T1,T0=min(abs(T0-t1);dist,T0=Rdetection(CF_in,t1,Ts,T0);T0=t1(T0);else% 自动搜索C=3e8;auto_slip_error,T0=search(CF_in,N0);dist=C*(t1(T0)/2;dist,T0=Rdetection(CF_in,t1,Ts,T0,N0);%进入跟踪T0=t1(T0);track_flog=1;end%方位-俯仰角测量V elocity_functionAzimuth_m,Elevation_m=angel_m(Azimuthchannal,Elevationchannal,Sot_Add_IQ, Trec,f0);%速度测量Velocity=Velocity_function(Add_IQ,T,Ts,t1,N0,R,wavelength);站长资讯 http:/ cjb13x
