speech_recognition.rar_speech recognition_speech-recognition

/* * endsil.c * * Code generation for function 'endsil' * * C source code generated on: Wed Mar 02 11:54:34 2016 * */ /* Include files */ #include "rt_nonfinite.h" #include "endsil.h" #include "endsil_emxutil.h" /* Function Definitions */ void endsil(const emxArray_real_T *input, real32_T fs, emxArray_real_T *out) { emxArray_real_T *bgSample; real32_T samplePerFrame; real32_T bgSampleCount; int32_T k; int32_T d; int32_T xs; real_T cUnVoiced; real_T meanVal; int32_T ix; real_T xbar; real_T r; int8_T voicedUnvoiced[30000]; int8_T voiced[30000]; real32_T f0; int32_T i0; emxInit_real_T(&bgSample, 2); samplePerFrame = (real32_T)floor(fs / 100.0F); bgSampleCount = (real32_T)floor(fs / 5.0F); /* according to formula, 1600 sample needed for 8 khz */ /* ---------- */ /* calculation of mean and std */ k = bgSample->size[0] * bgSample->size[1]; bgSample->size[0] = 1; bgSample->size[1] = 0; emxEnsureCapacity((emxArray__common *)bgSample, k, (int32_T)sizeof(real_T)); for (d = 0; d < (int32_T)bgSampleCount; d++) { xs = bgSample->size[1]; k = bgSample->size[0] * bgSample->size[1]; bgSample->size[1] = xs + 1; emxEnsureCapacity((emxArray__common *)bgSample, k, (int32_T)sizeof(real_T)); bgSample->data[xs] = input->data[(int32_T)(1.0F + (real32_T)d) - 1]; } if (bgSample->size[1] == 0) { cUnVoiced = 0.0; } else { cUnVoiced = bgSample->data[0]; for (k = 2; k <= bgSample->size[1]; k++) { cUnVoiced += bgSample->data[k - 1]; } } meanVal = cUnVoiced / (real_T)bgSample->size[1]; if (bgSample->size[1] > 1) { d = bgSample->size[1] - 1; } else { d = bgSample->size[1]; } if (bgSample->size[1] == 0) { cUnVoiced = rtNaN; } else { ix = 0; xbar = bgSample->data[0]; for (k = 0; k <= bgSample->size[1] - 2; k++) { ix++; xbar += bgSample->data[ix]; } xbar /= (real_T)bgSample->size[1]; ix = 0; r = bgSample->data[0] - xbar; cUnVoiced = r * r; for (k = 0; k <= bgSample->size[1] - 2; k++) { ix++; r = bgSample->data[ix] - xbar; cUnVoiced += r * r; } cUnVoiced /= (real_T)d; } emxFree_real_T(&bgSample); xbar = sqrt(cUnVoiced); /* ---------- */ /* identify voiced or not for each value */ for (k = 0; k < 30000; k++) { voicedUnvoiced[k] = 0; voiced[k] = 0; } for (d = 0; d < input->size[0]; d++) { cUnVoiced = input->data[d] - meanVal; if (fabs(cUnVoiced) / xbar > 0.3) { voiced[d] = 1; } else { voiced[d] = 0; } } /* identify voiced or not for each frame */ /* discard insufficient samples of last frame */ if (samplePerFrame == 0.0F) { bgSampleCount = (real32_T)input->size[0]; } else if (samplePerFrame == samplePerFrame) { bgSampleCount = (real32_T)input->size[0] - (real32_T)floor((real32_T) input->size[0] / samplePerFrame) * samplePerFrame; } else { bgSampleCount = (real32_T)input->size[0] / samplePerFrame; bgSampleCount = (bgSampleCount - (real32_T)floor(bgSampleCount)) * samplePerFrame; } bgSampleCount = ((real32_T)input->size[0] - bgSampleCount) / samplePerFrame; for (d = 0; d < (int32_T)bgSampleCount; d++) { xbar = 0.0; cUnVoiced = 0.0; f0 = ((1.0F + (real32_T)d) * samplePerFrame - samplePerFrame) + 1.0F; k = (int32_T)((1.0F + (real32_T)d) * samplePerFrame + (1.0F - f0)); for (ix = 0; ix < k; ix++) { if (voiced[(int32_T)(f0 + (real32_T)ix) - 1] == 1) { xbar++; } else { cUnVoiced++; } } /* mark frame for voiced/unvoiced */ if (xbar > cUnVoiced) { voicedUnvoiced[(int32_T)(real32_T)(1 + d) - 1] = 1; } else { voicedUnvoiced[(int32_T)(real32_T)(1 + d) - 1] = 0; } } k = out->size[0] * out->size[1]; out->size[0] = 1; out->size[1] = 0; emxEnsureCapacity((emxArray__common *)out, k, (int32_T)sizeof(real_T)); /* ----- */ for (d = 0; d < (int32_T)bgSampleCount; d++) { if (voicedUnvoiced[(int32_T)(real32_T)(1 + d) - 1] == 1) { f0 = ((1.0F + (real32_T)d) * samplePerFrame - samplePerFrame) + 1.0F; k = (int32_T)((1.0F + (real32_T)d) * samplePerFrame + (1.0F - f0)); for (ix = 0; ix < k; ix++) { xs = out->size[1]; i0 = out->size[0] * out->size[1]; out->size[1] = xs + 1; emxEnsureCapacity((emxArray__common *)out, i0, (int32_T)sizeof(real_T)); out->data[xs] = input->data[(int32_T)(f0 + (real32_T)ix) - 1]; } } } /* ---display plot and play both sounds */ /* figure; plot(y); */ /* figure; plot(silenceRemovedSignal); */ /* %%play */ /* sound(y,Fs); */ /* sound(silenceRemovedSignal,Fs); */ } /* End of code generation (endsil.c) */