VC++ 下FFT试验

2019年11月6日 159次阅读

#include

using namespace std;

#define MAX
2048

#define FFT_1024 

#define FFTNUM
1024

#ifdef
FFT_1024

const double sin_tab[] = {

0.0000,0.00614,0.01227,0.01841,0.02454,0.03067,0.03681,0.04294,0.04907,0.0552,0.06132,0.06744,0.07356,0.07968,0.0858,0.09191,

0.09802,0.10412,0.11022,0.11632,0.12241,0.1285,0.13458,0.14066,0.14673,0.1528,0.15886,0.16491,0.17096,0.177,0.18304,0.18907,

0.19509,0.2011,0.20711,0.21311,0.2191,0.22508,0.23106,0.23702,0.24298,0.24893,0.25487,0.26079,0.26671,0.27262,0.27852,0.28441,

0.29028,0.29615,0.30201,0.30785,0.31368,0.3195,0.32531,0.33111,0.33689,0.34266,0.34842,0.35416,0.3599,0.36561,0.37132,0.37701,

0.38268,0.38835,0.39399,0.39962,0.40524,0.41084,0.41643,0.422,0.42756,0.43309,0.43862,0.44412,0.44961,0.45508,0.46054,0.46598,

0.4714,0.4768,0.48218,0.48755,0.4929,0.49823,0.50354,0.50883,0.5141,0.51936,0.52459,0.5298,0.535,0.54017,0.54532,0.55046,

0.55557,0.56066,0.56573,0.57078,0.57581,0.58081,0.5858,0.59076,0.5957,0.60062,0.60551,0.61038,0.61523,0.62006,0.62486,0.62964,

0.63439,0.63912,0.64383,0.64851,0.65317,0.65781,0.66242,0.667,0.67156,0.67609,0.6806,0.68508,0.68954,0.69397,0.69838,0.70275,

0.70711,0.71143,0.71573,0.72,0.72425,0.72846,0.73265,0.73682,0.74095,0.74506,0.74914,0.75319,0.75721,0.7612,0.76517,0.7691,

0.77301,0.77689,0.78074,0.78456,0.78835,0.79211,0.79584,0.79954,0.80321,0.80685,0.81046,0.81404,0.81758,0.8211,0.82459,0.82805,

0.83147,0.83486,0.83822,0.84155,0.84485,0.84812,0.85136,0.85456,0.85773,0.86087,0.86397,0.86705,0.87009,0.87309,0.87607,0.87901,

0.88192,0.8848,0.88764,0.89045,0.89322,0.89597,0.89867,0.90135,0.90399,0.9066,0.90917,0.91171,0.91421,0.91668,0.91911,0.92151,

0.92388,0.92621,0.92851,0.93077,0.93299,0.93518,0.93734,0.93946,0.94154,0.94359,0.94561,0.94759,0.94953,0.95144,0.95331,0.95514,

0.95694,0.9587,0.96043,0.96212,0.96378,0.96539,0.96698,0.96852,0.97003,0.9715,0.97294,0.97434,0.9757,0.97703,0.97832,0.97957,

0.98079,0.98196,0.98311,0.98421,0.98528,0.98631,0.9873,0.98826,0.98918,0.99006,0.9909,0.99171,0.99248,0.99321,0.99391,0.99456,

0.99518,0.99577,0.99631,0.99682,0.99729,0.99772,0.99812,0.99848,0.9988,0.99908,0.99932,0.99953,0.9997,0.99983,0.99992,0.99998,

1,0.99998,0.99992,0.99983,0.9997,0.99953,0.99932,0.99908,0.9988,0.99848,0.99812,0.99772,0.99729,0.99682,0.99631,0.99577,

0.99518,0.99456,0.99391,0.99321,0.99248,0.99171,0.9909,0.99006,0.98918,0.98826,0.9873,0.98631,0.98528,0.98421,0.98311,0.98196,

0.98079,0.97957,0.97832,0.97703,0.9757,0.97434,0.97294,0.9715,0.97003,0.96852,0.96698,0.96539,0.96378,0.96212,0.96043,0.9587,

0.95694,0.95514,0.95331,0.95144,0.94953,0.94759,0.94561,0.94359,0.94154,0.93946,0.93734,0.93518,0.93299,0.93077,0.92851,0.92621,

0.92388,0.92151,0.91911,0.91668,0.91421,0.91171,0.90917,0.9066,0.90399,0.90135,0.89867,0.89597,0.89322,0.89045,0.88764,0.8848,

0.88192,0.87901,0.87607,0.87309,0.87009,0.86705,0.86397,0.86087,0.85773,0.85456,0.85136,0.84812,0.84485,0.84155,0.83822,0.83486,

0.83147,0.82805,0.82459,0.8211,0.81758,0.81404,0.81046,0.80685,0.80321,0.79954,0.79584,0.79211,0.78835,0.78456,0.78074,0.77689,

0.77301,0.7691,0.76517,0.7612,0.75721,0.75319,0.74914,0.74506,0.74095,0.73682,0.73265,0.72846,0.72425,0.72,0.71573,0.71143,

0.70711,0.70275,0.69838,0.69397,0.68954,0.68508,0.6806,0.67609,0.67156,0.667,0.66242,0.65781,0.65317,0.64851,0.64383,0.63912,

0.63439,0.62964,0.62486,0.62006,0.61523,0.61038,0.60551,0.60062,0.5957,0.59076,0.5858,0.58081,0.57581,0.57078,0.56573,0.56066,

0.55557,0.55046,0.54532,0.54017,0.535,0.5298,0.52459,0.51936,0.5141,0.50883,0.50354,0.49823,0.4929,0.48755,0.48218,0.4768,0.4714,

0.46598,0.46054,0.45508,0.44961,0.44412,0.43862,0.43309,0.42756,0.422,0.41643,0.41084,0.40524,0.39962,0.39399,0.38835,0.38268,

0.37701,0.37132,0.36561,0.3599,0.35416,0.34842,0.34266,0.33689,0.33111,0.32531,0.3195,0.31368,0.30785,0.30201,0.29615,0.29028,

0.28441,0.27852,0.27262,0.26671,0.26079,0.25487,0.24893,0.24298,0.23702,0.23106,0.22508,0.2191,0.21311,0.20711,0.2011,0.19509,

0.18907,0.18304,0.177,0.17096,0.16491,0.15886,0.1528,0.14673,0.14066,0.13458,0.1285,0.12241,0.11632,0.11022,0.10412,0.09802,

0.09191,0.0858,0.07968,0.07356,0.06744,0.06132,0.0552,0.04907,0.04294,0.03681,0.03067,0.02454,0.01841,0.01227,0.00614,0,

-0.00614,-0.01227,-0.01841,-0.02454,-0.03067,-0.03681,-0.04294,-0.04907,-0.0552,-0.06132,-0.06744,-0.07356,-0.07968,-0.0858,

-0.09191,-0.09802,-0.10412,-0.11022,-0.11632,-0.12241,-0.1285,-0.13458,-0.14066,-0.14673,-0.1528,-0.15886,-0.16491,-0.17096,

-0.177,-0.18304,-0.18907,-0.19509,-0.2011,-0.20711,-0.21311,-0.2191,-0.22508,-0.23106,-0.23702,-0.24298,-0.24893,-0.25487,

-0.26079,-0.26671,-0.27262,-0.27852,-0.28441,-0.29028,-0.29615,-0.30201,-0.30785,-0.31368,-0.3195,-0.32531,-0.33111,-0.33689,

-0.34266,-0.34842,-0.35416,-0.3599,-0.36561,-0.37132,-0.37701,-0.38268,-0.38835,-0.39399,-0.39962,-0.40524,-0.41084,-0.41643,

-0.422,-0.42756,-0.43309,-0.43862,-0.44412,-0.44961,-0.45508,-0.46054,-0.46598,-0.4714,-0.4768,-0.48218,-0.48755,-0.4929,

-0.49823,-0.50354,-0.50883,-0.5141,-0.51936,-0.52459,-0.5298,-0.535,-0.54017,-0.54532,-0.55046,-0.55557,-0.56066,-0.56573,

-0.57078,-0.57581,-0.58081,-0.5858,-0.59076,-0.5957,-0.60062,-0.60551,-0.61038,-0.61523,-0.62006,-0.62486,-0.62964,-0.63439,

-0.63912,-0.64383,-0.64851,-0.65317,-0.65781,-0.66242,-0.667,-0.67156,-0.67609,-0.6806,-0.68508,-0.68954,-0.69397,-0.69838,

-0.70275,-0.70711,-0.71143,-0.71573,-0.72,-0.72425,-0.72846,-0.73265,-0.73682,-0.74095,-0.74506,-0.74914,-0.75319,-0.75721,

-0.7612,-0.76517,-0.7691,-0.77301,-0.77689,-0.78074,-0.78456,-0.78835,-0.79211,-0.79584,-0.79954,-0.80321,-0.80685,-0.81046,

-0.81404,-0.81758,-0.8211,-0.82459,-0.82805,-0.83147,-0.83486,-0.83822,-0.84155,-0.84485,-0.84812,-0.85136,-0.85456,-0.85773,

-0.86087,-0.86397,-0.86705,-0.87009,-0.87309,-0.87607,-0.87901,-0.88192,-0.8848,-0.88764,-0.89045,-0.89322,-0.89597,-0.89867,

-0.90135,-0.90399,-0.9066,-0.90917,-0.91171,-0.91421,-0.91668,-0.91911,-0.92151,-0.92388,-0.92621,-0.92851,-0.93077,-0.93299,

-0.93518,-0.93734,-0.93946,-0.94154,-0.94359,-0.94561,-0.94759,-0.94953,-0.95144,-0.95331,-0.95514,-0.95694,-0.9587,-0.96043,

-0.96212,-0.96378,-0.96539,-0.96698,-0.96852,-0.97003,-0.9715,-0.97294,-0.97434,-0.9757,-0.97703,-0.97832,-0.97957,-0.98079,

-0.98196,-0.98311,-0.98421,-0.98528,-0.98631,-0.9873,-0.98826,-0.98918,-0.99006,-0.9909,-0.99171,-0.99248,-0.99321,-0.99391,

-0.99456,-0.99518,-0.99577,-0.99631,-0.99682,-0.99729,-0.99772,-0.99812,-0.99848,-0.9988,-0.99908,-0.99932,-0.99953,-0.9997,

-0.99983,-0.99992,-0.99998,-1,-0.99998,-0.99992,-0.99983,-0.9997,-0.99953,-0.99932,-0.99908,-0.9988,-0.99848,-0.99812,-0.99772,

-0.99729,-0.99682,-0.99631,-0.99577,-0.99518,-0.99456,-0.99391,-0.99321,-0.99248,-0.99171,-0.9909,-0.99006,-0.98918,-0.98826,

-0.9873,-0.98631,-0.98528,-0.98421,-0.98311,-0.98196,-0.98079,-0.97957,-0.97832,-0.97703,-0.9757,-0.97434,-0.97294,-0.9715,

-0.97003,-0.96852,-0.96698,-0.96539,-0.96378,-0.96212,-0.96043,-0.9587,-0.95694,-0.95514,-0.95331,-0.95144,-0.94953,-0.94759,

-0.94561,-0.94359,-0.94154,-0.93946,-0.93734,-0.93518,-0.93299,-0.93077,-0.92851,-0.92621,-0.92388,-0.92151,-0.91911,-0.91668,

-0.91421,-0.91171,-0.90917,-0.9066,-0.90399,-0.90135,-0.89867,-0.89597,-0.89322,-0.89045,-0.88764,-0.8848,-0.88192,-0.87901,

-0.87607,-0.87309,-0.87009,-0.86705,-0.86397,-0.86087,-0.85773,-0.85456,-0.85136,-0.84812,-0.84485,-0.84155,-0.83822,-0.83486,

-0.83147,-0.82805,-0.82459,-0.8211,-0.81758,-0.81404,-0.81046,-0.80685,-0.80321,-0.79954,-0.79584,-0.79211,-0.78835,-0.78456,

-0.78074,-0.77689,-0.77301,-0.7691,-0.76517,-0.7612,-0.75721,-0.75319,-0.74914,-0.74506,-0.74095,-0.73682,-0.73265,-0.72846,

-0.72425,-0.72,-0.71573,-0.71143,-0.70711,-0.70275,-0.69838,-0.69397,-0.68954,-0.68508,-0.6806,-0.67609,-0.67156,-0.667,

-0.66242,-0.65781,-0.65317,-0.64851,-0.64383,-0.63912,-0.63439,-0.62964,-0.62486,-0.62006,-0.61523,-0.61038,-0.60551,-0.60062,

-0.5957,-0.59076,-0.5858,-0.58081,-0.57581,-0.57078,-0.56573,-0.56066,-0.55557,-0.55046,-0.54532,-0.54017,-0.535,-0.5298,-0.52459,

-0.51936,-0.5141,-0.50883,-0.50354,-0.49823,-0.4929,-0.48755,-0.48218,-0.4768,-0.4714,-0.46598,-0.46054,-0.45508,-0.44961,-0.44412,

-0.43862,-0.43309,-0.42756,-0.422,-0.41643,-0.41084,-0.40524,-0.39962,-0.39399,-0.38835,-0.38268,-0.37701,-0.37132,-0.36561,-0.3599,

-0.35416,-0.34842,-0.34266,-0.33689,-0.33111,-0.32531,-0.3195,-0.31368,-0.30785,-0.30201,-0.29615,-0.29028,-0.28441,-0.27852,

-0.27262,-0.26671,-0.26079,-0.25487,-0.24893,-0.24298,-0.23702,-0.23106,-0.22508,-0.2191,-0.21311,-0.20711,-0.2011,-0.19509,

-0.18907,-0.18304,-0.177,-0.17096,-0.16491,-0.15886,-0.1528,-0.14673,-0.14066,-0.13458,-0.1285,-0.12241,-0.11632,-0.11022,

-0.10412,-0.09802,-0.09191,-0.0858,-0.07968,-0.07356,-0.06744,-0.06132,-0.0552,-0.04907,-0.04294,-0.03681,-0.03067,-0.02454,

-0.01841,-0.01227,-0.00614};

