有偿求助！！！vba做数据分析，计算两个4149*4149数组相乘，跑了很久没跑出来，见下述

vba做pls

      我做的是一个pls分析（就是一种数据分析），其中我用的是自己写的矩阵乘法运算，加了一点优化在里面（对含很多0元素才有效果），里面的各种函数都是自己编写的（如jian(x(),y())函数是两个数组对应元素相减），但都是用循环实现的。调试过程中，我发现有一步要计算两个4149*4149数组相乘（图中我标断点处，主过程名为test，调用pls函数，pls中调用了我自己编写的矩阵运算函数），pls函数中第一次循环时由于第一个矩阵chg是单位矩阵，我的优化有点效果（我猜的），我调试到第二次时chg不再是单位矩阵了，这下就是复杂度O（n^3)的矩阵相乘运算了，直接跑不出来了，请指教。（也尝试了内置函数mmult，同样跑死。）
我自己考虑了一些问题，但我是第一次使用vba，所以对这些也不清楚，只是感觉会影响到运行速度，不知道值不值得去改：
1、我定义的数组都没有指定类型（原因是我从单元格拉数据进我定义的single或者定义的double数组时报错类型不匹配）
2、 计算过程我用了很多中间变量，因为这些函数是自己写的，如果直接嵌套写会报错如a = mult(x,trans(y))（mult和trans都是我自己写的，因为在网上看到内置的函数如转置在数组维度很大时会出错，我尝试了一下的确转不了，后面就都自己写了）

希望能教我解决一下这个矩阵相乘跑死的问题，我所考虑的问题不知道对整个程序运行减少时间有没有帮助，目前整个程序的运行速度也比较慢，如能帮助解决有偿！！！！由于文件超过2M，无法附件上传，劳烦在百度云上下载一下。
链接：https://pan.baidu.com/s/1GQx2bEmf6VH7EL-oTho7kA
提取码：tz28

我上传的文件中pls里大循环只有一次，第二次循环时计算量很大跑死，我要解决的问题最多可能跑到10几次循环

约定的童话

看的瑟瑟发抖，你这其他语言转换来的吧？

zopey

A=4149*98
B=98*1

C=A*B=4149*1

vba做pls

约定的童话 发表于 2019-12-6 14:42
看的瑟瑟发抖，你这其他语言转换来的吧？

嗯，matlab先写

lsdongjh

要优化的地方太多，比如下面这个函数，简单处理一下，就不知道快多少了

1. Function eye(n)
   
2.     ReDim earr(1 To n, 1 To n) As Double 'As Double
   
3.     Dim i As Integer
   
4.     For i = 1 To n
   
5.         earr(i, i) = 1
   
6.      Next
   
7.     eye = earr
   
8. End Function

复制代码

vba做pls

lsdongjh 发表于 2019-12-9 09:33
要优化的地方太多，比如下面这个函数，简单处理一下，就不知道快多少了

您好，谢谢您回复。我一开始是这样您这样写的，指定数组类型double后默认为0值。
但是外面调用时是这样：
dim chg()

chg = eye(n)
如果指定了数据类型会报错，类型不匹配，也就是我提问中，我自己考虑的问题里第一点：
1、我定义的数组没有指定类型（原因是我从单元格拉数据进来时提示数据类型不匹配）
不知道您觉得有哪些地方可以改进

lsdongjh

从单元格中，用Variant 类型导入时，先进行转换，后续所有的传参、入参都指定类型！

lsdongjh

1. '数组转为Double， 也可行列转换
   
2. Function TransArrToDouble(arr As Variant, Optional blIsTrans As Boolean = True) As Double()
   
3.     Dim lngR As Long, lngC As Long
   
4.     Dim dblTemp() As Double
   
5. 
   
6.     If blIsTrans Then
   
7.         ReDim dblTemp(LBound(arr, 2) To UBound(arr, 2), LBound(arr) To UBound(arr)) As Double
   
8.         For lngR = LBound(arr) To UBound(arr)
   
9.             For lngC = LBound(arr, 2) To UBound(arr, 2)
   
10.                 dblTemp(lngC, lngR) = Val(arr(lngR, lngC))
    
11.             Next
    
12.         Next
    
13.     Else
    
14.         ReDim dblTemp(LBound(arr) To UBound(arr), LBound(arr, 2) To UBound(arr, 2)) As Double
    
15.         For lngR = LBound(arr) To UBound(arr)
    
16.             For lngC = LBound(arr, 2) To UBound(arr, 2)
    
17.                 dblTemp(lngR, lngC) = Val(arr(lngR, lngC))
    
18.             Next
    
19.         Next
    
20.     End If
    
21. 
    
22.     TransArrToDouble = dblTemp
    
23. End Function

复制代码

1. 'w到w*变换矩阵的初始化
   
2. Function EYE(lngCount As Long) As Double()
   
3.     Dim dblArr() As Double, lngID As Long
   
4. 
   
5.     ReDim dblArr(1 To lngCount, 1 To lngCount) As Double
   
6. 
   
7.     For lngID = 1 To lngCount
   
8.         dblArr(lngID, lngID) = 1
   
9.     Next
   
10. 
    
11.     EYE = dblArr
    
12. End Function

复制代码

1.     arrTemp = sh.Range("A1").Resize(lngSample, lngCountInd)
   
2.     varIndependent_Old = TransArrToDouble(arrTemp, False)
   
3.     varIndependent_New = TransArrToDouble(arrTemp)
   

复制代码

1.     Dim dblChg() As Double   
   
2.     dblChg = EYE(4149) 'w到w*变换矩阵的初始化
   

复制代码

vba做pls

lsdongjh 发表于 2019-12-10 11:38

十分感谢，我刚刚尝试了一下，的确我看有的教程上说指定类型会比用variant类型快，我打算后期优化时对这个方面进行改进。方便的话可以加我微信 17718151693（注明一下）。您的代码规范性十足

vba做pls

lsdongjh 发表于 2019-12-10 11:38

十分感谢，我刚刚尝试了一下，的确我看有的教程上说指定类型会比用variant类型快，我打算后期优化时对这个方面进行改进。方便的话可以加我微信 17718151693。您的代码规范性十足