void m_queue_TEvent(DataAndFreq t) { IQDataAndFreq nIQDataAndFreq = null; int dataLength = t.Ai_Data.Length / 20 * 2; //(重要)保持多线程数据一致行所做的必要处理 int mm_Resolution = m_Resolution, Resolution_length = mm_Resolution * 512; if (Resolution_length != IQDataResolution.Length) { ResolutionIndex = 0; IQDataResolution = new double[Resolution_length]; IQDataResolutionOrg = new double[Resolution_length]; } int BaseIndex = dataLength * ResolutionIndex; for (int i = 0; i < dataLength / 2; ++i) { IQDataResolutionOrg[BaseIndex + 2 * i] = (double)t.Ar_Data[i]; IQDataResolutionOrg[BaseIndex + 2 * i + 1] = (double)t.Ai_Data[i]; } if (++ResolutionIndex < mm_Resolution)//凑包个数计数 { return; } ResolutionIndex = 0; //对不同分辨率(长度)的数据进行加窗处理 int mm_Resolution_Index = (int)Math.Log(mm_Resolution, 2); for (int i = 0; i < IQDataResolution.Length / 2; i++) { IQDataResolution[2 * i] = IQDataResolutionOrg[2 * i] * Window[mm_Resolution_Index][i]; IQDataResolution[2 * i + 1] = IQDataResolutionOrg[2 * i + 1] * Window[mm_Resolution_Index][i]; } nIQDataAndFreq = new IQDataAndFreq(t.StartFreq, t.StopFreq, t.Type, t.RF_Gain, t.Digital_Gain);//构造结构体 nIQDataAndFreq.Data = IQDataResolution; ///分发至时域处理类 DataIQWaveform.PushData(IQDataResolutionOrg); //进入此处必然处理下一包 需要为下一包申请指向新空间 IQDataResolution = new double[Resolution_length]; IQDataResolutionOrg = new double[Resolution_length]; ///分发至频域处理类 DataFrequencySpectrum.PushData(nIQDataAndFreq); }
void m_queue_TEvent(IQDataAndFreq t) { double[] outData = new double[t.Data.Length]; int mm_SmoothNum = m_SmoothNum; double offset = 0; int m_RFGain; int m_DigitalGain; if (t.Type == 3) { Array.Copy(t.Data, 0, outData, 0, t.Data.Length); for (int k = 0; k < outData.Length; k++) { IQdata[k][SmoothIndex] = outData[k]; } offset = 133 - 11.85; } else { hdin = GCHandle.Alloc(t.Data, GCHandleType.Pinned); hdout = GCHandle.Alloc(outData, GCHandleType.Pinned); dplan = fftw.dft_1d(t.Data.Length / 2, hdin.AddrOfPinnedObject(), hdout.AddrOfPinnedObject(), fftw_direction.Forward, fftw_flags.Measure); fftw.execute(dplan); fftw.destroy_plan(dplan); hdin.Free(); hdout.Free(); for (int k = 0; k < outData.Length / 2; k++) { double dou1 = Math.Pow(outData[k * 2], 2.0); double dou2 = Math.Pow(outData[k * 2 + 1], 2.0); double dou3 = dou1 + dou2; IQdata[k][SmoothIndex] = Math.Pow(dou3, 0.5); } int Resolution_Offset = (int)(6 * (Math.Log(t.Data.Length, 2) - 10)); if (t.Type == 2) { offset = 133 - 7 + Resolution_Offset; } else if (t.Type == 1) { offset = 133 + 8.63 + Resolution_Offset + (m_DataProcessing.DigitialGain_24db ? 24 : 0) + (m_DataProcessing.LowNoise ? 12 : 0); } } double[] power = null; if (t.Type == 3) { power = new double[outData.Length]; for (int i = 0; i < outData.Length; i++) { double sum = 0; for (int j = 0; j < mm_SmoothNum; j++) { sum += IQdata[i][j]; } sum /= mm_SmoothNum; m_DigitalGain = Convert.ToInt32(t.Digital_Gain[i * 57 / outData.Length]); m_RFGain = Convert.ToInt32(t.RF_Gain[i / outData.Length]); power[i] = sum / 10 - offset + m_DigitalGain + m_RFGain; } } else if (t.Type == 2) { power = new double[outData.Length / 2]; for (int i = 0; i < outData.Length / 2; i++) { double sum = 0; for (int j = 0; j < mm_SmoothNum; j++) { sum += IQdata[i][j]; } sum /= mm_SmoothNum; m_DigitalGain = Convert.ToInt32(t.Digital_Gain[i * 40 / outData.Length]); m_RFGain = Convert.ToInt32(t.RF_Gain[0]); power[i] = Math.Log10(sum) * 20 - offset + m_DigitalGain + m_RFGain; } } else { power = new double[outData.Length / 2]; for (int i = 0; i < outData.Length / 2; i++) { double sum = 0; for (int j = 0; j < mm_SmoothNum; j++) { sum += IQdata[i][j]; } sum /= mm_SmoothNum; power[i] = Math.Log10(sum) * 20 - offset; } } //平滑处理 if (++SmoothIndex >= mm_SmoothNum) { SmoothIndex = 0; } //数据传递 PowerAndFreq nPowerAndFreq = new PowerAndFreq(); power.Reverse <double>();//倒谱处理 long powerLength = power.Length / 2; if (t.Type == 3) { nPowerAndFreq.Power = power; } else { nPowerAndFreq.Power = new double[power.Length]; Array.Copy(power, powerLength, nPowerAndFreq.Power, 0, powerLength); Array.Copy(power, 0, nPowerAndFreq.Power, powerLength, powerLength); } switch (t.Type) { case 1: { if (m_DataProcessing.nbddcbandwidth >= 1 && m_DataProcessing.nbddcbandwidth <= 6) { double off = 0.00425 - (m_DataProcessing.nbddcbandwidth - 1) * 10 * 0.00005; nPowerAndFreq.StartFreq = t.StartFreq - off; nPowerAndFreq.StopFreq = t.StopFreq + off; } else if (m_DataProcessing.nbddcbandwidth >= 9 && m_DataProcessing.nbddcbandwidth <= 15) { double off = 0.005 - (m_DataProcessing.nbddcbandwidth - 9) * 0.0005; nPowerAndFreq.StartFreq = t.StartFreq - off; nPowerAndFreq.StopFreq = t.StopFreq + off; } else if (m_DataProcessing.nbddcbandwidth == 30) { nPowerAndFreq.StartFreq = t.StartFreq - 0.0089; nPowerAndFreq.StopFreq = t.StopFreq + 0.0089; } else if (m_DataProcessing.nbddcbandwidth == 50) { nPowerAndFreq.StartFreq = t.StartFreq - 0.0228; nPowerAndFreq.StopFreq = t.StopFreq + 0.0228; } break; } case 2: { nPowerAndFreq.StartFreq = t.StartFreq - 0.07 - 0.640 / power.Length; nPowerAndFreq.StopFreq = t.StopFreq + 0.07; break; } case 3: { nPowerAndFreq.StartFreq = t.StartFreq; nPowerAndFreq.StopFreq = t.StopFreq; break; } } SaveFrequencySpectrum(nPowerAndFreq.Power); if (passPowerAndFreq != null) { passPowerAndFreq(nPowerAndFreq); } }
public void PushData(IQDataAndFreq t) { m_queue.Enqueue(t); }