Commit 8b91ab20 authored by Rafael Laboissiere's avatar Rafael Laboissiere

New upstream version 0.1.3

parent 937783f4
Name: lssa
Version: 0.1.2
Date: 2012-08-17
Version: 0.1.3
Date: 2018-07-23
Author: Ben Lewis <benjf5@gmail.com>
Maintainer: Ben Lewis <benjf5@gmail.com>
Maintainer: Ben Lewis <benjf5@gmail.com>, John Donoghue <john.donoghue@ieee.org>
Title: Least squares spectral analysis
Description: A package implementing tools to compute spectral decompositions of
irregularly-spaced time series. Currently includes functions based off the
......
## Copyright 2015-2016 Carnë Draug
## Copyright 2015-2016 Oliver Heimlich
## Copyright 2017 Julien Bect <jbect@users.sf.net>
## Copyright 2017 Olaf Till <i7tiol@t-online.de>
## Copyright 2018 John Donoghue <john.donoghue@ieee.org>
##
## Copying and distribution of this file, with or without modification,
## are permitted in any medium without royalty provided the copyright
## notice and this notice are preserved. This file is offered as-is,
## without any warranty.
## Some basic tools (can be overriden using environment variables)
SED ?= sed
TAR ?= tar
GREP ?= grep
CUT ?= cut
TR ?= tr
## Note the use of ':=' (immediate set) and not just '=' (lazy set).
## http://stackoverflow.com/a/448939/1609556
package := $(shell $(GREP) "^Name: " DESCRIPTION | $(CUT) -f2 -d" " | \
$(TR) '[:upper:]' '[:lower:]')
version := $(shell $(GREP) "^Version: " DESCRIPTION | $(CUT) -f2 -d" ")
## These are the paths that will be created for the releases.
target_dir := target
release_dir := $(target_dir)/$(package)-$(version)
release_tarball := $(target_dir)/$(package)-$(version).tar.gz
html_dir := $(target_dir)/$(package)-html
html_tarball := $(target_dir)/$(package)-html.tar.gz
## Using $(realpath ...) avoids problems with symlinks due to bug
## #50994 in Octaves scripts/pkg/private/install.m. But at least the
## release directory above is needed in the relative form, for 'git
## archive --format=tar --prefix=$(release_dir).
real_target_dir := $(realpath .)/$(target_dir)
installation_dir := $(real_target_dir)/.installation
package_list := $(installation_dir)/.octave_packages
install_stamp := $(installation_dir)/.install_stamp
## These can be set by environment variables which allow to easily
## test with different Octave versions.
ifndef OCTAVE
OCTAVE := octave
endif
OCTAVE := $(OCTAVE) --no-gui --silent --norc
MKOCTFILE ?= mkoctfile
## Command used to set permissions before creating tarballs
FIX_PERMISSIONS ?= chmod -R a+rX,u+w,go-w,ug-s
## Detect which VCS is used
vcs := $(if $(wildcard .hg),hg,$(if $(wildcard .git),git,unknown))
ifeq ($(vcs),hg)
release_dir_dep := .hg/dirstate
endif
ifeq ($(vcs),git)
release_dir_dep := .git/index
endif
## .PHONY indicates targets that are not filenames
## (https://www.gnu.org/software/make/manual/html_node/Phony-Targets.html)
.PHONY: help
## make will display the command before runnning them. Use @command
## to not display it (makes specially sense for echo).
help:
@echo "Targets:"
@echo " dist - Create $(release_tarball) for release."
@echo " html - Create $(html_tarball) for release."
@echo " release - Create both of the above and show md5sums."
@echo " install - Install the package in $(installation_dir), where it is not visible in a normal Octave session."
@echo " check - Execute package tests."
@echo " doctest - Test the help texts with the doctest package."
@echo " run - Run Octave with the package installed in $(installation_dir) in the path."
@echo " clean - Remove everything made with this Makefile."
##
## Recipes for release tarballs (package + html)
##
.PHONY: release dist html clean-tarballs clean-unpacked-release
## To make a release, build the distribution and html tarballs.
release: dist html
md5sum $(release_tarball) $(html_tarball)
@echo "Upload @ https://sourceforge.net/p/octave/package-releases/new/"
@echo " and note the changeset the release corresponds to"
## dist and html targets are only PHONY/alias targets to the release
## and html tarballs.
dist: $(release_tarball)
html: $(html_tarball)
## An implicit rule with a recipe to build the tarballs correctly.
%.tar.gz: %
$(TAR) -c -f - --posix -C "$(target_dir)/" "$(notdir $<)" | gzip -9n > "$@"
clean-tarballs:
@echo "## Cleaning release tarballs (package + html)..."
-$(RM) $(release_tarball) $(html_tarball)
@echo
## Create the unpacked package.
##
## Notes:
## * having ".hg/dirstate" (or ".git/index") as a prerequesite means it is
## only rebuilt if we are at a different commit.
## * the variable RM usually defaults to "rm -f"
## * having this recipe separate from the one that makes the tarball
## makes it easy to have packages in alternative formats (such as zip)
## * note that if a commands needs to be run in a specific directory,
## the command to "cd" needs to be on the same line. Each line restores
## the original working directory.
$(release_dir): $(release_dir_dep)
-$(RM) -r "$@"
ifeq (${vcs},hg)
hg archive --exclude ".hg*" --type files "$@"
endif
ifeq (${vcs},git)
git archive --format=tar --prefix="$@/" HEAD | $(TAR) -x
$(RM) "$@/.gitignore"
endif
## Don't fall back to run the supposed necessary contents of
## 'bootstrap' here. Users are better off if they provide
## 'bootstrap'. Administrators, checking build reproducibility, can
## put in the missing 'bootstrap' file if they feel they know its
## necessary contents.
ifneq (,$(wildcard src/bootstrap))
cd "$@/src" && ./bootstrap && $(RM) -r "autom4te.cache"
endif
## Uncomment this if your src/Makefile.in has these targets for
## pre-building something for the release (e.g. documentation).
# cd "$@/src" && ./configure && $(MAKE) prebuild && \
# $(MAKE) distclean && $(RM) Makefile
##
${FIX_PERMISSIONS} "$@"
run_in_place = $(OCTAVE) --eval ' pkg ("local_list", "$(package_list)"); ' \
--eval ' pkg ("load", "$(package)"); '
# html_options = --eval 'options = get_html_options ("octave-forge");'
## Uncomment this for package documentation.
##html_options = --eval 'options = get_html_options ("octave-forge");' \
## --eval 'options.package_doc = "$(package).texi";'
html_options = --eval 'options = get_html_options ("octave-forge");'
$(html_dir): $(install_stamp)
$(RM) -r "$@";
$(run_in_place) \
--eval ' pkg load generate_html; ' \
$(html_options) \
--eval ' generate_package_html ("$(package)", "$@", options); ';
$(FIX_PERMISSIONS) "$@";
clean-unpacked-release:
@echo "## Cleaning unpacked release tarballs (package + html)..."
-$(RM) -r $(release_dir) $(html_dir)
@echo
##
## Recipes for installing the package.
##
.PHONY: install clean-install
octave_install_commands = \
' llist_path = pkg ("local_list"); \
mkdir ("$(installation_dir)"); \
load (llist_path); \
local_packages(cellfun (@ (x) strcmp ("$(package)", x.name), local_packages)) = []; \
save ("$(package_list)", "local_packages"); \
pkg ("local_list", "$(package_list)"); \
pkg ("prefix", "$(installation_dir)", "$(installation_dir)"); \
pkg ("install", "-local", "-verbose", "$(release_tarball)"); '
## Install unconditionally. Maybe useful for testing installation with
## different versions of Octave.
install: $(release_tarball)
@echo "Installing package under $(installation_dir) ..."
$(OCTAVE) --eval $(octave_install_commands)
touch $(install_stamp)
## Install only if installation (under target/...) is not current.
$(install_stamp): $(release_tarball)
@echo "Installing package under $(installation_dir) ..."
$(OCTAVE) --eval $(octave_install_commands)
touch $(install_stamp)
clean-install:
@echo "## Cleaning installation under $(installation_dir) ..."
-$(RM) -r $(installation_dir)
@echo
##
## Recipes for testing purposes
##
.PHONY: run doctest check
## Start an Octave session with the package directories on the path for
## interactice test of development sources.
run: $(install_stamp)
$(run_in_place) --persist
## Test example blocks in the documentation. Needs doctest package
## https://octave.sourceforge.io/doctest/index.html
doctest: $(install_stamp)
$(run_in_place) --eval 'pkg load doctest;' \
--eval "targets = '$(shell (ls inst; ls src | $(GREP) .oct) | $(CUT) -f2 -d@ | $(CUT) -f1 -d.)';" \
--eval "targets = strsplit (targets, ' '); doctest (targets);"
## Test package.
octave_test_commands = \
' pkgs = pkg("list", "$(package)"); \
cd ("$(target_dir)/"); \
dirs = {sprintf(".installation/%s-%s", pkgs{1}.name, pkgs{1}.version)}; \
__run_test_suite__ (dirs, {}); '
## the following works, too, but provides no overall summary output as
## __run_test_suite__ does:
##
## else cellfun (@runtests, horzcat (cellfun (@ (dir) ostrsplit (([~, dirs] = system (sprintf ("find %s -type d", dir))), "\n\r", true), dirs, "UniformOutput", false){:})); endif '
check: $(install_stamp)
$(run_in_place) --eval $(octave_test_commands)
clean-check:
@echo "## Removing fntests.log..."
-$(RM) $(target_dir)/fntests.log
@echo
##
## CLEAN
##
.PHONY: clean
clean: clean-tarballs clean-unpacked-release clean-install clean-check
@echo "## Removing target directory (if empty)..."
-rmdir $(target_dir)
@echo
@echo "## Cleaning done"
@echo
Summary of changes in lssa 0.1.3:
** Mark fastlscomplex BIT as a known failure (Bug #53963)
** Update package with maintainer makefile
** Code updates to support install in Octave 4.4
Summary of changes in lssa 0.1.2:
** All functions now have input checks in place to return useful errors as
......
Welcome to the first release of lssa, 0.1.1
Current status:
** lscomplex and lsreal both produce accurate results; they can be slow for
very large datasets.
** fastlscomplex is accurate for the first octave of results; there is still an
error I need to pin down in the merging for additional octaves. fastlsreal
is disabled at the moment as I move to an implementation based on the new
fastlscomplex.
** lscorrcoeff works, although I'm still attempting to understand the initial
author's reasoning. Its generated results are relevant to any given data
set, but it does not appear to be normalized to any great extent.
** There are two wavelet functions under development, but they are not included
in this release as they are currently not functional. For all your wavelet
needs, the specific transformation used is available in the lswaveletcoeff
function, and will generate a single cosine/sine magnitude pair (as a
complex number) for a complex-valued series (this function may be joined by
a companion for real-valued series) and can be looped to simulate a full
wavelet transform.
** For all the working functions, tests have been written and formatted to
Octave coding standards. These tests should pass on any given architecture
(there was some question about that previously) and often provide examples
of how the function operates. For a few functions, there are demo scripts.
This diff is collapsed.
149.1,5679,2342,284.7
173.1,6828,3634,272.8
177.4,7043,3833,268.1
228.6,9523,6220,262.2
250.3,10579,7327,254.6
266,11334,8113,259.6
302.6,13449,10123,261.6
321.2,14538,11013,263.7
331.6,15208,11326,244.8
342.1,15922,11719,238.3
365.5,17747,13405,236.2
375.6,18580,13989,225.3
443.5,24315,17695,182.2
470.8,26578,19988,189.2
506.4,29630,22977,191.6
544.7,32844,26303,188.5
554.2,33645,27062,191.7
602.3,37421,31447,205.4
625.6,39310,33884,209.1
700.2,45242,39880,209.1
765,50610,44766,189.3
788.2,52446,47024,188.4
800,53436,48229,210.1
812.2,54474,49414,215.7
834.8,56300,51174,190.4
902.2,61783,57068,221.8
912,62689,57799,210.4
986.2,69618,63687,195.4
1011.3,71767,65701,191.4
1023.5,72797,66883,195
1087.2,78042,72849,227.4
1112.5,80057,75360,229.2
1162.1,83615,78995,217.1
1175,84515,80059,221.8
1209.9,86840,82858,231
1237.2,88808,84929,241.1
1251.5,89864,85727,236.4
1261.2,90609,86323,228.1
1274.2,91560,87180,214.2
1289.2,92632,88051,217
1309.2,94039,89363,208
1338.2,96047,91691,224.3
1349,96791,92460,228.4
1387.2,99498,95349,232.1
1451.5,103733,99842,225.9
1463.2,104566,100833,230.9
1476.1,105492,101829,236.9
1505,107599,103372,228.2
1526.3,109404,105213,236.9
1542.1,110674,106203,230.7
1575.2,113363,108308,238.2
1582.8,113952,108994,245.7
1598,115077,110253,251.3
1615,116228,111456,256.8
1627.9,117072,112577,266.3
1637.6,117671,113472,261.4
1644,118074,114082,274.6
1651,118499,114738,273.3
1669.2,119601,116175,262.5
1687.2,120680,117519,267.6
1700.9,121485,118396,273.8
1716,122345,119273,272
1726.8,122965,120002,265.2
1736.8,123535,120652,277.7
1758.2,124721,121961,272.2
1770,125380,122606,276.5
1789.2,126447,123815,268.7
1790,126491,123858,266.6
1799,126955,124306,266.3
1804,127210,124571,279.8
1810,127510,124876,277.2
1825.7,128293,125746,273.8
1830,128501,126023,267.1
1836,128804,126475,262.5
1841.6,129097,126809,262.6
1852.4,129674,127445,275.4
1869.3,130599,128300,274.1
1870.2,130653,128399,287.1
1875.9,130998,128652,286.8
1882.5,131406,129007,282.7
1890,131908,129411,264.1
1895,132264,129755,263.4
1903.5,132873,130167,259
1932,135308,131789,240.4
1954.5,137601,133334,224
1969.8,139285,134205,208.9
1980.2,140452,135003,204.6
1987.4,141257,135683,198.1
1990.6,141615,135976,201.8
1994.6,142064,136359,202.5
1998,142440,136659,195.9
2005.8,143315,137383,194.4
2009.5,143732,137694,193.4
2015,144346,138226,190.2
2025.7,145527,139445,192.3
2041.5,147292,141312,196.5
2050.3,148287,142357,190.4
2077.5,151234,145435,197
2116,155625,150303,191.9
2157,160293,154471,189
2164,161077,155299,185.5
2203,165646,160494,204.4
2225,168391,162996,191.6
2247,170922,165278,183.8
2280,174920,169870,197.9
2302,177644,172596,197.8
2325,180546,175440,190.3
2333,181502,176271,190.1
2348,183525,178550,207.7
2363,185640,180779,213.2
2372,186927,181617,217.7
2386,188892,183355,199.8
2399,190587,185063,203.5
2414,192382,187199,210.7
2425,193816,189335,231.4
2437,195298,191057,231.5
2451,197086,192632,218
2475,199904,195625,220.1
2499,202496,199025,242.6
2525,205458,202212,251
2533,206496,203191,239.1
2543,207803,204283,247.7
2552.01,208947,205148,244.4
2557.71,209648,205715,232.2
2560.91,210053,206119,228.7
2574.21,211628,207991,238.2
2584.71,212846,209414,242.2
2588.81,213315,210022,244.6
2595.11,214036,210830,247.3
2596.71,214219,211005,252
2606.61,215343,212281,257.4
2621.71,217099,214153,251.2
2629.41,217989,215041,241.4
2634.41,218602,215593,240.3
2636.71,218908,215879,242.7
2640.41,219393,216459,247.5
2644.41,219935,217009,251.7
2646.61,220255,217271,251.2
2650.41,220822,217676,245.4
2656.21,221718,218342,240.5
2666.71,223548,219680,212.2
2670.41,224240,220182,216.2
2674.61,224993,220760,207.2
2677.41,225469,221054,208.9
2682.61,226339,221612,205.7
2691.01,227885,222958,203.4
2693.61,228349,223446,215.7
2698.01,229167,224630,236.9
2701.41,229802,225299,234.5
2702.71,230039,225509,233.1
2705.61,230547,225888,224.5
2711.71,231601,226710,232.4
2715.41,232206,227384,233.9
2717.71,232571,227840,241.7
2732.71,234795,230703,245.2
2735.71,235232,231382,252.2
2738.71,235652,231990,241.4
2741.71,236077,232570,247.4
2744.61,236467,233102,243.1
2747.61,236866,233646,239.2
2751.11,237301,234126,245.7
2753.61,237602,234470,245.9
2756.21,237889,234781,247.4
2759.11,238206,235213,252.9
2765.21,238908,236236,259.8
2773.51,240006,237831,279
2782.71,241535,238935,263.8
2785.51,242092,239250,252.4
2788.51,242675,239545,249.9
2794.51,243813,240201,230.4
2797.51,244446,240577,219.4
2806.51,246379,242068,214.7
2815.61,248364,243653,200.2
2818.61,249046,244215,213.9
2821.51,249670,244863,195.4
2824.51,250309,245483,196.7
2833.81,252279,247447,199
2836.51,252879,248087,201.9
2839.51,253559,248980,204
2845.51,254873,250461,203.9
2851.51,256182,251521,209.7
2857.51,257604,252959,208.9
2860.51,258351,253880,214.7
2866.51,259882,255233,228.2
2870.51,260936,256053,199.9
2872.71,261526,256501,211.7
2876.21,262455,257247,188.7
2881.42,263844,258477,194.2
2884.51,264666,259228,198.9
2887.51,265457,259958,184.7
2890.51,266255,260754,190.4
2893.51,267039,261595,193.9
2896.51,267825,262411,194.2
2899.51,268558,263207,198.4
2902.51,269307,264046,193.2
2905.51,270052,264834,202.2
2911.46,271454,266492,211
2914.51,272171,267434,215.4
2919.41,273256,268679,223.7
2926.51,274802,270680,231.4
2935.51,276731,273012,226.4
2941.51,278017,274445,230.4
2944.51,278711,275218,231
2956.51,281740,277925,220.4
2959.51,282539,278602,217.2
2964.51,283810,279543,207.7
2973.81,286475,282301,212.7
2979.51,288013,283492,213.2
2988.61,290554,286217,224.4
2994.51,292121,287846,236.2
3003.51,294495,290571,240.2
3007.01,295477,291769,240.7
3009.51,296158,292474,250.2
3015.51,297810,293676,244.9
3018.51,298623,294615,225.9
3022.71,299716,295849,227.9
3027.51,300928,297131,233.2
3030.41,301665,298051,237.9
3033.51,302452,299020,239
3036.51,303201,299877,241.9
3039.51,303939,300646,251.7
3042.51,304745,301496,256.8
3045.51,305612,302456,257.2
3048.56,306561,303334,246.9
3051.51,307464,303953,272.7
3054.51,308358,304590,251.7
3057.71,309281,305306,244.7
3063.51,310905,307131,255.9
3066.51,311704,308101,249.2
3072.51,313292,310039,256.3
3075.41,314032,310930,260.4
3078.51,314822,311774,260.3
3084.51,316304,313493,266.3
3090.51,317775,315143,266.2
3093.51,318509,315940,270.2
3096.46,319231,316681,271.9
3099.51,319978,317445,275.2
3105.51,321448,318980,265
3109.01,322216,319754,271.8
3111.51,322746,320378,272.7
3114.81,323441,321386,273.2
3117.51,324021,322111,282.4
3119.51,324461,322582,289.2
3120.61,324711,322827,288.4
3123.51,325400,323485,298.7
3126.51,326200,324189,278.2
3129.91,327237,324991,285.8
3132.41,328058,325527,278.7
3135.51,329135,326239,270.5
3138.51,330278,327114,255.7
3141.51,331513,328097,241.9
3145.01,333111,329267,239.7
3147.51,334356,330208,234.2
3153.51,337304,332293,250.2
3156.51,338770,333627,200.7
3159.51,340242,335290,205.2
3162.81,341849,336972,204.9
3169.01,344844,340165,220.4
3174.51,347589,342998,221.2
3177.81,349298,344735,216.2
3183.41,352377,347610,209.2
3189.51,355795,350765,193
3192.51,357450,352412,186.2
3200.01,361445,356838,201.2
3204.71,363962,359688,206.4
3210.51,367048,362766,201.9
3216.45,370095,366221,214.7
3222.51,373629,369563,229.7
3228.91,377177,373014,227
3231.51,378469,374561,240
3237.51,381834,378194,246.9
3240.51,383395,379633,245.9
3249.51,388757,384909,264.7
3252.45,390641,386579,259.3
3258.51,394634,390589,255.2
3261.51,396423,392451,250.2
3264.51,398091,394628,266.3
3267.51,399733,396713,274.7
3273.81,403173,400390,278
3283.51,408236,405844,279.7
3289.45,411202,409022,283.7
3292.91,413010,410831,276.3
3299.01,416332,414085,285.5
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
......@@ -32,7 +32,7 @@
function a = cubicwgt (s)
if (nargin != 1)
print_usage ();
print_usage ();
endif
## s is the value/vector/matrix to be windowed
......
......@@ -25,16 +25,13 @@
function coeff = lombcoeff (T, X, o)
if (nargin != 3)
print_usage ();
endif
if (! all (size (T) == size (X)))
error ("lombcoeff: Time series vectors of uneven size.\n");
endif
if (! isscalar (o))
error ("lombcoeff: Supplied frequency is not a scalar.\n");
endif
if (o == 0)
error ("lombcoeff: Supplied frequency is not a frequency.\n");
print_usage ();
elseif (! all (size (T) == size (X)))
error ("lombcoeff: Time series vectors of uneven size");
elseif (! isscalar (o))
error ("lombcoeff: Supplied frequency is not a scalar");
elseif (o == 0)
error ("lombcoeff: Supplied frequency is not a frequency");
endif
oT = o .* T;
......
......@@ -2,7 +2,7 @@
##
## This program 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 of the License, or (at your option) any later
## Foundation; either version 3 of the License, or (at your option) any later
## version.
##
## This program is distributed in the hope that it will be useful, but WITHOUT
......@@ -27,28 +27,25 @@
function coeff = lombnormcoeff (T, X, omega)
if (nargin != 3)
print_usage ();
endif
if (! all (size (T) == size (X)))
error ("lombnormcoeff: Time series vectors of uneven size.\n");
endif
if (! isscalar (omega))
error ("lombnormcoeff: Supplied frequency is not a scalar.\n");
endif
if (omega == 0)
error ("lombnormcoeff: Supplied frequency is not a frequency.\n");
print_usage ();
elseif (! all (size (T) == size (X)))
error ("lombnormcoeff: Time series vectors of uneven size");
elseif (! isscalar (omega))
error ("lombnormcoeff: Supplied frequency is not a scalar");
elseif (omega == 0)
error ("lombnormcoeff: Supplied frequency is not a frequency");
endif
xmean = mean (X);
theta = atan2 (sum (sin (2 .* omega .*T)),
sum (cos (2 .* omega .* T))) / (2*omega);
sum (cos (2 .* omega .* T))) / (2*omega);
coeff = ((sum ((X-xmean) .* cos (omega .* T - theta)) .^ 2 /
sum (cos (omega .* T - theta) .^ 2) +
sum (cos (omega .* T - theta) .^ 2) +
sum ((X-xmean) .* sin (omega .* T - theta)) .^ 2 /
sum (sin (omega .* T - theta) .^ 2 )) /
sum (sin (omega .* T - theta) .^ 2 )) /
(2 * var(X)));
endfunction
......
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
......@@ -37,25 +37,19 @@
function coeff = lswaveletcoeff (x, y, t, o, wgt = @cubicwgt, wgtrad = 1)
if (! (nargin >= 4) && (nargin <= 6))
print_usage ();
endif
if (! isvector (x))
error ("lswaveletcoeff: Time values are not a vector.\n");
endif
if (! isvector (y))
error ("lswaveletcoeff: Magnitude values are not a vector.\n");
endif
if (! all (size (x) == size (y)))
error ("lswaveletcoeff: Time series vectors of uneven size.\n");
endif
if (! isscalar (t))
error ("lswaveletcoeff: Window centre specified is not scalar.\n");
endif
if (! isscalar (o))
error ("lswaveletcoeff: Frequency specified is not scalar.\n");
endif
if (! isscalar (wgtrad))
error ("lswaveletcoeff: Window radius specified is not scalar.\n");
print_usage ();
elseif (! isvector (x))
error ("lswaveletcoeff: Time values are not a vector");
elseif (! isvector (y))
error ("lswaveletcoeff: Magnitude values are not a vector");
elseif (! all (size (x) == size (y)))
error ("lswaveletcoeff: Time series vectors of uneven size");
elseif (! isscalar (t))
error ("lswaveletcoeff: Window centre specified is not scalar");
elseif (! isscalar (o))
error ("lswaveletcoeff: Frequency specified is not scalar");
elseif (! isscalar (wgtrad))
error ("lswaveletcoeff: Window radius specified is not scalar");
endif
so = 0.05 .* o;
......
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment