]>
Commit | Line | Data |
---|---|---|
1 | //----------------------------------------------------------------------------- | |
2 | // Copyright (C) 2009 Michael Gernoth <michael at gernoth.net> | |
3 | // Copyright (C) 2010 iZsh <izsh at fail0verflow.com> | |
4 | // | |
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 | |
7 | // the license. | |
8 | //----------------------------------------------------------------------------- | |
9 | // UI utilities | |
10 | //----------------------------------------------------------------------------- | |
11 | ||
12 | #include "ui.h" | |
13 | double CursorScaleFactor; | |
14 | int PlotGridX, PlotGridY, PlotGridXdefault= 64, PlotGridYdefault= 64; | |
15 | int offline; | |
16 | int flushAfterWrite = 0; | |
17 | extern pthread_mutex_t print_lock; | |
18 | ||
19 | static char *logfilename = "proxmark3.log"; | |
20 | ||
21 | void PrintAndLog(char *fmt, ...) | |
22 | { | |
23 | char *saved_line; | |
24 | int saved_point; | |
25 | va_list argptr, argptr2; | |
26 | static FILE *logfile = NULL; | |
27 | static int logging = 1; | |
28 | ||
29 | // lock this section to avoid interlacing prints from different threats | |
30 | pthread_mutex_lock(&print_lock); | |
31 | ||
32 | if (logging && !logfile) { | |
33 | logfile = fopen(logfilename, "a"); | |
34 | if (!logfile) { | |
35 | fprintf(stderr, "Can't open logfile, logging disabled!\n"); | |
36 | logging=0; | |
37 | } | |
38 | } | |
39 | ||
40 | int need_hack = (rl_readline_state & RL_STATE_READCMD) > 0; | |
41 | ||
42 | if (need_hack) { | |
43 | saved_point = rl_point; | |
44 | saved_line = rl_copy_text(0, rl_end); | |
45 | rl_save_prompt(); | |
46 | rl_replace_line("", 0); | |
47 | rl_redisplay(); | |
48 | } | |
49 | ||
50 | va_start(argptr, fmt); | |
51 | va_copy(argptr2, argptr); | |
52 | vprintf(fmt, argptr); | |
53 | printf(" "); // cleaning prompt | |
54 | va_end(argptr); | |
55 | printf("\n"); | |
56 | ||
57 | if (need_hack) { | |
58 | rl_restore_prompt(); | |
59 | rl_replace_line(saved_line, 0); | |
60 | rl_point = saved_point; | |
61 | rl_redisplay(); | |
62 | free(saved_line); | |
63 | } | |
64 | ||
65 | if (logging && logfile) { | |
66 | vfprintf(logfile, fmt, argptr2); | |
67 | fprintf(logfile,"\n"); | |
68 | fflush(logfile); | |
69 | } | |
70 | va_end(argptr2); | |
71 | ||
72 | if (flushAfterWrite == 1) { | |
73 | fflush(NULL); | |
74 | } | |
75 | //release lock | |
76 | pthread_mutex_unlock(&print_lock); | |
77 | } | |
78 | ||
79 | void SetLogFilename(char *fn) { | |
80 | logfilename = fn; | |
81 | } | |
82 | ||
83 | void iceIIR_Butterworth(int *data, const size_t len){ | |
84 | ||
85 | int i,j; | |
86 | ||
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; | |
91 | ||
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}; | |
96 | ||
97 | float sample = 0; // input sample read from array | |
98 | float complex x_prime = 1.0f; // save sample for estimating frequency | |
99 | float complex x; | |
100 | ||
101 | for (i=0; i<adjustedLen; ++i) { | |
102 | ||
103 | sample = data[i]; | |
104 | ||
105 | // remove DC offset and mix to complex baseband | |
106 | x = (sample - 127.5f) * cexpf( _Complex_I * 2 * M_PI * fc * i ); | |
107 | ||
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] + | |
113 | b[1]*iir_buf[1] + | |
114 | b[2]*iir_buf[2]; | |
115 | ||
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 | |
120 | ||
121 | output[i] =(freq > 0)? 10 : -10; | |
122 | } | |
123 | ||
124 | // show data | |
125 | for (j=0; j<adjustedLen; ++j) | |
126 | data[j] = output[j]; | |
127 | ||
128 | free(output); | |
129 | } | |
130 | ||
131 | float complex cexpf (float complex Z) | |
132 | { | |
133 | float complex Res; | |
134 | double rho = exp (__real__ Z); | |
135 | __real__ Res = rho * cosf(__imag__ Z); | |
136 | __imag__ Res = rho * sinf(__imag__ Z); | |
137 | return Res; | |
138 | } |