转载URL: https://blog.csdn.net/harriszhuang/article/details/94642872
/* 这是一个简单的信号发生器,用于产生符合CSMPTE标准的千周信号 */
/* 左声道:1kHz,每间隔3秒间断约0.4秒,即出现2.6秒,静音0.4秒 */
/* 右声道:1kHz,连续 */
/* 输出格式:PCM_S16LE,16比特,低位在前,高位在后 */
/* C-Free 4.0 mingw 2.95 */
#include <math.h>
#include <stdio.h>
#define SAMPLE_RATE 48000
#define LEN 60.0
#define GATE_OFF (0.4 * SAMPLE_RATE)
#define GATE_LOOP_LEN (3 * SAMPLE_RATE)
#define GATE_ON (GATE_LOOP_LEN - GATE_OFF)
#define FALSE 0
#define TRUE 1
#define INT16 1
#define GAIN_L_DB -20
#define GAIN_R_DB -20
short sp_to_short_word(double sp);
double dB_to_decimal(double x);
int cnt_max = (double)LEN * SAMPLE_RATE;
const char *Path = "D:\\work\\abc.pcm";
const double Pi = 3.1415926535897932384626433832795;
const double Left_freq = 1000;
const double Right_freq = 1000;
const double Left_gain = 0.1;
const double Right_gain = 0.1;
int main()
{
int sp_cnt = 0; /* Running on mingw2.95 */
int gate_cnt = 0;
int bytes_cnt = 0;
int tmp1 = 0;
double sp_left = 0, sp_right = 0;
double left_freq = Left_freq, right_freq = Right_freq;
double left_gain = dB_to_decimal(GAIN_L_DB);
double right_gain = dB_to_decimal(GAIN_R_DB);
double len = LEN;
signed short sp_out[2] = { 0, 0 };
FILE *fp1;
fp1 = fopen(Path, "wb+");
if (fp1 == NULL) {
printf("Error in opening output file.\n");
return -1;
}
printf("[Debug] cnt_max = %ld\n", cnt_max);
printf("[Debug] len = %f\n", len);
printf("[Debug] left_freq = %f\n", left_freq);
printf("[Debug] right_freq = %f\n", right_freq);
printf("[Debug] left_gain = %f\n", left_gain);
printf("[Debug] right_gain = %f\n", right_gain);
printf("[Debug] cos 30deg = %.16f\n", cos(Pi / 6));
printf("[Debug] cos 45deg = %.16f\n", cos(Pi / 4));
printf("[Debug] cos 60deg = %.16f\n", cos(Pi / 3));
printf("[Debug] sizeof(signed short) = %d\n", sizeof(signed short));
printf("[Debug] sizeof(signed int) = %d\n", sizeof(signed int));
printf("[Debug] writing out samples...\n");
for (sp_cnt = 0; sp_cnt < cnt_max; sp_cnt++)
{
double t;
t = (double)sp_cnt / 48000;
sp_left = cos(Pi * 2 * left_freq * t) * left_gain;
sp_right = cos(Pi * 2 * right_freq * t) * right_gain;
if (TRUE && gate_cnt >= GATE_ON)
sp_left = 0;
if (FALSE && gate_cnt >= GATE_ON)
sp_right = 0;
#ifdef INT16
sp_out[0] = sp_to_short_word(sp_left);
sp_out[1] = sp_to_short_word(sp_right);
tmp1 = fwrite(sp_out, sizeof(signed short), 2, fp1);
bytes_cnt += tmp1;
#endif
#ifdef IEEE64
tmp1 = fwrite(&sp_left, sizeof(double), 1, fp1);
bytes_cnt += tmp1;
tmp1 = fwrite(&sp_right, sizeof(double), 1, fp1);
bytes_cnt += tmp1;
#endif
gate_cnt = (gate_cnt >= GATE_LOOP_LEN ? 0 : gate_cnt + 1);
} /* for */
printf("[Debug] bytes_cnt = %ld\n", bytes_cnt);
fclose(fp1);
} /* main */
short sp_to_short_word(double sp)
{
signed short ret = 0;
ret = sp * 32767;
return ret;
}
double dB_to_decimal(double x)
{
double y = pow(10, x / 20);
return y;
}
/*
long sp_to_int_word(double sp)
{
signed int ret = 0;
ret = sp * 2147483647L;
}
*/