【运动学】脉冲多普勒测距测速【含Matlab源码 1947期】.zip

clear all a=1; Bi=[ -0.000046451962915, -0.00075518384827 ,-0.000135308415921, 0.002841555023619 ,0.00084437987058 ,-0.007634010375115, -0.002839935656942 ,0.017363950908455 ,0.007712529489581 ,-0.038285191153445 ,-0.021562253256757 ,0.115631928738445 ,0.24000074917771 ,0.192654053061175 ,0.036063683235075 ,-0.045213871437468 ,-0.015108046955929 ,0.019192080579196 ,0.007479138501009 ,-0.008227322716728 ,-0.003467948663973 ,0.003044762311011 ,0.001333679741865 ,-0.000823931106908 ,-0.000362828687586 ,0.000096763709473]; Bq=[0.000096763709473, -0.000362828687586 ,-0.000823931106908 ,0.001333679741865 ,0.003044762311011, -0.003467948663973, -0.008227322716728 ,0.007479138501009 ,0.019192080579196 ,-0.015108046955929 ,-0.045213871437468 ,0.036063683235075 ,0.192654053061175 ,0.24000074917771 ,0.115631928738445 ,-0.021562253256757 ,-0.038285191153445 ,0.007712529489581 ,0.017363950908455 ,-0.002839935656942 ,-0.007634010375115 ,0.00084437987058 ,0.002841555023619 ,-0.000135308415921 ,-0.00075518384827 ,-0.000046451962915]; f0=60*10^6-390.625*40; fs=48*10^6; n=240*8192*4; signal=zeros(1,80*12*8192); %设置回波脉冲位置 for i=1:8192 for k=1:12 signal((i-1)*960+12*66+8+k)=cos(2*pi*f0*((i-1)*960+12*66+8+k)*(1/fs)); end end %高斯白噪声 signal=signal+20*rand(1,n); %奇偶抽取 num_i=1:2:(240*8192*4-1); num_q=2:2:(240*8192*4); signal_i_temp=signal(num_i); signal_q_temp=signal(num_q); %符号变换 k=1:(240*8192*2); signal_i=((-1).^(k+1).*signal_i_temp(k)); signal_q=((-1).^(k).*signal_q_temp(k)); %多相滤波 Yi=filter(Bi,a,signal_i); Yq=filter(Bq,a,signal_q); %修正延迟 Yi(1:12)=[]; Yi=[Yi,zeros(1,12)]; Yq(1:12)=[]; Yq=[Yq,zeros(1,12)]; %2倍抽取 kk=0:2:(240*8192*2-2); Yii=Yi(kk+1); Yqq=Yq(kk+1); i=1:(80*8192); Yii_1=Yii(3*i-2)+Yii(3*i-1)+Yii(3*i); A=reshape(Yii_1,80,8192); i=1:(80*8192); Yqq_1=Yqq(3*i-2)+Yqq(3*i-1)+Yqq(3*i); A1=reshape(Yqq_1,80,8192); sepctrum=zeros(80,128); %光谱 for p=1:80 B=A(p,:); C=reshape(B,128,64); D=C'; B1=A1(p,:); C1=reshape(B1,128,64); D1=C1'; E=D+j*D1; for i=1:64 sepctrum(p,:)=sepctrum(p,:)+abs(fft((E(i,:).*(hann(128)')))); end end figure(1); mesh(sepctrum); xlabel('速度对应频率点数N：1-128'); ylabel('距离对应距离门数n：1-80'); zlabel('视频信号幅度');