ÆäÖD￡¬s1oís2¾ùÎa80*80μÄêy×é¡￡Dèòa¼ÆËãμÄ1«ê½êÇ￡o
s2[p][q]=max(max(w[p-i,q-j]*s1[i][j]))
ÆäÖDp,q,i,j¶¼êÇ′ó0μ½79
w[m,n]=exp(-beita^2(m^2+n^2))

3ìDòèçÏÂ￡¬òѾ-êÔ׿õéù¼ÆËãá¿áË¡￡


NN=40;
for (i=0;i<inputN1;i=i+1)
{
for (j=0;j<inputN1;j=j+1)
{
nb1=j-NN;
if (nb1<0)
nb1=0;
nb2=j+NN;
if (nb2>inputN1)
nb2=inputN1;

nb3=i-NN;
if (nb3<0)
nb3=0;
nb4=i+NN;
if (nb4>inputN1)
nb4=inputN1;

aa=0;
for (i1=nb3;i1<nb4;i1=i1+1)
{
for (j1=nb1;j1<nb2;j1=j1+1)
{
j2=abs(i-i1);
j3=abs(j-j1);
fai=exp(-beita*(j2*j2+j3*j3));
bb=fai*s1[i1][j1];
if(aa<bb) aa=bb;
}
}
s2[i][j]=aa;
}
}