2015-02-10 70 views
0

p(i,x)=a0,i+a1,icos(2π/T*x)+b1,isin(2π/T*x)+c1,ix如何,使用R

在哪裏寫這個方程函數:

T=365a0,ia1,ib1,ic1,i是參數和x是唯一的價值,我們的輸入。

我很困惑,如何寫功能,而不當參數未知

回答

1

你剛纔寫的函數一般

p_gernal <- function(x, a0, a1, b1, c1, t){ 
    a0 + a1 * cos(2*pi/t * x) + b1 * sin(2*pi/t * x) + c1 * x 
} 

[RX,A0 ,a1,b1,c1參數應該是相同長度的向量 或長度1

head(p_gernal(x = x, a0 = a0, a1 = a1, b1 = b1, c1 = c1, 365)) 
[1] 0.0000000000 0.0205558610 0.0000614598 0.0866477574 0.0137448690 

因此,對於一個固定的指數

p_i <- function(x){ 
    p_gernal(x = x, a0 = 0, a1 = 1, b1 = 2, c1 = 3, 365) 
} 
> p_i(1:5) 
[1] 4.034279 7.068251 10.101906 13.135235 16.168228 
+0

你的參數被命名爲't',但你把它作爲'T'。 – 2015-02-10 21:12:55

+0

THX,只是修好了 – Rentrop 2015-02-10 21:13:44