octave: scripts/polynomial/polyaffine.m annotate

annotate scripts/polynomial/polyaffine.m @ 14237:11949c9795a0

Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Add clf() to all demos using plot features to get reproducibility. Use 64 as input to all colormaps (jet (64)) to get reproducibility. * bicubic.m, cell2mat.m, celldisp.m, cplxpair.m, interp1.m, interp2.m, interpft.m, interpn.m, profile.m, profshow.m, convhull.m, delaunay.m, griddata.m, inpolygon.m, voronoi.m, autumn.m, bone.m, contrast.m, cool.m, copper.m, flag.m, gmap40.m, gray.m, hot.m, hsv.m, image.m, imshow.m, jet.m, ocean.m, pink.m, prism.m, rainbow.m, spring.m, summer.m, white.m, winter.m, condest.m, onenormest.m, axis.m, clabel.m, colorbar.m, comet.m, comet3.m, compass.m, contour.m, contour3.m, contourf.m, cylinder.m, daspect.m, ellipsoid.m, errorbar.m, ezcontour.m, ezcontourf.m, ezmesh.m, ezmeshc.m, ezplot.m, ezplot3.m, ezpolar.m, ezsurf.m, ezsurfc.m, feather.m, fill.m, fplot.m, grid.m, hold.m, isosurface.m, legend.m, loglog.m, loglogerr.m, pareto.m, patch.m, pbaspect.m, pcolor.m, pie.m, pie3.m, plot3.m, plotmatrix.m, plotyy.m, polar.m, quiver.m, quiver3.m, rectangle.m, refreshdata.m, ribbon.m, rose.m, scatter.m, scatter3.m, semilogx.m, semilogxerr.m, semilogy.m, semilogyerr.m, shading.m, slice.m, sombrero.m, stairs.m, stem.m, stem3.m, subplot.m, surf.m, surfc.m, surfl.m, surfnorm.m, text.m, title.m, trimesh.m, triplot.m, trisurf.m, uigetdir.m, uigetfile.m, uimenu.m, uiputfile.m, waitbar.m, xlim.m, ylim.m, zlim.m, mkpp.m, pchip.m, polyaffine.m, spline.m, bicgstab.m, cgs.m, gplot.m, pcg.m, pcr.m, treeplot.m, strtok.m, demo.m, example.m, rundemos.m, speed.m, test.m, calendar.m, datestr.m, datetick.m, weekday.m: Revamp %!demos to use Octave coding conventions on spacing, etc.

author	Rik <octave@nomad.inbox5.com>
date	Fri, 20 Jan 2012 12:59:53 -0800
parents	72c96de7a403
children	ce2b59a6d0e5

