private void waveformGraph1_Paint(object sender, PaintEventArgs e)
{
waveformPlot1.PlotY(waveform);
// Calculate the single-sided, scaled auto power spectrum of signal
// Set the sample period
double dt, df;
dt = 1.0 / (actualPts);
//UnscaledWindow.Blackman(waveform);
//UnscaledWindow.BlackmanHarris(waveform);
//UnscaledWindow.ExactBlackman(waveform);
//UnscaledWindow.FlatTop(waveform);
//UnscaledWindow.Hamming(waveform);
UnscaledWindow.Hanning(waveform);
autoPowerSpectrum = Measurements.AutoPowerSpectrum(waveform, dt, out df);
double maximum = ArrayOperation.GetMax(autoPowerSpectrum);
//maximum = Math.Sqrt(maximum);
for (int i = 0; i < actualPts / 2; i++)
{
autoPowerSpectrum[i] = 10.0 * Math.Log10(Math.Max(autoPowerSpectrum[i] / maximum, 1.0e-15));
}
waveformPlot2.PlotY(autoPowerSpectrum);
}
public override void ToDo(object Arg)
{
b = new BoxController();
b.Init("USB0::0x0957::0x1718::TW54334510::INSTR");
var _ch = new AI_ChannelConfig[4]
{
new AI_ChannelConfig()
{
ChannelName = AnalogInChannelsEnum.AIn1, Enabled = true, Mode = ChannelModeEnum.DC, Polarity = PolarityEnum.Polarity_Bipolar, Range = RangesEnum.Range_1_25
},
new AI_ChannelConfig()
{
ChannelName = AnalogInChannelsEnum.AIn2, Enabled = false, Mode = ChannelModeEnum.DC, Polarity = PolarityEnum.Polarity_Bipolar, Range = RangesEnum.Range_1_25
},
new AI_ChannelConfig()
{
ChannelName = AnalogInChannelsEnum.AIn3, Enabled = false, Mode = ChannelModeEnum.DC, Polarity = PolarityEnum.Polarity_Bipolar, Range = RangesEnum.Range_1_25
},
new AI_ChannelConfig()
{
ChannelName = AnalogInChannelsEnum.AIn4, Enabled = false, Mode = ChannelModeEnum.DC, Polarity = PolarityEnum.Polarity_Bipolar, Range = RangesEnum.Range_1_25
}
};
b.ConfigureAI_Channels(_ch);
var freq = 500000;
var updNumber = 1;
var avgNumber = 1000;
double[] autoPSDLowFreq;
double[] autoPSDHighFreq;
Point[] noisePSD = new Point[] { };
if (freq % 2 != 0)
{
throw new ArgumentException("The frequency should be an even number!");
}
b.AcquisitionInProgress = true;
b.AI_ChannelCollection[AnalogInChannelsEnum.AIn1].DataReady += DefResistanceNoise_DataReady;
var sb = new StringBuilder();
double dtLowFreq = 0.0, dtHighFreq = 0.0;
double dfLowFreq = 1.0, dfHighFreq = 0.0;
double equivalentNoiseBandwidthLowFreq, equivalentNoiseBandwidthHighFreq;
double coherentGainLowFreq, coherentGainHighFreq;
Parallel.Invoke(
() =>
{
b.StartAnalogAcquisition(freq);
IsRunning = false;
},
() =>
{
while (true)
{
if (!IsRunning)
{
b.AcquisitionInProgress = false;
break;
}
if (averagingCounter >= avgNumber)
{
b.AcquisitionInProgress = false;
break;
}
Point[] timeTrace;
var dataReadingSuccess = b.AI_ChannelCollection[AnalogInChannelsEnum.AIn1].ChannelData.TryDequeue(out timeTrace);
if (dataReadingSuccess)
{
var query = from val in timeTrace
select val.Y;
var counter = 0;
var traceData = new double[timeTrace.Length];
foreach (var item in query)
{
traceData[counter] = item;
++counter;
}
var unit = new System.Text.StringBuilder("V", 256);
var sw = ScaledWindow.CreateRectangularWindow();
// Calculation of the low-frequency part of the spectrum
sw.Apply(traceData, out equivalentNoiseBandwidthLowFreq, out coherentGainLowFreq);
dtLowFreq = 1.0 / (double)freq;
autoPSDLowFreq = Measurements.AutoPowerSpectrum(traceData, dtLowFreq, out dfLowFreq);
var singlePSDLowFreq = Measurements.SpectrumUnitConversion(autoPSDLowFreq, SpectrumType.Power, ScalingMode.Linear, DisplayUnits.VoltsPeakSquaredPerHZ, dfLowFreq, equivalentNoiseBandwidthLowFreq, coherentGainLowFreq, unit);
// Calculation of the hugh-frequency part of the spectrum
var selection64Hz = PointSelector.SelectPoints(ref traceData, 64);
sw.Apply(selection64Hz, out equivalentNoiseBandwidthHighFreq, out coherentGainHighFreq);
dtHighFreq = 64.0 * 1.0 / (double)freq;
autoPSDHighFreq = Measurements.AutoPowerSpectrum(selection64Hz, dtHighFreq, out dfHighFreq);
var singlePSDHighFreq = Measurements.SpectrumUnitConversion(autoPSDHighFreq, SpectrumType.Power, ScalingMode.Linear, DisplayUnits.VoltsPeakSquaredPerHZ, dfHighFreq, equivalentNoiseBandwidthHighFreq, coherentGainHighFreq, unit);
var lowFreqSpectrum = singlePSDLowFreq.Select((value, index) => new Point((index + 1) * dfLowFreq, value)).Where(value => value.X <= 1064);
var highFreqSpectrum = singlePSDLowFreq.Select((value, index) => new Point((index + 1) * dfHighFreq, value)).Where(value => value.X > 1064);
noisePSD = new Point[lowFreqSpectrum.Count() + highFreqSpectrum.Count()];
counter = 0;
foreach (var item in lowFreqSpectrum)
{
noisePSD[counter].X = item.X;
noisePSD[counter].Y += item.Y;
++counter;
}
foreach (var item in highFreqSpectrum)
{
noisePSD[counter].X = item.X;
noisePSD[counter].Y += item.Y;
++counter;
}
//for (int i = 0; i < singlePSDLowFreq.Length; i++)
// noisePSD[i] += singlePSDLowFreq[i];
if (averagingCounter % updNumber == 0)
{
sb = new StringBuilder();
for (int i = 0; i < noisePSD.Length; i++)
{
sb.AppendFormat("{0}\t{1}\r\n", (noisePSD[i].X).ToString(NumberFormatInfo.InvariantInfo), (noisePSD[i].Y / (double)averagingCounter).ToString(NumberFormatInfo.InvariantInfo));
}
onDataArrived(new ExpDataArrivedEventArgs(sb.ToString()));
}
}
}
//sb = new StringBuilder();
//for (int i = 0; i < noisePSD.Length; i++)
// sb.AppendFormat("{0}\t{1}\r\n", (noisePSD[i].X).ToString(NumberFormatInfo.InvariantInfo), (noisePSD[i].Y / (double)averagingCounter).ToString(NumberFormatInfo.InvariantInfo));
//onDataArrived(new ExpDataArrivedEventArgs(sb.ToString()));
});
b.Close();
}
public Point[] GetTwoPartsFFT(double[] timeTrace, int samplingFrequency = 500000, int nDataSamples = 1, double kAmpl = 1.0, double lowFreqStartFreq = 1.0, double cutOffLowFreq = 1600, double cutOffHighFreq = 102400, int filterOrder = 8, double filterFrequency = 6400, int lowFreqPeriod = 10, int highFreqPeriod = 10)
{
Point[] autoPSDLowFreq = new Point[] { };
Point[] autoPSDHighFreq = new Point[] { };
if (filterFrequency == -1)
{
filterFrequency = cutOffLowFreq;
}
var sb = new StringBuilder();
double dtLowFreq = 0.0, dtHighFreq = 0.0;
double dfLowFreq = 0.0, dfHighFreq = 0.0;
double equivalentNoiseBandwidthLowFreq, equivalentNoiseBandwidthHighFreq;
double coherentGainLowFreq, coherentGainHighFreq;
// Subsetting samples from the entire trace
var timeTraceSelectionList = new LinkedList <double[]>();
var range = (int)((timeTrace.Length) / nDataSamples);
for (int i = 0; i != nDataSamples;)
{
var selection = timeTrace.Where((value, index) => index >= i * range && index < (i + 1) * range).Select(val => val);
timeTraceSelectionList.AddLast(selection.ToArray());
++i;
}
var filter = new NationalInstruments.Analysis.Dsp.Filters.EllipticLowpassFilter(filterOrder, samplingFrequency, filterFrequency, 0.1, 100.0);
var unit = new System.Text.StringBuilder("V", 256);
var noisePSD = new Point[] { };
// Calculating FFT in each sample
foreach (var trace in timeTraceSelectionList)
{
// Calculation of the LOW-FREQUENCY part of the spectrum
// Filtering data for low frequency selection
var filteredData = filter.FilterData(trace);
var sw = ScaledWindow.CreateRectangularWindow();
// Selecting lower amount of data points to reduce the FFT noise
var selectionLowFreq = filteredData.Where((value, index) => (index) % lowFreqPeriod == 0).ToArray();
sw.Apply(selectionLowFreq, out equivalentNoiseBandwidthLowFreq, out coherentGainLowFreq);
dtLowFreq = lowFreqPeriod * 1.0 / (double)samplingFrequency;
var singlePSD_LOW_Freq = Measurements.AutoPowerSpectrum(selectionLowFreq, dtLowFreq, out dfLowFreq);
autoPSDLowFreq = (singlePSD_LOW_Freq.Select((value, index) => new Point(index * dfLowFreq, value)).Where(p => p.X >= 1 && p.X <= cutOffLowFreq)).ToArray();
// Calculation of the HIGH-FREQUENCY part of the spectrum
dtHighFreq = 1.0 / (double)samplingFrequency;
var highFreqSelectionRange = (int)((timeTrace.Length) / highFreqPeriod);
LinkedList <double[]> selectionList = new LinkedList <double[]>();
for (int i = 0; i != highFreqPeriod;)
{
var arr = new double[highFreqSelectionRange];
Array.Copy(timeTrace, i * highFreqSelectionRange, arr, 0, highFreqSelectionRange);
selectionList.AddLast(arr);
++i;
}
var cumulativePSD_HIGH_Freq = new double[] { };
foreach (var selection in selectionList)
{
sw.Apply(selection, out equivalentNoiseBandwidthHighFreq, out coherentGainHighFreq);
var singlePSD_HIGH_Freq = Measurements.AutoPowerSpectrum(selection, dtHighFreq, out dfHighFreq);
if (cumulativePSD_HIGH_Freq.Length == 0)
{
cumulativePSD_HIGH_Freq = Enumerable.Repeat(0.0, singlePSD_HIGH_Freq.Length).ToArray();
}
for (int i = 0; i != singlePSD_HIGH_Freq.Length;)
{
cumulativePSD_HIGH_Freq[i] += singlePSD_HIGH_Freq[i];
++i;
}
}
autoPSDHighFreq = (cumulativePSD_HIGH_Freq
.Select((value, index) => new Point(index * dfHighFreq, value)))
.Where(p => p.X > cutOffLowFreq && p.X <= cutOffHighFreq).ToArray();
if (noisePSD == null || noisePSD.Length == 0)
{
noisePSD = new Point[autoPSDLowFreq.Length + autoPSDHighFreq.Length];
}
var counter = 0;
for (int i = 0; i != autoPSDLowFreq.Length;)
{
var item = autoPSDLowFreq[i];
noisePSD[counter].X = item.X;
noisePSD[counter].Y += item.Y;
++counter;
++i;
}
for (int i = 0; i != autoPSDHighFreq.Length;)
{
var item = autoPSDHighFreq[i];
noisePSD[counter].X = item.X;
noisePSD[counter].Y += item.Y / ((double)(highFreqPeriod * highFreqPeriod));
++counter;
++i;
}
}
return((from item in noisePSD
select new Point(item.X, item.Y / (kAmpl * kAmpl))).ToArray());
}
