1、做一仿真

2、是基于相位编码的雷达回波仿真

3、具体编码采用P4码、回波分析采用小波变换

4、下面给一个资料，关于相位编码(P4码)的波形参数设计以及一些现成的仿真图像

再给一个线性调频连续波回波差频信号处理的代码(别人的  我不太懂)：

%构造一个线性调频信号作为差频信号进行2D FFT处理

clear all;

clc;

%仿真参数===================================================================

B=30e3;%调频带宽

f0=10e6;%调频起始频率

T=12e-3;%调频周期

N2=128;%单个调频周期内的采样点数

t=0:T/N2:T-T/N2;%采样序列

R0=1e5;%目标距离

v=250;%目标速度

c=3e8;%光速c

N1=512;%相参积累周期数N1=512;

n=0:N1-1;

k=B/T;%调频速率

Rn=R0+v*T.*n;%Rn为一个256维的行向量

%构造多普勒频率=============================================================

fdv=2*f0*v/c; %fdv为速度多普勒频率

fRn=2*k.*Rn/c; %fRn为距离多普勒频率

fdvRn=4*k*v.*Rn/(c^2); % fdvRn为距离速度多普勒频率，

%生成信号===================================================================

phase0=4*pi*f0.*Rn/c-4*pi*k.*(Rn.*Rn)/(c^2);

phaseMat=zeros(N2,N1);

for count=1:N2

phaseMat(count,:)=phase0;

end

phasen=2*pi.*(t'*(fdv+fRn+fdvRn))+phaseMat;

si=cos(phasen);%差频信号

%第一维傅里叶变换===========================================================

window1=window(@kaiser,N2,2);

fftre1=zeros(N2,N1);

for countercol=1:N1

si(:,countercol)=si(:,countercol).*window1;

fftre1(:,countercol)=fft(si(:,countercol));

end

%figure(1);

%mesh(1:512,(1:64)*5e3,db(fftre1(1:N2/2,:)));axis xy;

%第二维傅里叶变换===========================================================

window2=window(@kaiser,N1);

fftre2=zeros(N2,N1);

fftre1w=fftre1;

%for countercol=1:N1

%   fftre1(:,countercol)=fftre1(:,countercol).*window2;

%end

for counterrow=1:N2

fftre1w(counterrow,:)=fftre1w(counterrow,:).*window2.';

fftre2(counterrow,:)=fft(fftre1w(counterrow,:));

end

%fftre2=abs(fft(fftre1')).';

figure(1);

%mesh(db(fftre2(:,1:N1/2)));axis xy;

mesh((1:N1/2)*2.4316,(1:N2/4)*5e3,db(fftre2(1:N2/4,1:N1/2)));axis xy;

%colormap gray;

title('2DFFT');

xlabel('Velocity--m/s');

ylabel('Distance--m');

5、说下我要的最终结果：是将材料中的速度图和距离图  放到一个三维图形中 。

还有一个要求就是代码要有详细注释。

6、奖励：加我QQ：515519872    重金交易至少150 决不食言

[本帖最后由 cwjy 于 2010-5-25 17:34 编辑]