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

���������� MATLAB

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

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

 

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

���������� ������������ ������������ ������ �� ������� �����������, ����� �������� ��������� �� ���. 17.12, �. ��� ������������� �������������� ����� ������������ �� ���. 17.12, �.

a

���. 17.12. �������������� (�) � �������������(6) ����� ��������� �� ������� �����������

���� ���� ����������� � ������� ������������� ��������� ����������� �������. ��������� ������ ��� ����� plots:

> restart:with(plots):

Warning, the name changecoords has been redefined

�� ������� ������� ��������, ��� ����� �����, ��������� � ������ ���� � ���������� �� ���� ����� 0. �������������, ��� ����� ������������� ����� ���. 17.12 ����� �������� ��������� ������� ��������� � ����������� �����:

 

���������� ���������� �� ����� ����� ��������� �� �������������� ������� ������ �������. ��� ���� ����������� ����� �� ������������� ��������, ��������� �� ����� ���������. ��� �� ����� �� ������ ���������� �� ���������� �������, �������� ���� ����� � ������� ������ ����� ��������� � ����������� ���� ����������:

> solve({eql,eq2,eq3.eq4}�{Vl,V2.V3,Vo}):

��������� ���������� ���������� Vo, VI, V2 � V3 ��������� �� ������� ������� ��������� ��������:

> assign(%):

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

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

��� ��������� ����� � ��� ������� �� ������� f ����� � ������������� ����:

��� ���� �������� ���������� ���������, � ��� ���������� ��� «������» ������� ����������� �� �� ��� ������ �����������������. ����� ��� Maple 7 ��������� «� ��� �����» ���������� �� ���� ����������-��������� (AVM) � ����- ��������� (PhaseAV) �������������� ��������� ��� ������� �������:

> AVM=-evalc(abs(H)):

> PhaseAV:=evalc(argument(H)):

����������� AVD � ��������������� ��������������, ���������� �������� � ��������� (dB):

> AVdB:=20*1og10(AVM):

����� �������������� ����� �������� ��� ������������ � ����������������. �������������� ����-��������� �������������� ������� � ��������:

> R2D:=evalf(360/(2*Pi));R2D := 57.29577950 

> AVdeg:=R2D*PhaseAV:

������ ����� ������� � ������� ��������� ��������. ������� ���������� �������� ��������� ������������� ����� ���������:

> Rl:=100: R2:=100000: R3:=1000: R4:=10000: Cl:=1.*10^(-6): �2:=5*10^(-12): ��:=1*10^(-6): mu:=50:

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

> gaindata:-NULL:

 phasedata:=NULL:

 for a from 0 to 8 do:

 for i from 2*10^a to l(T(a+l) by 10^a do

 gaindata:=gaindata,  [1. evalf(subs(f=i,AVdB))];

  phasedata:=phasedata, [i, eva1f(subs(f=i,AVdeg))]:

  od: od: 

> 1oglogp1ot([gaindata]. thickness»2, color=black, style=1ine, axes=boxed,

title=`����������� �������� K(f)`,1abels=['������� (Hz)VK(d8)']):

��� �������� �� ���. 17.13.

���. 17.13. ����������-��������� �������������� ���������

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

> 1og1ogplot([phasedata], thickness=2, color=b1ue, style=line, axes=boxed, title='������� ����� (� ��������)`, labels=['������� (Hz)','����']);

��� ������������ �� ���. 17.14.

���. 17.14. ����-��������� �������������� ���������

������ ����������� ����������� �������� �� ������� f=1000 (��):

> AVmid:=eva1f(subs(f=1000, AVdB)):

AVmid=33.12074854

���� ������������� ��������� ��� ����������-��������� ��������������, ����� ��������� ��������� ��� ���������� ��������� ������ (�� ����� �������� �� -dAV � dB):

> dAV:=3:   #���������� (� dB �� ��������� ��������) 

> eq5:=AVmid-dAV=20*log10(AVM):

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

> flow:=fsolve(eq5,f. f-10..2000):flow:= 23.61659476

> fhigh:=fsolve(eqS,f, f-2000..100*10*6);

fliigh := .5737800225 107

�� ����� ��������� � ����� ��������� ����������-��������� �������������� � �������, ���������������� ��������� ��������:

> with(plottools) :h:=log10(AVnvid-dAV):

aplot:= Loglogplot([gaindata], thickness=2, color=b1ack. style=line, axes=boxed,

title='������� flow � fhigh �����', labels=['������� (Hz)VK(dB)']):

bplot:=line([0.1,h], [7.1,h], color=black, linestyle=3):

cplot:=line([log10(flow),0.58],[logHK flow). 1.6], color=blue, linestyle=3):

dplot:=line([log10(fh1gh).0.58],. [log10(fhigh).1.6],. color=red,. 1inestyle=3):

display([aplot.bplot,cplotJ,dplot]):

��� �������������� �������� �� ���. 17.15.

�� ��� ����������� ����� � ������� ���������� ���������, ��������������� ��������� ��������� �������� flowfhigh, � ����� ���������� �����������, ��������������� ������������ �������� �� ���� ��������. ��� ��������� �������� ������� ��������� �������� ������ ���������.

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

���. 17.15. ����������-��������� �������������� � ����������� ������� ��������� ������

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