Variance des résistances en parallèle

10

Supposons que vous ayez un ensemble de résistances R, qui sont toutes distribuées avec la moyenne μ et la variance σ.

Considérons une section d'un circuit avec la disposition suivante: (r) || (r + r) || (r + r + r). La résistance équivalente de chaque partie est r, 2r et 3r. La variance de chaque section serait alors σ2 , 2σ2 , 3σ2 .

Quelle est la variance de la résistance de l'ensemble du circuit?

Après avoir échantillonné plusieurs millions de points, nous avons constaté que la variance est d'environ .10286σ2 .

Comment arriverions-nous à cette conclusion analytiquement?

Edit: Les valeurs de résistance sont supposées être normalement distribuées avec une résistance moyenne r et une variance σ2 .

lrAndroid
la source
1
Je ne suis pas convaincu que ce soit un modèle approprié pour commencer. Connaissez-vous la théorie de Nyquist-Johnson sur le bruit des circuits thermiques? Si vous faites délibérément quelque chose de différent, il serait intéressant de voir la motivation. Sinon, il pourrait être utile d'envisager un modèle plus standard. :)
cardinal
Oui, alors que j'écrivais ma tentative de réponse, j'ai également réalisé que le modèle n'était apparemment pas maniable comme il était posé. Cependant, je pensais que cela ressemblait plus à un problème académique qu'à un problème pratique (ils font des simulations, après tout).
Néstor
Mes excuses pour avoir sigma comme variance, j'ai utilisé à l'origine VAR et quelqu'un l'a édité en sigma.
lrAndroid
Merci pour la mise à jour. Je suis toujours intéressé par la motivation derrière cette question, si vous êtes prêt à ajouter un peu à ce sujet à votre question. :)
cardinal

Réponses:

9

La résistance équivalente de l'ensemble du circuit résout On suppose que , pour certaines variables aléatoires indépendantes , centrées et de variance .R

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

Sans autre indication, on ne peut pas calculer la variance de , donc, pour aller plus loin, on considère le régime où Ensuite, donc où On voit que De plus, Ainsi, dans la limiteR

σμ.
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, et Ces asymptotiques de et peuvent être généralisés à n'importe quel nombre de résistances en parallèle, chacune étant le résultat de résistances élémentaires en série, les résistances élémentaires étant indépendantes et chacune ayant une moyenne et une variance . Ensuite, lorsque , où
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.
Fait
la source
8

Je ne pense pas que la réponse exacte dépende uniquement de et . Lorsque vous avez échantillonné, je suppose que vous devez avoir utilisé une distribution concrète - probablement une distribution normale? Dans tous les cas, nous pouvons calculer la moyenne et la variance de la résistance du circuit en approximation linéaire, puis la forme exacte de la distribution est sans importance.μσ2

La résistance du circuit est . En approximation linéaire, la moyenne et la variance de l'inverse d'une variable aléatoire avec la moyenne et la variance sont respectivement et . Nous avons donc une somme de termes avec les moyennes , et et les variances , et , respectivement, ce qui donne une moyenne de et une variance de(R11+R21+R31)1μσ21/μσ2/μ41/μ1/(2μ)1/(3μ)σ2/μ4σ2/(8μ4)σ2/(27μ4)116/μ251216σ2/μ4. La prise de l'inverse de cela donne une moyenne de et une variance de , en accord avec votre résultat.611μ(251216σ2/μ4)/(116/μ)4=150614641σ20.10286σ2

joriki
la source
Ceci suppose, bien entendu, que les résistances sont des variables aléatoires indépendantes.
@Robert: Oui (les résistances, plutôt). Cela était déjà supposé dans le calcul des variances , et dans la question, et cela a un sens physique (bien que si nous prenons toutes les résistances du même lot de production, leurs résistances seront quelque peu corrélées ). σ2σ3σ
joriki
Dans une conception réelle, bien sûr, les résistances sont loin d'être des RV indépendants. En fait, beaucoup de travail est consacré à la mise en page pour que certains groupes d'éléments se suivent (ce qui est sans surprise appelé «correspondance»).
1
Utilisez-vous ? Je suis plus habitué à voir ceci écrit comme . σ=E(XEX)2σ2
@ copper.hat: Vous avez tout à fait raison au sujet de , bien sûr - j'avais adopté la notation utilisée dans la question sans réfléchir. σ2
joriki
5

Cela dépend de la forme de la distribution de la résistance. Sans connaître la répartition, je ne peux même pas dire la résistance moyenne, même si je pense qu'il y a des contraintes.

Alors, prenons une distribution qui est traitable: Soit l'écart type de la résistance d'une résistance. Soit la résistance , chaque signe apparaissant avec une probabilité . Cela nous donne cas à considérer, ou si nous combinons certains cas. Bien sûr, nous supposerons que les résistances sont indépendantes.sμ±s1/226=642×3×4=24

Si nous choisissons et la moyenne est (légèrement inférieure à ), et la variance est de . Si nous choisissons et , alors la variance est de .μ=100s=154.543291100×6110.102864μ=5s=10.103693

