excel宏编程用最小二乘法计算两参数weibull分布的参数

以下是使用Excel VBA编写的最小二乘法计算两参数Weibull分布的参数的示例代码：

Sub WeibullLSM()

    '定义变量
    Dim n As Integer '样本个数
    Dim i As Integer '循环变量
    Dim x() As Double '自变量x
    Dim y() As Double '因变量y
    Dim lnX() As Double '自变量x的自然对数
    Dim lnY() As Double '因变量y的自然对数
    Dim sX As Double 'x的和
    Dim sY As Double 'y的和
    Dim sXX As Double 'x平方的和
    Dim sXY As Double 'x*y的和
    Dim sLnX As Double 'ln(x)的和
    Dim sLnY As Double 'ln(y)的和
    Dim b As Double '斜率
    Dim c As Double '截距
    Dim a As Double 'Weibull分布的参数a
    Dim b1 As Double 'Weibull分布的参数b

    '输入样本数据
    n = InputBox("请输入样本个数：")
    ReDim x(1 To n)
    ReDim y(1 To n)
    For i = 1 To n
        x(i) = InputBox("请输入第" & i & "个x值：")
        y(i) = InputBox("请输入第" & i & "个y值：")
    Next i

    '计算各项和
    For i = 1 To n
        sX = sX + x(i)
        sY = sY + y(i)
        sXX = sXX + x(i) ^ 2
        sXY = sXY + x(i) * y(i)
        sLnX = sLnX + Log(x(i))
        sLnY = sLnY + Log(y(i))
    Next i

    '计算斜率和截距
    b = (n * sXY - sX * sY) / (n * sXX - sX ^ 2)
    c = (sY - b * sX) / n

    '计算Weibull分布的参数a和b
    b1 = Exp(-c / b)
    a = b / Gamma(1 + 1 / b)

    '输出结果
    MsgBox "Weibull分布的参数a为：" & a & vbCrLf & "Weibull分布的参数b为：" & b1

End Sub

'计算Gamma函数的值
Function Gamma(x As Double) As Double
    If x = 1 Then
        Gamma = 1
    ElseIf x < 1 Then
        Gamma = Sqr(PI) / (Sin(x * PI) * Gamma(1 - x))
    Else
        Gamma = (x - 1) * Gamma(x - 1)
    End If
End Function

使用方法：

打开Excel，按下Alt + F11打开VBA窗口。
在VBA窗口中，选择插入 -> 模块，然后将上述代码复制粘贴到模块中。
点击运行按钮或按下F5键，然后按照提示输入样本数据。
程序将计算出Weibull分布的参数a和b，并在消息框中输出结果。

注意事项：

样本数据必须是正数。
样本数据的顺序不影响计算结果。
如果样本数据包含0或负数，则程序会出现错误