clear; clf; sita=0:pi/499.5:2*pi; %取 1000 个样本 xnoise=sqrt(0.03)*randn(1,1000);%产生 x 轴方向噪声 ynoise=sqrt(0.02)*randn(1,1000);%产生 y 轴方向噪声 x=cos(sita)+xnoise;%产生 x 轴方向观测信号 y=sin(sita)+ynoise;%产生 y 轴方向观测信号 %产生维纳滤波中 x 方向上观测信号的自相关矩阵 rxx=xcorr(x); for i=1:100 %滤波器的阶数设置为 100 for j=1:100 mrxx(i,j)=rxx(1000-i+j