Voici une extension de la série de puissances pour les rapports entre les variances lorsque la moyenne est et la variance est : . Lorsque est petit, le terme dominant est .1x150614641+360001771561x+21801619487171x2+O(x3)x150614641=0.102862

Bien que la question que vous posiez dépende techniquement de la distribution, vous êtes probablement intéressé par des situations où l'écart-type est faible par rapport à la moyenne, et je pense qu'il existe une limite bien définie qui ne dépend pas de la distribution. Linéariser la dépendance de la résistance du circuit en fonction des résistances de chaque pièce:

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

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

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

Avec ce circuit spécifique, les dérivées partielles mises à l'échelle sont et36121,9121,9121,4121,4121,4121

(36121)2+2(9121)2+3(4121)2=150614641=0.102862
Douglas Zare
la source
1
Cela me rappelle le théorème delta multivarié, c'est-à-dire que ont respectivement la moyenne et la variance , puis devrait avoir une variance asymptotique comme , où et . La réponse finale est la même que @Douglas Zare et OP, soit 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

Je préviens que, comme je l'ai raisonné, c'est une longue réponse , mais peut-être que quelqu'un peut trouver quelque chose de mieux à partir de ma tentative (qui n'est peut-être pas optimale). De plus, j'ai mal lu la question OP d'origine et j'ai pensé qu'elle disait que les résistances étaient normalement distribuées. Je vais quand même laisser la réponse, mais c'est une supposition sous-jacente.

1. Raisonnement physique du problème

Mon raisonnement est le suivant: rappelons que, pour les résistances qui sont en parallèle, la résistance équivalente est donnée par:Req

Req1=iN1Ri,

où sont les résistances de chaque partie du circuit. Dans votre cas, cela nous donneRi

Req=(1R1+1R2+1R3)1,   ()
où est la partie du circuit avec 1 résistance, et a donc une distribution normale avec la moyenne et la variance , et par le même raisonnement est le la résistance équivalente de la partie du circuit à deux résistances et, enfin, est la résistance équivalente de la partie du circuit à trois résistances. Vous devez trouver la distribution de et en obtenir la variance.R1μσ2R2N(2μ,2σ2)R3N(3μ,3σ2)Req

2. Obtention de la distribution deReq

Une façon de trouver la distribution est de noter que: À partir d'ici, nous notons également que nous pouvons écrire (qui a été obtenu via le théorème de Bayes), qui, en supposant l'indépendance entre , et (qui est physiquement plausible), peut s'écrire Remplacer cela dans et noter qu'une autre conséquence de l'indépendance entre les trois résistances est que

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)on obtient: Notre dernier problème est alors de trouver , c'est-à-dire la distribution de la rv . Ce problème est analogue à celui que nous avons trouvé ici, sauf que maintenant vous remplacez dans l'eq. par une constante, disons . En suivant les mêmes arguments que ci-dessus, vous pouvez trouver que Apparemment, le reste est remplacement des distributions connues, sauf pour un petit problème: la distribution de peut être obtenue à partir de en notant que
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 est gaussien, donc, vous devez essentiellement trouver la distribution de la variable aléatoire où et sont des constantes, et est gaussien avec une moyenne et une variance . Si mes calculs sont corrects, cette distribution est: où, donc la distribution de serait
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),
où et . Le fait est que je ne sais pas si cela est traitable analytiquement afin de résoudre l'intégrale dans l'équation , ce qui nous conduira alors à résoudre le problème en remplaçant son résultat dans l'équation . Du moins pour moi à cette heure de la nuit, ce n'est pas le cas.a=1/R2b=1/R3(3)(2)
Néstor
la source
Vous supposez une distribution normale, même si la résistance ne peut pas être négative? Je suppose que cela fera diverger la variance du circuit.
Douglas Zare
1
Je sais, cela m'a aussi gêné, mais en pratique cela dépend vraiment des valeurs de et . Si et , alors nous pouvons "enregistrer" le modèle. Dans des conditions normales, la dispersion d'une résistance n'est pas très élevée, donc la dernière hypothèse est clairement remplie. C'était quelque chose qui me dérangeait au départ aussi lorsque les gens modélisaient la taille comme une variable aléatoire normale, mais pour la même raison que je l'ai donnée ici, certaines personnes ici à Stack-exchange m'ont fait me sentir d'accord avec ça :-). μσ2μ>>0μ>>σ
Néstor
Hmm, je pense que la modélisation de la hauteur comme normale est si mauvaise que je l'utilise comme exemple d'une distribution qui n'est évidemment pas normale. Je suppose que ce ne serait pas terrible si vous avez une population d'hommes adultes en bonne santé de la même origine génétique. Cependant, je voudrais entendre un biologiste que ce n'est pas grave. Le raisonnement que j'ai trop souvent entendu selon lequel la taille de chaque os est indépendante est un non-sens total.
Douglas Zare
Je viens de réaliser que les résistances n'étaient pas normalement distribuées (je pourrais jurer que je les ai lues sur la réponse originale des OP, mais je pense que c'était juste mon imagination).
Néstor