���������� 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.
�� ��� ����������� ����� � ������� ���������� ���������, ��������������� ��������� ��������� �������� flow � fhigh, � ����� ���������� �����������, ��������������� ������������ �������� �� ���� ��������. ��� ��������� �������� ������� ��������� �������� ������ ���������.
����� �������, ������ ������� ��������� � �������������� ������ ��������� ������. �� �������� �������� ������������ ������������ ��������, ���������� ������ � ������� ��������� �������, �������� ������������� ��������� ��� ����������-��������� � ����-��������� ������������� ��������� � ��������� �� ��������� �������.
���. 17.15. ����������-��������� �������������� � ����������� ������� ��������� ������