1 //-----------------------------------------------------------------------------
2 // Copyright (C) 2009 Michael Gernoth <michael at gernoth.net>
3 // Copyright (C) 2010 iZsh <izsh at fail0verflow.com>
5 // This code is licensed to you under the terms of the GNU GPL, version 2 or,
6 // at your option, any later version. See the LICENSE.txt file for the text of
8 //-----------------------------------------------------------------------------
10 //-----------------------------------------------------------------------------
13 double CursorScaleFactor
;
14 int PlotGridX
, PlotGridY
, PlotGridXdefault
= 64, PlotGridYdefault
= 64;
16 int flushAfterWrite
= 0;
17 extern pthread_mutex_t print_lock
;
19 static char *logfilename
= "proxmark3.log";
21 void PrintAndLog(char *fmt
, ...)
25 va_list argptr
, argptr2
;
26 static FILE *logfile
= NULL
;
27 static int logging
= 1;
29 // lock this section to avoid interlacing prints from different threats
30 pthread_mutex_lock(&print_lock
);
32 if (logging
&& !logfile
) {
33 logfile
= fopen(logfilename
, "a");
35 fprintf(stderr
, "Can't open logfile, logging disabled!\n");
40 int need_hack
= (rl_readline_state
& RL_STATE_READCMD
) > 0;
43 saved_point
= rl_point
;
44 saved_line
= rl_copy_text(0, rl_end
);
46 rl_replace_line("", 0);
50 va_start(argptr
, fmt
);
51 va_copy(argptr2
, argptr
);
53 printf(" "); // cleaning prompt
59 rl_replace_line(saved_line
, 0);
60 rl_point
= saved_point
;
65 if (logging
&& logfile
) {
66 vfprintf(logfile
, fmt
, argptr2
);
67 fprintf(logfile
,"\n");
72 if (flushAfterWrite
== 1) {
76 pthread_mutex_unlock(&print_lock
);
79 void SetLogFilename(char *fn
) {
83 void iceFsk3(int * data
, const size_t len
){
87 int * output
= (int* ) malloc(sizeof(int) * len
);
88 memset(output
, 0x00, len
);
89 float fc
= 0.1125f
; // center frequency
90 size_t adjustedLen
= len
;
92 // create very simple low-pass filter to remove images (2nd-order Butterworth)
93 float complex iir_buf
[3] = {0,0,0};
94 float b
[3] = {0.003621681514929, 0.007243363029857, 0.003621681514929};
95 float a
[3] = {1.000000000000000, -1.822694925196308, 0.837181651256023};
97 float sample
= 0; // input sample read from file
98 float complex x_prime
= 1.0f
; // save sample for estimating frequency
101 for (i
=0; i
<adjustedLen
; ++i
) {
103 sample
= data
[i
]+128;
105 // remove DC offset and mix to complex baseband
106 x
= (sample
- 127.5f
) * cexpf( _Complex_I
* 2 * M_PI
* fc
* i
);
108 // apply low-pass filter, removing spectral image (IIR using direct-form II)
109 iir_buf
[2] = iir_buf
[1];
110 iir_buf
[1] = iir_buf
[0];
111 iir_buf
[0] = x
- a
[1]*iir_buf
[1] - a
[2]*iir_buf
[2];
112 x
= b
[0]*iir_buf
[0] +
116 // compute instantaneous frequency by looking at phase difference
117 // between adjacent samples
118 float freq
= cargf(x
*conjf(x_prime
));
119 x_prime
= x
; // retain this sample for next iteration
121 output
[i
] =(freq
> 0)? 10 : -10;
125 for (j
=0; j
<adjustedLen
; ++j
)
132 for (j
=0; j
<adjustedLen
; ++j
){
133 if ( data
[j
] == 10) break;
137 for (;j
<adjustedLen
; ++j
){
138 if ( data
[j
] == -10 ) break;
142 int fieldlen
= stopOne
-startOne
;
144 fieldlen
= (fieldlen
== 39 || fieldlen
== 41)? 40 : fieldlen
;
145 fieldlen
= (fieldlen
== 59 || fieldlen
== 51)? 50 : fieldlen
;
146 if ( fieldlen
!= 40 && fieldlen
!= 50){
147 printf("Detected field Length: %d \n", fieldlen
);
148 printf("Can only handle 40 or 50. Aborting...\n");
153 // FSK sequence start == 000111
155 for (i
=0; i
<adjustedLen
; ++i
){
157 for ( j
= 0; j
< 6*fieldlen
; ++j
){
166 printf("000111 position: %d \n", startPos
);
168 startPos
+= 6*fieldlen
+5;
173 for (i
=startPos
; i
< adjustedLen
; i
+= 40){
174 bit
= data
[i
]>0 ? 1:0;
180 for (i
=startPos
; i
< adjustedLen
; i
+= 50){
181 bit
= data
[i
]>0 ? 1:0;
182 printf("%d", bit
); }
188 float complex cexpf (float complex Z
)
191 double rho
= exp (__real__ Z
);
192 __real__ Res
= rho
* cosf(__imag__ Z
);
193 __imag__ Res
= rho
* sinf(__imag__ Z
);