blob: 9611161e685385151a257dc37b112d4b5f8d9da2 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
// Copyright (C) 2018 - IIT Bombay - FOSSEE
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
/// Author : Abinash Singh Under FOSSEE Internship
// Modifieded by: Abinash Singh Under FOSSEE Internship
// Last Modified on : 3 Feb 2024
// Organization: FOSSEE, IIT Bombay
// Email: toolbox@scilab.in
function [ dSided, minH, minD, minW, maxW ] = parser ( varargin )
// Default values
// This is an helper function for findpeaks
// It parses the input arguments and returns the values of the options
dSided = %f ;
minH = %eps ;
minD = 1 ;
minW = %eps;
maxW = %inf ;
idx = 1 ;
while idx <= length(varargin)
select lower(varargin(idx))
case 'doublesided'
dSided = %t ;
idx = idx + 1 ;
case 'minpeakheight'
if idx+1 > length(varargin) || ~( isscalar(varargin(idx+1)) && varargin(idx+1) >=0 ) then
error('findpeaks: MinPeakHeight must be a postive scalar') ;
end
minH = varargin(idx+1) ;
idx = idx + 2 ;
case 'minpeakdistance'
if idx+1 > length(varargin) || ~( isscalar(varargin(idx+1)) && varargin(idx+1) >=0 ) then
error('findpeaks: MinPeakDistance must be a postive scalar') ;
end
minD = varargin(idx+1) ;
idx = idx + 2 ;
case 'minpeakwidth'
if idx+1 > length(varargin) || ~( isscalar(varargin(idx+1)) && varargin(idx+1) >=0 ) then
error('findpeaks: MinPeakWidth must be a postive scalar') ;
end
minW = varargin(idx+1) ;
idx = idx + 2 ;
case 'maxpeakwidth'
if idx+1 > length(varargin) || ~( isscalar(varargin(idx+1)) && varargin(idx+1) >=0 ) then
error('findpeaks: MaxPeakWidth must be a postive scalar') ;
end
maxW = varargin(idx+1) ;
idx = idx + 2 ;
else
warning("findpeaks: Ignoring unknown option ") ;
idx = idx + 1 ;
end
end
endfunction
function y = lower (y)
// This function converts the input string to lower case
// Returns the input without modification if its not a string
if type(y) == 10 then
y = convstr(y, 'l') ;
else
y = y ;
end
endfunction
function out = bsminuseq(A)
// This function returns the difference between the elements of the input array
// This is only useful for the findpeaks function
A = A(:).' ;
atemp = [];
btemp=[] ;
for i=1:length(A)
atemp = [atemp ; A ]
btemp = [btemp A.' ]
end
out = atemp - btemp
endfunction
|
