Mercurial > hg > octave-nkf
diff scripts/signal/sinc.m @ 2303:5cffc4b8de57
[project @ 1996-06-24 09:15:24 by jwe]
author | jwe |
---|---|
date | Mon, 24 Jun 1996 09:15:24 +0000 |
parents | 3470f1e25a79 |
children | 2b5788792cad |
line wrap: on
line diff
--- a/scripts/signal/sinc.m +++ b/scripts/signal/sinc.m @@ -1,13 +1,32 @@ +### Copyright (C) 1996 John W. Eaton +### +### This file is part of Octave. +### +### Octave is free software; you can redistribute it and/or modify it +### under the terms of the GNU General Public License as published by +### the Free Software Foundation; either version 2, or (at your option) +### any later version. +### +### Octave is distributed in the hope that it will be useful, but +### WITHOUT ANY WARRANTY; without even the implied warranty of +### MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +### General Public License for more details. +### +### You should have received a copy of the GNU General Public License +### along with Octave; see the file COPYING. If not, write to the Free +### Software Foundation, 59 Temple Place - Suite 330, Boston, MA +### 02111-1307, USA. + function result = sinc (x) -# usage: sinc(x) -# -# Returns sin(pi*x)/(pi*x). + ## usage: sinc(x) + ## + ## Returns sin(pi*x)/(pi*x). -# We either need to set the do_fortran_indexing variable to "true" -# or use reshape to convert the input matrix to a vector, so that -# we can use find to determine the elements of x that equal zero. -# I prefer reshaping. + ## We either need to set the do_fortran_indexing variable to "true" + ## or use reshape to convert the input matrix to a vector, so that + ## we can use find to determine the elements of x that equal zero. + ## I prefer reshaping. [nr, nc] = size(x); @@ -15,10 +34,10 @@ x = reshape(x,nels,1); -# Set result to all ones initially. + ## Set result to all ones initially. result = ones(nels,1); -# Find non-zero elements in the input matrix. + ## Find non-zero elements in the input matrix. i = find(x); if (!isempty(i))