diff

差分と近似微分

構文

Y = diff(X)

Y = diff(X,n)

Y = diff(X,n,dim)

説明

Y = diff(X) は、X の隣接する要素間の差分を計算します。既定では、diff はサイズが 1 に等しくない最初の配列次元に沿って動作します。

• X が、長さ m のベクトルの場合、Y = diff(X) は長さが m-1 のベクトルを返します。Y の要素は、X の隣接する要素間の差分です。

  Y = [X(2)-X(1) X(3)-X(2) ... X(m)-X(m-1)]

• X が空ではなく p 行 m 列の非ベクトル行列である場合、Y = diff(X) はサイズが (p-1) 行 m 列の行列を返します。この行列の要素は X の行間の差分になります。

  Y = [X(2,:)-X(1,:); X(3,:)-X(2,:); ... X(p,:)-X(p-1,:)]

• X が 0 行 0 列の空の行列である場合、Y = diff(X) は 0 行 0 列の空の行列を返します。

• X が p 行 m 列の table または timetable である場合、Y = diff(X) はサイズが (p-1) 行 m 列の table または timetable を返します。この行列の要素は X の行間の差分になります。X が 1 行 m 列の table または timetable の場合、Y のサイズは 0 行 m 列です。 (R2023a 以降)

Y = diff(X,n) は、サイズが 1 に等しくない最初の配列次元に沿って diff(X) 演算子を再帰的に n 回適用することで n 階差分を計算します。

Y = diff(X,n,dim) は、dim で指定された次元に沿って計算された n 階差分です。dim 入力は、正の整数のスカラーです。

例

ベクトル要素間の差分

ベクトルを作成し、要素間の差分を計算します。

X = [1 1 2 3 5 8 13 21];
Y = diff(X)

Y = 1×7

     0     1     1     2     3     5     8

Y の要素数が X より 1 つ少なくなることに注意してください。

行列の行間の差分

3 行 3 列の行列を作成し、行間の 1 階差分を計算します。

X = [1 1 1; 5 5 5; 25 25 25];
Y = diff(X)

Y = 2×3

     4     4     4
    20    20    20

Y は 2 行 3 列の行列です。

複数の差分

ベクトルを作成し、要素間の 2 階差分を計算します。

X = [0 5 15 30 50 75 105];
Y = diff(X,2)

Y = 1×5

     5     5     5     5     5

行列の列間の差分

3 行 3 列の行列を作成し、列間の 1 次差分を計算します。

X = [1 3 5;7 11 13;17 19 23];
Y = diff(X,1,2)

Y = 3×2

     2     2
     4     2
     2     4

Y は 3 行 2 列の行列です。

diff による導関数の近似

関数 diff を構文 Y = diff(f)/h で使用して偏導関数を近似します。f は、領域 X で評価される関数値のベクトルで、h は適切なステップ サイズです。

たとえば、x に関する sin(x) の 1 次導関数は cos(x) で、x に関する 2 次導関数は -sin(x) になります。diff を使用して、これらの導関数を近似できます。

h = 0.001;       % step size
X = -pi:h:pi;    % domain
f = sin(X);      % range
Y = diff(f)/h;   % first derivative
Z = diff(Y)/h;   % second derivative
plot(X(:,1:length(Y)),Y,'r',X,f,'b', X(:,1:length(Z)),Z,'k')

このプロットで、青のラインは元の関数 sin に対応します。赤のラインは計算された 1 次導関数 cos に対応し、黒のラインは計算された 2 次導関数 -sin に対応します。

datetime 値の差分

等間隔の datetime 値のシーケンスを作成し、それらの時間差を求めます。

t1 = datetime('now');
t2 = t1 + minutes(5);
t = t1:minutes(1.5):t2

t = 1×4 datetime
   09-Aug-2025 12:09:38   09-Aug-2025 12:11:08   09-Aug-2025 12:12:38   09-Aug-2025 12:14:08

dt = diff(t)

dt = 1×3 duration
   00:01:30   00:01:30   00:01:30

diff は duration 配列を返します。

入力引数

X — 入力配列

ベクトル | 行列 | 多次元配列 | table | timetable

入力配列。ベクトル、行列、多次元配列、table、または timetable として指定します。X は数値配列、logical 配列、datetime 配列、duration 配列、またはそれらのデータ型のいずれかが含まれた変数をもつ table または timetable にすることができます。

複素数のサポート: あり

n — 微分の次数

1 (既定値) | 正の整数スカラー | []

微分の次数。正の整数スカラーまたは [] として指定します。n の既定値は 1 です。n が操作次元より大きい場合、diff は空の配列を返します。

データ型: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

dim — 演算の対象の次元

正の整数スカラー

