Дисперсия резисторов параллельно


10

Предположим, у вас есть набор резисторов R, все из которых имеют среднее значение μ и дисперсию σ.

Рассмотрим часть схемы со следующим расположением: (r) || (r + r) || (г + г + г). Эквивалентное сопротивление каждой части составляет r, 2r и 3r. Тогда дисперсия каждого сечения будет σ2 , 2σ2 , 3σ2 .

Какова разница в сопротивлении всей цепи?

После отбора нескольких миллионов точек мы обнаружили, что дисперсия составляет приблизительно .10286σ2 .

Как бы мы пришли к такому выводу аналитически?

Редактировать: Предполагается, что значения сопротивления обычно распределены с некоторым средним сопротивлением r и дисперсией σ2 .


1
Я не уверен, что это подходящая модель для начала. Вы знакомы с теорией теплового шума Найквиста-Джонсона ? Если вы целенаправленно делаете что-то другое, было бы интересно увидеть мотивацию. В противном случае, возможно, стоит рассмотреть более стандартную модель. :)
кардинал

Да, когда я писал свою попытку ответа, я также понял, что модель, по-видимому, не отслеживается, как она была задана. Однако я подумал, что это скорее академическая проблема, нежели практическая (в конце концов, они делают симуляции).
Нестор

Мои извинения за то, что у меня есть сигма как отклонение, я первоначально использовал VAR, и кто-то отредактировал это к сигме.
lrAndroid

Спасибо за обновления. Я по-прежнему заинтересован в мотивации этого вопроса, если вы хотите немного добавить к этому свой вопрос. :)
кардинал

Ответы:


9

Эквивалентное сопротивление всей схемы решает Предполагается, что для некоторых независимых случайных величин центром и с дисперсией .1RRi=iμ+σ

1R=i=131Ri.
Zя1Ri=iμ+σiZiZi1

Без дополнительных указаний невозможно вычислить дисперсию , поэтому, чтобы пойти дальше, мы рассмотрим режим, в котором Тогда следовательно, где Видно, что Кроме того, таким образом, в пределеσ μ . 1R

σμ.
1
1Ri=1iμσμ2Ziii+higher order terms,
a= 3 i=11
1R=aμσμ2Z+higher order terms,
E(Z)=0,
a=i=131i=116,Z=i=13Ziii.
R=μ
E(Z)=0,E(Z2)=b,b=i=131i3=251216.
σ0E(R)μ
R=μaσa2Z+higher order terms,
σ0, и Эти асимптотики и можно обобщить на любое количество параллельных сопротивлений, каждое из которых является результатом последовательных элементарных сопротивлений, причем элементарные сопротивления являются независимыми и каждое из них имеет среднее значение и дисперсию . Затем, когда , где Var(R)σ2b
E(R)μa=611μ,
E(R)Var(R)niμσ2σ0E(R)μ
Var(R)σ2ba4=σ2(611)4251216=σ20.10286
E(R)Var(R)niμσ2σ0а=я1
E(R)μa,σ2Var(R)ba4,
a=i1ni,b=i1ni3.

8

Я не думаю, что точный ответ зависит только от и . Когда вы делали выборку, я полагаю, вы использовали какой-то конкретный дистрибутив - возможно, нормальный дистрибутив? В любом случае мы можем рассчитать среднее значение и дисперсию сопротивления цепи в линейном приближении, и тогда точный вид распределения не имеет значения.σ 2μσ2

Сопротивление цепи . В линейном приближении среднее и дисперсия обратной величины случайной величины со средним значением и дисперсией равны и соответственно. Таким образом, мы имеем сумму слагаемых со средними , и и дисперсиями , и соответственно, что в сумме составляет среднее значение и дисперсию(R11+R21+R31)1μσ21/μσ2/μ41/μ1/(2μ)1/(3μ)σ2/μ4σ2/(8μ4)σ2/(27μ4)116/μ251216σ2/μ4, Тогда взятие обратной дает среднее значение и дисперсию , в соответствии с вашим результатом.611μ(251216σ2/μ4)/(116/μ)4=150614641σ20.10286σ2


