Міністерство освіти і науки України при Національному університеті "Львівська політехніка" РОЗРАХУНКОВА РОБОТА З ПРЕДМЕТУ «ТЕОРІЯ ЕЛЕКТРИЧНИХ KІЛ» Задача Дано: R=2 Ом; l=1,5 м; N=4;
м; f=30+2N=38; w=13000/f=342; Математична модель


Будуємо дискретну модель за явним методом Ейлера:

Алгоритм прямого методу Ейлера: t=0; Задаємося початковими значеннями
Обчислюємо статичну індуктивність
Обчислюємо струм
Запам’ятовуємо значення
Перевіряємо умову чи час t досяг свого кінцевого значення tEND. Якщо так то виходимо з алгоритму якщо ні то йдемо далі. Обчислюємо значення потокозчеплення згідно із явною формулою Ейлера на наступному кроці.
Змінюємо крок
Йдем на пункт 3. Апроксимація кривої намагнічування методом вибраних точок Нам дана характеристика намагнічування сталі 3411 В, Тл 1,0 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9 2,0 2,1  Н, А/м 170 220 280 370 500 700 1000 2000 4000 9200 40000 90000   Графічно вона виглядає так:
Для того щоб апроксимувати ми виберемо 2 точки з яких нарисуємо прямі. Нашим рядом апроксимації являється вираз:
Рівняннями прямих будуть:
і
Запишемо систему рівнянь:

На графіку це виглядатиме так: де вибрані точки будуть в: В1=1,3; В2=1,8; В3=1,6, відповідно. Обрахуємо коефіцієнти при прямих:



Тепер знайдемо коефіцієнти а4, а5, а6:
розв’язавши цю систему ми отримаємо:


Остаточне рівняння:
Перерахована крива:
Аналіз перехідних процесів з розв’язанням диференціальних рівнянь явними методом Ейлера. Текст програми написаний мовою Сі++ void __fastcall TForm1::Button1Click(TObject *Sender) { double const l=1.5,S=25E-4,N=4; double B,f,w,Af,t,psi,h,i, ,T; f=30+2*N; w=13000/f;h=1/(f*360);T=1/f; B=l/(w*S); psi=0;i=0;t=0; fstream file; file.open("viv(t-psi-i).txt",ios_base::out|ios_base::trunc); for (double t=0;t<=3.25*T;t=t+h) { file<<t<<" "<<psi<<" "<<i<<endl; psi=psi+h*(310*sin(6.28*f*t+30)-2*i); i=5076.455*B*psi+4965.48*B*B*B*psi*psi*psi+120.37*pow(B,4)*pow(psi,4); } file.close(); } Графік залежностей і від ψ:
Графік заложнестей t від i:
Графік заложестей t від ψ: