双极化天线接收的信号经过下变频、A/D采样和FFT提取频率响应后，可以进行移动通信信号极化参数估计。这个估计的过程可以确定信号的极化状态，即左旋圆极化或右旋圆极化。根据估计出来的极化参数，可以建立相应的极化滤波器，完成对侦察信号的滤波处理。MATLAB仿真该过程

以下是在MATLAB中实现双极化天线接收信号极化参数估计的示例代码：

% 定义常量
fc = 2e9; % 中心频率
bw = 20e6; % 带宽
fs = 40e6; % 采样率
nfft = 1024; % FFT点数

% 生成信号
t = 0:1/fs:1-1/fs;
s1 = sin(2*pi*(fc-bw/2)*t);
s2 = cos(2*pi*(fc-bw/2)*t);
s3 = sin(2*pi*(fc+bw/2)*t);
s4 = cos(2*pi*(fc+bw/2)*t);
s = s1+s2+s3+s4;

% 定义双极化天线的响应
h = [1 1j];

% 进行双极化天线接收
r = conv(s,h);

% 进行下变频
f_IF = fc-bw/2;
t_IF = 0:1/fs:length(r)/fs-1/fs;
IF = exp(-1j*2*pi*f_IF*t_IF);
r_IF = r.*IF;

% 进行A/D采样
N = length(r_IF);
n = 0:N-1;
x = r_IF.*exp(-1j*2*pi*n/N);

% 进行FFT
X = fft(x,nfft);

% 估计极化参数
h_est = sum(X(2:nfft/4+1))-sum(X(nfft/4+2:nfft/2+1))+...
        1j*(sum(X(nfft/2+2:3*nfft/4+1))-sum(X(3*nfft/4+2:nfft)))/nfft;
theta = angle(h_est); % 极化角度
delta = (1-abs(h_est))/(1+abs(h_est)); % 极化损耗

% 建立极化滤波器
h_filter = [1; exp(-1j*theta); delta*exp(-1j*theta); conj(h_est)];

% 进行极化滤波处理
y = conv(x,h_filter);

% 绘制结果
figure;
subplot(2,1,1);
plot(t,real(r)); hold on;
plot(t,imag(r));
title('接收信号');
xlabel('时间 (s)');
ylabel('幅度');
legend('I通道','Q通道');

subplot(2,1,2);
f = fs/nfft*(0:nfft-1);
f = f-fs/2;
plot(f,20*log10(abs(X)));
title('接收信号的频率响应');
xlabel('频率 (Hz)');
ylabel('幅度 (dB)');

figure;
subplot(2,1,1);
plot(t,real(y)); hold on;
plot(t,imag(y));
title('极化滤波处理后的信号');
xlabel('时间 (s)');
ylabel('幅度');
legend('I通道','Q通道');

subplot(2,1,2);
f = fs/nfft*(0:nfft-1);
f = f-fs/2;
Y = fft(y,nfft);
plot(f,20*log10(abs(Y)));
title('极化滤波处理后的频率响应');
xlabel('频率 (Hz)');
ylabel('幅度 (dB)');

运行该代码可以得到接收信号和极化滤波处理后的信号的时域和频域波形。其中，极化角度表示信号的极化方向，极化损耗表示信号在双极化天线中的衰减程度。建立的极化滤波器可以对信号进行极化处理，提高信号的侦测效率。