< M A T L A B > Copyright 1984-2005 The MathWorks, Inc. Version 7.1.0.183 (R14) Service Pack 3 August 02, 2005 To get started, select MATLAB Help or Demos from the Help menu. >> % confronto tra due interpolazioni interp1 e interp2 >> % definiamo la funzione f >> f=@(x) 1./(1+25*x.*x) f = @(x) 1./(1+25*x.*x) >> c=interp2(20,f) c = 1.0e+03 * 0.0010 0.0000 -0.0157 -0.0000 0.1480 0.0000 -0.8050 -0.0000 2.6548 0.0000 -5.4789 -0.0000 7.1122 0.0000 -5.6385 -0.0000 2.4936 0.0000 -0.4714 -0.0000 >> % controlliamo qualche punto... >> phorn(1,c) ans = 0.0385 >> f(1) ans = 0.0385 >> % facciamo un grafico accurato >> x=-1:0.001:1; >> gg=f(x); >> size(x) ans = 1 2001 >> for k=1:2001 pp(k)=phorn(x(k),c); end >> plot(x,pp,x,gg) >>% proviamo l'altra interpolazione >> a=-1;b=1; >> c=interp1(20,a,b,f); >> for k=1:2001 pp(k)=phorn(x(k),c); end >> plot(x,pp,x,gg) >>% si vede che questa ha una oscillazione maggiore