並列の抵抗の分散


10

抵抗Rのセットがあり、そのすべてが平均μと分散σで分布しているとします。

次のレイアウトの回路のセクションを考えてみましょう:(r)|| (r + r)|| (r + r + r)。各部品の等価抵抗は、r、2r、および3rです。各セクションの分散は次のようになりσ22σ23σ2

回路全体の抵抗の変動はどのくらいですか?

数百万点をサンプリングした後、分散は約.10286 \ sigma ^ 2であることがわかりました.10286σ2

この結論に分析的にどのように到達するのでしょうか?

編集:抵抗値は、いくつかの平均抵抗rと分散σ^ 2で正規分布していると想定されていますσ2


1
これが最初から適切なモデルであるとは思いません。熱回路ノイズのナイキストジョンソン理論を知っていますか?意図的に何か違うことをしているのであれば、その動機を見るのは興味深いでしょう。それ以外の場合は、より標準的なモデルを検討する価値があります。:)
枢機卿

ええ、私が答えへの私の試みを書いている間に、私はモデルがポーズをとられたので明らかに扱いにくいことにも気づきました。しかし、これは実際的な問題というよりはむしろ学術的な問題に近いと思いました(結局、彼らはシミュレーションを行っています)。
ネストル2012

分散としてsigmaを使用したことに対する私の謝罪、私は元々VARを使用し、誰かがそれをsigmaに編集しました。
lrAndroid 2012

更新していただきありがとうございます。あなたが質問に少しだけ追加したいのであれば、私はこの質問の背後にある動機にまだ興味があります。:)
枢機卿

回答:


9

回路全体の等価抵抗は、解き 一つは、前提と、いくつかの独立したランダム変数の、中心と分散を有する。1RRは、私は=μ+σ

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

さらに指示がないと、の分散を計算できないため、さらに進むために、のレジームを検討し 次に、 したがって 一つは、見ている さらに、 したがって、限界σ « μ1R

σμ.
1
1Ri=1iμσμ2Ziii+higher order terms,
= 3 Σ iは=11
1R=aμσμ2Z+higher order terms,
EZ=0
a=i=131i=116,Z=i=13Ziii.
R=μ
E(Z)=0,E(Z2)=b,b=i=131i3=251216.
σ0ERμ
R=μaσa2Z+higher order terms,
σ0、 および これらの漸近形とは、並列の任意の数の抵抗に一般化できます。それぞれは、直列の基本抵抗の結果であり、基本抵抗は独立しており、それぞれ平均および分散。次に、、 ここで ヴァーRσ2B
E(R)μa=611μ,
ERヴァーRNIμσ2σ0ERμ
Var(R)σ2ba4=σ2(611)4251216=σ20.10286
E(R)Var(R)niμσ2σ0a=i1
E(R)μa,σ2Var(R)ba4,
a=i1ni,b=i1ni3.

8

正確な答えはとだけに依存するとは思いません。あなたがサンプリングしたとき、私はあなたがいくつかの具体的な分布-おそらく正規分布を使用したに違いないと思いますか?いずれの場合でも、回路の抵抗の平均と分散を線形近似で計算できます。その場合、分布の正確な形式は関係ありません。σ 2μσ2

回路の抵抗はです。線形近似では、平均および分散を持つ確率変数の逆数の平均および分散は、それぞれおよびです。したがって、平均が、および分散と分散、とはそれぞれ、平均してと分散 μσ21/μσ2/μ41/μ1/2μ1/3μσ2/μ4σ2/8μ4σ2/(R11+R21+R31)1μσ21/μσ2/μ41/μ1/(2μ)1/(3μ)σ2/μ4σ2/(8μ4)11σ2/(27μ4)251116/μ6251216σ2/μ4。次に、その逆数を取ると、の平均と分散が得られます、結果と一致しています。251611μ(251216σ2/μ4)/(116/μ)4=150614641σ20.10286σ2


もちろん、これは抵抗が独立したランダム変数であることを前提としています。

@ロバート:はい(むしろ抵抗)。これは、質問の分散、、およびの計算ですでに想定されており、物理的に理にかなっています(ただし、同じ生産バッチからすべての抵抗を取得すると、それらの抵抗はある程度相関します) )。2 σ 3 σσ2σ3σ
joriki

