Re: [問題] 矩陣生成

[celestialgod (天)]

要寫很多行QQ

# 第一種
nr = nc = 4
m = matrix(0, nr, nc)
m[row(m) == col(m)] = 2
m[1, 1] = 1
m[nrow(m), nrow(m)] = 1
m[row(m) == col(m) - 1] = -1
m[row(m) == col(m) + 1] = -1
m[1, 2] = 0
m[nrow(m), nrow(m)-1] = 0


# 第二種
nr = nc = 4
m = matrix(0, nr, nc)
m[row(m) == col(m)] = c(1, rep(2, nrow(m)-2), 1)
m[row(m) == col(m) - 1] = c(0, rep(-1, nrow(m)-2))
m[row(m) == col(m) + 1] = c(rep(-1, nrow(m)-2), 0)


※ 引述《hankckdy (漢克ㄦ)》之銘言：
: [問題類型]:
: 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
: [軟體熟悉度]:
: 新手(沒寫過程式，R 是我的第一次)
: [問題敘述]:
: Create a square matrix, say A, with entries aii = 2, i = 2, . . . , n 1,
: a11 = ann = 1, ai,i+1 = ai,i-1 = -1 , and aij = 0 elsewhere.
: 若寫成4x4的矩陣的話 看起來應該是像這樣
:  1  0  0  0
: -1  2 -1  0
:  0 -1  2 -1
:  0  0  0  1
: 請問要如何以程式碼呈現呢?
: 謝謝
: 程式範例]:

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 123.205.27.107
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1442362526.A.17E.html
推 hankckdy: 感謝回答 目前是R的新手 手中只有一本R for everyone可
→ hankckdy: 參考 想請問板大有其他推薦的參考書籍或資源可以利用的
→ hankckdy: 嗎? 再次謝謝你!!
置底很多XD

我自己只有讀過R in nutshell
※ 編輯: celestialgod (123.205.27.107), 09/16/2015 12:18:40
推 thea: 版大好強ＱＱ