[C++] Синусоида убывающей частоты

[ZenBooster]

Пишу сирену на с++ Столкнулся с тем, что на отрезке от верхней частоты к нижней, синусоида замедляется до нулевой частоты, а затем частота вновь растет до некоторого значения! С чем это может быть связано, и как это побороть?

#define FREQUENCY 44100
#define MAX_SAMPLE_VAL 0xff

#define FMIN 400.0
#define FMAX 1400.0

cWail::cWail(void)
{
	double i, j;

	//int dt_i = 16;
	//int dt_j = 4;
	double dt_i = 0.04;
	int max_j = 1;

	sz_pdata += 2 * (max_j * (FMAX - FMIN + dt_i) / dt_i);
	p_data = new char[sz_pdata];
	ZeroMemory(p_data, sz_pdata);
	char *p = p_data;

	double t = 0;
	/*for (i = FMIN; i < FMAX; i+= dt_i)
	{
		for (int j = 0; j < max_j; j++)
		{
			double v = sin(2.0 * PI * t / ((double)FREQUENCY / i));
			//*p++ = (v > 0) ? MAX_SAMPLE_VAL : 0; // meander
			*p++ = (1 + v) * 0.5 * MAX_SAMPLE_VAL; // sine
			t += 1.0;
		}
	}*/

	for (i = FMAX; i > FMIN; i-= dt_i)
	{
		for (int j = 0; j < max_j; j++)
		{
			double v = sin(2.0 * PI * t / ((double)FREQUENCY / i));
			//*p++ = (v > 0) ? MAX_SAMPLE_VAL : 0; // meander
			*p++ = (1+v) * 0.5 * MAX_SAMPLE_VAL; // sine
			t += 1.0;
		}
	}
	/*FILE *f;
	errno_t e = fopen_s(&f, "sine.bin", "wb");
	if (!e)
	{
		fwrite(p_data, 1, sz_pdata, f);
		fclose(f);
	}*/
	sound_data.SetSource(p_data, 0, sz_pdata);
}

[digitalis]

Я не вижу здесь на осциллограмме нулевой частоты. Уменьшилась до минимальной и затем в точке, близкой к минимуму, начала зеркально увеличиваться (если спроецировать на механику, то колесо в нижней точке замерло на миг и начало повораитваться в обратную сторону).

[ZenBooster]

В итоге вот так все работает:

cWail::cWail(void)
{
	double i, j;

	double dt_i = 0.02;
	double phase = 0.0;
	int t = 0;

	//sz_pdata += 2 * ((FMAX - FMIN + dt_i) / dt_i);
	for (i = FMIN; i <= FMAX; i += dt_i)
	{
		double v = sin(phase);
		phase = fmod((phase + 2.0 * PI * i / (double)FREQUENCY), (2 * PI));
		t++;
	}
	for (i = FMAX; i >= FMIN || phase > 0.01; i -= dt_i)
	{
		double v = sin(phase);
		phase = fmod((phase + 2.0 * PI * i / (double)FREQUENCY), (2 * PI));
		t++;
	}
	sz_pdata = t;

	p_data = new char[sz_pdata];
	ZeroMemory(p_data, sz_pdata);
	char *p = p_data;
	phase = 0.0;
	t = 0;

	for (i = FMIN; i <= FMAX; i+= dt_i)
	{
		double v = sin(phase);
		phase = fmod((phase + 2.0 * PI * i / (double)FREQUENCY), (2 * PI));
		//*p++ = (v > 0) ? MAX_SAMPLE_VAL : 0; // meander
		*p++ = (1 + v) * 0.5 * MAX_SAMPLE_VAL; // sine
		t++;
	}
	for (i = FMAX; i >= FMIN || phase > 0.01; i-= dt_i)
	{
		double v = sin(phase);
		phase = fmod((phase + 2.0 * PI * i / (double)FREQUENCY), (2 * PI));
		//*p++ = (v > 0) ? MAX_SAMPLE_VAL : 0; // meander
		*p++ = (1 + v) * 0.5 * MAX_SAMPLE_VAL; // sine
		t++;
	}

	/*FILE *f;
	errno_t e = fopen_s(&f, "sine.bin", "wb");
	if (!e)
	{
		fwrite(p_data, 1, sz_pdata, f);
		fclose(f);
	}*/
	sound_data.SetSource(p_data, 0, sz_pdata);
}

Спектр: