| 25 | |
| 26 | template <typename T> |
| 27 | void square_wave2() |
| 28 | { |
| 29 | typedef exprtk::symbol_table<T> symbol_table_t; |
| 30 | typedef exprtk::expression<T> expression_t; |
| 31 | typedef exprtk::parser<T> parser_t; |
| 32 | |
| 33 | const std::string wave_program = |
| 34 | " var r := 0; " |
| 35 | " " |
| 36 | " for (var i := 0; i < 1000; i += 1) " |
| 37 | " { " |
| 38 | " r += (1 / (2i + 1)) * sin((4i + 2) * pi * f * t); " |
| 39 | " }; " |
| 40 | " " |
| 41 | " r *= a * (4 / pi); "; |
| 42 | |
| 43 | static const T pi = T(3.141592653589793238462643383279502); |
| 44 | |
| 45 | T f = pi / T(10); |
| 46 | T t = T(0); |
| 47 | T a = T(10); |
| 48 | |
| 49 | symbol_table_t symbol_table; |
| 50 | symbol_table.add_variable("f",f); |
| 51 | symbol_table.add_variable("t",t); |
| 52 | symbol_table.add_variable("a",a); |
| 53 | symbol_table.add_constants(); |
| 54 | |
| 55 | expression_t expression; |
| 56 | expression.register_symbol_table(symbol_table); |
| 57 | |
| 58 | parser_t parser; |
| 59 | parser.compile(wave_program,expression); |
| 60 | |
| 61 | const T delta = (T(4) * pi) / T(1000); |
| 62 | |
| 63 | for (t = (T(-2) * pi); t <= (T(+2) * pi); t += delta) |
| 64 | { |
| 65 | const T result = expression.value(); |
| 66 | printf("%19.15f\t%19.15f\n", t, result); |
| 67 | } |
| 68 | } |
| 69 | |
| 70 | int main() |
| 71 | { |