ВУЗ:
Составители:
x
1
1:=y2
1
y2
0
dt a2− y2
0
⋅
()
⋅+:=y1
1
y1
0
dt a1− y1
0
⋅
()
⋅+:=
x
0
0:=y2
0
0:=y1
0
0:=dt 0.25:=
4. Задаем шаг дискретизации dt и вычисляем значения функций
y
k
,
y2
k
в точках
t
k
, k=1, 2
ORIGIN 0:=
t
y2 t()
d
d
a2 y2 t()⋅+ a2 y1 t()⋅
t
y1 t()
d
d
a1 y1 t()⋅+ Ka1⋅ xt()⋅
a2 s2−:=a1 s1−:=
3. Представляем исходное дифференциальное уравнение в виде
системы дифференциальных уравнений первого порядка.
gt( ) 37.5 1 2 exp s1 t⋅()⋅+ 3 exp s2 t⋅()⋅−()⋅:=
2. С применением встроенной функции системы MathCAD находим
аналитическое решение уравнения при x(t) = 1(t).
s2 0.133−:=s1 0.2−:=
hs( ) 37.5 s
2
⋅ 12.5 s⋅+ 1+:=
1. С применением встроенной функции системы MathCAD находим
корни характеристического уравнения.
Yp 0:=Yo 0:=K 37.5:=
Дано: 37.5y"(t) + 12.5y'(t) + y(t) = Kx(t). Начальное значение искомой
функции y(0) = Yo, производной - y'(0) = Yp; x(t) - возмущающая
функция (входной сигнал); dt - шаг дискретизации (с); K - коэффициент.
ПРИЛОЖЕНИЕ 3
Решение дифференциального уравнения
методом Адамса
44
ПРИЛОЖЕНИЕ 3 Решение дифференциального уравнения методом Адамса Дано: 37.5y"(t) + 12.5y'(t) + y(t) = Kx(t). Начальное значение искомой функции y(0) = Yo, производной - y'(0) = Yp; x(t) - возмущающая функция (входной сигнал); dt - шаг дискретизации (с); K - коэффициент. K := 37.5 Yo := 0 Yp := 0 1. С применением встроенной функции системы MathCAD находим корни характеристического уравнения. 2 h ( s) := 37.5 ⋅ s + 12.5 ⋅ s + 1 s1 := −0.2 s2 := −0.133 2. С применением встроенной функции системы MathCAD находим аналитическое решение уравнения при x(t) = 1(t). g( t) := 37.5 ⋅ ( 1 + 2 ⋅ exp( s1 ⋅ t) − 3 ⋅ exp( s2 ⋅ t) ) 3. Представляем исходное дифференциальное уравнение в виде системы дифференциальных уравнений первого порядка. a1 := −s1 a2 := −s2 d y1 ( t) + a1 ⋅ y1 ( t) K ⋅ a1 ⋅ x( t) dt d y2 ( t) + a2 ⋅ y2 ( t) a2 ⋅ y1 ( t) dt ORIGIN := 0 4. Задаем шаг дискретизации dt и вычисляем значения функций yk, y2k в точках tk, k=1, 2 dt := 0.25 y10 := 0 y20 := 0 x0 := 0 y11 := y10 + dt ⋅ ( −a1 ⋅ y10) y21 := y20 + dt ⋅ ( −a2 ⋅ y20) x1 := 1 44