������� ������� ������-����������� �����

���������� 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) (�). ��� ������� ����� �������� ��������� �������������� �������, ���������� ����� ���������, ���������� ��� �������������. ��������� ������� odeplotdispl 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

� ���, ��� ��������� �������������� ����� ������ �� ����, ��� ��� ������ ���� ��������� ����������� � ��������� ������ � ���������� ����, ����� �������� ������� ���������� �� �������������.

����, �� �����, ��� ������ ���� ��������� ������� ���������� ������������ �������������� ���������. ���� ������������ ������ ������������� ���� �� ���������� ����� ������� ������, � ���� �� ������� �� ����������� � ��������� ����������� ������ ���� �� ������� ��� ������� ������� ���������������� ���������, ����������� ������ ����. ��� ��������� ������ ���� ����� ��������� ����� �������������� �������.

������, ���������� ������ �������� �������� �������� � ���������