演算の対象の次元。正の整数のスカラーとして指定します。次元を指定しない場合、既定値はサイズが 1 でない最初の配列次元です。

2 次元の p 行 m 列の入力配列 A について考えます。

• diff(A,1,1) は、A の列の連続する要素に対して有効で、(p-1) 行 m 列の差分行列を返します。

• diff(A,1,2) は、A の行の連続する要素に対して有効で、p 行 (m-1) 列の差分行列を返します。

データ型: double | single | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

出力引数

Y — 差分配列

スカラー | ベクトル | 行列 | 多次元配列 | table | timetable

差分配列。スカラー、ベクトル、行列、多次元配列、table、または timetable として返されます。diff により動作する X の次元のサイズが n より大きい場合、この次元は出力では n 減らされます。この次元が n より小さい場合はゼロまで減らされ、出力は空になります。

拡張機能

tall 配列

メモリの許容量を超えるような多数の行を含む配列を計算します。

この関数は tall 配列を制限付きでサポートしています。

3 つの入力の構文 Y = diff(X,N,dim) を使用しなければなりません。

詳細については、tall 配列を参照してください。

C/C++ コード生成

MATLAB® Coder™ を使用して C および C++ コードを生成します。

使用上の注意および制限:

• dim 引数を使用する場合、コード生成時に定数でなければなりません。

• 入力引数 X では、次のようになります。

  • 次元を指定しない場合、コード ジェネレーターは、可変サイズであるかサイズが 1 と等しくない入力配列の最初の次元に沿って演算を行います。この次元がコード生成時に可変サイズで実行時に 1 になる場合、ランタイム エラーが発生することがあります。このエラーを回避するには、次元を指定します。

• コード生成では、この関数のスパース行列入力はサポートされません。

スレッドベースの環境

MATLAB® の backgroundPool を使用してバックグラウンドでコードを実行するか、Parallel Computing Toolbox™ の ThreadPool を使用してコードを高速化します。

この関数はスレッドベースの環境を完全にサポートしています。詳細については、スレッドベースの環境での MATLAB 関数の実行を参照してください。

GPU 配列

Parallel Computing Toolbox™ を使用してグラフィックス処理装置 (GPU) 上で実行することにより、コードを高速化します。

diff 関数は、GPU 配列を完全にサポートします。GPU 上で関数を実行するには、入力データを gpuArray (Parallel Computing Toolbox) として指定します。詳細については、GPU での MATLAB 関数の実行 (Parallel Computing Toolbox)を参照してください。

分散配列

Parallel Computing Toolbox™ を使用して、クラスターの結合メモリ上で大きなアレイを分割します。

この関数は分散配列を完全にサポートしています。詳細については、分散配列を使用した MATLAB 関数の実行 (Parallel Computing Toolbox)を参照してください。

バージョン履歴

R2006a より前に導入

R2025a: 2 入力の構文でサイズが 1 より大きい最初の次元のみに沿って差分を計算

diff 関数の 2 入力の構文 diff(X,n) で、n のすべての値について、サイズが 1 より大きい最初の次元のみに沿って演算が行われるようになりました。その次元よりも n が大きい場合、diff は空の配列を返します。以前は、サイズが 1 より大きい最初の次元よりも n が大きい場合、diff(X,n) はその次元のサイズを 1 に減らしてから、サイズが 1 より大きい次の次元に沿って差分を計算していました。

R2023a: table および timetable で直接計算を実行

関数 diff は、table または timetable 内のすべての変数に対して、それらの変数にアクセスするためのインデックス付けを行うことなく計算できます。すべての変数のデータ型で計算がサポートされている必要があります。詳細については、table および timetable での直接計算を参照してください。

R2022a: 要素数が多い場合のパフォーマンスが向上

関数 diff は、10⁵ 個以上の要素をもつベクトルの演算を行う場合や、5 x 10⁵ 個以上の要素をもつ行列または多次元配列の 1 番目または 2 番目の次元に沿って演算を行う場合のパフォーマンスが向上しました。

たとえば、次のコードは、2.5 x 10⁷ 個の要素をもつ double を構築し、隣接する要素間の差を計算します。直前のリリースと比較して、次のコードは約 2.4 倍速くなっています。

function timingDiff
rng default
N = 5000;
A = rand(N);

tic
for k = 1:40
   D = diff(A);
end
toc
end

おおよその実行時間は以下のとおりです。

R2021b: 2.43 秒

R2022a: 1.00 秒

コードは、Windows^® 10、Intel^® Xeon^® CPU E5-1650 v4 (3.60 GHz) のテスト システム上で関数 timingDiff を呼び出して時間測定されました。

参考

gradient | prod | cumsum | sum