Question

不确定如何填充DenseMatrix类型：

let rows = [|for line in File.ReadAllLines("Z:\\mypath.csv") 
            |> Seq.skip 1 do yield line.Split(',') |> Array.map float|]
let data = DenseMatrix.ofRowArrays rows
let mutable data_logdiff = DenseMatrix.zero<float> (data.RowCount-1) (data.ColumnCount)

for i in [0 .. data.ColumnCount-1] do
    for j in [1 .. data.RowCount-1] do
        data_logdiff.At(j-1, i) <- data.At(j, i) / data.At(j-1, i) |> log

最后一行产生错误“无效分配”。

还想知道是否有一种表达上述逻辑的功能方式。

感谢。

Answer 1

我假设您使用的是MathNet，不知道哪个版本。

首先，我不认为你的data_logdiff需要是可变的，它已经是一个具有你可以变异的属性的对象，我的意思是你要改变它的属性，而不是对象本身。

然后你的代码中的问题是你如何改变这些属性，你应该使用索引器而不是.At，它只返回值：

#r @"packages\MathNet.Numerics.3.8.0\lib\net40\MathNet.Numerics.dll"
#r @"packages\MathNet.Numerics.FSharp.3.8.0\lib\net40\MathNet.Numerics.FSharp.dll"

open System.IO
open MathNet.Numerics.LinearAlgebra

let rows = [|for line in File.ReadAllLines("Z:\\mypath.csv") 
        |> Seq.skip 1 do yield line.Split(',') |> Array.map float|]
let data = DenseMatrix.ofRowArrays rows
let data_logdiff = DenseMatrix.zero<float> (data.RowCount-1) (data.ColumnCount)

for i in [0 .. data.ColumnCount-1] do
    for j in [1 .. data.RowCount-1] do
        data_logdiff.[j-1, i] <- data.At(j, i) / data.At(j-1, i) |> log

要使其功能更强大，请使用DenseMatrix.init代替DenseMatrix.zero，然后循环播放：

let data_logdiff = 
    DenseMatrix.init 
        (data.RowCount-1) 
        (data.ColumnCount) 
        (fun j i -> if j = 0 then 0. else data.At(j, i) / data.At(j-1, i) |> log)

一般来说，每次使用.zero并使用循环来初始化矩阵或向量时，请考虑使用.init函数，这需要一个额外的参数，就像循环的主体一样。

DenseMatrix无效分配

1 个答案: