Re: A new version of bicubic in Octave

[Paul Kienzle]

Fixed bug in line:

   t = 1 + (YI-Y(1))./((Y(rz)-Y(1))/(rz-1));

One of the Y was an X.


Doesn't work for X non-uniform.  E.g.:

 A=[13,-1,12;5,4,3;1,6,2];
 x=[0,1,4]+10; y=[-10,-9,-8];
 xi=linspace(min(x),max(x),17);
 yi=linspace(min(y),max(y),26);
 mesh(xi,yi,interp2(x,y,A,xi,yi,'cubic'));
 [x,y] = meshgrid(x,y);
 __gnuplot_raw__ ("set nohidden3d;\n")
 hold on; plot3(x(:),y(:),A(:),"b*"); hold off;

- Paul


On Aug 28, 2005, at 9:43 AM, Hoxide Ma wrote:

> It's really a new version , not based on any commerce
> codes.
>
> It's based on my note:
>
> http://wiki.woodpecker.org.cn/moin/Octave/Bicubic? 
> action=AttachFile&do=get&target=bicubic.pdf
>
> You can test my code by run intt1.m and intt2.m .
>
> intt2.m is a example to interpolation  by a
> anomalistic surface.
>
> BTW: I'm sorry for my poor English :P
>
> ==========================================================
>
> A student of Math in Soochow University in China.
> Intrested in Math, Python, Octave, Lisp, Maxima, Prolog and .NET/MONO.
>
> Blog: http://blog.sina.com.cn/blog/1142604745
> Wiki: http://wiki.woodpecker.org.cn/moin/Hoxide
>
>
>         
>
>         
>                 
> ___________________________________________________________
> ~{QE;"Cb7Q~}G~{SJOd#-VP9z5ZR;>xN^@,;xSJ<~I'HE3,4sSJOd~}
> http://cn.mail.yahoo.com
> <interp2.m><intt1.m><intt2.m><bicubic.m>



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------