matlab中仿真自适应信号处理LMSNewton算法

p=0.01;N=16;u=0.05;%p表示随机信号的功率，N表示一个周期的采样点数，u表示权调整的步进

%****************************画出投影在平面的椭圆********************************%
syms w0 w1 
 y=（0.5+p）*（w0.^2+w1.^2）+w0*w1*cos（2*pi/N）+2*w1*sin（2*pi/N）+2;
 ezplot（（0.5+p）*（w0.^2+w1.^2）+w0*w1*cos（2*pi/N）+2*w1*sin（2*pi/N）+2-0.49）
 hold on
 ezplot（（0.5+p）*（w0^2+w1^2）+w0*w1*cos（2*pi/N）+2*w1*sin（2*pi/N）+2-0.78）
 hold on
 ezplot（（0.5+p）*（w0^2+w1^2）+w0*w1*cos（2*pi/N）+2*w1*sin（2*pi/N）+2-1.9[-8 20 -20 6]）
 hold on
 ezplot（（0.5+p）*（w0^2+w1^2）+w0*w1*cos（2*pi/N）+2*w1*sin（2*pi/N）+2-6.3[-8 20 -20 6]）
 hold on
 
 %****************************定义W的初始值及R矩阵********************************%
 w=[6 -9 zeros（1200）]‘; 
 R=0.5*[1