public double[] getAnalyzeDataC()
{
g_OffSet = 0;//初始化
//使用互相关算法分析
int len = 262144;
double[] A = new double[len];
double[] B = new double[len];
for (int i = 0; i < len; i++)
{
A[i] = g_DataA[i];
B[i] = g_DataB[i];
}
CrossCorrelation ccorr = new CrossCorrelation();
double[] R = ccorr.Rxy(A, B);
int R_length = R.Length;
Position p = new Position();
int location = p.max(R);
int d = location - R_length / 2;
g_OffSet = d;//偏移值
return R;//经过互相关计算后的数组
}
private void button3_Click(object sender, EventArgs e)
{
EraseGraph();
CreateGraph();
panel3.Controls.Clear();
/*Label label13 = new Label();
* Label label14 = new Label();
* label13.Font = label14.Font = new System.Drawing.Font("宋体", 15.75F, System.Drawing.FontStyle.Bold);
* label13.ForeColor = label14.ForeColor = Color.Red;
* label13.Location = new System.Drawing.Point(100, 30);
* label14.Location = new System.Drawing.Point(900, 30);
* label13.AutoSize = label14.AutoSize = true;
* panel3.Controls.Add(label13);
* panel3.Controls.Add(label14);*/
ZedGraphControl zed3 = new ZedGraphControl();
zed3.Width = panel3.Width;
zed3.Height = panel3.Height;
//读取数据
int datalength = WaterLeak.get_DataLength();
double[] dataA = WaterLeak.get_DataA();
double[] dataB = WaterLeak.get_DataB();
double dataA_avg = WaterLeak.get_DataA_Avg();
double dataB_avg = WaterLeak.get_DataB_Avg();
if (比较ToolStripMenuItem.Checked == true)
{
byte[] A = data_A;
byte[] B = data_B;
int[] m = new int[2 * datalength];
int j;
int k = 0;
for (int i = 0; i < 2 * datalength; i++)
{
if (A[i] != B[i])
{
j = i;
m[k] = j;
k++;
}
}
int index = 0;
for (int i = 0; i < 2 * datalength; i++)
{
if (m[i] == 0)
{
index = i;
}
}
int[] n = new int[index];
for (int i = 0; i < index; i++)
{
n[i] = m[i];
}
MessageBox.Show("" + index);
}
//处理数据，采用均方差算法
if (算术平均ToolStripMenuItem.Checked == true)
{
progressBar1.Value = 0;
int j = 0;
double MINC = 0;
int bias = 5000;
int w = 5000;
double[] C = new double[2 * bias];
for (int i = -bias; i < bias; i++)
{
C[i + bias] = 0;
for (int k = w; k < datalength - w; k++)
{
C[i + bias] += Math.Abs(dataA[k + i] - dataB[k]); //曼哈顿距离
}
if (i == -bias)
{
MINC = C[0];
}
else if (C[bias + i] < MINC)
{
j = i;
MINC = C[bias + i];
}
Application.DoEvents();
progressBar1.Value = (i + bias) * 50 / bias;
}
progressBar1.Value += 1; //进度条
double w_avg;
w_avg = 0;
for (int i = 0; i < 2 * bias; i++)
{
w_avg += C[i] / (2 * bias);
}
for (int i = 0; i < 2 * bias; i++)
{
C[i] = (C[i] - w_avg) / w_avg;
}
panel3.Controls.Add(zed3);
GraphPane Pane = zed3.GraphPane;
/*Pane.Title = "偏差结果";
* Pane.XAxis.Title = "";
* Pane.YAxis.Title = "";*/
PointPairList list1 = new PointPairList();
for (int i = 0; i < 2 * bias; i++)
{
list1.Add(i, C[i]);
}
LineItem myCurve1 = Pane.AddCurve("Porsche", list1, Color.Blue, SymbolType.None);
zed3.AxisChange();
//计算漏水点位置
double positon = (WaterLeak.get_PipelineLength() + WaterLeak.get_WaveNumber() * j / (1000 * WaterLeak.get_CaptureRate())) / 2;
/*label14.Text = "漏水点距离A探头" + positon.ToString() + "米!";
*
* //显示标注数据
* if (j > 0)
* label13.Text = "A探头滞后B探头" + j.ToString("G") + "个基点!";
* if (j < 0)
* label13.Text = "A探头超前B探头" + (-j).ToString() + "个基点!";
* if (j == 0)
* label13.Text = "A探头与B探头无偏移!";*/
//显示基点
textBoxOffSet.Text = j.ToString();
//保存参数
WaterLeak.set_dataC_factor(500);
WaterLeak.set_dataC_shift(100);
WaterLeak.set_Bias(bias);
WaterLeak.set_Min_Position(j);
WaterLeak.set_Min_Value(MINC);
WaterLeak.set_DataC(C);
}
//if (LMSToolStripMenuItem.Checked == true)
//{
// progressBar1.Value = 0;
// MatlabClass mat = new MatlabClass();
// MWNumericArray dataA_m = new MWNumericArray(1, datalength, dataA);
// MWNumericArray dataB_m = new MWNumericArray(1, datalength, dataB);
// MWArray M_m = 100;
// MWArray mu_m = 0.00001;
// MWArray W_m = mat.my_LMS(dataA_m, dataB_m, M_m, mu_m);
// Array W = ((MWNumericArray)W_m).ToArray();
// int M = ((MWNumericArray)M_m).ToScalarInteger();
// int W_length = W.Length;
// int d = 0;
// double temp = Convert.ToDouble(W.GetValue(0, 0));
// for (int i = 1; i < W_length; i++)
// {
// if (temp < Convert.ToDouble(W.GetValue(i, 0)))
// {
// temp = Convert.ToDouble(W.GetValue(i, 0));
// d = i - 2 * M;
// }
// }
// panel3.Controls.Add(zed3);
// GraphPane Pane = zed3.GraphPane;
// Pane.Title = "偏差结果";
// Pane.XAxis.Title = "";
// Pane.YAxis.Title = "";
// PointPairList list = new PointPairList();
// for (int i = 0; i < W_length; i++)
// {
// list.Add(i - 2 * M, Convert.ToDouble(W.GetValue(i, 0)));
// }
// LineItem myCurve = Pane.AddCurve("Porsche", list, Color.Blue, SymbolType.None);
// zed3.AxisChange();
// //计算漏水点位置
// double positon = (WaterLeak.get_PipelineLength() + WaterLeak.get_WaveNumber() * d / (1000 * WaterLeak.get_CaptureRate())) / 2;
// label14.Text = "漏水点距离A探头" + positon.ToString() + "米!";
// //显示标注数据
// if (d > 0)
// label13.Text = "A探头滞后B探头" + d.ToString("G") + "个基点!";
// if (d < 0)
// label13.Text = "A探头超前B探头" + (-d).ToString() + "个基点!";
// if (d == 0)
// label13.Text = "A探头与B探头无偏移!";
//}
if (广义互相关ToolStripMenuItem.Checked == true)
{
DateTime startDT = System.DateTime.Now;
progressBar1.Value = 0;
int len = 262144;
double[] A = new double[len];
double[] B = new double[len];
for (int i = 0; i < len; i++)
{
A[i] = dataA[i];
B[i] = dataB[i];
}
CrossCorrelation ccorr = new CrossCorrelation();
double[] R = ccorr.Rxy(A, B);
int R_length = R.Length;
Position p = new Position();
int location = p.max(R);
int d = location - R_length / 2;
panel3.Controls.Add(zed3);
GraphPane Pane = zed3.GraphPane;
/*Pane.Title = "互相关系数";
* Pane.XAxis.Title = "";
* Pane.YAxis.Title = "";*/
PointPairList list = new PointPairList();
for (int i = 0; i < R_length; i++)
{
list.Add(i, R[i]);
}
LineItem myCurve = Pane.AddCurve("Porsche", list, Color.Blue, SymbolType.None);
zed3.AxisChange();
//计算漏水点位置
double positon = (WaterLeak.get_PipelineLength() + WaterLeak.get_WaveNumber() * d / (1000 * WaterLeak.get_CaptureRate())) / 2;
/*label14.Text = "漏水点距离A探头" + positon.ToString() + "米!";
* //显示标注数据
* if (d > 0)
* label13.Text = "A探头滞后B探头" + d.ToString("G") + "个基点!";
* if (d < 0)
* label13.Text = "A探头超前B探头" + (-d).ToString() + "个基点!";
* if (d == 0)
* label13.Text = "A探头与B探头无偏移!";*/
//显示基点
textBoxOffSet.Text = d.ToString();
DateTime endDT = System.DateTime.Now;
TimeSpan time = endDT.Subtract(startDT);
//MessageBox.Show(time.TotalMilliseconds.ToString() + "ms");
}
//if (互功率谱ToolStripMenuItem.Checked == true)
//{
// progressBar1.Value = 0;
// MatlabClass mat = new MatlabClass();
// MWNumericArray dataA_m = new MWNumericArray(1, datalength, dataA);
// MWNumericArray dataB_m = new MWNumericArray(1, datalength, dataB);
// MWArray datalength_m = datalength;
// MWArray R_m = mat.my_xcorrPow(dataA_m, dataB_m);
// Array R = ((MWNumericArray)R_m).ToArray();
// int R_length = R.Length;
// int d = 0;
// double temp = Convert.ToDouble(R.GetValue(0, 0));
// for (int i = 1; i < R_length; i++)
// {
// if (temp < Convert.ToDouble(R.GetValue(0, i)))
// {
// temp = Convert.ToDouble(R.GetValue(0, i));
// d = i - datalength + 1;
// }
// }
// panel3.Controls.Add(zed3);
// GraphPane Pane = zed3.GraphPane;
// Pane.Title = "偏差结果";
// Pane.XAxis.Title = "";
// Pane.YAxis.Title = "";
// PointPairList list = new PointPairList();
// for (int i = 0; i < R_length; i++)
// {
// list.Add(i - datalength + 1, Convert.ToDouble(R.GetValue(0, i)));
// }
// LineItem myCurve = Pane.AddCurve("Porsche", list, Color.Blue, SymbolType.None);
// zed3.AxisChange();
// //计算漏水点位置
// double positon = (WaterLeak.get_PipelineLength() + WaterLeak.get_WaveNumber() * d / (1000 * WaterLeak.get_CaptureRate())) / 2;
// label14.Text = "漏水点距离A探头" + positon.ToString() + "米!";
// //显示标注数据
// if (d > 0)
// label13.Text = "A探头滞后B探头" + d.ToString("G") + "个基点!";
// if (d < 0)
// label13.Text = "A探头超前B探头" + (-d).ToString() + "个基点!";
// if (d == 0)
// label13.Text = "A探头与B探头无偏移!";
//}
//int bias = 5000;
//int w = 2500;
//double[] C = new double[bias];
//double minc = 0;
//int d = 0;
//int n = 5;//计算次数
//int sp = datalength / n;//每次计算的间隔
//int[] T = new int[n];//时延点数
//double[] S = new double[n];//离A点距离
//DataTable dt = new DataTable();
//dt.Columns.Add("A探头超前点数", typeof(int));
//dt.Columns.Add("漏水点离A探头距离", typeof(double));
//for (int m = 0; m < n; m++)
//{
// for (int i = 0 + m * sp; i < bias + m * sp; i++)
// {
// C[i - m * sp] = 0;
// for (int j = w + m * sp; j < w + bias + m * sp; j++)
// {
// C[i - m * sp] += Math.Abs(dataA[i + j - w - m * sp] - dataB[j]);
// }
// if (i == 0) minc = C[0];
// else if (C[i - m * sp] < minc)
// {
// d = i - m * sp;
// minc = C[i - m * sp];
// }
// }
// T[m] = w - d;
// S[m] = 500 - (double)T[m] / 10;
// DataRow dr = dt.NewRow();
// dr[0] = T[m].ToString();
// dr[1] = S[m].ToString();
// dt.Rows.Add(dr);
// Application.DoEvents();
// progressBar1.Value = (m - 1) * 100 / n;
//}
//dataGridView1.DataSource = dt;
//dataGridView1.AutoSizeColumnsMode = DataGridViewAutoSizeColumnsMode.Fill;
/*
*
* //处理数据，采用LMS算法
* //初始化
* int M=1000;
* float mu = (float)0.1;
* float[] en = new float[datalength];
* for (int i = 0; i < datalength; i++)
* en[i] = 0;
* double[] w = new double[M];
* for (int i = 0; i < M; i++)
* w[i] = 0;
* float y;
* //迭代运算
* for (int k = M - 1; k < datalength; k++)
* {
* y=0;
* for (int i = 0; i < M; i++)
* y += (float)(w[i]) * dataB[k - i];
* if (double.IsNaN(y)) y = 0;
* en[k] =dataA[k] - y;
*
* for (int i = 0; i < M; i++)
* {
* w[i] = w[i] + 2 * mu * en[k] * dataB[k - i];
* if (double.IsNaN(w[i])) w[i] = 0;
* // w[i] = Math.Round(w[i], 6);
* }
*
* progressBar1.Value = (k-M) * 100 / (datalength - M-1);
*
* Application.DoEvents();
*
* }
*
*
* float w_avg;
* w_avg = 0;
* for (int i = 0; i < M; i++)
* w_avg += (float)w[i];
*
* w_avg /= M;
* // w_avg = Math.Round(w_avg, 6);
*
* // w[0] = Math.Round(w[0], 10);
* textBox1.Text = w[0].ToString("G");
*/
}