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
|
function [psi,x] = morlet(lb, ub, n)
//Compute the Morlet wavelet.
//Calling sequence:
//[psi,x]= morlet(lb,ub,n)
//Parameters:
//lb: Real or complex valued vector/scalar
//ub: Real or complex valued vector/scalar
//n: Real positive scalar number
//Description:
//This function returns values of the Morlet wavelet in the specified interval for all the sample points.
//Example:
//[a,b] = morlet([1 2 3], [1 2 3], 1)
//a = [0.1720498; -0.1135560; -0.0084394]
//b = [1; 2; 3]
funcprot(0);
rhs = argn(2);
if (rhs ~= 3) then
error("morlet: wrong number of input arguments");
end
if (length(lb) ~= length(ub)) then
error("morlet: arg1 and arg2 msut have same dimension");
end
if (~isreal(n) | n <= 0) then
error("morlet: n must be a strictly positive real number");
end
x = [];
for i=1:length(lb)
x(i, :) = linspace(lb(i), ub(i), n);
end
psi = cos(5.*x) .* exp(-x.^2/2);
endfunction
//input validation:
//assert_checkerror("morlet()", "morlet: wrong number of input arguments");
//assert_checkerror("morlet(1, 2)", "morlet: wrong number of input arguments");
//assert_checkerror("morlet(1, 2, 3, 4)", "Wrong number of input arguments.");
//assert_checkerror("morlet(1, 2, -1)", "morlet: n must be a strictly positive real number");
//assert_checkerror("morlet(1, 2, 2+3*%i)", "morlet: n must be a strictly positive real number");
//assert_checkerror(" morlet([5, 2, 7], [1, 3], 3)", "morlet: arg1 and arg2 msut have same dimension");
//test basic input:
//[a, b] = morlet(1, 2, 3);
//assert_checkalmostequal(a, [0.17205, 0.11254, -0.11356], 5e-5);
//assert_checkalmostequal(b, [1, 1.5, 2], %eps);
//test complex input:
//[a, b] = morlet(3+2*%i, 4+8*%i, 2);
//assert_checkalmostequal(a, [-495.15886-756.35443*%i, -5.08135E26-3.07588E27*%i], 5e-5);
//assert_checkalmostequal(b, [3+2*%i, 4+8*%i], %eps);
//test real vector:
//[a, b] = morlet([1, 2], [3, 4], 2);
//A = [0.1720498, -0.0084394; -0.113556, 0.0001369];
//B = [1, 3; 2, 4];
//assert_checkalmostequal(A, a, 5e-5);
//assert_checkalmostequal(B, b, %eps);
//test complex vector:
//[a, b] = morlet([1 + 2*%i, 3*%i], [2*%i, 2 + 3*%i], 1);
//A = [81377.39587; -19069291.16508 + 5732843.75676*%i];
//B = [2*%i; 2+3*%i];
//assert_checkalmostequal(a, A, 5e-5);
//assert_checkalmostequal(b, B, %eps);
|
