xwjb.asv

有关雷达系统的仿真

%
%相位检波；
%

function [Y1,Y2]=xwjb(x)
N=length(x);
f0=10e6;                %中频 10M
fs=100e6;               %采样频率100M
Ts=1/fs;
T=0.00003;              %脉冲重复周期30微秒
t=0:Ts:T-Ts;
               

x1=x2=x;
x1=x1.*cos(2*pi*f0*t);
x2=x2.*sin(2*pi*f0*t);

n=50;
Wn=0.2; 
b=fir1(n,Wn);           %低通滤波器系数
x1=[x1,zeros(1,25)];
x2=[x2,zeros(1,25)];
Y1=filter(b,1,x1);
Y2=filter(b,1,x2);
Y1=Y1(26:length(Y1));
Y2=Y2(26:length(Y2));