樣條插值中的R

-1

h = c[157 144 80 106 124 46 207 188 190 208 143 170 162 178 155 163 162 149 135 160 149 147 133 146 126 120 151 74 122 145 160 155 173 126 172 93]

一些數據，然後我得到的局部最大值與差異功能

max = [1 5 7 10 12 14 16 20 24 27 31 33 35 36]

我有簡單的MATLAB代碼，找出樣條插值

maxenv = spline(max,h(max),1:N);

此代碼將顯示結果

maxenv = 157 86.5643152828762 67.5352696350679 84.9885891697257 124 169.645228239041 207 224.396380746179 223.191793341491 208 185.421032390413 170 172.173624690130 178 172.759468849065 163 158.147870987344 157.874968589889 159.414581897490 160 157.622863516083 153.308219179638 148.839465253375 146 146.051320982064 148.167322961480 151 153.474200222188 155.606188003845 157.685081783579 160 163.653263154551 173 186.027639098700 172 93

現在，我有R努力卻得到了一些錯誤

maxenv <- spline(max, h(max), seq(N))

，則錯誤在xy.coords

錯誤（X，Y）： 'X' 和 'Y' 的長度不同

來源

2015-11-27 public_html

這個錯誤本身就很明顯。 'max'和'h（max）'長度不同。閱讀'？spline'。 – shekeine

你在用什麼，Matlab或R？ – nico

我想試試R –

你的代碼和問題真的是不明確的，所以我不是100％肯定這是你在找什麼：

# vector of values 

h <- c(157, 144, 80, 106, 124, 46, 207, 188, 190, 208, 143, 170, 162, 178, 155, 
     163, 162, 149, 135, 160, 149, 147, 133, 146, 126, 120, 151, 74, 122, 145, 
     160, 155, 173, 126, 172, 93) 

# local maxima 
lmax <- h[c(1, which(diff(sign(diff(h)))==-2)+1, length(h))] 

# spline calculation 
spl <- spline(1:length(lmax), lmax) 

# visual inspection 
plot(spl) 
lines(spl)

來源

2015-11-27 15:34:19 hrbrmstr

樣條插值中的R

回答

相關問題