internal Disturbance(int start, int finish, int detailsLength, int signalLength)
{
Start = start;
Finish = finish;
SignalStart = (int)WaveMath.Scale(start, 0, detailsLength, 0, signalLength) + 1;
SignalFinish = (int)WaveMath.Scale(finish, 0, detailsLength, 0, signalLength) + 1;
if (start % 2 != 0 && finish % 2 != 0)
{
SignalFinish--;
}
}
/// <summary>
/// Executes the block
/// </summary>
public override void Execute()
{
var inputNode = InputNodes[0].ConnectingNode as BlockOutputNode;
if (inputNode == null || inputNode.Object == null)
{
return;
}
OutputNodes[0].Object.Clear();
foreach (var signal in inputNode.Object)
{
var output = WaveMath.Scale(signal, _amplitudeScalingFactor, _timeScalingFactor);
OutputNodes[0].Object.Add(output);
}
if (Cascade && OutputNodes[0].ConnectingNode != null)
{
OutputNodes[0].ConnectingNode.Root.Execute();
}
}
/// <summary>
/// Estimates the disturbances in the signal based on the normal distribution of the details coefficients
/// </summary>
/// <param name="threshold">The higher the threshold, the higher the tolerance in flutuations on energy of the details</param>
/// <param name="minimunDistance">Minimun distance between disturbances to consider a new one</param>
/// <returns></returns>
public List <Disturbance> GetDisturbances(double threshold = 0.1, int minimunDistance = 3)
{
threshold = 1 - threshold;
var disturbances = new List <Disturbance>();
var mean = WaveMath.Mean(Details);
var deviation = WaveMath.StandardDeviation(Details);
var samples = new List <KeyValuePair <int, double> >();
var min = double.MaxValue;
var max = double.MinValue;
for (var i = 0; i < Details.Length; i++)
{
var sample = Details[i];
var norm = WaveMath.ProbabilityDensityFunction(sample, mean, deviation);
if (double.IsNaN(norm))
{
continue;
}
if (norm < min)
{
min = norm;
}
if (norm > max)
{
max = norm;
}
samples.Add(new KeyValuePair <int, double>(i, norm));
}
//ajusta a escala da distribuição normal para 0..1, removendo os valores maiores que threshold
for (var i = samples.Count - 1; i >= 0; i--)
{
var scaledNorm = WaveMath.Scale(samples[i].Value, min, max, 0, 1);
samples[i] = new KeyValuePair <int, double>(i, scaledNorm);
if (scaledNorm > threshold)
{
samples.RemoveAt(i);
}
else
{
samples[i] = new KeyValuePair <int, double>(i, scaledNorm);
}
}
int?start = null;
var startIndex = 0;
for (var i = 0; i < samples.Count; i++)
{
if (start == null)
{
start = i;
startIndex = samples[i].Key;
}
else if (samples[i].Key - samples[i - 1].Key >= minimunDistance || i == samples.Count - 1)
{
disturbances.Add(new Disturbance(startIndex, samples[i - 1].Key, Details.Length, Signal.Samples.Length));
if (i < samples.Count - 1)
{
start = i;
startIndex = samples[i].Key;
}
else
{
disturbances.Add(new Disturbance(samples[i].Key, samples[i].Key, Details.Length, Signal.Samples.Length));
}
}
}
if (disturbances.Count > 1 && (disturbances[disturbances.Count - 1].Finish - disturbances[disturbances.Count - 2].Finish) < minimunDistance)
{
disturbances[disturbances.Count - 2] = new Disturbance(disturbances[disturbances.Count - 2].Start, disturbances[disturbances.Count - 1].Finish, Details.Length, Signal.Samples.Length);
disturbances.RemoveAt(disturbances.Count - 1);
}
return(disturbances);
}