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
83
84
85
86
87
88
89
90
91
92
93
94
95
|
function [Y] = peak2peak(X, dim)
funcprot(0);
narginchk(1,2,argn(2));
// peak2peak finds the maximum-to-minimum difference
// Calling Sequence
// Y=peak2peak(X)
// Y=peak2peak(X,dim)
// Parameters
// X: Real or complex valued input vector or matrix.
// dim: Dimension to compute maximum-to-minimum difference (Row/Column)
// Description
// Y= peak2peak(X) returns the difference between maximum and minimum values in X, operating along the first non-singleton
// dimension of X by default.
// Y= peak2peak(X,dim) returns the maximum-to-minimum difference in X, along dimension dim.
// Examples
// 1. t = 0:0.001:1-0.001;
// x = cos(2*pi*100*t);
// y = peak2peak(x)
// 2. t = 0:0.001:1-0.001;
// x = (1:4)'*cos(2*pi*100*t);
// y = peak2peak(x,2)
// Authors
// Rahul Dalmia
//function only accepts real values due to limitations of the 'max' function in Scilab
//Modifications made by Debdeep Dey
[nr, nc] = size (X); // Dimensions of Input calculated
if(type(X)==10) then //if i/p is a char type vector or matrix or a string
w=X;
[nr,nc]=size(X);
if(nr==1 & nc==1) then
X=ascii(X);
X=matrix(X,length(w));
X=X';
else
X=ascii(X);
X=matrix(X,size(w));
end
end
if (~exists('dim','local')) then
if (nr==1) then
Y = zeros(nr, 1); // preset all output fields to 0
for i= 1:nr
maxim=max(X(i,:)); // maximum and minimum values are found
minim=min(X(i,:));
Y(i,1)=maxim-minim; // Peak to peak value is calculated from the difference of max and min
end
else
Y = zeros(1, nc); // preset all output fields to 0
for i= 1:nc
maxim=max(X(:,i)); // maximum and minimum values are found
minim=min(X(:,i));
Y(1,i)=maxim-minim; // Peak to peak value is calculated from the difference of max and min
end
end
elseif (exists('dim','local')) then
if (dim<1) then
error("Dimension argument must be a positive integer scalar within indexing range.");
end
if (dim==1) then
Y = zeros(1, nc); // preset all output fields to 0
for i= 1:nc
maxim=max(X(:,i)); // maximum and minimum values are found
minim=min(X(:,i));
Y(1,i)=maxim-minim; // Peak to peak value is calculated from the difference of max and min
end
elseif (dim==2) then
Y = zeros(nr, 1); // preset all output fields to 0
for i= 1:nr
maxim=max(X(i,:)); // maximum and minimum values are found
minim=min(X(i,:));
Y(i,1)=maxim-minim; // Peak to peak value is calculated from the difference of max and min
end
else// for cases when dim >2
Y=zeros(nr,nc);
end
end
endfunction
function narginchk(l,h,ni)
if(ni<l) then //ensure that the no. of input args is either 1 or 2
error("Not enough input arguments");
end
if (ni>h) then
error("Too many input arguments");
end
endfunction
|
