���������� MATLAB
������ ������� � ���������� ����������
���������� ���������� � ��������� ��������, �. �. ������ � �������� �-�����, ������������� �������� �������. ������� help ���� ������� �����������, ����������� ������ � �������� �������, ����� ��� ������� type name ������� ���� ������� m-�����. ��� ��� �������� � ��������� m-����� ���������� ������ ������������ �� � ��������� ������ ������, �� � ������ �-������. ��� ��������� � ������� ������������� ������� MATLAB ������ ����� ������������� m-���� �� ������� ��������� ����������. ���� ��� ����������, �� �������� ��� ��������� �������. ��������� ���������� �������� ���������� �� ����� ����� ��������� � ������� �������� ������� �������. ���� � ������� � ����������� ������ �������������� ����� ����������, �� ���� �������� ����������� ��� � �������, ��� � � �� �����������.
������������� ���� �� ���������� �����
� ������ �������� �������������� ���� ���������� ����. �������� ��� ��� ����, ������� ������� �� ��������������� ���������� ��������� ���������� Es, ��������� Rs, ������������� L � ����������� �����, �������� N-�������� ������������� �������������� (���). ���������� ���� �������� �������� �, ��� ����������� ������������� �, ������������ ����������� ����������� �����. ����� ��� ��������� ����������� ����� ������ ����������:
> restart:
> A:=.3t: �:=10: �:=1*10^(-8): b:=20:
> Id:=Ud->A*Ud*exp(-a*Ud)+B*(exp(b*Ud-D):
Id:=Ud->AUde(-aUd)+Be(bUd-1)
�������� ������ ���:
> plot(Id(Ud), Ud=-.02..0.76,color=black):
���� ������ ����������� �� ���. 17.25. �������� ��������, ��� ��� ����������� ����� �� ������ ����� ���������, �� � �������� ����������� ������� ������������� ���������������� ������������, �� ������� ��� ������ � ������ ����������. ��� �������� ��������� ����, ��� ����� ���� �������� �� ���������� ���� �������� ������� �� ������� ���� � ��������� � ������������� ��������� � ��� ���������� ����.
������ ���� ����������� �������� �� ���� ���������������� ���������:
di/dt=(Es-i(t)*Rs-u(t))/L
du/dt=(i(t)-Id(u(t))/C
���. 17.25. ��� ����������� �����
����� ������ Es = 0,35 �, Rs= 15 ��, � = 10*10-12, L = 30*10-9 � ������������ ����� ������������� tm=10*10-9. ����, ������ �������� ������:
> Es:=.35:Rs:=15:C:=10*10^(-12):L:=30*10^(-6):tm:=10*10^(-9):
�������� ������� ���������������� ��������� ���� � �������� �� ������� � ������� ������� dsolve:
��������� �������� ��������, ��� ����� ����� ����� �������� L � �, �� ������ � ������� ���������� ���������� ����� ��� ������� ��� ������� dsolve � stepsize=l(T(-11) (�). ��� ������� ����� �������� ��������� �������������� �������, ���������� ����� ���������, ���������� ��� �������������. ��������� ������� odeplot � displ ay ������ plots, �������� ������� ������� � ���� ��������� ������������ u(t) � 10*i (t) � �����, ��������������� ���������� Es ��������� �������:
> gu:=odeplot(F,[t,u(t)],0,tm,color=black,
labels=['tVu(t),10*i(tr]):
> gi:=odeplot(F,[t,10*i(t)],0..tm.color-black):
> ge:=odeplot(F,[t,Es].0..tm.color=red): .
> display(gu.gi,ge);
��� ����������� ������������ �� ���. 17.26. �� ��� ������ �����, ��� ���� ������� ������������� ��������������� ����. �� ����� ������ ���������� �� ��������������.
���. 17.26. ��������� ����������� ���������� �� ���������� ����� � ����
������� ����� ����������� ����� � ���� �������� ��������, ������������ �� ���� ����������� ��� � ����� �������� ��������� Rs:
> gv:=plot({Id(Ud),(Es-Ud)/Rs},Ud=-.05..0.75,color=black,
labels=[Ud,Id]):
> gpp:=odeplot(F.[u(t),i(t)],0..tm,color=blue):
> display(gv,gpp);
������� ������� ��������� ������� �� ���. 17.27.
���. 17.27. ������� ������� ��������� �� ���� ��� ����������� ����� � ����� �������� ��������� Rs
� ���, ��� ��������� �������������� ����� ������ �� ����, ��� ��� ������ ���� ��������� ����������� � ��������� ������ � ���������� ����, ����� �������� ������� ���������� �� �������������.
����, �� �����, ��� ������ ���� ��������� ������� ���������� ������������ �������������� ���������. ���� ������������ ������ ������������� ���� �� ���������� ����� ������� ������, � ���� �� ������� �� ����������� � ��������� ����������� ������ ���� �� ������� ��� ������� ������� ���������������� ���������, ����������� ������ ����. ��� ��������� ������ ���� ����� ��������� ����� �������������� �������.