const double cos_tab[] = {

1,0.99998,0.99992,0.99983,0.9997,0.99953,0.99932,0.99908,0.9988,0.99848,0.99812,0.99772,0.99729,0.99682,0.99631,0.99577,0.99518,

0.99456,0.99391,0.99321,0.99248,0.99171,0.9909,0.99006,0.98918,0.98826,0.9873,0.98631,0.98528,0.98421,0.98311,0.98196,0.98079,

0.97957,0.97832,0.97703,0.9757,0.97434,0.97294,0.9715,0.97003,0.96852,0.96698,0.96539,0.96378,0.96212,0.96043,0.9587,0.95694,

0.95514,0.95331,0.95144,0.94953,0.94759,0.94561,0.94359,0.94154,0.93946,0.93734,0.93518,0.93299,0.93077,0.92851,0.92621,0.92388,

0.92151,0.91911,0.91668,0.91421,0.91171,0.90917,0.9066,0.90399,0.90135,0.89867,0.89597,0.89322,0.89045,0.88764,0.8848,0.88192,

0.87901,0.87607,0.87309,0.87009,0.86705,0.86397,0.86087,0.85773,0.85456,0.85136,0.84812,0.84485,0.84155,0.83822,0.83486,0.83147,

0.82805,0.82459,0.8211,0.81758,0.81404,0.81046,0.80685,0.80321,0.79954,0.79584,0.79211,0.78835,0.78456,0.78074,0.77689,0.77301,

0.7691,0.76517,0.7612,0.75721,0.75319,0.74914,0.74506,0.74095,0.73682,0.73265,0.72846,0.72425,0.72,0.71573,0.71143,0.70711,0.70275,

0.69838,0.69397,0.68954,0.68508,0.6806,0.67609,0.67156,0.667,0.66242,0.65781,0.65317,0.64851,0.64383,0.63912,0.63439,0.62964,

0.62486,0.62006,0.61523,0.61038,0.60551,0.60062,0.5957,0.59076,0.5858,0.58081,0.57581,0.57078,0.56573,0.56066,0.55557,0.55046,

0.54532,0.54017,0.535,0.5298,0.52459,0.51936,0.5141,0.50883,0.50354,0.49823,0.4929,0.48755,0.48218,0.4768,0.4714,0.46598,0.46054,

0.45508,0.44961,0.44412,0.43862,0.43309,0.42756,0.422,0.41643,0.41084,0.40524,0.39962,0.39399,0.38835,0.38268,0.37701,0.37132,

0.36561,0.3599,0.35416,0.34842,0.34266,0.33689,0.33111,0.32531,0.3195,0.31368,0.30785,0.30201,0.29615,0.29028,0.28441,0.27852,

0.27262,0.26671,0.26079,0.25487,0.24893,0.24298,0.23702,0.23106,0.22508,0.2191,0.21311,0.20711,0.2011,0.19509,0.18907,0.18304,

0.177,0.17096,0.16491,0.15886,0.1528,0.14673,0.14066,0.13458,0.1285,0.12241,0.11632,0.11022,0.10412,0.09802,0.09191,0.0858,0.07968,

0.07356,0.06744,0.06132,0.0552,0.04907,0.04294,0.03681,0.03067,0.02454,0.01841,0.01227,0.00614,0,-0.00614,-0.01227,-0.01841,

-0.02454,-0.03067,-0.03681,-0.04294,-0.04907,-0.0552,-0.06132,-0.06744,-0.07356,-0.07968,-0.0858,-0.09191,-0.09802,-0.10412,

-0.11022,-0.11632,-0.12241,-0.1285,-0.13458,-0.14066,-0.14673,-0.1528,-0.15886,-0.16491,-0.17096,-0.177,-0.18304,-0.18907,

-0.19509,-0.2011,-0.20711,-0.21311,-0.2191,-0.22508,-0.23106,-0.23702,-0.24298,-0.24893,-0.25487,-0.26079,-0.26671,-0.27262,

-0.27852,-0.28441,-0.29028,-0.29615,-0.30201,-0.30785,-0.31368,-0.3195,-0.32531,-0.33111,-0.33689,-0.34266,-0.34842,-0.35416,

-0.3599,-0.36561,-0.37132,-0.37701,-0.38268,-0.38835,-0.39399,-0.39962,-0.40524,-0.41084,-0.41643,-0.422,-0.42756,-0.43309,

-0.43862,-0.44412,-0.44961,-0.45508,-0.46054,-0.46598,-0.4714,-0.4768,-0.48218,-0.48755,-0.4929,-0.49823,-0.50354,-0.50883,

-0.5141,-0.51936,-0.52459,-0.5298,-0.535,-0.54017,-0.54532,-0.55046,-0.55557,-0.56066,-0.56573,-0.57078,-0.57581,-0.58081,-0.5858,

-0.59076,-0.5957,-0.60062,-0.60551,-0.61038,-0.61523,-0.62006,-0.62486,-0.62964,-0.63439,-0.63912,-0.64383,-0.64851,-0.65317,

-0.65781,-0.66242,-0.667,-0.67156,-0.67609,-0.6806,-0.68508,-0.68954,-0.69397,-0.69838,-0.70275,-0.70711,-0.71143,-0.71573,-0.72,

-0.72425,-0.72846,-0.73265,-0.73682,-0.74095,-0.74506,-0.74914,-0.75319,-0.75721,-0.7612,-0.76517,-0.7691,-0.77301,-0.77689,

-0.78074,-0.78456,-0.78835,-0.79211,-0.79584,-0.79954,-0.80321,-0.80685,-0.81046,-0.81404,-0.81758,-0.8211,-0.82459,-0.82805,

-0.83147,-0.83486,-0.83822,-0.84155,-0.84485,-0.84812,-0.85136,-0.85456,-0.85773,-0.86087,-0.86397,-0.86705,-0.87009,-0.87309,

-0.87607,-0.87901,-0.88192,-0.8848,-0.88764,-0.89045,-0.89322,-0.89597,-0.89867,-0.90135,-0.90399,-0.9066,-0.90917,-0.91171,

-0.91421,-0.91668,-0.91911,-0.92151,-0.92388,-0.92621,-0.92851,-0.93077,-0.93299,-0.93518,-0.93734,-0.93946,-0.94154,-0.94359,

-0.94561,-0.94759,-0.94953,-0.95144,-0.95331,-0.95514,-0.95694,-0.9587,-0.96043,-0.96212,-0.96378,-0.96539,-0.96698,-0.96852,

-0.97003,-0.9715,-0.97294,-0.97434,-0.9757,-0.97703,-0.97832,-0.97957,-0.98079,-0.98196,-0.98311,-0.98421,-0.98528,-0.98631,

-0.9873,-0.98826,-0.98918,-0.99006,-0.9909,-0.99171,-0.99248,-0.99321,-0.99391,-0.99456,-0.99518,-0.99577,-0.99631,-0.99682,

-0.99729,-0.99772,-0.99812,-0.99848,-0.9988,-0.99908,-0.99932,-0.99953,-0.9997,-0.99983,-0.99992,-0.99998,-1,-0.99998,-0.99992,

-0.99983,-0.9997,-0.99953,-0.99932,-0.99908,-0.9988,-0.99848,-0.99812,-0.99772,-0.99729,-0.99682,-0.99631,-0.99577,-0.99518,

-0.99456,-0.99391,-0.99321,-0.99248,-0.99171,-0.9909,-0.99006,-0.98918,-0.98826,-0.9873,-0.98631,-0.98528,-0.98421,-0.98311,

-0.98196,-0.98079,-0.97957,-0.97832,-0.97703,-0.9757,-0.97434,-0.97294,-0.9715,-0.97003,-0.96852,-0.96698,-0.96539,-0.96378,

-0.96212,-0.96043,-0.9587,-0.95694,-0.95514,-0.95331,-0.95144,-0.94953,-0.94759,-0.94561,-0.94359,-0.94154,-0.93946,-0.93734,

-0.93518,-0.93299,-0.93077,-0.92851,-0.92621,-0.92388,-0.92151,-0.91911,-0.91668,-0.91421,-0.91171,-0.90917,-0.9066,-0.90399,

-0.90135,-0.89867,-0.89597,-0.89322,-0.89045,-0.88764,-0.8848,-0.88192,-0.87901,-0.87607,-0.87309,-0.87009,-0.86705,-0.86397,

-0.86087,-0.85773,-0.85456,-0.85136,-0.84812,-0.84485,-0.84155,-0.83822,-0.83486,-0.83147,-0.82805,-0.82459,-0.8211,-0.81758,

-0.81404,-0.81046,-0.80685,-0.80321,-0.79954,-0.79584,-0.79211,-0.78835,-0.78456,-0.78074,-0.77689,-0.77301,-0.7691,-0.76517,

-0.7612,-0.75721,-0.75319,-0.74914,-0.74506,-0.74095,-0.73682,-0.73265,-0.72846,-0.72425,-0.72,-0.71573,-0.71143,-0.70711,

-0.70275,-0.69838,-0.69397,-0.68954,-0.68508,-0.6806,-0.67609,-0.67156,-0.667,-0.66242,-0.65781,-0.65317,-0.64851,-0.64383,

-0.63912,-0.63439,-0.62964,-0.62486,-0.62006,-0.61523,-0.61038,-0.60551,-0.60062,-0.5957,-0.59076,-0.5858,-0.58081,-0.57581,

-0.57078,-0.56573,-0.56066,-0.55557,-0.55046,-0.54532,-0.54017,-0.535,-0.5298,-0.52459,-0.51936,-0.5141,-0.50883,-0.50354,

-0.49823,-0.4929,-0.48755,-0.48218,-0.4768,-0.4714,-0.46598,-0.46054,-0.45508,-0.44961,-0.44412,-0.43862,-0.43309,-0.42756,

-0.422,-0.41643,-0.41084,-0.40524,-0.39962,-0.39399,-0.38835,-0.38268,-0.37701,-0.37132,-0.36561,-0.3599,-0.35416,-0.34842,

-0.34266,-0.33689,-0.33111,-0.32531,-0.3195,-0.31368,-0.30785,-0.30201,-0.29615,-0.29028,-0.28441,-0.27852,-0.27262,-0.26671,

-0.26079,-0.25487,-0.24893,-0.24298,-0.23702,-0.23106,-0.22508,-0.2191,-0.21311,-0.20711,-0.2011,-0.19509,-0.18907,-0.18304,

-0.177,-0.17096,-0.16491,-0.15886,-0.1528,-0.14673,-0.14066,-0.13458,-0.1285,-0.12241,-0.11632,-0.11022,-0.10412,-0.09802,

-0.09191,-0.0858,-0.07968,-0.07356,-0.06744,-0.06132,-0.0552,-0.04907,-0.04294,-0.03681,-0.03067,-0.02454,-0.01841,-0.01227,

-0.00614,0,0.00614,0.01227,0.01841,0.02454,0.03067,0.03681,0.04294,0.04907,0.0552,0.06132,0.06744,0.07356,0.07968,0.0858,

0.09191,0.09802,0.10412,0.11022,0.11632,0.12241,0.1285,0.13458,0.14066,0.14673,0.1528,0.15886,0.16491,0.17096,0.177,0.18304,

0.18907,0.19509,0.2011,0.20711,0.21311,0.2191,0.22508,0.23106,0.23702,0.24298,0.24893,0.25487,0.26079,0.26671,0.27262,0.27852,

0.28441,0.29028,0.29615,0.30201,0.30785,0.31368,0.3195,0.32531,0.33111,0.33689,0.34266,0.34842,0.35416,0.3599,0.36561,0.37132,

0.37701,0.38268,0.38835,0.39399,0.39962,0.40524,0.41084,0.41643,0.422,0.42756,0.43309,0.43862,0.44412,0.44961,0.45508,0.46054,

0.46598,0.4714,0.4768,0.48218,0.48755,0.4929,0.49823,0.50354,0.50883,0.5141,0.51936,0.52459,0.5298,0.535,0.54017,0.54532,0.55046,

0.55557,0.56066,0.56573,0.57078,0.57581,0.58081,0.5858,0.59076,0.5957,0.60062,0.60551,0.61038,0.61523,0.62006,0.62486,0.62964,

0.63439,0.63912,0.64383,0.64851,0.65317,0.65781,0.66242,0.667,0.67156,0.67609,0.6806,0.68508,0.68954,0.69397,0.69838,0.70275,

0.70711,0.71143,0.71573,0.72,0.72425,0.72846,0.73265,0.73682,0.74095,0.74506,0.74914,0.75319,0.75721,0.7612,0.76517,0.7691,

0.77301,0.77689,0.78074,0.78456,0.78835,0.79211,0.79584,0.79954,0.80321,0.80685,0.81046,0.81404,0.81758,0.8211,0.82459,0.82805,

0.83147,0.83486,0.83822,0.84155,0.84485,0.84812,0.85136,0.85456,0.85773,0.86087,0.86397,0.86705,0.87009,0.87309,0.87607,0.87901,

0.88192,0.8848,0.88764,0.89045,0.89322,0.89597,0.89867,0.90135,0.90399,0.9066,0.90917,0.91171,0.91421,0.91668,0.91911,0.92151,

0.92388,0.92621,0.92851,0.93077,0.93299,0.93518,0.93734,0.93946,0.94154,0.94359,0.94561,0.94759,0.94953,0.95144,0.95331,0.95514,

0.95694,0.9587,0.96043,0.96212,0.96378,0.96539,0.96698,0.96852,0.97003,0.9715,0.97294,0.97434,0.9757,0.97703,0.97832,0.97957,

0.98079,0.98196,0.98311,0.98421,0.98528,0.98631,0.9873,0.98826,0.98918,0.99006,0.9909,0.99171,0.99248,0.99321,0.99391,0.99456,

0.99518,0.99577,0.99631,0.99682,0.99729,0.99772,0.99812,0.99848,0.9988,0.99908,0.99932,0.99953,0.9997,0.99983,0.99992,0.99998};

#endif

int FFTdotnum = 0;

int M = 0;//M 级

float data[MAX];

float brevdataRe[MAX];

float brevdataIm[MAX];

float magdata[MAX];

//获取FFT点数 和 待处理数据

void GetFFT_Dots()

{

unsigned int i = 0;

cout<<“Please enter the dot num of FFT:”;

cin>>FFTdotnum;

cout<<“Please enter the “<<FFTdotnum <<” datas”<<endl;

while(i<FFTdotnum)

{

cin>>data[i++];

}

}

//将输入数组 data[] 倒位排序存放至 bredata[]

//如输入数组下标为bit3 bit2 bit1 bit0 将其存放至 bit0 bit1 bit2 bit3 对应的新地址中

void FFT_brev(float* data,float* bredata,unsigned int Dotnums)

{

unsigned int i = 0,j = 0,k = 0,Dotnum = 1;

//log2N , 求N

for(i=0;Dotnum<Dotnums;i++)Dotnum = Dotnum<<1;

M = i;

for(i=0;i<Dotnums;i++)

{

k = 0;

for(j=0;j<M;j++)

{

if(i&(1<<j))k |= (1<<(M-1))>>j;

}

*(bredata+k) = *(data+i);

}

}

//FFT计算

void FFT_butterfly()

