Технический форум

Технический форум (http://www.tehnari.ru/)
-   Помощь студентам (http://www.tehnari.ru/f41/)
-   -   Помогите с модифицированным методом Эйлера (http://www.tehnari.ru/f41/t88274/)

Alexei18 15.05.2013 15:34

Помогите с модифицированным методом Эйлера
 
помогите пожалуйста, уже не знаю что делать, препод дал курсач и мне достался модифицированный метод Эйлера, а в примере только усовершенствованный, а учитель такой, что ни в чем не поможет и ни чего не объяснит. Прошу пожалуйста хотябы примерчик какой нито скиньте, а я по нему попробую сделать. Внизу усовершенствованный метод.

Program Kursach;
uses crt;
const n=200;
f=90;
fi=72;
L=0.00557;
C=0.00002;
W=2*PI*f;
PI=3.14;
E0=15;
R1=30;
R2=25;
R3=50;
R4=1.88;
R5=15;
R6=50;
t1=0.003;
t2=0.009;
var A,B,G,D,U0,I0,t0,t,I,U,hx,k,j,E,KU1,KU2,KI1,KI2:re al;
m:integer;
O:TEXT;
function E1(ft:real):real;
begin
E:=E0+E0*sin(W*ft+fi);
If (t<t1)Then
E1:=E
Else
E1:=0;
end;
begin
assign(O,'D:\result.txt');
rewrite(O);
k:=0;
J:=5;
t0:=0;
U0:=0;
I0:=0;
hx:=(t2-t0)/n;
Writeln(O,' t I U ');
Writeln(O,t:9:4,I:11:7,U:11:6);
A:=R2/(R1+R2);
B:=R1*R2/(R1+R2);
G:=(R5+R6)/(R3+R5+R6);
D:=1/(R3+R5+R6);
For m:=1 to n do
begin
k:=k+1;
KU1:=hx*(1/C*(I0*G-U0*D));
t:=t0+hx;
KU2:=hx*(1/C*((I0+KI1)*G-(U0+KU1)*D));
KI2:=hx*(1/L*(E1(t0+hx)*A-(I0+KI1)*(R4+B+G*R3)-(U0+KU1)*G));
U:=U0+((KU1+KU2)/2);
I:=I0+((KI1+KI2)/2);
I0:=I;
U0:=U;
t0:=t;
if k=j then
begin
Writeln(O,' ',t:6:4,' ',I:10:7,' ',U:10:6);
j:=j+5;
end;
end;
close(O);
end.

Vladimir_S 15.05.2013 16:07

Цитата:

Сообщение от Alexei18 (Сообщение 907989)
помогите пожалуйста, уже не знаю что делать, препод дал курсач и мне достался модифицированный метод Эйлера, а в примере только усовершенствованный, а учитель такой, что ни в чем не поможет и ни чего не объяснит. Прошу пожалуйста хотябы примерчик какой нито скиньте, а я по нему попробую сделать.

Примерчик хотите? А пожалуйста: http://www.tehnari.ru/f41/t87015/

Alexei18 15.05.2013 17:31

спасибо!!!надеюсь все получится.

Alexei18 15.05.2013 17:46

не подскажите как мне переделать итерационные формулы под мои, что то мне не совсем все понятно...

Vladimir_S 15.05.2013 20:26

Цитата:

Сообщение от Alexei18 (Сообщение 908017)
не подскажите как мне переделать итерационные формулы под мои, что то мне не совсем все понятно...

Ну так может быть Вы всё-таки снизойдёте до формулировки задачи? Вид уравнения, опорная точка, диапазон и т.п.

Alexei18 16.05.2013 13:12

• гармонический источник тока, B E(t)=E0(t)•sin(ωt+φ)
• амплитуда колебаний, В E0=15
• катушка индуктивности, мГн L=5.57
• конденсатор, мкФ C=20
• циклическая частота ω=2πf
• текущее время, c T1=0,003
T2=0.009
• линейная частота, Гц 90
• фаза, рад 2*пи/5
• резисторы, Ом R1=30, R2=25, R3=50, R4=1.88, R5=15, R6=50

Решить Дифференциальное уравнение Модифицированным методом Эйлера.

Vladimir_S 16.05.2013 14:01

Цитата:

Сообщение от Alexei18 (Сообщение 908180)
Решить Дифференциальное уравнение Модифицированным методом Эйлера.

Эх, еще бы самО уравнение увидеть... Ну или по крайней мере схему. Мечты, мечты...

Alexei18 16.05.2013 14:30

Вложений: 2
вот надеюсь эти формулы.

Vladimir_S 16.05.2013 15:33

Цитата:

Сообщение от Alexei18 (Сообщение 908199)
вот надеюсь эти формулы.

Черт, да тут система... Тогда не знаю, сумею ли быстро помочь. Надо вспоминать, разбираться. Попробую, конечно, но ничего не обещаю.

Alexei18 16.05.2013 16:05

большое вам спасибо за вашу помощь, я очень вам благодарен,надеюсь у вас получится.


Часовой пояс GMT +4, время: 19:32.

Powered by vBulletin® Version 4.5.3
Copyright ©2000 - 2024, Jelsoft Enterprises Ltd.