/* * Copyright (C) 2017 Thales Lima Oliveira * * 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 . */ #ifndef TRANSFERFUNCTION_H #define TRANSFERFUNCTION_H #include "ControlElement.h" #include #include "OpenGLText.h" class TransferFunctionForm; /** * @class TransferFunction * @author Thales Lima Oliveira * @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 > A; std::vector B; std::vector 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 GetNumerator() const { return m_numerator; } virtual std::vector GetDenominator() const { return m_denominator; } virtual void SetNumerator(std::vector numerator) { m_numerator = numerator; } virtual void SetDenominator(std::vector 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: *
\f$ y_{n+1}=y_n+\frac{1}{2}h(y'_n+y'_{n+1}) \f$
* @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 m_numerator; std::vector m_denominator; SpaceState m_ss; std::vector m_x; std::vector m_dx; double m_error = 1e-3; int m_maxIteration = 100; }; #endif // TRANSFERFUNCTION_H