私はCコードを持っています:
mel[i] = malloc(sizeof(double)*mellength[i]);
memcpy(mel[i],&temp[melstart[i]],mellength[i]*sizeof(double));
これを C++ に変換し、配列の代わりにベクトルを使用します。
vector< vector< double > > mel; // 2D vector
mel[i].resize(mellength[i]);
「temp」(1D) という配列があり、「temp」のすべての要素をベクトル mel にコピーしたいと考えています。私はこれを試しました:
for(int k=0; (k < mellength[i]); k++)
{
mel[i][k] = temp[i];
cout << mel[i][k] << endl;
}
しかし、それは機能しません。存在しないはずの値を出力するだけです。私がやろうとしていることの詳細なビューが必要な場合は、C 関数と以下の C++ 関数を含めました。
C 関数:
void Setup_Mel(int fft_size, int sample_rate) {
int i,j,k,tap;
double fmax;
double dphi;
double fsample;
double freq;
double temp[fft_size/2];
fmax=2595*log10(8000.0f/700+1);
dphi = fmax/17;
freq = (double)sample_rate/fft_size;
for (i=0; i<16; i++) {
melstart[i]=fft_size/2;
mellength[i]=0;
memset(temp,0,sizeof(double)*fft_size/2);
for (j=0; j<fft_size/2; j++) {
fsample = 2595*log10(freq*j/700 + 1);
if ((dphi*i <= fsample) && (fsample < dphi*(i+1))) temp[j] = (fsample-dphi*i)/(dphi*(i+1)-dphi*i);
if ((dphi*(i+1) <= fsample) && (fsample < dphi*(i+2))) temp[j] = (fsample-dphi*(i+2))/(dphi*(i+1)-dphi*(i+2));
if ((temp[j] != 0) && (melstart[i] > j)) melstart[i] = j;
if (temp[j] != 0) mellength[i]++;
}
mel[i] = malloc(sizeof(double)*mellength[i]);
memcpy(mel[i],&temp[melstart[i]],mellength[i]*sizeof(double));
// for (k=0; k<mellength[i]; tap++,k++) printf("mel filter: %d, %d, %d, %f, %f\n",i,melstart[i]+k,tap,mel[i][k],(melstart[i]+k)*freq);
}
}
C++ コード:
void setUp_Mel(int fft_size, int sample_rate, vector< vector< double > > &mel, int *melstart, int *mellength)
{
double fmax;
double dphi;
double fsample;
double freq;
double temp[fft_size/2];
fmax = 2595*log10(8000.0f/700+1);
dphi = fmax / 17;
freq = (double)sample_rate/fft_size;
for(int i=0; (i < 16); i++)
{
melstart[i] = fft_size/2;
mellength[i] = 0;
double temp[fft_size/2];
//double *temp = new double[fft_size/2];
for(int j=0; (j < fft_size/2); j++)
{
fsample = 2595*log10(freq*j/700 + 1);
if ((dphi*i <= fsample) && (fsample < dphi*(i+1)))
{
temp[j] = (fsample-dphi*i)/(dphi*(i+1)-dphi*i);
}
if ((dphi*(i+1) <= fsample) && (fsample < dphi*(i+2))) {
temp[j] = (fsample-dphi*(i+2))/(dphi*(i+1)-dphi*(i+2));
}
if ((temp[j] != 0) && (melstart[i] > j))
{
melstart[i] = j;
}
if (temp[j] != 0) mellength[i]++;
}
mel[i].resize(mellength[i]);
for(int k=0; (k < mellength[i]); k++)
{
mel[i] = temp[i];
}
}
}