peparabimp.m

来自「分析了人工免疫系统的原理,简单建立了一种基于否定选择算法的网络入侵检测系统, 」· M 代码 · 共 30 行

function u = peParabImp(c,dt,n,minx,maxx,lbu,rbu,M)
format long;
h = (maxx-minx)/(n-1);
u0(1) = lbu;
u0(n) = rbu;
for j=2:n-1
    u0(j) = PrIniU(minx+(j-1)*h);
end 
u1 = u0;

for k=1:M
    A = zeros(n-2,n-2);
    cb =  - transpose(u0(2:(n-1)));
    cb(1) = cb(1) - dt*c*lbu/h/h;
    cb(n-2) = cb(n-2) - dt*c*rbu/h/h;
    A(1,1) = -2*dt*c/h/h -1;
    A(1,2) = dt*c/h/h ;
    for i=2:n-3
        A(i,i-1) = dt*c/h/h ;
        A(i,i) = - 2*dt*c/h/h -1 ;
        A(i,i+1) = dt*c/h/h  ;
    end
    A(n-2,n-2) = -2*dt*c/h/h -1;
    A(n-2,n-3) = dt*c/h/h;
    u1(2:(n-1)) = A\cb;
    u0 = u1;
end

u = u1;
format short;