matlab上机实验

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);fori=1:100%滤波器的阶数设置为100forj=1:100mrxx(i,j)=rxx(1000-i+j);00为原始信号的长度endendxd=cos(sita