summaryrefslogtreecommitdiffstats
path: root/Project/TransferFunction.h
blob: e51de5aceb2ecf949324997174eb24d9f061445b (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
/*
 *  Copyright (C) 2017  Thales Lima Oliveira <thales@ufu.br>
 *
 *  This program is free software; you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation; either version 2 of the License, or
 *  any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <https://www.gnu.org/licenses/>.
 */

#ifndef TRANSFERFUNCTION_H
#define TRANSFERFUNCTION_H

#include "ControlElement.h"

#include <wx/dcscreen.h>
#include "OpenGLText.h"

class TransferFunctionForm;

/**
 * @class TransferFunction
 * @author Thales Lima Oliveira <thales@ufu.br>
 * @date 06/10/2017
 * @brief Calculates the time response by a frequency domain transfer function.
 * @file TransferFunction.h
 */
class TransferFunction : public ControlElement
{
   public:
    struct SpaceState {
        std::vector<std::vector<double> > A;
        std::vector<double> B;
        std::vector<double> C;
        double D;
    };

    TransferFunction(int id);
    ~TransferFunction();

    virtual void Draw(wxPoint2DDouble translation, double scale) const;
    virtual void DrawDC(wxPoint2DDouble translation, double scale, wxGraphicsContext* gc) const;
    virtual bool Contains(wxPoint2DDouble position) const { return m_rect.Contains(position); }
    virtual bool Intersects(wxRect2DDouble rect) const { return m_rect.Intersects(rect); }
    virtual bool ShowForm(wxWindow* parent, Element* element);
    virtual void Rotate(bool clockwise = true);

    virtual std::vector<double> GetNumerator() const { return m_numerator; }
    virtual std::vector<double> GetDenominator() const { return m_denominator; }
    virtual void SetNumerator(std::vector<double> numerator) { m_numerator = numerator; }
    virtual void SetDenominator(std::vector<double> denominator) { m_denominator = denominator; }
    virtual void UpdateTFText();
    virtual bool UpdateText();
    virtual SpaceState GetSpaceState() { return m_ss; }
    /**
     * @brief Convert the transfer function to space state on controllable canonical form (CCF).
     * @param maxIteration Max number of solution iteration.
     * @param error Tolerance for calculation.
     */
    virtual void CalculateSpaceState(int maxIteration = 100, double error = 1e-3);
    /**
     * @brief Calculates the time response by the space state form of transfer function.
     *
     * Uses the implicit trapezoidal rule to solve:
     * <center>\f$ y_{n+1}=y_n+\frac{1}{2}h(y'_n+y'_{n+1}) \f$</center>
     * @param input Input value.
     * @param timeStep Time step.
     * @return true if the calculation converges, false otherwise.
     */
    virtual bool Solve(double* input, double timeStep);
    
    virtual rapidxml::xml_node<>* SaveElement(rapidxml::xml_document<>& doc, rapidxml::xml_node<>* elementListNode);
    virtual bool OpenElement(rapidxml::xml_node<>* elementNode);

    virtual Element* GetCopy();

   protected:
    virtual void SetText(wxString numerator, wxString denominator);
    virtual wxString GetSuperscriptNumber(int number);
    virtual void GetTFString(wxString& numerator, wxString& denominator);

    wchar_t m_supNumber[10];

    OpenGLText* m_glTextNum = nullptr;
    OpenGLText* m_glTextDen = nullptr;
    int m_fontSize = 10;

    std::vector<double> m_numerator;
    std::vector<double> m_denominator;
    SpaceState m_ss;

    std::vector<double> m_x;
    std::vector<double> m_dx;
    double m_error = 1e-3;
    int m_maxIteration = 100;
};

#endif  // TRANSFERFUNCTION_H