lattice_ser.rar_格型 SER

clear; clc; N=16; u=0.02; %u=0.1; k=200; alpha=0.9; pha=0.01; p0=0.5+pha; p1=0.5*cos(2*pi/N); p2=0.5*cos(4*pi/N); k0=[-2:0.1:0.5]; k1=[-2:0.01:4]; [k0,k1]=meshgrid(k0,k1); z=(k1.^2+1).*((k0.^2+1)*p0+2*k0*p1)+2*k1.*(p0*k0.^2+2*k0*p1+p2); %求性能表面 figure; contour(k0,k1,z,[0.048 0.060 0.095 0.048 0.060 0.095 0.172 0.600]);%画性能表面的等高线 rk=randn(1,k)*sqrt(pha); k0=zeros(1,k); k1=zeros(1,k); for j=1:k x(j)=sin(2*pi*j/N)+rk(j); s0(j)=x(j); s_0(j)=x(j); end s1k=zeros(1,k); %P0(1)=p0.^2-p1.^2; P0(1)=p0; P1(1)=p0.^2-p1.^2; for i=2:k %x(i)=sin(2*pi*i/N)+rk(i) ; s1k(i)=s0(i)+k0(i)*s_0(i-1); s_1(i)=k0(i)*s0(i)+s_0(i-1); P0(i)=((1-alpha)*s0(i).^2+alpha*(1-alpha.^i)*P0(i-1))/(1-alpha.^(k+1)); k0(i+1)=k0(i)-u* s1k(i)*s_0(i-1)/P0(i); s2k(i)=s1k(i)+k1(i)*s_1(i-1); P1(i)=((1-alpha)*s1k(i).^2+alpha*(1-alpha.^i)*P1(i-1))/(1-alpha.^(k+1)); k1(i+1)=k1(i)-u*s2k(i)*s_1(i-1)/P1(i); end x1=-0.906; y1=0.708; hold on; plot(k0,k1,'r'); %画迭代时权值的变化 xlabel('weight，k0');ylabel('weight，k1'); hold on;plot(x1,y1,'*'); %标注最佳权值的位置