MATLAB 双通道PCNN算法

function Fused=Dual_PCNN(Ahigh,Bhigh,Abeta,Bbeta)

SF_Ahigh =SF(abs(Ahigh));
SF_Bhigh = SF(abs(Bhigh));
[p1,q1]=size(Ahigh);
[p2,q2]=size(Bhigh);
if(p1~=p2||q1~=q2)
    error('The size of image A and image B must be the same......');
end

%% set the parameter for dual channel pcnn
alphaL = 0.01;
alphaTheta = 0.1;
Vl = 0.2;
con=1;

%% Dual channel pcnn process
L = zeros(p1,q1);
U = zeros(p1,q1);
Theta = zeros(p1,q1);
Y = zeros(p1,q1);
fuse=zeros(p1,q1);
VTheta = 25;
%beta=2.6;
W=[0.707,1,0.707;1,0,1;0.707,1,0.707];
N1 = 200;%iteration times
for num = 1:N1;
    K = conv2(Y,W,'same');
    L = exp(-alphaL)*L + Vl*K;
    Theta = exp(-alphaTheta)*Theta + VTheta*Y;
    U = max(SF_Ahigh.*(con+Abeta.*L),SF_Bhigh.*(con+Bbeta.*L));
    Y = im2double(U>Theta);
    %Y0 = Y0 +Y;
end
for i = 1:p1
    for j = 1:q1
        if(SF_Ahigh(i,j)*(con + Abeta(i,j)*L(i,j)) == U(i,j))
            fuse(i,j) = Ahigh(i,j);
        else
            if(SF_Bhigh(i,j)*(con + Bbeta(i,j)*L(i,j)) == U(i,j))
                fuse(i,j) = Bhigh(i,j);
            end
        end
    end
end
Fused=fuse;
end