【求助】那位高手帮帮忙？

[feifei]

clc;
clear;
T=10e-6;
B=30e6;
K=B/T;
fs=10*B;
N=T*fs;
A0=1;
R0=100000;
f0=3e9;
v=200;
Kt=2;
c=3.0e8;
%发射信号
t=-T/2:1/fs:T/2-1/fs;
st=A0.*exp(j*pi*K*t.^2+j*2*pi*f0*t);
st_fft=fft(st,1024);
st_fft=fftshift(st_fft);
st_fft_result=abs(st_fft);
st_fft_result=st_fft_result/max(st_fft_result);
St=20*log10(st_fft_result+1e-6);
figure(1);
plot([1:1024]/1024*fs,St);
title('发射信号');
grid;
%回波信号
for i=1:N
t=-T/2+i/fs;
tao(i)=2*(R0-v*t)/c;
sr(i)=Kt*A0*exp(j*pi*K*(t-tao(i)).^2+j*2*pi*f0*(t-tao(i)));
end
sr_fft=fft(sr,1024);
sr_fft_result=fftshift(sr_fft);
sr_fft_result=abs(sr_fft_result);
sr_fft_result=sr_fft_result/max(sr_fft_result);
Sr=20*log10(sr_fft_result+1e-6);
figure(2);
plot([1:1024]/1024*fs,Sr);
title('回驳信号');
grid;
%上扫频段差拍信号
ss=sr.*st;
ss_fft=fft(ss,1024);
ss_fft_result=fftshift(ss_fft);
ss_fft_result=abs(ss_fft_result);
ss_fft_result=ss_fft_result/max(ss_fft_result);
Ss=20*log10(ss_fft_result+1e-6);
figure(3);
plot([1:1024]/N*fs,Ss,'r.');
title('上扫频差拍信号');
xlabel('时间\times\itB(s)');
ylabel('幅度');
grid;
%当目标固定时
t=-T/2:1/fs:T/2-1/fs;
sm=Kt*A0*exp(j*pi*K*(t-2*R0/c).^2+j*pi*2*f0*t);
sm_fft=fft(sm,1024);
sm_fft_result=fftshift(sm_fft);
sm_fft_result=abs(sm_fft_result);
sm_fft_result=sm_fft_result/max(sm_fft_result);
Sm=20*log10(sm_fft_result+1e-6);
figure(4);
plot([1:1024]/1024*fs,Sm,[1:1024]/1024*fs,Ss,'r.');
title('上扫频差拍信号与固定目标回波信号的距离差值图');
legend('\fontsize{9}固定目标回波信号','上扫频差拍信号');
xlabel('\fontsize{9}时间\times\itB(s)');
ylabel('\fontsize{9}幅度');
grid;


t=-T/2:1/fs:T/2-1/fs;
sm1=Kt*A0*exp(-j*pi*K*(t-2*R0/c).^2+j*pi*2*f0*t);
sm1_fft=fft(sm1,1024);
sm1_fft_result=fftshift(sm1_fft);
sm1_fft_result=abs(sm1_fft_result);
sm_1fft_result=sm1_fft_result/max(sm1_fft_result);
Sm1=20*log10(sm1_fft_result+1e-6);
%同上，下扫频信号
t=-T/2:1/fs:T/2-1/fs;
st1=A0.*exp(-j*pi*K*t.^2+j*2*pi*f0*t);
for i=1:N
t=-T/2+i/fs;
tao(i)=2*(R0-v*t)/c;
sr1(i)=Kt*A0*exp(-j*pi*K*(t-tao(i)).^2+j*2*pi*f0*(t-tao(i)));
end
ss1=sr1.*st1;
ss1_fft=fft(ss1,1024);
ss1_fft_result=fftshift(ss1_fft);
ss1_fft_result=abs(ss1_fft_result);
ss1_fft_result=ss1_fft_result/max(ss1_fft_result);
Ss1=20*log10(ss1_fft_result+1e-6);
figure(5);
plot([1:1024]/1024*fs,Sm1,[1:1024]/1024*fs,Ss1,'gr.');
title('下扫频差拍信号与固定目标回波信号的距离差值图');
legend('\fontsize{9}固定目标回波信号','下扫频差拍信号');
xlabel('\fontsize{9}时间\times\itB(s)');
ylabel('\fontsize{9}幅度');
grid;

这个程序结果应该是上下扫频差拍信号对称
可结果不对
那位高手帮忙看看