[问答]

此函数为“直线-三次抛物线拟合”算法( 见西安交通大学《电机测试技术》）

问答对人有帮助，内容完整，我也想知道答案 0 这段多上的带码如何实现，求大神 Public Sub LinearAndCubicFit(X() As Single, Y() As Single, n As Integer, n1 As Integer, n2 As Integer) '** 此函数为“直线-三次抛物线拟合”算法( 见西安交通大学《电机测试技术》）** Dim n∑Xi As Single Dim n∑XiE2 As Single Dim n∑Yi As Single Dim n∑XiYi As Single Dim n1∑Xi As Single Dim n1∑Yi As Single Dim n1∑XiE2 As Single Dim n1∑XiYi As Single Dim n2∑XiE2 As Single Dim n2∑XiE3 As Single Dim n2∑XiE4 As Single Dim n2∑XiE5 As Single Dim n2∑XiE6 As Single Dim n2∑XiE2Yi As Single Dim n2∑XiE3Yi As Single Dim A(1 To 3, 1 To 3) As Single '方程组“系数矩阵” Dim C(1 To 3) As Single '方程组“常数项向量” Dim B() As Single ReDim B(1 To 3) As Single '方程组“解向量” Dim Q_13(1 To 20) As Single '直线-三次抛物线拟合之“残差平方和” Dim n1∑ViE2 As Single Dim n2∑ViE2 As Single Dim i, j As Integer Dim X0 As Single For i = 1 To n n∑Xi = n∑Xi + X(i) n∑Yi = n∑Yi + Y(i) n∑XiYi = n∑XiYi + X(i) * Y(i) n∑XiE2 = n∑XiE2 + X(i) ^ 2 Next i DoEvents For i = n2 + 1 To n n1∑Xi = n1∑Xi + X(i) n1∑Yi = n1∑Yi + Y(i) n1∑XiE2 = n1∑XiE2 + X(i) ^ 2 n1∑XiYi = n1∑XiYi + X(i) * Y(i) Next i DoEvents For i = 1 To n2 n2∑XiE2 = n2∑XiE2 + X(i) ^ 2 n2∑XiE3 = n2∑XiE3 + X(i) ^ 3 n2∑XiE4 = n2∑XiE4 + X(i) ^ 4 n2∑XiE5 = n2∑XiE5 + X(i) ^ 5 n2∑XiE6 = n2∑XiE6 + X(i) ^ 6 n2∑XiE2Yi = n2∑XiE2Yi + X(i) ^ 2 * Y(i) n2∑XiE3Yi = n2∑XiE3Yi + X(i) ^ 3 * Y(i) Next i DoEvents j = 1 For X0 = 0.3 To 1.1 Step 0.05 A(1, 1) = n A(1, 2) = n∑Xi: A(2, 1) = A(1, 2) A(1, 3) = n2∑XiE3 - 3 * X0 * n2∑XiE2 - 3 * X0 ^ 2 * n1∑Xi + n1 * X0 ^ 3: A(3, 1) = A(1, 3) A(2, 2) = n∑XiE2 A(2, 3) = n2∑XiE4 - 3 * X0 * n2∑XiE3 - 3 * X0 ^ 2 * n1∑XiE2 + X0 ^ 3 * n1∑Xi: A(3, 2) = A(2, 3) A(3, 3) = n2∑XiE6 - 6 * X0 * n2∑XiE5 + 9 * X0 ^ 2 * n2∑XiE4 + 9 * X0 ^ 4 * n1∑XiE2 - 6 * X0 ^ 5 * n1∑Xi + n1 * X0 ^ 6 C(1) = n∑Yi C(2) = n∑XiYi C(3) = n2∑XiE3Yi - 3 * X0 * n2∑XiE2Yi - 3 * X0 ^ 2 * n1∑XiYi + X0 ^ 3 * n1∑Yi B = SolvingEquations(A(), C(), 3) '调解方程组函数 Dim ValueTemp As Single For i = n2 + 1 To n ValueTemp = Y(i) - (B(1) + B(2) * X(i) + B(3) * (X0 ^ 3 - 3 * X0 ^ 2 * X(i))) ValueTemp = ValueTemp ^ 2 n1∑ViE2 = n1∑ViE2 + ValueTemp Next i For i = 1 To n2 ValueTemp = Y(i) - (B(1) + B(2) * X(i) + B(3) * (X(i) ^ 3 - 3 * X0 * X(i) ^ 2)) ValueTemp = ValueTemp ^ 2 n2∑ViE2 = n2∑ViE2 + ValueTemp Next i Q_13(j) = n1∑ViE2 + n2∑ViE2 If Q_13(1) >= Q_13(j) Then Q_13(1) = Q_13(j) P_X0 = X0 AA(0) = B(1) + B(3) * P_X0 ^ 3 AA(1) = B(2) - 3 * B(3) * P_X0 ^ 2 BB(0) = B(1) BB(1) = B(2) BB(2) = -3 * B(3) * P_X0 BB(3) = B(3) End If j = j + 1 n1∑ViE2 = 0: n2∑ViE2 = 0 Next X0 P_QLineCurve = Q_13(1) End Sub 2 已退回5积分
2016-7-16 15:17:13　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × jihuahua 该类别下有 140 个回答，其中被选为最佳答案 4 次。邀请回答羊驼啊该类别下有 77 个回答，其中被选为最佳答案 2 次。邀请回答唐少华该类别下有 61 个回答。邀请回答酸牛奶甜该类别下有 55 个回答。邀请回答 btlook 该类别下有 48 个回答。邀请回答 qiujingshui 该类别下有 47 个回答。邀请回答世人皆不同该类别下有 42 个回答。邀请回答 hjc6063209 该类别下有 42 个回答。邀请回答 907510265 该类别下有 41 个回答。邀请回答 cc栗子该类别下有 41 个回答。邀请回答 woshisu 该类别下有 39 个回答。邀请回答哇哈哈哈6 该类别下有 34 个回答。邀请回答周建远该类别下有 34 个回答，其中被选为最佳答案 87 次。邀请回答留白该类别下有 34 个回答。邀请回答 ww 该类别下有 31 个回答，其中被选为最佳答案 2 次。邀请回答志尊宝该类别下有 31 个回答。邀请回答 wenT 该类别下有 30 个回答。邀请回答 yingtn 该类别下有 30 个回答。邀请回答舞出精彩该类别下有 29 个回答。邀请回答阳光1下该类别下有 29 个回答。邀请回答举报 loge2012 相关推荐 • 野外无线传感器网络有什么特点？ 957 • 学习模拟电路的苦恼？请老师们指点迷津！ 3496 • 步进电机的加、减速是怎么加的？ 4285 • 关于Q3D产生3D模型RLC参数后，输出.cir文件后怎么在PSPICE中进行仿真 8735 • 如何用labview中的波形图画除抛物线；还有我有两个数组，... 3538 • 求助-双三次B样条拟合背景的matlab实现 2222 • 什么是直线电机 2169 • 使用matlab对数据进行拟合 4279 • 抛物线时基发生器是什么？ 1552 • 如何对一波形所有极大（小）值点用三次样条插值函数拟... 2808 2个回答

答案对人有帮助，有参考价值 0 使用ＶＢ写的吗？

2016-7-17 09:40:35 评论举报 chenwei6991627

答案对人有帮助，有参考价值 0 本帖最后由 loge2012 于 2016-7-18 09:56 编辑 chenwei6991627 发表于 2016-7-17 09:40 使用ＶＢ写的吗？对，VB写的网上也可以搜到的，准备用labview来做软件，好多这样的代码无从下手

2016-7-18 09:55:49 评论举报 loge2012