selfitm.m

ｍａｔｌａｂ程序

function [w,f,v]=selfitm(t,n,k);
w=[0.5 0.5 0.5 0.5 0.5 0.5;0.5 0.5 0.5 0.5 0.5 0.5];
num=size(t);
for j=1:n
  for i=2:9,
    t1(i+1,:)=t(i,:).*w(1,:)+w(2,:).*t(i-1,:);
    e(i+1,:)=t(i+1,:)-t1(i+1,:);
    w(2,:)=w(2,:)+2*k*e(i+1,:).*t(i-1,:);
    w(1,:)=w(1,:)+2*k*e(i+1,:).*t(i,:);
  end;
end;
f=sum(sum(abs(e)));
v=t*100;