rev	line source
14138 72c96de7a403 maint: update copyright notices for 2012 John W. Eaton <jwe@octave.org> parents: 14104 diff changeset	1 ## Copyright (C) 2009-2012 Tony Richardson, Jaroslav Hajek
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	2 ##
11104 2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	3 ## This file is part of Octave.
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	4 ##
11104 2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	5 ## Octave is free software; you can redistribute it and/or modify it
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	6 ## under the terms of the GNU General Public License as published by
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	7 ## the Free Software Foundation; either version 3 of the License, or (at
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	8 ## your option) any later version.
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	9 ##
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	10 ## Octave is distributed in the hope that it will be useful, but
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	13 ## General Public License for more details.
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	14 ##
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	15 ## You should have received a copy of the GNU General Public License
11104 2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	16 ## along with Octave; see the file COPYING. If not, see
2c356a35d7f5 fix copyright notices John W. Eaton <jwe@octave.org> parents: 10821 diff changeset	17 ## <http://www.gnu.org/licenses/>.
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	18
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	19 ## -- texinfo --
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	20 ## @deftypefn {Function File} {} polyaffine (@var{f}, @var{mu})
11538 6eded7561d9d Add undocumented polyaffine function to documentation. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	21 ## Return the coefficients of the polynomial vector @var{f} after an affine
11575 d6619410e79c Spellcheck documentation before 3.4 release. Rik <octave@nomad.inbox5.com> parents: 11538 diff changeset	22 ## transformation. If @var{f} is the vector representing the polynomial f(x),
11538 6eded7561d9d Add undocumented polyaffine function to documentation. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	23 ## then @code{@var{g} = polyaffine (@var{f}, @var{mu})} is the vector
6eded7561d9d Add undocumented polyaffine function to documentation. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	24 ## representing:
10821 693e22af08ae Grammarcheck documentation of m-files Rik <octave@nomad.inbox5.com> parents: 9758 diff changeset	25 ##
9116 a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	26 ## @example
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	27 ## g(x) = f((x-@var{mu}(1))/@var{mu}(2)).
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	28 ## @end example
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11575 diff changeset	29 ##
14104 614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	30 ## @seealso{polyval, polyfit}
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	31 ## @end deftypefn
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	32
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	33
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	34 function g = polyaffine (f, mu)
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	35
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	36 if (nargin != 2)
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	37 print_usage ();
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	38 endif
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	39
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	40 if (! isvector (f))
11588 d5bd2766c640 style fixes for warning and error messages in script files John W. Eaton <jwe@octave.org> parents: 11587 diff changeset	41 error ("polyaffine: F must be a vector");
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	42 endif
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	43
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	44 if (! isvector (mu) \|\| length (mu) != 2)
11588 d5bd2766c640 style fixes for warning and error messages in script files John W. Eaton <jwe@octave.org> parents: 11587 diff changeset	45 error ("polyaffine: MU must be a two-element vector");
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	46 endif
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	47
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	48 lf = length (f);
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	49
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	50 ## Ensure that f is a row vector
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	51 if (rows (f) > 1)
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	52 f = f.';
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	53 endif
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	54
9116 a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	55 g = f;
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	56
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	57 ## Scale.
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	58 if (mu(2) != 1)
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	59 g = g ./ (mu(2) .^ (lf-1:-1:0));
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	60 endif
79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	61
9116 a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	62 ## Translate.
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	63 if (mu(1) != 0)
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	64 w = (-mu(1)) .^ (0:lf-1);
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	65 ii = lf:-1:1;
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	66 g = g(ii) * (toeplitz (w) .* pascal (lf, -1));
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	67 g = g(ii);
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	68 endif
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	69
9114 79e90381b402 merge polytrans & polyscale -> polyaffine, make compatible with polyval Jaroslav Hajek <highegg@gmail.com> parents: diff changeset	70 endfunction
9116 a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	71
11538 6eded7561d9d Add undocumented polyaffine function to documentation. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	72
14104 614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	73 %!demo
614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	74 %! f = [1/5 4/5 -7/5 -2];
614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	75 %! g = polyaffine (f, [1, 1.2]);
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	76 %! x = linspace (-4,4,100);
11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	77 %! plot (x,polyval(f, x), x,polyval(g, x));
14104 614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	78 %! legend ("original", "affine");
614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	79 %! axis ([-4 4 -3 5]);
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	80 %! grid on;
14104 614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	81
9116 a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	82 %!test
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	83 %! f = [1/5 4/5 -7/5 -2];
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	84 %! mu = [1, 1.2];
a8be2f7c81ee improve docs for polyaffine, fix bugs Jaroslav Hajek <highegg@gmail.com> parents: 9114 diff changeset	85 %! g = polyaffine (f, mu);
14237 11949c9795a0 Revamp %!demos in m-files to use Octave coding conventions on spacing, etc. Rik <octave@nomad.inbox5.com> parents: 14138 diff changeset	86 %! x = linspace (-4,4,100);
14104 614505385171 doc: Overhaul docstrings for polynomial functions. Rik <octave@nomad.inbox5.com> parents: 11588 diff changeset	87 %! assert (polyval (f, x, [], mu), polyval (g, x), 1e-10);
11538 6eded7561d9d Add undocumented polyaffine function to documentation. Rik <octave@nomad.inbox5.com> parents: 11523 diff changeset	88

Mercurial > octave

annotate scripts/polynomial/polyaffine.m @ 14237:11949c9795a0