Это, конечно, при условии, что резисторы являются независимыми случайными величинами.

@ Роберт: Да (скорее, сопротивления). Это уже предполагалось при расчете дисперсий , и в вопросе, и это имеет физический смысл (хотя, если мы возьмем все резисторы из одной партии продукции, их сопротивления будут несколько коррелированными ). σ2σ3σ
Йорики

Конечно, в реальном дизайне сопротивления далеки от независимых rvs. Фактически, в макете много работы, чтобы некоторые группы элементов отслеживали друг друга (что неудивительно).

1
Вы используете ? Я более привык видеть это как . σ=E(XEX)2σ2

@ copper.hat: Вы, конечно, правы относительно - я, не задумываясь, принял обозначения, использованные в вопросе. σ2
Йорики

5

Это зависит от формы распределения сопротивления. Не зная распределения, я не могу даже сказать среднее сопротивление, хотя я думаю, что есть ограничения.

Итак, давайте выберем распределение, которое можно отследить: пусть будет стандартным отклонением сопротивления одного резистора. Пусть сопротивление будет , причем каждый знак встречается с вероятностью . Это дает нам случая для рассмотрения или если мы объединяем некоторые случаи. Конечно, мы предположим, что сопротивления независимы.sμ±s1/226=642×3×4=24

Если мы выберем и то среднее значение будет (немного меньше, чем ), а дисперсия - . Если мы выберем и , тогда дисперсия будет .μ=100s=154.543291100×6110.102864μ=5s=10.103693

Вот разложение в степенной ряд для отношений между дисперсиями, когда среднее значение равно а дисперсия равна : . Когда мало, доминирующим членом является .1x150614641+360001771561x+21801619487171x2+O(x3)x150614641=0.102862

Хотя вопрос, который вы задаете, технически зависит от распределения, вас, вероятно, интересуют ситуации, когда стандартное отклонение мало по сравнению со средним, и я думаю, что есть четко определенный предел, который не зависит от распределения. Линеаризуйте зависимость сопротивления цепи как функцию сопротивлений каждого элемента:

C=11/R1+1/(R2+R3)+1/(R4+R5+R6)

611μ+i=16(Riμ)CRi(μ,μ,μ,μ,μ,μ)

Var(C)i=16Var(Ri)(CRi(μ,μ,μ,μ,μ,μ))2

В этой конкретной схеме масштабные частные производные имеют следующие значения: и36121,9121,9121,4121,4121,4121

(36121)2+2(9121)2+3(4121)2=150614641=0.102862

1
Это напоминает мне о многомерной дельта-теореме, то есть имеет среднее значение и дисперсию соответственно, затем должны иметь асимптотическую дисперсию как , где и . Окончательный ответ такой же, как у @Douglas Zare и OP, то есть 0.1028 . R1,R2,R3μ,2μ,3μσ2,2σ2,3σ2g(R1,R2,R3)=((1/R1)+(1/R2)+(1/R3))1g(μ)Σg(μ)g(μ)=(36121,9121,4121)Σ=\[(.σ20002σ20003σ2)\]σ2
VitalStatistix

1

Я предупреждаю, что, как я рассуждал, это длинный ответ , но, возможно, кто-то может придумать что-то лучшее, начиная с моей попытки (что может быть не оптимальным). Кроме того, я неправильно прочитал исходный вопрос ОП и подумал, что это говорит о том, что сопротивления обычно распределяются. Я все равно оставлю ответ, но это основополагающее предположение.

1. Физическое обоснование проблемы

Я рассуждаю следующим образом: напомним, что для резисторов, которые находятся в паралеле, эквивалентное сопротивление определяется как:Req

Req1=iN1Ri,

где - сопротивления каждой части цепи. В вашем случае это дает намRi

