gradient with irregular grid
显示 更早的评论
I have as a variable the temperature at a given depth and longitude, so it is a vector (317,1) where the 317 are all latitudes.
How do I get the gradient of this variable (which I call T) knowing that the latitude grid is irregular?
I can try
gradient_T=gradient(T,e1v) where e1v is the irregular vector?
采纳的回答
Star Strider
2023-5-2
0 个投票
The gradient function assumes a fixed step size for the second argument.
The way I calculate the numerical derivative using an irregular grid for the reference (assuming vectors here) is:
gradient_T = gradient(T) ./ gradient(e1v);
That will essentially calculate 
and generally produces the result I want.
and generally produces the result I want. .
4 个评论
Joshua Port
2025-9-10
This would probably work if the grid were irregular in both directions. If one dimension has fixed spacing and the other has irregular spacing, though, the gradient in the direction with the fixed spacing will be Inf since gradient([fixed spacing dimension]) = 0 everywhere.
@Joshua Port --
i may be missing something in your Comment, however it seems to work correctly here. (The gradient function has changed, and now permits the inclusion of the independent variable vectors as arguments, so doing element-wise division as a separate operation is no longer necessary.)
Lonv = sort(rand(1,20))*10 + 117;
Latv = linspace(35, 45, 25);
[Lat,Lon] = ndgrid(Latv, Lonv);
T = exp(-((Lon-122).^2/2.5 + (Lat - 40).^2/3))*10 + 10;
figure
surfc(Lat, Lon, T, EdgeColor='interp')
colormap(turbo)
colorbar
[gradient_T_Lon, gradient_T_Lat] = gradient(T, Lonv, Latv);
figure
surfc(Lat, Lon, gradient_T_Lon)
colormap(turbo)
colorbar
figure
surfc(Lat, Lon, gradient_T_Lat)
colormap(turbo)
colorbar
.
Joshua Port
2025-9-10
编辑：Joshua Port
2025-9-10
I should have realized what this was getting at. If your data looks like this:
x_1d = linspace(0, 100)';
x = repmat(x_1d, 1, 100);
And you do:
gradient(x)
You'll get 0 everywhere. To compute df/dx you would need to do:
[~, dx] = gradient(x);
dfdx = gradient(f) ./ dx;
This method doesn't work when one of your dimensions is curvy, though. I wish gradient could take in an array of points, even scattered ones, and a function defined at those points, and compute the gradient. I have a method for dealing with this, but it requires a lot of headache.
Star Strider
2025-9-10
Of course, the derivative of a constant is zero. Beyond that, I am not following your use of the gradient function.
Stopping here.
更多回答（1 个）
You can give the positions of the corresponding values as the second function argument. In your case the lattitude for each temperature measurement.
temperature = [30 29 28 27 26 25 24 23];
latitude = [1 2 3 4 5 10 11 20];
T = gradient(temperature, latitude)
类别
在 帮助中心 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
