[問題] 微分非鄰近數

[xj04a83 (shin)]

目前有一筆十萬個data資料

希望的微分形式是 x(5001)-x(1)/5001-1, x(5002)-x(2)/5002-2 ...

這種情況下的微分應該如何寫呢?


之前我是令另一個變數 (y)  (用for loop)

在y(5000)前 y(i) = x(i)，第5001開始 y(5001) = x(1)

之後再另 dy = x - y; dt = 5000; dydt = dy./dt


可是這樣跑時間會拉很長，所以想請問還有沒有其他的方法呢?

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 24.87.35.10
推 tml:M=5000;d=(x(M+1:end)-x(1:end-M))/M;