もちろん、実際の設計では、抵抗は独立したrvとはかけ離れています。実際、いくつかの要素グループを互いに追跡させるために、多くの作業がレイアウトに組み込まれています(当然、「マッチング」と呼ばれます)。

1
あなたは使用している?私はこれを書くことに慣れています。σ 2σ=E(XEX)2σ2

@ copper.hat:もちろん、あなたはについてはまったく正しいです-私は質問で使用された表記法を考えずに採用しました。σ2
joriki '19

5

これは、抵抗の分布の形状に依存します。制約があるとは思いますが、分布を知らなければ平均抵抗すら言えません。

それで、扱いやすい分布を選びましょう:を1つの抵抗の抵抗の標準偏差とします。抵抗をとし、各符号は確率発生します。これにより、考慮すべきケース、またはいくつかのケースを組み合わせる場合はます。もちろん、抵抗は独立していると仮定します。μ ± S 1 / 2 2 6 = 64 2 × 3 × 4 = 24sμ±s1/226=642×3×4=24

およびを選択した場合、平均は(よりわずかに低い)で、分散はです。とを選択すると、分散はます。s = 1 54.543291 100 × 6μ=100s=154.543291 0.102864μ=5、S=10.103693100×6110.102864μ=5s=10.103693

以下は、平均がで分散が場合の分散間の比率のべき級数展開です。。場合小さい場合、支配的な用語である。x 15061xx1506150614641+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

私がそれを推論したように、これは長い答えですが、誰かが私の試みから始めるより良い何かを思いつくかもしれないと警告します(これは最適ではないかもしれません)。また、私は元のOPの質問を誤解し、抵抗が正常に分布していると言っていたと思いました。とにかく答えは残しておきますが、それは根本的な仮定です。

1.問題の物理的な理由

私の推論は次のとおりです。並列の抵抗の場合、等価抵抗は次ので与えられることを思い出してください。Req

Req1=iN1Ri,

ここで、は回路の各部分の抵抗です。あなたの場合、これは私たちに与えますRi

Req=(1R1+1R2+1R3)1,   ()
ここでは抵抗が1つの回路の部分であり、したがって平均および分散正規分布を持ち、同じ理由では2つの抵抗を持つ回路の部分の等価抵抗、そして最後に、は、3つの抵抗を持つ回路の部分の等価抵抗です。分布を見つけ、そこからその分散を取得する必要があります。R1μσ2R2N(2μ,2σ2)R3N(3μ,3σ2)Req

2.分布を取得するReq

分布を求めるための一つの方法は、注目することによってである: ここから、 (これはベイズの定理によって得られた)、、、および(物理的にもっともらしい)の間の独立性は、 これをで置き換え、3つの抵抗間の独立性の別の結果は、ことに注意してください

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,R3=1/R2B=1/R332
p(Req|R2,R3)=1[1Req(a+b)]212πσ2exp(X(Req)μ2σ2),
ここでおよびです。問題は、これが式積分を解くために分析的に扱いやすいかどうかわからないということです。これにより、式の結果を置き換えることによって問題を解決することができます。少なくとも私にとっては、この夜はそうではありません。a=1/R2b=1/R3(3)(2)

抵抗が負にならない場合でも、正規分布を想定していますか?私の推測では、これにより回路の分散が発散します。
ダグラスザレ2012

1
私も知っています。私も困惑しましたが、実際にはと値に本当に依存しています。場合はと、我々はモデルを「保存」することができます。通常の状態では、抵抗のばらつきはそれほど大きくないため、最後の仮定は明らかに満たされます。これは、人々が高さを通常のランダム変数としてモデル化したときに最初に私を悩ませたものでもありましたが、私がここで与えたのと同じ理由で、Stack-exchangeの何人かの人々が私にそれを受け入れました:-) σ 2 μ > > 0 μ > > σμσ2μ>>0μ>>σ
ネストル2012

うーん、通常の高さのモデリングは非常に悪いので、明らかに通常ではない分布の例として使用します。同じ遺伝的背景を持つ健康な成人男性の集団があれば、それは恐ろしいことではないでしょう。しかし、これは大丈夫だと生物学者に聞きたいです。各骨のサイズが独立していると私が頻繁に聞いた理由はまったくナンセンスです。
ダグラスザレ2012

抵抗が通常は分散されていないことに気づきました(元のOPの回答のどこで抵抗を読んだかを確信できましたが、それは私の想像でした)。
ネストル
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.