DenseMatrix无效赋值

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

最后一行生成错误"invalid assignment"

还想知道是否有一种函数式的方式来表达上面的逻辑。

谢谢。

Answer 1

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

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

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

#r @"packagesMathNet.Numerics.3.8.0libnet40MathNet.Numerics.dll"
#r @"packagesMathNet.Numerics.FSharp.3.8.0libnet40MathNet.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函数，它接受一个额外的参数，类似于循环体。