{

unsigned int m,n,k;//m级，每级下分n个Wn，每个Wn下有几次蝶形运算

unsigned int N,p;

float temR,temI,tempRp,tempIp;

for(m=0;m<M;m++)

{

N = 1<<m;//N = 2^m

for(n=0;n<N;n++)//每级下W因子的个数

{

p = (1<<(M-m-1));//p = 2^(M-m-1)同级不同旋转因子之间的下标间距

p = p*n;

p = p*(FFTNUM>>M);

for(k=n;k<FFTdotnum;k += 2*N)//计算相同W因子的所有蝶形

{

tempRp = brevdataRe[k];

tempIp = brevdataIm[k];

temR = brevdataRe[k+N]*cos_tab[p] + brevdataIm[k+N]*sin_tab[p];

temI = – brevdataRe[k+N]*sin_tab[p] + brevdataIm[k+N]*cos_tab[p];

brevdataRe[k+N] = tempRp – temR;

brevdataIm[k+N] = tempIp – temI;

brevdataRe[k] = tempRp + temR;

brevdataIm[k] = tempIp + temI;

}

}

}

}

int main()

{

unsigned int i = 0;

GetFFT_Dots();

FFT_brev(data,brevdataRe,FFTdotnum);

FFT_butterfly();

for(i=0;i<FFTdotnum;i++)cout<<i<<“:  “<<brevdataRe[i]<<”     “<<brevdataIm[i]<<endl;

while(1);

return 1;

}

点赞