function?L=EAV（f）
f=double（f）;
[mn]=size（f）;
f=[zeros（m1）?f?zeros（m1）];
f=[zeros（1n+2）;f;zeros（1n+2）];

for?i=2:m??%i代表行
????for?j=2:n???%j代表列
????????H=（4+2*sqrt（2））*f（ij）