Gabor:求综合窗函数的对偶函数

花瓣雨

下面是我求求综合窗函数的对偶函数的过程，为什么图形不对，请各位指导，谢谢
%%%%%%%%%%%%%%%%%%%%%%%
L=128;
delM = 16;
delN=8;
N=L/delN;
M=L/delM;
hi=[];
aa=(delM*L)/(delN*2*pi);
temp1=(pi*aa).^(-0.25);
%%hi
i=0:L-1
hi=temp1.*exp(-((i-0.5*(L-1)).^2)/(2*aa));
plot(i,hi,'b');
hold on;
%%H
h=[];
aa=(delM/delN)*(L/(2*pi));
for i=0:L-1
    for m=0:delN-1
        for n=0:delM-1
            H(m*delM+n+1,i+1)=temp1.*exp(-((i+m*N-0.5*(L-1)).^2)/(2*aa))*exp(-2*pi*n*M*i/L);
        end
    end
end
%%%%%%%%%%%%%%
temp1=H*H';
temp2=inv(temp1);
ropt=H'*temp2;
p=delM*delN;
i=1:p-1;
u=[1,i-i];
u=u';
%temp2=H'*u;
%ropt=temp1\temp2;
ropt=ropt*u;
i=0:L-1;
plot(i,ropt,'r');
y=norm(ropt);