Smalltalk Vandermonde矩阵

时间:2016-02-23 22:01:59

标签: smalltalk gnu-smalltalk

简而言之,这是一个Vandermonde矩阵,我有一个问题是在数组的第二维中运行for。

'add meg M-et majd N-et (enter kozotte)(az 1. sor az 1-es szam hatvanyai)' displayNl.
M := stdin nextLine asInteger.
N := stdin nextLine asInteger.
|tomb|
tomb := Array new: M.
x := 1.
y := 1.
a := M + 1.
b := N + 1.
x to: a do: [ :i|
  tomb at:x put: (Array new: N) y to: b do: [ :j |
    x at: y put: (x raisedTo: y - 1) ] ].
tomb printNl.

1 个答案:

答案 0 :(得分:1)

这是创建矩阵的好方法,我们有一个通用条目aij的表达式:

Matrix class >> fromBlock: aBlock rows: n columns: m
  | matrix |
  matrix := self rows: n columns: m.
  matrix indicesDo: [:i :j | | aij |
    aij := aBlock value: i value: j.
    matrix at: i at: j put: aij].
  ^matrix

使用上述方法,您现在可以实现

Matrix class >> vandermonde: anArray degree: anInteger
  ^self
    fromBlock: [:i :j | (anArray at: i) raisedTo: j - 1]
    rows: anArray size
    columns: anInteger + 1

修改

我刚刚意识到在Pharo中有一种方法可以从aij的表达式创建一个矩阵,它被命名为rows:columns:tabulate:,所以我的答案简化为:

Matrix class >> vandermonde: anArray degree: anInteger
  ^self
    rows: anArray size
    columns: anInteger + 1
    tabulate: [:i :j | (anArray at: i) raisedTo: j - 1]