discrete second order derivative operator for unequally spaced data

4 次查看（过去 30 天）
David Kusnirak
David Kusnirak 2013-2-14
Hi,
I need help to obtain a discrete second order derivative operator, that could be applied on a matrix. The best way, I guess, is to reshape the matrix into vector form and create a diagonal matrix with the operator, which is easy for equally spaced data in the matrix,but how can I incorporate the distances (width and heigth of the cells) for unequally spaced data?
Thanks
  2 个评论
Matt J
Matt J 2013-2-14
I don't see how a diagonal matrix could ever give you a derivative operator, even with equally spaced data. A diagonal matrix can only weight each element of the operand independently, not take differences between them. Please describe the organization of your data more clearly. Does each column of your matrix contain signal samples to be differentiated? Is the spacing between these samples the same for all columns?
David Kusnirak
David Kusnirak 2013-2-15
sorry, I wrote it unclear, I meant diagonal matrix with 4 off-diagonals.
example:
if ma data matrix is
M = magic(3)
and the width and height of each cell is 1, the operator looks like:
L=[ -4 1 0 0 1 0 0 0 0;
1 -4 1 0 0 1 0 0 0;
0 1 -4 1 0 0 1 0 0;
0 0 1 -4 1 0 0 1 0;
1 0 0 1 -4 1 0 0 1;
0 1 0 0 1 -4 1 0 0;
0 0 1 0 0 1 -4 1 0;
0 0 0 1 0 0 1 -4 1;
0 0 0 0 1 0 0 1 -4];
so the product of the operator (in vector form) is
dz = L * reshape(m,1,[]);
and back to matrix
dz = reshape(dz,3,3)
My question is how to implement proper geometry of the cells, for example width = 2*heigth

请先登录，再进行评论。

请先登录，再回答此问题。

回答（4 个）

Teja Muppirala
Teja Muppirala 2013-2-15
Could you possibly use the GRADIENT command? This allows you to pass in unevenly spaced values for X and Y.
help gradient
Youssef Khmou
Youssef Khmou 2013-2-15
Discrete Laplacian :
del2(matrix)
Example :
[x,y]=meshgrid(-4:.1:4);
z=exp(-x.^2-y.^2);
dz2=del2(z);
Jan
Jan 2013-2-15
A fast gradient with 2nd order method for unevenly distributed data: FEX: DGradient.
David Kusnirak
David Kusnirak 2013-2-15
Thanks for hints. However I'm not as interested in the solution as in the operator itself, I need to know how does it look and work.
thanks for any help

类别

在 Help Center 和 File Exchange 中查找有关 Linear Algebra 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by