dewt_glg.m

人工免疫是对生物免疫系统的工作机理及所具有各种优良特性的模拟。在邮件系统中 引入免疫的思想,构建了一种基于人工免疫的反垃圾邮件模型。

function yy = DEWT_glg(f,h,a,b,gama,y0,order,varvec)
%一阶常微分方程的一般表达式的右端函数：f
%积分步长：h
%自变量取值下限：a
%自变量取值上限：b
%外推参数，参考外推公式：gama
%函数初值：y0
%外推阶数：order
%常微分方程的变量组：varvec

ArrayH = [1;2;4;6;8;12;16;24;32;48;64;96];
N = (b-a)/h;
yy = zeros(N+1,1);
for i = 2:N+1
    dh = h;
    s = zeros(order,1);
    for j=1:order
        dh = h/ArrayH(j);
        NI = (i-1)* ArrayH(j);
        tmpY = zeros(NI+2,1);
        tmpY(1) = y0;
        tmpY(2) = y0 + dh*Funval(f,varvec,[a+dh y0]);
        for m=3:NI+2
            tmpY(m) = tmpY(m-2)+ 2*dh*Funval(f,varvec,[a+(m-1)*dh tmpY(m-1)]);
        end
        s(j) = tmpY(NI)/4+ tmpY(NI+1)/2+ tmpY(NI+2)/4;
    end
    tmpS = zeros(order,1);
    for j=1:order-1
        for k=(j+1):order
            tmpS(k) = s(k)+(s(k)-s(k-1))/((ArrayH(k)/ArrayH(j))^gama-1);
        end
        s(1:(order-j)) = tmpS((j+1):order);
    end
    yy(i) = tmpS(order);
end