�������������� ������ ����������������� ������ ������� ������� �������� �������� ��������� �������� ���������������� ���� �������� �� ������� ������ ��������� ������ �������������� ������������� ������ ���������������� ����������� ������

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

���������� MATLAB

������ ������� � ���������� ����������

���������� ���������� � ��������� ��������, �. �. ������ � �������� �-�����, ������������� �������� �������. ������� help ���� ������� �����������, ����������� ������ � �������� �������, ����� ��� ������� type name ������� ���� ������� m-�����. ��� ��� �������� � ��������� m-����� ���������� ������ ������������ �� � ��������� ������ ������, �� � ������ �-������. ��� ��������� � ������� ������������� ������� MATLAB ������ ����� ������������� m-���� �� ������� ��������� ����������. ���� ��� ����������, �� �������� ��� ��������� �������. ��������� ���������� �������� ���������� �� ����� ����� ��������� � ������� �������� ������� �������. ���� � ������� � ����������� ������ �������������� ����� ����������, �� ���� �������� ����������� ��� � �������, ��� � � �� �����������.

������ ����������� ������� �� ������������ ���������

������ ���������� �������������� ����������� ���������� ������� �� ������������ ���������, ����� �������� ��������� �� ���. 17.16.

���. 17.16. ����� ���������� ������� �� ������������ ������������ ���������

������������ � ������� �������:

> restart:

������� �������� ���������, ����������� ������ ������� �� ����� �������:

������ �������� �������:

> omega := 2*Pi*f;

W := 2�f

������ ����������� �������� ������� � ��� ����-��������� �������������� ��� ������� �� �������:

> gain := abs(eva1c(Vo/Vi)):

> phase := evalc(op(2,convert(Vo/Vi.polar))):

��� ��������� ���������� ������������� ��������� ��� ���� ���������� �������� ����� ��������� � ��������� ��� gain � phase �� ���� ����� � �������. ����� ������ ���������� �������� ������ ��� �������:

> R3 :=1000:

> R4 := 3000:

> �� :=0.08*10^(-6):

> �4 := 0.01*10^(-6):

�������� ��� ������� ��� ����������� ������������ �������� � ��������� (dB) �� ������� f � ��:

> plot(DogWf), 20*log10(gain), f=[10..50000], color=black, title='����������� �������� dB ��� ������� �� ������� f � ��'):

��� �������������� ������������ �� ���. 17.17. ����� ������� �������� �������� �� ��, ��� ���� �������� �� ������ � ������� �������� ���������� �������� �������� ��-�� ������ ������� �������.

���. 17.17. ��� ������� �� ������������ ���������

����� �������� ����-��������� �������������� ������� ��� ����������� ���� � �������� �� ������� f � ��:

> plot ([log10(f),phase, f=10..50000], color=black, title=*����-��������� �������������� �������*);

����-��������� �������������� (���) ������� �������� �� ���. 17.18

�� ��� ������� ����� �������� ����������� ������, ��������� � ����������� ������� ����� ���������� �������� �. ����� ������ ������������� �������� ������ ����������, ��������� ��� ��������� ��������� ������� � �������� �� -� �� �.

���. 17.18. ��� ������� �� ������������ ���������

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