* Данный текст распознан в автоматическом режиме, поэтому может содержать ошибки

1166 Аналитическое и спектральное моделирование > se:=diff(i(t),t)=(Es-i(t)*Rs-u(t))/L, diff(u(t),t)=(i(t)-Id(u(t)))/C; > F:=dsolve({se,i(0)=0,u(0)=0},{i(t),u(t)},type=numeric, method=classical, stepsize=10^(-11), output=listprocedure); F := [t = (proc(t) … end proc), u(t) = (proc(t) … end proc), i(t) = (proc(t) … end proc)] Поскольку заведомо известно, что схема имеет малые значения L и N, мы зада ли с помощью параметров достаточно малый шаг решения для функции dsolve – stepsize=10^(-11) c. При больших шагах возможна численная не устойчивость решения, искажающая форму колебаний, получаемую при модели ровании. Используя функции odeplot и display пакета plots, построим графи ки решения в виде временных зависимостей u(t) и 10*i(t) и линии, соответствующей напряжению Es источника питания: > gu:=odeplot(F,[t,u(t)],0..tm,color=black, labels=[`t`,`u(t),10*i(t)`]): > gi:=odeplot(F,[t,10*i(t)],0..tm,color=black): > ge:=odeplot(F,[t,Es],0..tm,color=red): > display(gu,gi,ge); Эти зависимости представлены на рис. 15.40. Из них хорошо видно, что цепь создает автоколебания релаксационного типа. Их форма сильно отличается от синусоидальной. Рис. 15.40. Временные зависимости напряжения на туннельном диоде и тока