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