static int power_surge_detector_file_test(const char *file) { SNDFILE *inhandle; SNDFILE *outhandle; int inframes; power_surge_detector_state_t *sig; int i; int16_t amp[8000]; int16_t amp_out[2*8000]; int sample; float signal_power; int32_t signal_level; int signal_present; int prev_signal_present; if ((inhandle = sf_open_telephony_read(file, 1)) == NULL) { printf(" Cannot open speech file '%s'\n", file); exit(2); } if ((outhandle = sf_open_telephony_write(OUT_FILE_NAME, 1)) == NULL) { fprintf(stderr, " Cannot create audio file '%s'\n", OUT_FILE_NAME); exit(2); } sig = power_surge_detector_init(NULL, -50.0f, 6.0f); prev_signal_present = FALSE; sample = 0; while ((inframes = sf_readf_short(inhandle, amp, 8000))) { for (i = 0; i < inframes; i++) { signal_level = power_surge_detector(sig, amp[i]); signal_present = (signal_level != 0); if (prev_signal_present != signal_present) { signal_power = power_surge_detector_current_dbm0(sig); if (signal_present) printf("On at %f (%fdBm0)\n", (sample + i)/8000.0, signal_power); else printf("Off at %f (%fdBm0)\n", (sample + i)/8000.0, signal_power); prev_signal_present = signal_present; } amp_out[2*i] = amp[i]; amp_out[2*i + 1] = signal_present*5000; } sf_writef_short(outhandle, amp_out, inframes); sample += inframes; } sf_close(inhandle); sf_close(outhandle); return 0; }
static int power_surge_detector_tests(void) { SNDFILE *outhandle; power_surge_detector_state_t *sig; int i; int sample; int16_t amp[8000]; int16_t amp_out[2*8000]; awgn_state_t *awgnx; int32_t phase_rate; uint32_t phase_acc; int16_t phase_scale; float signal_power; int32_t signal_level; int signal_present; int prev_signal_present; int ok; int extremes[4]; if ((outhandle = sf_open_telephony_write(OUT_FILE_NAME, 2)) == NULL) { fprintf(stderr, " Cannot create audio file '%s'\n", OUT_FILE_NAME); exit(2); } sig = power_surge_detector_init(NULL, -50.0f, 5.0f); prev_signal_present = false; phase_rate = dds_phase_rate(450.0f); phase_acc = 0; phase_scale = dds_scaling_dbm0(-33.0f); awgnx = awgn_init_dbm0(NULL, 1234567, -45.0f); extremes[0] = 8001; extremes[1] = -1; extremes[2] = 8001; extremes[3] = -1; for (sample = 0; sample < 800000; sample += 8000) { ok = 0; for (i = 0; i < 8000; i++) { amp[i] = awgn(awgnx); if (i < 4000) amp[i] += dds_mod(&phase_acc, phase_rate, phase_scale, 0); signal_level = power_surge_detector(sig, amp[i]); signal_present = (signal_level != 0); if (prev_signal_present != signal_present) { signal_power = power_surge_detector_current_dbm0(sig); if (signal_present) { if (ok == 0 && i >= 0 && i < 25) ok = 1; if (extremes[0] > i) extremes[0] = i; if (extremes[1] < i) extremes[1] = i; printf("On at %f (%fdBm0)\n", (sample + i)/8000.0, signal_power); } else { if (ok == 1 && i >= 4000 + 0 && i < 4000 + 35) ok = 2; if (extremes[2] > i) extremes[2] = i; if (extremes[3] < i) extremes[3] = i; printf("Off at %f (%fdBm0)\n", (sample + i)/8000.0, signal_power); } prev_signal_present = signal_present; } amp_out[2*i] = amp[i]; amp_out[2*i + 1] = signal_present*5000; } sf_writef_short(outhandle, amp_out, 8000); if (ok != 2 || extremes[0] < 1 || extremes[1] > 30 || extremes[2] < 4001 || extremes[3] > 4030) { printf(" Surge not detected correctly (%d)\n", ok); exit(2); } } if (sf_close_telephony(outhandle)) { fprintf(stderr, " Cannot close audio file '%s'\n", OUT_FILE_NAME); exit(2); } printf("Min on %d, max on %d, min off %d, max off %d\n", extremes[0], extremes[1], extremes[2], extremes[3]); power_surge_detector_free(sig); awgn_free(awgnx); return 0; }