Function convolve(arr1, arr2)
'arr1、arr2一维数组,卷积运算返回一维数组(数组从1开始计数)
Dim i&, j&, n1&, n2&, max_n&
'转为从1开始计数
If LBound(arr1) = 0 Then arr1 = WorksheetFunction.Transpose(WorksheetFunction.Transpose(arr1))
If LBound(arr2) = 0 Then arr2 = WorksheetFunction.Transpose(WorksheetFunction.Transpose(arr2))
n1 = UBound(arr1): n2 = UBound(arr2)
max_n = n1 + n2 - 1: ReDim result(1 To max_n) '返回数组大小
For i = 1 To max_n
For j = 1 To n1
If i - j >= 0 And i - j < n2 Then result(i) = result(i) + arr1(j) * arr2(i - j + 1)
Next
Next
convolve = result
End Function
举例
用复利理解卷积,假设每年存入100元,年利率5%,按复利计算每年结束时的累计本息总额
第1年本息=100*1.05=105
第2年本息=100*1.05^2+100*1.05=215.25
第3年本息=100*1.05^3+100*1.05^2+100*1.05=331.0125
第4年本息=100*1.05^4+100*1.05^3+100*1.05^2+100*1.05=452.563125
第5年本息=100*1.05^5+100*1.05^4+100*1.05^3+100*1.05^2+100*1.05=580.19128125
卷积运算
Sub 卷积测试()
a = Array(100, 100, 100, 100, 100)
b = Array(1.05, 1.05 ^ 2, 1.05 ^ 3, 1.05 ^ 4, 1.05 ^ 5)
Debug.Print Join(convolve(a, b), ", ")
End Sub
105, 215.25, 331.0125, 452.563125, 580.19128125, 475.19128125, 364.94128125, 249.17878125, 127.62815625
得到9个结果,前5个就是5年每年累计本息总额