Req=(1R1+1R2+1R3)1,   ()
где является частью схемы с 1 сопротивлением и поэтому имеет нормальное распределение со средним и дисперсией , и по тем же соображениям является эквивалентное сопротивление части схемы с двумя сопротивлениями и, наконец, - это эквивалентное сопротивление части схемы с тремя сопротивлениями. Вы должны найти распределение и оттуда получить его дисперсию.R1μσ2R2N(2μ,2σ2)R3N(3μ,3σ2)Req

2. Получение распределенияReq

Один из способов найти распределение - это отметить: Здесь также отметим, что мы можем записать (который был получен с помощью теоремы Байеса), который, предполагая Независимость между , и (которая является физически правдоподобной) может быть записана как Заменив это в и отметив, что другим следствием независимости между тремя сопротивлениями является то, что

p(Req)=p(Req,R1,R2,R3)dR1dR2dR3=p(R1|Req,R2,R3)p(Req,R2,R3)dR1dR2dR3.   (1)
p(Req,R2,R3)=p(R2|Req,R3)p(Req,R3)=p(R2|Req,R3)p(Req|R3)p(R3)
R1R2R3
p(Req,R2,R3)=p(R2|Req)p(Req|R3)p(R3).
(1)p(R1|Req,R2,R3)=p(R1|Req)получаем: Наша последняя задача - найти , т. Распределение rv , Эта проблема аналогична той, которую мы нашли здесь, за исключением того, что теперь вы заменяете в уравнении. по константе, скажем, . Следуя тем же аргументам, что и выше, вы можете обнаружить, что Видимо, остальное заменяя известные распределения, за исключением небольшой проблемы: распределение можно получить из , отметив, что
p(Req)=p(R1|Req)p(R2|Req)p(Req|R3)p(R3)dR1dR2dR3=p(Req|R3)p(R3)dR3.   (2)
p(Req|R3)Req|R3R3()r3
p(Req|R3)=p(Req|R2,R3)p(R2)dR2.   (3)
Req|R2,R3()X1 является гауссовским, поэтому вам необходимо найти распределение случайной величины где и - постоянные, и является гауссовым со средним значением и дисперсией . Если мои расчеты верны, это распределение: где поэтому будет
W=(1X+a+b)1,
abXμσ2
p(W)=1[1W(a+b)]212πσ2exp(X(W)μ2σ2),
X(W)=1W1ab,
Req|R2,R3a=1/R2b=1/R3(3)(2)
p(Req|R2,R3)=1[1Req(a+b)]212πσ2exp(X(Req)μ2σ2),
где и . Дело в том, что я не знаю, является ли это аналитически управляемым для решения интеграла в уравнении , что затем приведет нас к решению задачи, заменив ее результат в уравнении . По крайней мере, для меня в это время ночи это не так.a=1/R2b=1/R3(3)(2)

Вы предполагаете нормальное распределение, хотя сопротивление не может быть отрицательным? Я предполагаю, что это приведет к расхождению дисперсии схемы.
Дуглас Заре

1
Я знаю, это меня тоже беспокоило, но на практике это действительно зависит от значений и . Если и , то мы можем «сохранить» модель. В нормальных условиях разброс сопротивления не очень высок, поэтому последнее предположение явно выполняется. Это было то, что поначалу меня тоже беспокоило, когда люди моделировали рост как нормальную случайную величину, но по той же причине, что я привел здесь, некоторые люди здесь, на Stack-exchange, заставили меня чувствовать себя хорошо :-). сг 2 ц > > 0 ц > > сгμσ2μ>>0μ>>σ
Нестор

Хм, я думаю, что высота моделирования как нормальная настолько плоха, что я использую ее как пример распределения, которое явно не является нормальным. Я предполагаю, что это может быть не страшно, если у вас есть популяция здоровых взрослых мужчин с таким же генетическим фоном. Тем не менее, я хотел бы услышать от биолога, что это нормально. Я слишком часто слышал рассуждения о том, что размер каждой кости независимый, это полная чепуха.
Дуглас Заре

Я только что понял, что сопротивления обычно не распределяются (могу поклясться, что прочитал их там, где на оригинальном ответе ОП, но я думаю, что это было только мое воображение).
Нестор
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.