2017-02-13 5 views

Répondre

1

En effet, sur mon Scilab (5.5.2), "interp1" juste bandes la partie imaginaire.

Essayez de remplacer:

res = interp1(x,y,xx) 

par:

c = size(y)*[0;1] 
D = eye(c,c) 
res = interp1(x,[real(y),imag(y)],xx)*[D;%i*D] 

Il devrait fonctionner avec toutes les dimensions pour y. Vous pouvez en faire un seul revêtement si vous le souhaitez, mais il sera moins lisible.

ou (plus simple, mais peut-être plus lent):

res = interp1(x,real(y),xx) + %i * interp1(x,imag(y),xx)