【回波】雷达FSK 回波处理的matlab程序.rar

close all; clear ; clc %%%%%%%发射波形2-FSK%%%%%%% f1=3e3;%调制频率1 f2=1e4;%调制频率2 ph0=0;%初相位 T=0.2; f=15*f2;%左右采样频率 dt=1/f;%采样周期即间隔 t=0:dt:0.5-dt; %N=T/dt; fsk1=0.*t; fsk2=0.*t; i=0; % j=sqrt(-1); % fsk1=fsk1+1/2*(exp(j*(2*pi*f1*t+ph0))+exp(-j(2*pi*f1*t+ph0))).*rectpuls((t-T/2-2*i*T)/T); % fsk2=fsk2+1/2*(exp(j*(2*pi*f2*t+ph0))+exp(-j(2*pi*f2*t+ph0))).*rectpuls((t-3*T/2-2*i*T)/T); fsk1=fsk1+cos(2*pi*f1*t+ph0).*rectpuls((t-T/2-2*i*T)/T); fsk2=fsk2+cos(2*pi*f2*t+ph0).*rectpuls((t-3*T/2-2*i*T)/T); fsk=fsk1+fsk2; figure subplot(3,1,1);plot(t,fsk);title("发射波形2-FSK");xlabel('时间/s');ylabel('幅度'); axis([0 0.4 -4 4]);grid on; subplot(3,1,2);plot(t,fsk1);title("发射波形FSK1");xlabel('时间/s');ylabel('幅度'); axis([0 0.4 -4 4]);grid on; subplot(3,1,3);plot(t,fsk2);title("发射波形FSK2");xlabel('时间/s');ylabel('幅度'); axis([0 0.4 -4 4]);grid on; %%%%%%%理想回波波形%%%%%%% v=300;%目标速度为30 R0=ones(1,length(t))*10000;%初始距离为100 R=R0-v.*t;%即时距离 c=3e+05; td=2*R/c; fd1=2*v*f1/c; fd2=2*v*f2/c; rt1=0.*t; rt2=0.*t; vt1=0.*t; vt2=0.*t; i=0; rt1=rt1+cos(2*pi*(f1+fd1)*(t-td)+ph0).*rectpuls(((t-td)-T/2-2*i*T)/T); rt2=rt2+cos(2*pi*(f2+fd2)*(t-td)+ph0).*rectpuls(((t-td)-3*T/2-2*i*T)/T); rt=rt1+rt2; vt1=vt1+cos(2*pi*fd1*(t-td)+ph0).*rectpuls(((t-td)-T/2-2*i*T)/T); vt2=vt2+cos(2*pi*fd2*(t-td)+ph0).*rectpuls(((t-td)-T/2-2*i*T)/T); vt=vt1+vt2; figure subplot(3,1,1);plot(t,rt);title("回波波形2-FSK");xlabel('时间/s');ylabel('幅度'); axis([0 0.5 -4 4]);grid on; subplot(3,1,2);plot(t,rt1);title("回波波形FSK1");xlabel('时间/s');ylabel('幅度'); axis([0 0.5 -4 4]);grid on; subplot(3,1,3);plot(t,rt2);title("回波波形FSK2");xlabel('时间/s');ylabel('幅度'); axis([0 0.5 -4 4]);grid on; N=2^25; n=0:N-1; t=n/f; %时间序列 y1=fft(fsk1,N); %对信号进行快速Fourier变换 mag1=abs(y1); %求得Fourier变换后的振幅 f11=n*f/N; %进行对应的频率转换 figure plot(f11,mag1,'b'); %绘出随频率变化的振幅 xlabel('频率/Hz'); ylabel('振幅');title('图1：发射FSK1 FFT','color','b');grid on; y2=fft(rt1,N); mag2=abs(y2); % t=td+n/f; %时间序列 f22=n*f/N; %进行对应的频率转换 figure plot(f22,mag2,'r'); %绘出随频率变化的振幅 xlabel('频率/Hz'); ylabel('振幅');title('图2：回波FSK1 FFT','color','r');grid on; y3=fft(vt1,N); mag3=abs(y3); % t=td+n/f; %时间序列 f33=n*f/N; %进行对应的频率转换 figure plot(f33,mag3,'r'); %绘出随频率变化的振幅 xlabel('频率/Hz'); ylabel('振幅');title('图3：FSK1多普勒频谱','color','r');grid on; figure v2=v*ones(1,length(f33)); plot(v2,mag3,'r'); xlabel('速度m/s'); ylabel('振幅');title('图4：FSK1处理结果','color','r');grid on; % % i=0; % alvt1=0.*t; % alfd1=fd1*c/2/f1; % alvt1=alvt1+cos(2*pi*alfd1*(t-td)+ph0).*rectpuls(((t-td)-T/2-2*i*T)/T); % y4=fft(alvt1,N); % mag4=abs(y4); % t=td+n/f; %时间序列 % f33=n*f/N; %进行对应的频率转换 % figure % plot(f44,mag4,'b'); %绘出随频率变化的振幅 % xlabel('速度/m/s'); % ylabel('振幅');title('图4：速度频谱','color','r');grid on;