自适应滤波matlab代码

clear;
R = [1 0.4045;0.4045 1];
la = max（eig（R））;
p = [0;0.2939];
d = 0.5;
w0 =  R\p;
syms w1 w2 s
w = [w1;w2];
s = d-2*w.‘*p+w.‘*R*w;
expand（s）
w1d = -2:0.1:2;
w2d = -2:0.1:2;
sd = zeros（length（w1d）length（w2d））;
for k = 1:length（w2d）
    sd（:k） = subs（s{w1w2}{w1dw2d（k）}）;
end
figure;contour（w1dw2dsd‘10）;hold on 
plot（w0（1）w0（2）‘*‘）;xlabel（‘w1‘）;ylabel（‘w2‘）;

u = 0.1/la;
N = 50;
I = ey