float Functions::AM(int t)
{



	return (1 + kAM*First(t))*sin(2 * M_PI*fm*t);
}

float Functions::PM(int t)
{



	return sin(2 * M_PI*fm*t + kPM*First(t));

}


////////////////////////////////w innej funkcji
case 4: //AM
			for (int t = 0; t < N; t++)
			{
				xAxis[t] = AM(t);
			}
			break;
		case 5: //PM
			for (int t = 0; t < N; t++)
			{
				xAxis[t] = PM(t);
			}
			break;
		}
}

///////////////////////////////////////
void Functions::Power_Spectrum()
{
	float SumA = 0, SumB = 0, a, b;
	float max = 0;
	float *temp = new float[N];
	for (int n = 0; n < N; n++)
	{
		temp[n] = xAxis[n];
		if (temp[n]>max)
			max = temp[n];
	}
	yAxis = new float[N];
	for (int k = 0; k < N; k++)
	{
		SumA = SumB = 0;
		for (int n = 0; n < N; n++)
		{
			a = cos((2 * M_PI*n*k) / N);
			b = sin((2 * M_PI*n*k) / N);
			SumA += temp[n] * a;
			SumB += temp[n] * b;
		}
		xAxis[k] = (SumA*SumA) + (SumB*SumB);
		xAxis[k] = 10 * log10f(xAxis[k]);
		xAxis[k] = xAxis[k] - max;
		yAxis[k] = (float)k*((float)fs / (float)N);
	}
	delete(temp);
}