diff --git a/MemMapr_DZ_2.sln b/MemMapr_DZ_2.sln new file mode 100644 index 0000000..49be185 --- /dev/null +++ b/MemMapr_DZ_2.sln @@ -0,0 +1,31 @@ + +Microsoft Visual Studio Solution File, Format Version 12.00 +# Visual Studio Version 17 +VisualStudioVersion = 17.14.36616.10 d17.14 +MinimumVisualStudioVersion = 10.0.40219.1 +Project("{8BC9CEB8-8B4A-11D0-8D11-00A0C91BC942}") = "MemMapr_DZ_2", "MemMapr_DZ_2.vcxproj", "{D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}" +EndProject +Global + GlobalSection(SolutionConfigurationPlatforms) = preSolution + Debug|x64 = Debug|x64 + Debug|x86 = Debug|x86 + Release|x64 = Release|x64 + Release|x86 = Release|x86 + EndGlobalSection + GlobalSection(ProjectConfigurationPlatforms) = postSolution + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Debug|x64.ActiveCfg = Debug|x64 + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Debug|x64.Build.0 = Debug|x64 + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Debug|x86.ActiveCfg = Debug|Win32 + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Debug|x86.Build.0 = Debug|Win32 + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Release|x64.ActiveCfg = Release|x64 + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Release|x64.Build.0 = Release|x64 + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Release|x86.ActiveCfg = Release|Win32 + {D677E0E2-76AB-4F8E-87E6-6B8FFC9E15FA}.Release|x86.Build.0 = Release|Win32 + EndGlobalSection + GlobalSection(SolutionProperties) = preSolution + HideSolutionNode = FALSE + EndGlobalSection + GlobalSection(ExtensibilityGlobals) = postSolution + SolutionGuid = {469D3EE6-5DB4-468D-875B-43DF886C506B} + EndGlobalSection +EndGlobal diff --git a/MemMapr_DZ_2.vcxproj b/MemMapr_DZ_2.vcxproj new file mode 100644 index 0000000..6901591 --- /dev/null +++ b/MemMapr_DZ_2.vcxproj @@ -0,0 +1,144 @@ + + + + + Debug + Win32 + + + Release + Win32 + + + Debug + x64 + + + Release + x64 + + + + 17.0 + Win32Proj + {d677e0e2-76ab-4f8e-87e6-6b8ffc9e15fa} + MemMaprDZ2 + 10.0 + + + + Application + true + v143 + Unicode + + + Application + false + v143 + true + Unicode + + + Application + true + v143 + Unicode + + + Application + false + v143 + true + Unicode + + + + + + + + + + + + + + + + + + + + + + Level3 + true + WIN32;_DEBUG;_CONSOLE;%(PreprocessorDefinitions) + true + + + Console + true + + + + + Level3 + true + true + true + WIN32;NDEBUG;_CONSOLE;%(PreprocessorDefinitions) + true + + + Console + true + + + + + Level3 + true + _DEBUG;_CONSOLE;%(PreprocessorDefinitions) + true + + + Console + true + + + + + Level3 + true + true + true + NDEBUG;_CONSOLE;%(PreprocessorDefinitions) + true + + + Console + true + + + + + + + + + + + + + + + + + + Данный проект ссылается на пакеты NuGet, отсутствующие на этом компьютере. Используйте восстановление пакетов NuGet, чтобы скачать их. Дополнительную информацию см. по адресу: http://go.microsoft.com/fwlink/?LinkID=322105. Отсутствует следующий файл: {0}. + + + + \ No newline at end of file diff --git a/MemMapr_DZ_2.vcxproj.filters b/MemMapr_DZ_2.vcxproj.filters new file mode 100644 index 0000000..e09ba90 --- /dev/null +++ b/MemMapr_DZ_2.vcxproj.filters @@ -0,0 +1,30 @@ + + + + + {4FC737F1-C7A5-4376-A066-2A32D752A2FF} + cpp;c;cc;cxx;c++;cppm;ixx;def;odl;idl;hpj;bat;asm;asmx + + + {93995380-89BD-4b04-88EB-625FBE52EBFB} + h;hh;hpp;hxx;h++;hm;inl;inc;ipp;xsd + + + {67DA6AB6-F800-4c08-8B7A-83BB121AAD01} + rc;ico;cur;bmp;dlg;rc2;rct;bin;rgs;gif;jpg;jpeg;jpe;resx;tiff;tif;png;wav;mfcribbon-ms + + + + + Исходные файлы + + + + + + + + Файлы заголовков + + + \ No newline at end of file diff --git a/Solver.hpp b/Solver.hpp new file mode 100644 index 0000000..3f59c93 --- /dev/null +++ b/Solver.hpp @@ -0,0 +1,2 @@ +#pragma once + diff --git a/main.cpp b/main.cpp new file mode 100644 index 0000000..0e10110 --- /dev/null +++ b/main.cpp @@ -0,0 +1,76 @@ +#include +#include +#define M_PI 3.14159265358979323846 // pi + +using namespace Eigen; + +#define STATE_VARS 3*2 // +#define PHI 5 // +#define EDS 2 // +#define MATRIX_DIM STATE_VARS+PHI+EDS // + +const double TK = 1e-3; // +const double T_Start = 1e-10; // +const double SMN = 1e-15; // +const double ACR = 1e-2; // + +// +double TAU = 5e-9; +double MFt = 0.026; + +double C1 = 0.1; +double Cb = 2e-12; + +double L = 0.001; + +double R1 = 1000; +double R2 = 10000; +double Ru = 1e+6; +double Rb = 20; + +double I = 1; +double It = 1e-12; + +double IL_prev = 0; +double UC1_prev = 0; +double UCb_prev = 0; + +double E2 = 5.; +// +double E1(double t) { + const double amplitude = 10.0; + const double period = 1e-4; + const double frequency = 2 * M_PI / period; + const double phase = 0.0; + return amplitude * sin(frequency * t + phase); +} + +// +MatrixXd Y = MatrixXd::Zero(MATRIX_DIM, MATRIX_DIM); + +// +VectorXd x = VectorXd::Zero(MATRIX_DIM); + +void Jacobian(double dt) { + Y(0, 0) = 1, Y(0, 3) = -1. / dt; + Y(1, 1) = 1, Y(1, 4) = -1. / dt; + Y(2, 2) = 1, Y(2, 5) = -1. / dt; + + Y(3, 3) = 1., Y(3, 0) = -1.; + Y(4, 4) = 1., Y(4, 8) = -1., Y(4, 9) = 1.; + Y(5, 2) = L, Y(5, 7) = -1., Y(5, 10) = 1.; + + Y(6, 6) = 1. / R1, Y(6, 7) - 1. / R1, Y(6, 11) = -1.; + Y(7, 5) = 1., Y(7, 6) = -1. / R1, Y(7, 7) = 1. / R1 + 1. / Rb, Y(7, 8) = -1. / Rb; + Y(8, 1) = Cb, Y(8, 7) = -1. / Rb, Y(8, 8) = 1. / Rb + 1. / Ru, Y(8, 9) = -1. / Ru; + Y(9, 1) = -Cb, Y(9, 8) = -1. / Ru, Y(9, 9) = 1. / Ru, Y(9, 11) = 1.; + Y(10, 0) = C1, Y(10, 5) = -1., Y(10, 10) = 1. / R2; + + Y(11, 6) = -1.; + Y(12, 9) = -1.; +} + +int main() { + Jacobian(0.1); + std::cout << Y; +} \ No newline at end of file diff --git a/packages.config b/packages.config new file mode 100644 index 0000000..5e4558c --- /dev/null +++ b/packages.config @@ -0,0 +1,4 @@ + + + + \ No newline at end of file