删除复制文件

ningshuxia

2025-04-01 f2feed624a3186e020bc0283b41e36eb89164a6c

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
using MathNet.Numerics;
 
namespace IStation.Curve
{
    /// <summary>
    /// 样条曲线（二次拟合）：使用Math.Net
    /// </summary>
    public class TwinRelateCurve2M : IFitPoint
    {
        /// <summary>
        /// 
        /// </summary>
        public TwinRelateCurve2M(CurveExpress express)
        {
            _express = express;
        }
 
        /// <summary>
        /// 
        /// </summary>
        public TwinRelateCurve2M(List<CurvePoint> points)
        {
            if (points == null || points.Count < 3)
            {
                throw new Exception("出现数值计算问题, 无法解出二次拟合曲线！");
            }
 
 
            _express = BuildCurveExpress(points);
        }
 
        //曲线表达式
        private CurveExpress _express = null;
 
        /// <summary>
        /// 获取拟合点Y
        /// </summary>
        public double GetFitPointY(double x)
        {
            if (_express == null)
                return default;
            return _express.Index2 * x * x + _express.Index1 * x + _express.Index0;
        }
 
        /// <summary>
        /// 获取拟合点列表
        /// </summary>
        public List<CurvePoint> GetFitPoints(int pointNumber)
        {
            if (_express == null)
            {
                return default;
            }
            return GetFitPointsByXRange(_express.Min, _express.Max, pointNumber);
        }
 
        /// <summary>
        /// 通过X区间获取拟合点列表
        /// </summary>
        public List<CurvePoint> GetFitPointsByXRange(double minX, double maxX, int pointNumber)
        {
            if (_express == null)
            {
                return default;
            }
            if (pointNumber < 1)
            {
                pointNumber = 12;
            }
            if (pointNumber > 10000)
            {
                pointNumber = 10000;
            }
 
            double space = (maxX - minX) / (pointNumber - 1);
            if (space < 0.0001)
            {
                return default;
            }
            var points = new List<CurvePoint>();
            for (int i = 0; i < pointNumber; i++)
            {
                double x = space * i + minX;
                double y = GetFitPointY(x);
                points.Add(new CurvePoint(x, y));
            }
            return points;
        }
 
        /// <summary>
        /// 创建曲线表达式
        /// </summary>
        public static CurveExpress BuildCurveExpress(List<CurvePoint> points)
        {
            if (points == null || points.Count < 3)
            {
                return default;
            }
            points = points.OrderBy(x => x.X).ToList();
            CurveExpress express = new CurveExpress
            {
                FitType = eFitType.ConicCurve,
                DefinePoints = points.Select(x => new CurvePoint(x)).ToList(),
                Min = points.First().X,
                Max = points.Last().X
            };
 
            var X = (from x in points select x.X).ToArray();
            var Y = (from x in points select x.Y).ToArray();
            var res = Fit.Polynomial(X, Y, 2);
 
            //0次方
            express.Index0 = res[0];
            //1次方
            express.Index1 = res[1];
            //2次方
            express.Index2 = res[2];
            //3次方
            express.Index3 = 0;// res[3];
 
            return express;
        }
 
 
 
 
    }
}