]> cvs.zerfleddert.de Git - proxmark3-svn/blame - armsrc/lfsampling.c
Refactor parity functions
[proxmark3-svn] / armsrc / lfsampling.c
CommitLineData
31abe49f
MHS
1//-----------------------------------------------------------------------------
2// This code is licensed to you under the terms of the GNU GPL, version 2 or,
3// at your option, any later version. See the LICENSE.txt file for the text of
4// the license.
5//-----------------------------------------------------------------------------
6// Miscellaneous routines for low frequency sampling.
7//-----------------------------------------------------------------------------
8
9#include "proxmark3.h"
10#include "apps.h"
11#include "util.h"
12#include "string.h"
31abe49f 13#include "lfsampling.h"
e04475c4 14#include "usb_cdc.h" // for usb_poll_validate_length
15//#include "ticks.h" // for StartTicks
10a8875c 16
0cd2a41a 17sample_config config = { 1, 8, 1, 95, 0 } ;
31abe49f
MHS
18
19void printConfig()
20{
e2012d1b 21 Dbprintf("LF Sampling config: ");
31abe49f
MHS
22 Dbprintf(" [q] divisor: %d ", config.divisor);
23 Dbprintf(" [b] bps: %d ", config.bits_per_sample);
24 Dbprintf(" [d] decimation: %d ", config.decimation);
25 Dbprintf(" [a] averaging: %d ", config.averaging);
26 Dbprintf(" [t] trigger threshold: %d ", config.trigger_threshold);
27}
28
29
30/**
31 * Called from the USB-handler to set the sampling configuration
32 * The sampling config is used for std reading and snooping.
33 *
34 * Other functions may read samples and ignore the sampling config,
35 * such as functions to read the UID from a prox tag or similar.
36 *
37 * Values set to '0' implies no change (except for averaging)
38 * @brief setSamplingConfig
39 * @param sc
40 */
41void setSamplingConfig(sample_config *sc)
42{
43 if(sc->divisor != 0) config.divisor = sc->divisor;
44 if(sc->bits_per_sample!= 0) config.bits_per_sample= sc->bits_per_sample;
45 if(sc->decimation!= 0) config.decimation= sc->decimation;
46 if(sc->trigger_threshold != -1) config.trigger_threshold= sc->trigger_threshold;
47
48 config.averaging= sc->averaging;
49 if(config.bits_per_sample > 8) config.bits_per_sample = 8;
50 if(config.decimation < 1) config.decimation = 1;
51
52 printConfig();
53}
54
55sample_config* getSamplingConfig()
56{
57 return &config;
58}
10a8875c 59
31abe49f
MHS
60typedef struct {
61 uint8_t * buffer;
62 uint32_t numbits;
63 uint32_t position;
64} BitstreamOut;
65
31abe49f
MHS
66/**
67 * @brief Pushes bit onto the stream
68 * @param stream
69 * @param bit
70 */
10a8875c 71void pushBit( BitstreamOut* stream, uint8_t bit)
31abe49f
MHS
72{
73 int bytepos = stream->position >> 3; // divide by 8
74 int bitpos = stream->position & 7;
75 *(stream->buffer+bytepos) |= (bit > 0) << (7 - bitpos);
76 stream->position++;
77 stream->numbits++;
78}
10a8875c 79
31abe49f
MHS
80/**
81* Setup the FPGA to listen for samples. This method downloads the FPGA bitstream
82* if not already loaded, sets divisor and starts up the antenna.
83* @param divisor : 1, 88> 255 or negative ==> 134.8 KHz
84* 0 or 95 ==> 125 KHz
85*
86**/
87void LFSetupFPGAForADC(int divisor, bool lf_field)
88{
89 FpgaDownloadAndGo(FPGA_BITSTREAM_LF);
90 if ( (divisor == 1) || (divisor < 0) || (divisor > 255) )
91 FpgaSendCommand(FPGA_CMD_SET_DIVISOR, 88); //134.8Khz
92 else if (divisor == 0)
93 FpgaSendCommand(FPGA_CMD_SET_DIVISOR, 95); //125Khz
94 else
95 FpgaSendCommand(FPGA_CMD_SET_DIVISOR, divisor);
96
97 FpgaWriteConfWord(FPGA_MAJOR_MODE_LF_ADC | (lf_field ? FPGA_LF_ADC_READER_FIELD : 0));
98
99 // Connect the A/D to the peak-detected low-frequency path.
100 SetAdcMuxFor(GPIO_MUXSEL_LOPKD);
101 // Give it a bit of time for the resonant antenna to settle.
102 SpinDelay(50);
103 // Now set up the SSC to get the ADC samples that are now streaming at us.
104 FpgaSetupSsc();
105}
106
31abe49f
MHS
107/**
108 * Does the sample acquisition. If threshold is specified, the actual sampling
109 * is not commenced until the threshold has been reached.
110 * This method implements decimation and quantization in order to
111 * be able to provide longer sample traces.
112 * Uses the following global settings:
113 * @param decimation - how much should the signal be decimated. A decimation of N means we keep 1 in N samples, etc.
114 * @param bits_per_sample - bits per sample. Max 8, min 1 bit per sample.
115 * @param averaging If set to true, decimation will use averaging, so that if e.g. decimation is 3, the sample
116 * value that will be used is the average value of the three samples.
117 * @param trigger_threshold - a threshold. The sampling won't commence until this threshold has been reached. Set
118 * to -1 to ignore threshold.
119 * @param silent - is true, now outputs are made. If false, dbprints the status
120 * @return the number of bits occupied by the samples.
121 */
a37228c8 122uint32_t DoAcquisition(uint8_t decimation, uint32_t bits_per_sample, bool averaging, int trigger_threshold, bool silent, int bufsize)
31abe49f
MHS
123{
124 //.
0644d5e3 125 uint8_t *dest = BigBuf_get_addr();
a37228c8 126 bufsize = (bufsize > 0 && bufsize < BigBuf_max_traceLen()) ? bufsize : BigBuf_max_traceLen();
0644d5e3 127
3cec7061 128 //memset(dest, 0, bufsize); //creates issues with cmdread (marshmellow)
31abe49f
MHS
129
130 if(bits_per_sample < 1) bits_per_sample = 1;
131 if(bits_per_sample > 8) bits_per_sample = 8;
132
133 if(decimation < 1) decimation = 1;
134
135 // Use a bit stream to handle the output
136 BitstreamOut data = { dest , 0, 0};
137 int sample_counter = 0;
138 uint8_t sample = 0;
139 //If we want to do averaging
140 uint32_t sample_sum =0 ;
141 uint32_t sample_total_numbers =0 ;
142 uint32_t sample_total_saved =0 ;
143
e04475c4 144 while(!BUTTON_PRESS() && !usb_poll_validate_length() ) {
31abe49f
MHS
145 WDT_HIT();
146 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_TXRDY) {
147 AT91C_BASE_SSC->SSC_THR = 0x43;
148 LED_D_ON();
149 }
150 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_RXRDY) {
151 sample = (uint8_t)AT91C_BASE_SSC->SSC_RHR;
152 LED_D_OFF();
b29d55f2 153 // threshold either high or low values 128 = center 0. if trigger = 178
154 if ((trigger_threshold > 0) && (sample < (trigger_threshold+128)) && (sample > (128-trigger_threshold))) //
31abe49f 155 continue;
b29d55f2 156
31abe49f
MHS
157 trigger_threshold = 0;
158 sample_total_numbers++;
159
160 if(averaging)
161 {
162 sample_sum += sample;
163 }
164 //Check decimation
165 if(decimation > 1)
166 {
167 sample_counter++;
168 if(sample_counter < decimation) continue;
169 sample_counter = 0;
170 }
171 //Averaging
172 if(averaging && decimation > 1) {
173 sample = sample_sum / decimation;
174 sample_sum =0;
175 }
176 //Store the sample
177 sample_total_saved ++;
178 if(bits_per_sample == 8){
179 dest[sample_total_saved-1] = sample;
180 data.numbits = sample_total_saved << 3;//Get the return value correct
181 if(sample_total_saved >= bufsize) break;
182 }
183 else{
184 pushBit(&data, sample & 0x80);
185 if(bits_per_sample > 1) pushBit(&data, sample & 0x40);
186 if(bits_per_sample > 2) pushBit(&data, sample & 0x20);
187 if(bits_per_sample > 3) pushBit(&data, sample & 0x10);
188 if(bits_per_sample > 4) pushBit(&data, sample & 0x08);
189 if(bits_per_sample > 5) pushBit(&data, sample & 0x04);
190 if(bits_per_sample > 6) pushBit(&data, sample & 0x02);
191 //Not needed, 8bps is covered above
192 //if(bits_per_sample > 7) pushBit(&data, sample & 0x01);
193 if((data.numbits >> 3) +1 >= bufsize) break;
194 }
195 }
196 }
197
198 if(!silent)
199 {
200 Dbprintf("Done, saved %d out of %d seen samples at %d bits/sample",sample_total_saved, sample_total_numbers,bits_per_sample);
201 Dbprintf("buffer samples: %02x %02x %02x %02x %02x %02x %02x %02x ...",
202 dest[0], dest[1], dest[2], dest[3], dest[4], dest[5], dest[6], dest[7]);
203 }
204 return data.numbits;
205}
206/**
207 * @brief Does sample acquisition, ignoring the config values set in the sample_config.
208 * This method is typically used by tag-specific readers who just wants to read the samples
209 * the normal way
210 * @param trigger_threshold
211 * @param silent
212 * @return number of bits sampled
213 */
214uint32_t DoAcquisition_default(int trigger_threshold, bool silent)
215{
a37228c8 216 return DoAcquisition(1,8,0,trigger_threshold,silent,0);
31abe49f
MHS
217}
218uint32_t DoAcquisition_config( bool silent)
219{
220 return DoAcquisition(config.decimation
221 ,config.bits_per_sample
222 ,config.averaging
223 ,config.trigger_threshold
a37228c8 224 ,silent
225 ,0);
226}
227
228uint32_t DoPartialAcquisition(int trigger_threshold, bool silent, int sample_size) {
229 return DoAcquisition(1,8,0,trigger_threshold,silent,sample_size);
31abe49f
MHS
230}
231
1fbf8956 232uint32_t ReadLF(bool activeField, bool silent)
31abe49f 233{
1fbf8956 234 if (!silent) printConfig();
31abe49f
MHS
235 LFSetupFPGAForADC(config.divisor, activeField);
236 // Now call the acquisition routine
1fbf8956 237 return DoAcquisition_config(silent);
31abe49f
MHS
238}
239
240/**
241* Initializes the FPGA for reader-mode (field on), and acquires the samples.
242* @return number of bits sampled
243**/
1fbf8956 244uint32_t SampleLF(bool printCfg)
31abe49f 245{
e04475c4 246 uint32_t ret = ReadLF(true, printCfg);
247 FpgaWriteConfWord(FPGA_MAJOR_MODE_OFF);
248 return ret;
31abe49f
MHS
249}
250/**
251* Initializes the FPGA for snoop-mode (field off), and acquires the samples.
252* @return number of bits sampled
253**/
254
255uint32_t SnoopLF()
256{
e04475c4 257 uint32_t ret = ReadLF(false, true);
258 FpgaWriteConfWord(FPGA_MAJOR_MODE_OFF);
259 return ret;
31abe49f 260}
7cfc777b 261
262/**
263* acquisition of T55x7 LF signal. Similart to other LF, but adjusted with @marshmellows thresholds
264* the data is collected in BigBuf.
265**/
66837a03 266void doT55x7Acquisition(size_t sample_size) {
7cfc777b 267
6fe5c94b 268 #define T55xx_READ_UPPER_THRESHOLD 128+60 // 60 grph
269 #define T55xx_READ_LOWER_THRESHOLD 128-60 // -60 grph
7cfc777b 270 #define T55xx_READ_TOL 5
271
272 uint8_t *dest = BigBuf_get_addr();
273 uint16_t bufsize = BigBuf_max_traceLen();
274
66837a03 275 if ( bufsize > sample_size )
276 bufsize = sample_size;
7cfc777b 277
7cfc777b 278 uint16_t i = 0;
279 bool startFound = false;
280 bool highFound = false;
6fe5c94b 281 bool lowFound = false;
7cfc777b 282 uint8_t curSample = 0;
db829602 283 uint8_t lastSample = 0;
be2d41b7 284 uint16_t skipCnt = 0;
506672c4 285 while(!BUTTON_PRESS() && !usb_poll_validate_length() && skipCnt<1000 && i<bufsize ) {
be2d41b7 286 WDT_HIT();
7cfc777b 287 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_TXRDY) {
288 AT91C_BASE_SSC->SSC_THR = 0x43;
289 LED_D_ON();
290 }
291 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_RXRDY) {
292 curSample = (uint8_t)AT91C_BASE_SSC->SSC_RHR;
be2d41b7 293 LED_D_OFF();
294
9f669cb2 295 // skip until the first high sample above threshold
7cfc777b 296 if (!startFound && curSample > T55xx_READ_UPPER_THRESHOLD) {
6fe5c94b 297 //if (curSample > lastSample)
298 // lastSample = curSample;
7cfc777b 299 highFound = true;
300 } else if (!highFound) {
be2d41b7 301 skipCnt++;
7cfc777b 302 continue;
303 }
6fe5c94b 304 // skip until the first Low sample below threshold
305 if (!startFound && curSample < T55xx_READ_LOWER_THRESHOLD) {
306 //if (curSample > lastSample)
307 lastSample = curSample;
308 lowFound = true;
309 } else if (!lowFound) {
310 skipCnt++;
311 continue;
312 }
313
314
9f669cb2 315 // skip until first high samples begin to change
6fe5c94b 316 if (startFound || curSample > T55xx_READ_LOWER_THRESHOLD+T55xx_READ_TOL){
9f669cb2 317 // if just found start - recover last sample
318 if (!startFound) {
db829602 319 dest[i++] = lastSample;
9f669cb2 320 startFound = true;
321 }
322 // collect samples
7cfc777b 323 dest[i++] = curSample;
7cfc777b 324 }
325 }
326 }
327}
e04475c4 328
329/**
330* acquisition of Cotag LF signal. Similart to other LF, since the Cotag has such long datarate RF/384
331* and is Manchester?, we directly gather the manchester data into bigbuff
332**/
333#define COTAG_T1 384
334#define COTAG_T2 (COTAG_T1>>1)
335#define COTAG_ONE_THRESHOLD 128+30
336#define COTAG_ZERO_THRESHOLD 128-30
337#ifndef COTAG_BITS
338#define COTAG_BITS 264
339#endif
340void doCotagAcquisition(size_t sample_size) {
341
342 uint8_t *dest = BigBuf_get_addr();
343 uint16_t bufsize = BigBuf_max_traceLen();
344
345 if ( bufsize > sample_size )
346 bufsize = sample_size;
347
348 dest[0] = 0;
349 uint8_t sample = 0, firsthigh = 0, firstlow = 0;
350 uint16_t i = 0;
351
352 while (!BUTTON_PRESS() && !usb_poll_validate_length() && (i < bufsize) ) {
353 WDT_HIT();
354 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_TXRDY) {
355 AT91C_BASE_SSC->SSC_THR = 0x43;
356 LED_D_ON();
357 }
358
359 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_RXRDY) {
360 sample = (uint8_t)AT91C_BASE_SSC->SSC_RHR;
361 LED_D_OFF();
362
363 // find first peak
364 if ( !firsthigh ) {
365 if (sample < COTAG_ONE_THRESHOLD)
366 continue;
367 firsthigh = 1;
368 }
369 if ( !firstlow ){
370 if (sample > COTAG_ZERO_THRESHOLD )
371 continue;
372 firstlow = 1;
373 }
374
375 ++i;
376
377 if ( sample > COTAG_ONE_THRESHOLD)
378 dest[i] = 255;
379 else if ( sample < COTAG_ZERO_THRESHOLD)
380 dest[i] = 0;
381 else
382 dest[i] = dest[i-1];
383 }
384 }
385}
386
387uint32_t doCotagAcquisitionManchester() {
388
389 uint8_t *dest = BigBuf_get_addr();
390 uint16_t bufsize = BigBuf_max_traceLen();
391
392 if ( bufsize > COTAG_BITS )
393 bufsize = COTAG_BITS;
394
395 dest[0] = 0;
396 uint8_t sample = 0, firsthigh = 0, firstlow = 0;
397 uint16_t sample_counter = 0, period = 0;
398 uint8_t curr = 0, prev = 0;
399
400 while (!BUTTON_PRESS() && !usb_poll_validate_length() && (sample_counter < bufsize) ) {
401 WDT_HIT();
402 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_TXRDY) {
403 AT91C_BASE_SSC->SSC_THR = 0x43;
404 LED_D_ON();
405 }
406
407 if (AT91C_BASE_SSC->SSC_SR & AT91C_SSC_RXRDY) {
408 sample = (uint8_t)AT91C_BASE_SSC->SSC_RHR;
409 LED_D_OFF();
410
411 // find first peak
412 if ( !firsthigh ) {
413 if (sample < COTAG_ONE_THRESHOLD)
414 continue;
415 firsthigh = 1;
416 }
417
418 if ( !firstlow ){
419 if (sample > COTAG_ZERO_THRESHOLD )
420 continue;
421 firstlow = 1;
422 }
423
424 // set sample 255, 0, or previous
425 if ( sample > COTAG_ONE_THRESHOLD){
426 prev = curr;
427 curr = 1;
428 }
429 else if ( sample < COTAG_ZERO_THRESHOLD) {
430 prev = curr;
431 curr = 0;
432 }
433 else {
434 curr = prev;
435 }
436
437 // full T1 periods,
438 if ( period > 0 ) {
439 --period;
440 continue;
441 }
442
443 dest[sample_counter] = curr;
444 ++sample_counter;
445 period = COTAG_T1;
446 }
447 }
448 return sample_counter;
449}
Impressum, Datenschutz