[問題] 微分非鄰近數

看板MATLAB作者 (shin)時間12年前 (2013/08/01 04:13), 編輯推噓1(100)
留言1則, 1人參與, 最新討論串1/1
目前有一筆十萬個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

08/01 05:34, , 1F
M=5000;d=(x(M+1:end)-x(1:end-M))/M;
08/01 05:34, 1F
文章代碼(AID): #1H-M_qcY (MATLAB)
文章代碼(AID): #1H-M_qcY (MATLAB)