并联电阻的变化


10

假设您有一组电阻R,所有电阻均以均值μ和方差σ分布。

考虑具有以下布局的电路的一部分:(r)|| (r + r)|| (r + r + r)。每个部分的等效电阻分别为r,2r和3r。然后每个部分的方差是σ22σ23σ2

整个电路的电阻变化是多少?

在对数百万个点进行采样之后,我们发现方差约为.10286σ2

我们将如何分析得出这个结论?

编辑:假设电阻值是正态分布的,具有一些平均电阻r和方差σ2


1
我不认为这是一个合适的模型。您是否熟悉热电路噪声的奈奎斯特-约翰逊理论?如果您有目的地做一些不同的事情,那么看看动机会很有趣。否则,考虑使用更标准的模型可能是值得的。:)
红衣主教2012年

是的,当我尝试回答一个问题时,我也意识到该模型显然不易处理。但是,我认为这更像是一个学术问题,而不是实际问题(毕竟,他们正在进行模拟)。
内斯托尔·

我很抱歉以sigma作为方差,我最初使用VAR,有人将其编辑为sigma。
lrAndroid

感谢更新。如果您愿意在这个问题上加一点点,我仍然对这个问题背后的动机感兴趣。:)
主教

Answers:


9

整个电路 的等效电阻 有人假设,对于一些独立的随机变量,居中且具有方差。R

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

如果没有进一步的指示,则无法计算的方差,因此,为了进一步讲解,我们考虑了所在的体制 然后, 因此 其中 可以看到 此外, 因此,在极限R

σμ.
1Ri=1iμσμ2Ziii+higher order terms,
1R=aμσμ2Z+higher order terms,
a=i=131i=116,Z=i=13Ziii.
E(Z)=0,E(Z2)=b,b=i=131i3=251216.
R=μaσa2Z+higher order terms,
σ0, 和 这些和推广到任意数量的并联电阻,每个电阻都是基本电阻串联的结果,基本电阻是独立的,并且每个均值均具有和方差。然后,当, 其中
E(R)μa=611μ,
Var(R)σ2ba4=σ2(611)4251216=σ20.10286
E(R)Var(R)niμσ2σ0
E(R)μa,σ2Var(R)ba4,
a=i1ni,b=i1ni3.

8

我不认为确切的答案仅取决于和。当您采样时,我想您必须使用一些具体的分布-可能是正态分布?在任何情况下,我们都可以通过线性近似来计算电路电阻的均值和方差,则分布的确切形式无关紧要。μσ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σ
joriki

当然,在实际设计中,电阻远非独立的rvs。实际上,要进行大量工作才能使一些元素组相互跟踪(毫不奇怪,这称为“匹配”)。

1
您是否使用?我更习惯于将其写为。σ=E(XEX)2σ2

@ copper.hat:当然,您对是正确的-我不加思索地采用了问题中使用的符号。σ2
joriki

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'7

1

我警告,按照我的推理,这是一个很长的答案,但也许有人可以从我的尝试中提出更好的选择(这可能不是最佳选择)。另外,我误解了最初的OP问题,并认为它表示阻力在正态分布。无论如何,我都会留下答案,但这是一个潜在的假设。

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。这个问题类似于我们在这里找到的问题,除了现在您在eq中替换。以常数。按照与上述相同的参数,您可以发现 显然,其余的是替换已知的分布,除了一个小问题:可以从获得,注意
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,R3
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μ>>σ
内斯托尔·

嗯,我认为像正常情况那样建模高度太糟糕了,以至于我将其用作显然不正常的分布示例。我想,如果您有一群具有相同遗传背景的健康成年男性,那可能并不可怕。但是,我想从生物学家那里得知这是可以的。我经常听到的理由是,每个骨头的大小都是独立的,这完全是胡说八道。
道格拉斯·扎里

我只是意识到阻力并不是正态分布的(我可以发誓我在原始OP答案中看到了它们,但我认为那只是我的想象)。
内斯托尔·
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.