- #1
brainpushups
- 453
- 197
I'm working through John Taylor's An Introduction to Error Analysis and so far this is the only problem I haven't been able to solve. I was hoping someone could lend me some insight.
The problem asks you to use error propagation to verify that the uncertainties in A and B for a line of the form y = A + Bx is given by
[tex]\sigma[/tex]A = [tex]\sigma[/tex]y Sqrt[[tex]\Sigma[/tex]x^2/∆]
[tex]\sigma[/tex]B = [tex]\sigma[/tex]y Sqrt[N/∆]
with the assumptions that there are negligible uncertainties in x and all uncertainties in y have the same magnitude.
Here we go:
Recall that the constants for a linear regression are found by
A = (([tex]\Sigma[/tex]x^2[tex]\Sigma[/tex]y-([tex]\Sigma[/tex]x[tex]\Sigma[/tex]xy))/∆
B=((N[tex]\Sigma[/tex]xy - [tex]\Sigma[/tex]x[tex]\Sigma[/tex]y))/∆
where ∆ is dependent only on the values of x and therefore does not come into play for the partial derivatives with the assumptions.
The general formula for propagation of error is
[tex]\delta[/tex]q=Sqrt[(([tex]\partial[/tex]q/[tex]\partial[/tex]x)[tex]\delta[/tex]x)^2+...+[tex]\delta[/tex]q=Sqrt[(([tex]\partial[/tex]q/[tex]\partial[/tex]z)[tex]\delta[/tex]z)^2)
Applying this formula to find the uncertainty in A yeids
[tex]\delta[/tex]A=[tex]\sigma[/tex]y/Sqrt[∆]Sqrt(([tex]\Sigma[/tex]x^2 - x1[tex]\Sigma[/tex]x)^2+...+([tex]\Sigma[/tex]x^2 - xN[tex]\Sigma[/tex]x)^2
where I have factored out ∆ and the uncertainty in y squared that are shared by every term under the radical.
And here's where I'm stuck. I don't see how this can reduce to what it is supposed to. Somehow the stuff under the radical has to reduce to the sum of x squared. I've tried writing out some terms and it doesn't seem like it does. This makes me think I've made some other error. Any ideas?
Finding the uncertainty in the slope (B) has caused similar problems. Granted, I haven't spent as much time on it, but after applying the method above to the expression for B I end up with
[tex]\delta[/tex]B=[tex]\sigma[/tex]y/Sqrt[∆] Sqrt(N(x1-[tex]\Sigma[/tex]x))^2+...+(N(xN-[tex]\Sigma[/tex]x))^2
where the stuff under the radical has to reduce to N.
Thanks for your time
The problem asks you to use error propagation to verify that the uncertainties in A and B for a line of the form y = A + Bx is given by
[tex]\sigma[/tex]A = [tex]\sigma[/tex]y Sqrt[[tex]\Sigma[/tex]x^2/∆]
[tex]\sigma[/tex]B = [tex]\sigma[/tex]y Sqrt[N/∆]
with the assumptions that there are negligible uncertainties in x and all uncertainties in y have the same magnitude.
Here we go:
Recall that the constants for a linear regression are found by
A = (([tex]\Sigma[/tex]x^2[tex]\Sigma[/tex]y-([tex]\Sigma[/tex]x[tex]\Sigma[/tex]xy))/∆
B=((N[tex]\Sigma[/tex]xy - [tex]\Sigma[/tex]x[tex]\Sigma[/tex]y))/∆
where ∆ is dependent only on the values of x and therefore does not come into play for the partial derivatives with the assumptions.
The general formula for propagation of error is
[tex]\delta[/tex]q=Sqrt[(([tex]\partial[/tex]q/[tex]\partial[/tex]x)[tex]\delta[/tex]x)^2+...+[tex]\delta[/tex]q=Sqrt[(([tex]\partial[/tex]q/[tex]\partial[/tex]z)[tex]\delta[/tex]z)^2)
Applying this formula to find the uncertainty in A yeids
[tex]\delta[/tex]A=[tex]\sigma[/tex]y/Sqrt[∆]Sqrt(([tex]\Sigma[/tex]x^2 - x1[tex]\Sigma[/tex]x)^2+...+([tex]\Sigma[/tex]x^2 - xN[tex]\Sigma[/tex]x)^2
where I have factored out ∆ and the uncertainty in y squared that are shared by every term under the radical.
And here's where I'm stuck. I don't see how this can reduce to what it is supposed to. Somehow the stuff under the radical has to reduce to the sum of x squared. I've tried writing out some terms and it doesn't seem like it does. This makes me think I've made some other error. Any ideas?
Finding the uncertainty in the slope (B) has caused similar problems. Granted, I haven't spent as much time on it, but after applying the method above to the expression for B I end up with
[tex]\delta[/tex]B=[tex]\sigma[/tex]y/Sqrt[∆] Sqrt(N(x1-[tex]\Sigma[/tex]x))^2+...+(N(xN-[tex]\Sigma[/tex]x))^2
where the stuff under the radical has to reduce to N.
Thanks for your time