Matlab Интерполяция матрицы 4D

Question

Я хотел бы интерполировать 4D-матрицу. Размер матрицы - 400X200X400X200 комплексных чисел. Im используя функцию: interpn

И для этого я должен установить установить свою новую сетку с помощью: ndgrid

Но ndgrid занимает часы. Не могли бы вы помочь?

Спасибо!

Пример кода:

[Y X V U]=ndgrid (ky_vec,kx_vec,kv_vec,ku_vec); 
F=f(Y,X,V,U); 
[Yt Xt Vt Ut]=ndgrid (ky_vec,kx_vec,ky_vec,kx_vec,'single'); 
A=interpn(Y,X,V,U,F,0.9*Yt,0.9*Xt,(1-0.9)*Vt,(1-0.9)*Ut,'linear',0.0); 

Answer 1

'single' аргумент не принадлежит там. Вместо этого он интерпретируется как ones(1,length('single')) и, таким образом, дает даже матрицы huger.

Вы также можете сгенерировать griddedInterpolant вместо, который бы сохранить память:

F_interp = griddedInterpolant({ky_vec,kx_vec,kv_vec,ku_vec},F)