主要内容图形平台实现最小二乘数据拟合.docVIP

主要内容：图形平台实现最小二乘数据拟合 目的：对比通用语言平台和图形平台进行数据拟合的优缺点；自己设计算法和调用MATLAB引擎求解的优缺点；初步建立多个方案中择优的思想 重点：掌握通用语言平台调用MATLAB //图形平台直线拟合 Public Sub lineapproximate() Dim n As Integer 点数 Dim c() As Double 保存各点的动态数组 Dim zx As AcadLWPolyline 折线图 Dim nhqx As AcadLine 拟和曲线图 Dim pt As Variant Dim i As Integer 通过交互方式取得n的值 n = ThisDrawing.Utility.GetInteger(input a integer:) ReDim c(2 * n - 1) 读取各点数据 For i = 0 To 2 * n - 1 Step 2 pt = ThisDrawing.Utility.GetPoint(, input a point:) c(i) = pt(0) c(i + 1) = pt(1) Next 生成折线图 Set zx = ThisDrawing.ModelSpace.AddLightWeightPolyline(c) zx.Color = acRed 计算拟和函数参数a0、a1 Dim s1 As Double Dim s2 As Double Dim s3 As Double Dim s4 As Double Dim a0 As Double Dim a1 As Double Dim p1(2) As Double Dim p2(2) As Double Dim s As Double s1 = 0: s2 = 0: s3 = 0: s4 = 0 For i = 1 To n s1 = s1 + c(2 * i - 2) s2 = s2 + c(2 * i - 1) s3 = s3 + c(2 * i - 2) * c(2 * i - 2) s4 = s4 + c(2 * i - 2) * c(2 * i - 1) Next s = n * s3 - s1 * s1 If s 0# Then a0 = (n * s4 - s1 * s2) / (n * s3 - s1 * s1) a1 = (s3 * s2 - s1 * s4) / (n * s3 - s1 * s1) p1(0) = c(0) p1(1) = a0 * c(0) + a1 p1(2) = 0 p2(0) = c(2 * n - 2) p2(1) = a0 * p2(0) + a1 p2(2) = 0 Set nhqx = ThisDrawing.ModelSpace.AddLine(p1, p2) nhqx.Color = acBlue nhqx.Update End If End Sub //绘图自动化 Public Sub creatsinx() Dim p(61) As Double 轻量多义线点坐标；是等份段数的二倍 Dim i As Integer 循环控制变量 Dim n As Double 图形放大倍数 Dim startp(2) As Double 图形起点坐标 Dim ssin As AcadLWPolyline Dim delta As Double x坐标偏移量（步长） Dim x As Double 每个点的x坐标 Const pi = 3.1416 delta = pi / 16 一个周期等份段数为32段 startp(0) = 50 startp(1) = 50 startp(2) = 0 x = 0 For i = 0 To 60 Step 2 p(i) = x + startp(0) p(i + 1) = Sin(x) + startp(1) x = x + delta Next Set ssin = ThisDrawing.ModelSpace.AddLightWeightPolyline(p) n = 10 ssin.ScaleEntity startp, n End Sub