私はlibsndfileライブラリ(http://www.mega-nerd.com/libsndfile/api.html)を使用してwavオーディオファイルを読み込むCで簡単なソフトウェアを書いています。信号(2次バタワースローパスフィルタ、転置された直接形式IIで適用)。その後、結果を新しいwavファイルに書き出します。 フィルタの代わりに単純な操作(定数のmuliplyサンプルなど)を適用すると問題なく動作しますが、フィルタを適用するとノイズが大きくなります。 フィルタの後のサンプルの値を出力しようとしましたが、それらが新しいファイルに書き込まれる前にMatlabから取得したものと同じ値を得ましたが(出力は完全です)、値とは異なります私はライブラリによって書かれた出力を読んでもらえます。オーディオ信号にフィルタを適用したときの出力が異なる
static void processAudio (double *buffer, int length)
{
//arrays a and b are the coefficients
double a[] = { 1,
-1.799096409484668,
0.817512403384758};
double b[] = { 0.004603998475022,
0.009207996950045,
0.004603998475022};
//arrays a_ and b_ are the coefficients normalized
double a_[] = {1,a[1]/b[0],a[2]/b[0]};
double b_[] = {1,b[1]/b[0],b[2]/b[0]};
double gain = b[0]/a[0];
double reg[] ={0,0}; //memory registers
for(int i = 0; i<length; i++)
{
if(i%2==0) //just left channel is changed
{
//TRANSPOSED DIRECT FORM II
double input = *(buffer+i);
double output =(input + reg[0])*gain;
reg[0] = reg[1]+b_[1]*input-a_[1]*output;
reg[1] = b_[2]*input - a_[2]*output;
*(buffer+i) = output;
}
}
}
//that function is called in the main method inside this cycle
int main(void)
{
...
while ((framesRead = sf_read_double(inputFile, buffer, BUFFER_LENGTH)))
{
processAudio(buffer, framesRead);
sf_write_double(outputFile, buffer, framesRead);
}
...
}
私はちょうど私が手フィルタ後の結果を印刷する場合:
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
0.00000000000000
-0.00000014050288
-0.00000025277823
0.00000050310774
0.00000209530885
0.00000420138311
0.00000725078112
0.00001115570320
0.00001554761090
0.00002053775981
0.00002550357112
0.00002796948679
0.00000727086200
-0.00006401853354
...
私は、出力ファイルを読み込む場合の結果は、彼らは非常に異なっている
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
3.05175781250000e-05
3.05175781250000e-05
3.05175781250000e-05
3.05175781250000e-05
0
-6.10351562500000e-05
-0.000366210937500000
-0.00183105468750000
-0.00601196289062500
-0.0140075683593750
-0.0261840820312500
-0.0421142578125000
-0.0610656738281250
-0.0820617675781250
-0.104156494140625
-0.126342773437500
-0.147766113281250
-0.167663574218750
-0.185363769531250
-0.200378417968750
-0.212280273437500
...
です。私は何が起こっているのか本当に知りません。アイデアがあれば教えてください。前もって感謝します!
あなたのコードを実際の[mcve]にしてください。 'a'と' b'はどのような型ですか? –
私は使用している係数を追加しました... – thebesttony
このような全体の質問は1つの大きな赤いニシンです。質問はフィルタに焦点を当てていますが、フィルタの出力には明らかに満足しています。問題と思われるファイルの書き込み/読み込みのみです。言い換えれば、問題は 'sf_write_double'と何か関係があります。ですから、私は他の迷惑メールを取り除き、 'sf_write_double'を使ってどのように動作するかを調べるだけです。 – user3386109