微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

逆离散傅立叶变换

如何解决逆离散傅立叶变换

我一直在尝试实现逆 DFT(不使用 FFT),但在数学上有点挣扎。从我阅读的实现和伪代码中,我无法弄清楚我的实现存在的问题。

'''

spectrum.clear();
spectrum.add(new Complexnumber(36.0,0.0));
spectrum.add(new Complexnumber(-4.000000000000002,9.65685424949238));
spectrum.add(new Complexnumber(-4.000000000000002,3.9999999999999987));
spectrum.add(new Complexnumber(-4.0,1.6568542494923788));
spectrum.add(new Complexnumber(-4.0,0.0));
spectrum.add(new Complexnumber(-3.999999999999999,-1.656854249492381));
spectrum.add(new Complexnumber(-3.999999999999998,-4.000000000000001));
spectrum.add(new Complexnumber(-3.999999999999995,-9.65685424949238));
        
int N = spectrum.size();
double twoPiOnN = 2 * Math.PI / N;

for(Complexnumber c : spectrum) {
    waveform.add(0.0);
}
double twoPikOnN,twoPijkOnN;

for(int k=0; k<N; k++) {
    twoPikOnN = twoPiOnN *k;
    for(int j = 0; j < N; j++) {
        twoPijkOnN = twoPikOnN * j;
        double d = waveform.get(k);
        d += spectrum.get(j).getRe() * Math.cos(twoPijkOnN ) - spectrum.get(j).getIm() * Math.sin(twoPijkOnN );
        waveform.set(k,d);
        }
}
for(Double d : waveform) {
    System.out.println(d);
}

''' 这种情况应该导致 [1,2,3,4,5,6,7,8] 的波形数组列表。

非常感谢任何帮助。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。