|
7017
|
1 ## Copyright (C) 2006, 2007 John W. Eaton |
|
7016
|
2 ## |
|
|
3 ## This file is part of Octave. |
|
|
4 ## |
|
|
5 ## Octave is free software; you can redistribute it and/or modify it |
|
|
6 ## under the terms of the GNU General Public License as published by |
|
|
7 ## the Free Software Foundation; either version 3 of the License, or (at |
|
|
8 ## your option) any later version. |
|
|
9 ## |
|
|
10 ## Octave is distributed in the hope that it will be useful, but |
|
|
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
|
|
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
|
|
13 ## General Public License for more details. |
|
|
14 ## |
|
|
15 ## You should have received a copy of the GNU General Public License |
|
|
16 ## along with Octave; see the file COPYING. If not, see |
|
|
17 ## <http://www.gnu.org/licenses/>. |
|
|
18 |
|
5590
|
19 %% Automatically generated from DejaGNU files |
|
|
20 |
|
|
21 %% test/octave.test/signal/detrend-1.m |
|
|
22 %% detrend-1.m |
|
|
23 %% |
|
|
24 %% Author: David Billinghurst (David.Billinghurst@riotinto.com.au) |
|
|
25 %% Comalco Research and Technology |
|
|
26 %% 02 May 2000 |
|
|
27 %!test |
|
|
28 %! N=32; |
|
|
29 %! x = (0:1:N-1)/N + 2; |
|
|
30 %! y = detrend(x); |
|
|
31 %! assert(all (all (abs (y) < 10*eps))); |
|
|
32 |
|
|
33 %% test/octave.test/signal/detrend-2.m |
|
|
34 %% detrend-2.m |
|
|
35 %% |
|
|
36 %% Author: David Billinghurst (David.Billinghurst@riotinto.com.au) |
|
|
37 %% Comalco Research and Technology |
|
|
38 %% 02 May 2000 |
|
|
39 %!test |
|
|
40 %! N=32; |
|
|
41 %! t = (0:1:N-1)/N; |
|
|
42 %! x = t .* t + 2; |
|
|
43 %! y = detrend(x,2); |
|
|
44 %! assert(all (all (abs (y) < 30*eps))); |
|
|
45 |
|
|
46 %% test/octave.test/signal/detrend-3.m |
|
|
47 %% detrend-3.m |
|
|
48 %% |
|
|
49 %% Author: David Billinghurst (David.Billinghurst@riotinto.com.au) |
|
|
50 %% Comalco Research and Technology |
|
|
51 %% 02 May 2000 |
|
|
52 %!test |
|
|
53 %! N=32; |
|
|
54 %! t = (0:1:N-1)/N; |
|
|
55 %! x = [t;4*t-3]'; |
|
|
56 %! y = detrend(x); |
|
|
57 %! assert(all (all (abs (y) < 20*eps))); |
|
|
58 |
|
|
59 %% test/octave.test/signal/fft-1.m |
|
|
60 %% fft-1.m |
|
|
61 %% |
|
|
62 %% Author: David Billinghurst (David.Billinghurst@riotinto.com.au) |
|
|
63 %% Comalco Research and Technology |
|
|
64 %% 02 May 2000 |
|
|
65 %!test |
|
|
66 %! N=64; |
|
|
67 %! n=4; |
|
|
68 %! t = 2*pi*(0:1:N-1)/N; |
|
|
69 %! s = cos(n*t); |
|
|
70 %! S = fft(s); |
|
|
71 %! |
|
|
72 %! answer = 0*t; |
|
|
73 %! answer(n+1) = N/2; |
|
|
74 %! answer(N-n+1) = N/2; |
|
|
75 %! |
|
|
76 %! assert(all( abs(S-answer) < 4*N*eps )); |
|
|
77 |
|
|
78 %% test/octave.test/signal/ifft-1.m |
|
|
79 %% ifft-1.m |
|
|
80 %% |
|
|
81 %% Author: David Billinghurst (David.Billinghurst@riotinto.com.au) |
|
|
82 %% Comalco Research and Technology |
|
|
83 %% 02 May 2000 |
|
|
84 %!test |
|
|
85 %! N=64; |
|
|
86 %! n=7; |
|
|
87 %! t = 2*pi*(0:1:N-1)/N; |
|
|
88 %! s = cos(n*t); |
|
|
89 %! |
|
|
90 %! S = 0*t; |
|
|
91 %! S(n+1) = N/2; |
|
|
92 %! S(N-n+1) = N/2; |
|
|
93 %! |
|
|
94 %! assert(all( abs(ifft(S)-s) < 4*N*eps )); |
|
|
95 |
|
|
96 %% test/octave.test/signal/fft2-1.m |
|
|
97 %% fft2-1.m |
|
|
98 %% |
|
|
99 %% Author: David Billinghurst (David.Billinghurst@riotinto.com.au) |
|
|
100 %% Comalco Research and Technology |
|
|
101 %% 02 May 2000 |
|
|
102 %!test |
|
|
103 %! M=16; |
|
|
104 %! N=8; |
|
|
105 %! |
|
|
106 %! m=5; |
|
|
107 %! n=3; |
|
|
108 %! |
|
|
109 %! x = 2*pi*(0:1:M-1)/M; |
|
|
110 %! y = 2*pi*(0:1:N-1)/N; |
|
|
111 %! sx = cos(m*x); |
|
|
112 %! sy = sin(n*y); |
|
|
113 %! s=kron(sx',sy); |
|
|
114 %! S = fft2(s); |
|
|
115 %! answer = kron(fft(sx)',fft(sy)); |
|
|
116 %! assert(all( all( abs(S-answer) < 4*M*N*eps ) )); |
|
|
117 |
|
|
118 %% test/octave.test/signal/ifft2-1.m |
|
|
119 %% ifft2-1.m |
|
|
120 %% |
|
|
121 %% Author: David Billinghurst (David.Billinghurst@riotinto.com.au) |
|
|
122 %% Comalco Research and Technology |
|
|
123 %% 02 May 2000 |
|
|
124 %!test |
|
|
125 %! M=12; |
|
|
126 %! N=7; |
|
|
127 %! |
|
|
128 %! m=3; |
|
|
129 %! n=2; |
|
|
130 %! |
|
|
131 %! x = 2*pi*(0:1:M-1)/M; |
|
|
132 %! y = 2*pi*(0:1:N-1)/N; |
|
|
133 %! |
|
|
134 %! sx = cos(m*x); |
|
|
135 %! sy = cos(n*y); |
|
|
136 %! |
|
|
137 %! S = kron(fft(sx)',fft(sy)); |
|
|
138 %! answer=kron(sx',sy); |
|
|
139 %! s = ifft2(S); |
|
|
140 %! |
|
|
141 %! assert(all( all( abs(s-answer) < 30*eps ) )); |
|
|
142 |
|
|
143 %% test/octave.test/signal/unwrap-1.m |
|
6699
|
144 %!function t = xassert(a,b,tol) |
|
5590
|
145 %! if (nargin == 1) |
|
|
146 %! t = all(a(:)); |
|
|
147 %! else |
|
|
148 %! if (nargin == 2) |
|
|
149 %! tol = 0; |
|
|
150 %! endif |
|
|
151 %! if (any (size(a) != size(b))) |
|
|
152 %! t = 0; |
|
|
153 %! elseif (any (abs(a(:) - b(:)) > tol)) |
|
|
154 %! t = 0; |
|
|
155 %! else |
|
|
156 %! t = 1; |
|
|
157 %! endif |
|
|
158 %! endif |
|
|
159 %! |
|
|
160 %!test |
|
|
161 %! |
|
|
162 %! i = 0; |
|
|
163 %! t = []; |
|
|
164 %! |
|
|
165 %! r = [0:100]; # original vector |
|
|
166 %! w = r - 2*pi*floor((r+pi)/(2*pi)); # wrapped into [-pi,pi] |
|
|
167 %! tol = 1e3*eps; # maximum expected deviation |
|
|
168 %! |
|
6699
|
169 %! t(++i) = xassert(r, unwrap(w), tol); #unwrap single row |
|
|
170 %! t(++i) = xassert(r', unwrap(w'), tol); #unwrap single column |
|
|
171 %! t(++i) = xassert([r',r'], unwrap([w',w']), tol); #unwrap 2 columns |
|
|
172 %! t(++i) = xassert([r;r], unwrap([w;w],[],2), tol); #verify that dim works |
|
|
173 %! t(++i) = xassert(r+10, unwrap(10+w), tol); #verify that r(1)>pi works |
|
5590
|
174 %! |
|
6699
|
175 %! t(++i) = xassert(w', unwrap(w',[],2)); #unwrap col by rows should not change it |
|
|
176 %! t(++i) = xassert(w, unwrap(w,[],1)); #unwrap row by cols should not change it |
|
|
177 %! t(++i) = xassert([w;w], unwrap([w;w])); #unwrap 2 rows by cols should not change them |
|
5590
|
178 %! |
|
|
179 %! ## verify that setting tolerance too low will cause bad results. |
|
6699
|
180 %! t(++i) = xassert(any(abs(r - unwrap(w,0.8)) > 100)); |
|
5590
|
181 %! |
|
|
182 %! assert(all(t)); |
|
|
183 |
