前面的公式实际叫做离散傅里叶变换(DFT),其最初实现如下(百度到一篇,其点醒了我,很简单,相遇了):
int N = inputArr.GetLength(0);//输入声音波形byte[]数组,在128上下区间震动(0,255)
double PI = 3.1415926;
double[] real = new double[N];//实部
double[] imag = new double[N];//虚部
double[] resultN = new double[N];
byte[] x = new byte[N];
byte[] histN = new byte[N];
if (inputArr != null)
{
x = inputArr;
}
else return;
for (int k = 0; k < N; k++)
{
for (int n = 0; n < N; n++)
{
real[k] = real[k] + (x[n] - 128) * Math.Cos(2 * PI * k * n / N);//128与振幅相关
imag[k] = imag[k] - (x[n] - 128) * Math.Sin(2 * PI * k * n / N);
}
resultN[k] = Math.Sqrt(real[k] * real[k] + imag[k] * imag[k]);//很简单吧,一个二重循环便实现了。
histN[k] = (byte)(255 - resultN[k] / 10);//此处计算与显示相关
}
(待续,慢慢来。。。。。。。。。。。)每天一点小改变☺