📄 e.m
字号:
function Enum=E(WI1,WI2,width,height);
%计算能量E;
Enum=0;
for x1=1:width
for x2=1:height
%计算JI,JM;
for v=1:5
for u=1:8
JI((v-1)*8+u)=WI1(v,u,x1,x2);
JM((v-1)*8+u)=WI2(v,u,x1,x2);
end
end
%Enum=Enum-JI*rot90(JM)/sqrt(JI*rot90(JI))/sqrt(JM*rot90(JM));
JIM=0;JMM=0;
for k=1:40
JIM=JIM+JI(k)*JI(k);
end
for k=1:40
JMM=JMM+JM(k)*JM(k);
end
JIM=sqrt(JIM);JMM=sqrt(JMM);
for k=1:40
Enum=Enum-JI(k)*JM(k)/JIM/JMM;%12个中错了2个;如果用平方差,错4个(1/3);
end
end
end
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -