http://pspicelib.narod.ru

ElectronicsSolver считает эту схему за 5 минут

          На рис. 1 представлена достаточносложная схема, которая решена аналитически. В реальной практике аналитических расчёто её следовало бы упростить, избавившись от транзисторных каскадов. Однако это уже относится к технике получения решения, которая у каждого индивидуальна. Упрощения сложных схем, без потери существенной информации следует отнести к искусству . Очевидно, что эта операция должна быть самой первой и хорошо продуманной, и уже на этом этапе у вас должен созреть план решения.
задачи в целом.

Рис. 1. Схема  фильтра

[Maple Bitmap]

>    restart:with(MSpice):Приборы:=Oдинаковые:
ESolve(Q,`KKBJT.CKT`):

`Cистема Кирхгофа`

-V4/R15-V4*s*C3-(V4-V0_2)/s/(L1[2]-K2*sqrt(L1[2]*L2[2])) = 0

(V1-V0_1)/s/(L2[1]+K1*sqrt(L1[1]*L2[1]))+V0_1/s/K1/sqrt(L1[1]*L2[1])-(V0_1-V2)/s/(L1[1]+K1*sqrt(L1[1]*L2[1])) = 0

(V4-V0_2)/s/(L1[2]-K2*sqrt(L1[2]*L2[2]))+(V3-V0_2)/s/(L2[2]-K2*sqrt(L1[2]*L2[2]))-V0_2/s/K2/sqrt(L1[2]*L2[2]) = 0

-V5/R11+(V6-V5)/`Rэб`-(V5-VIN)*s*C7-V5/R9 = 0

(V3-V2)*s*C5-V2/R7+(V0_1-V2)/s/(L1[1]+K1*sqrt(L1[1]*L2[1]))-V2*s*C2 = 0

-V1*s*C1-(V5-V6)*beta/`Rэб`-V1/R13-(V1-V0_1)/s/(L2[1]+K1*sqrt(L1[1]*L2[1])) = 0

-V3/R8-V3*s*C4-(V3-V2)*s*C5-(V3-V0_2)/s/(L2[2]-K2*sqrt(L1[2]*L2[2])) = 0

-V6/R10-V6*s*C6+(V5-V6)*beta/`Rэб`-(V6-V5)/`Rэб` = 0

-VOUT/R14+(V4-VOUT)*S[J] = 0

{V3, V4, V0_2, V1, V2, V5, V6, VOUT, V0_1}

MSpice v8.96:   http://pspicelib.narod.ru

Заданы источники: [Jэ, VIN]

Заданы узлы: {VINP}

Получены решения:

V_NET:=[V3, V4, V0_2, V1, V2, V5, V6, VOUT, V0_1]:

J_NET:=[JR11, JRэб, JC5, JR10, JC1, JC2, JC3, JC4, JC6, JC7, JR7, JR8, JR9, JR15, JR13, JR14, JVIN, JL1[1], JL2[1], JL1[2], JL2[2]]:

>    K2:=K1: OUT:=simplify(VOUT):

>    Values(AC,RLCVI,[K1,K2]); Digits:=3:

Номиналы компонентов:   

C1:=150e-12:  [150P]

C2:=146e-12:  [146P]

C3:=150e-12:  [150P]

C4:=146e-12:  [146P]

C5:=4e-12:  [4P]

C6:=10e-9:  [10N]

C7:=1e-9:  [1N]

R7:=500e3:  [500K]

R8:=500e3:  [500K]

R9:=5e3:  [5K]

R10:=300:  [300]

R11:=1e3:  [1K]

R13:=100e3:  [100K]

R14:=1e3:  [1K]

R15:=100e3:  [100K]

S[J]:=3.5e-3:

beta:=100:  [100]

Jэ:=1e-3:

Rэб:=(k*Tkelvin/e*(1+beta)/Jэ):

L1[1]:=0.0008:  [0.0008]

L2[1]:=0.0008:  [0.0008]

L1[2]:=0.0008:  [0.0008]

L2[2]:=0.0008:  [0.0008]

K2:=0.025:  [0.025]

e:=.1602176462e-18:

k:=.1380650277e-22:

Tkelvin:=300:  [300]

AC источник: VIN:=0:  AC: VIN:=1e-3:   Pfase(degrees):=0:

>    VOUT:=simplify(VOUT);

VOUT := -.280e33*s^4*K1*(.100e7+3.*s)/(-.216e47*s+233.*K1^2*s^8-.940e21*s^5-.972e29*s^4+.144e38*K1^2*s^2-.112e55+.130e34*K1^2*s^3+.389e16*K1^2*s^6-.540e42*s^2-.777e34*s^3-233.*s^8-.778e16*s^6-.375e8*s^...
VOUT := -.280e33*s^4*K1*(.100e7+3.*s)/(-.216e47*s+233.*K1^2*s^8-.940e21*s^5-.972e29*s^4+.144e38*K1^2*s^2-.112e55+.130e34*K1^2*s^3+.389e16*K1^2*s^6-.540e42*s^2-.777e34*s^3-233.*s^8-.778e16*s^6-.375e8*s^...
VOUT := -.280e33*s^4*K1*(.100e7+3.*s)/(-.216e47*s+233.*K1^2*s^8-.940e21*s^5-.972e29*s^4+.144e38*K1^2*s^2-.112e55+.130e34*K1^2*s^3+.389e16*K1^2*s^6-.540e42*s^2-.777e34*s^3-233.*s^8-.778e16*s^6-.375e8*s^...

>    Digits:=10: OUT:=evalf(abs(subs(s=I*2*Pi*f,OUT))):

>    plot([eval(OUT,K1=0.01),eval(OUT,K1=0.015),eval(OUT,K1=0.025),eval(OUT,K1=0.035)],
f=4.3e5..4.9e5,
thickness=3,
legend=["K=0.01","K=0.015","K=0.025","K=0.035"],
color=[black,red,blue,gold],
labels=[`Частота f (Гц)`,"OUT"],
title=`Рис. 2. АЧХ фильтра`);

[Maple Plot]

Hosted by uCoz