M-K突变点检验源代码（C#版）

guxing-345

想用C#计算突变点了，网上各种版本的都有就是没C#版的，于是自己写了个，供大家一起用吧。

1. class Program
   
2.     {
   
3.         static void Main(string[] args)
   
4.         {
   
5.             try
   
6.             {
   
7.                 //输入文件及输出文件定义
   
8.                 string in_file = @"E:\year.txt";
   
9.                 string out_file = @"E:\out_year.txt";
   
10.                 //导入数据
    
11.                 List<mkdata> list = new List<mkdata>();
    
12.                 StreamReader sr = new StreamReader(in_file, Encoding.Default);
    
13.                 sr.ReadLine();
    
14.                 while (sr.Peek() > -1)
    
15.                 {
    
16.                     string[] line = sr.ReadLine().Trim().Split(new char[] { ' ', ',', '\t' }, StringSplitOptions.RemoveEmptyEntries);
    
17.                     mkdata mk = new mkdata();
    
18.                     mk.year = Int32.Parse(line[0]);
    
19.                     mk.value = double.Parse(line[1]);
    
20.                     list.Add(mk);
    
21.                 }
    
22.                 sr.Close();
    
23.                 sr.Dispose();
    
24.                 //计算MK突变点
    
25.                 int n = list.Count;
    
26.                 //正序列计算
    
27.                 double[] UFk = new double[n];
    
28.                 double[] y = new double[n];
    
29.                 for (int i = 0; i < n; i++)
    
30.                 {
    
31.                     y[i] = list[i].value;
    
32.                 }
    
33.                 UFk = GetUFk(y);
    
34.                 //逆序列计算
    
35.                 double[] UBk = new double[n];
    
36.                 double[] y2 = new double[n];
    
37.                 for (int i = 0; i < n; i++)
    
38.                 {
    
39.                     y2[i] = list[n - i - 1].value;
    
40.                 }
    
41.                 UBk = GetUFk(y2);
    
42.                 double[] UBk2 = new double[n];
    
43.                 for (int i = 0; i < n; i++)
    
44.                 {
    
45.                     UBk2[i] = UBk[n - i - 1];
    
46.                 }
    
47.                 StreamWriter sw = new StreamWriter(out_file, false, Encoding.Default);
    
48.                 for (int i = 0; i < n; i++)
    
49.                 {
    
50.                     sw.WriteLine("{0}\t{1}\t{2}\t{3}\t{4}", list[i].year, Math.Round(UFk[i], 3), Math.Round(UBk2[i], 3), -1.96, 1.96);
    
51.                 }
    
52.                 sw.Flush();
    
53.                 sw.Close();
    
54.                 sw.Dispose();
    
55.                 Console.WriteLine("计算完成!");
    
56.                 Console.ReadKey();
    
57.             }
    
58.             catch (Exception ex)
    
59.             {
    
60.                 Console.WriteLine(ex.Message);
    
61.                 Console.ReadKey();
    
62.             }
    
63.         }
    
64.         /// <summary>
    
65.         /// 计算UF或者UB
    
66.         /// </summary>
    
67.         /// <param name="y"></param>
    
68.         /// <returns></returns>
    
69.         static double[] GetUFk(double[] y)
    
70.         {
    
71.             try
    
72.             {
    
73.                 int n = y.Length;
    
74.                 double[] Sk = new double[n];
    
75.                 double[] UFk = new double[n];
    
76.                 double s = 0;
    
77.                 for (int i = 1; i < n; i++)
    
78.                 {
    
79.                     for (int j = 0; j < i; j++)
    
80.                     {
    
81.                         if (y[i] > y[j])
    
82.                         {
    
83.                             s += 1;
    
84.                         }
    
85.                     }
    
86.                     Sk[i] = s;
    
87.                     double E = ((i + 1) * i) / 4.0;
    
88.                     double var = ((i + 1) * i * (2 * (i + 1) + 5)) / 72.0;
    
89.                     UFk[i] = (Sk[i] - E) / Math.Sqrt(var);
    
90.                 }
    
91.                 return UFk;
    
92.             }
    
93.             catch (Exception ex)
    
94.             {
    
95.                 throw ex;
    
96.             }
    
97.         }
    
98.     }
    
99. 
    
100.     public struct mkdata
     
101.     {
     
102.         public int year;
     
103.         public double value;
     
104.     }

复制代码

hunter2019

有没有算出来的结果图？上一个来看看？

hunter2019

这个高大上