Error propogation in trig functions

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
10 replies · 4K views
richyw
Messages
179
Reaction score
0

Homework Statement



I have to use this in my calculation [tex]\theta=\tan^{-1}\left(\frac{19 \pm 1}{47 \pm 1}\right)[/tex] where both are in mm. How would I get this into [itex]\theta\pm \text{error}[/itex]?

Homework Equations



shown above

The Attempt at a Solution



looked through my lab manual, it wasn't there.
 
Physics news on Phys.org
ok I found this online

We can express the uncertainty in R for general functions of one or more observables. If R is a function of X and Y, written as R(X,Y), then the uncertainty in R is obtained by taking the partial derivatives of R with respect to each variable, multiplication with the uncertainty in that variable, and addition of these individual terms in quadrature.

[tex]R=R(X,Y,...)[/tex]
[tex]\delta R=\sqrt{\left(\frac{\partial R}{\partial X} \delta X\right)^2+\left(\frac{\partial R}{\partial Y}\delta Y\right)^2+...}[/tex]
 
so if I say that the numerator there is x and the denominator is L, then I wind up with
[tex]\delta\theta=\frac{L\delta x + x\delta L}{L^2+x^2}[/tex]

would someone be able to let me know if this is correct?
 
oops I made a mistake in the algebra. the formula I am getting is

[tex]\sqrt{ \left(\frac{L}{x^2+L^2}\delta x\right)^2+\left(-\frac{x}{x^2+L^2}\delta L\right)^2}[/tex]

I'm getting a very small error though. less than the number of significant digits. this seems a bit off...
 
richyw said:
oops I made a mistake in the algebra. the formula I am getting is

[tex]\sqrt{ \left(\frac{L}{x^2+L^2}\delta x\right)^2+\left(-\frac{x}{x^2+L^2}\delta L\right)^2}[/tex]

I'm getting a very small error though. less than the number of significant digits. this seems a bit off...

That looks good to me. Are you expressing θ in radians?
 
nope, it needs to be expressed in degrees. That simplfies to [tex]\frac{\sqrt{(L\delta x)^2 +(x\delta L)^2}}{L^2+x^2}[/tex] but it's giving me a value lower than the number of sig figs I have. what do I do now? haha
 
is this my relative error? it comes out unitless...
 
bumping this. I'm really confused as to what this error means physically.
 
First, a 'soapbox' about the quadrature method. It's perfectly fine and appropriate if the input uncertainties represent some number of standard deviations (the same number in each case) of roughly normally and independently distributed measurements. The answer then represents the same number of s.d. of the output.
But it's not unusual in engineering and in lab work that the input uncertainties would be better modeled as hard limits of uniform distributions. If there are many such inputs, the best approach is probably to turn each of those into a specific number of standard deviations, use quadrature, and understand that whatever number of s.d. you chose applies to the answer. But with only two inputs, you might be better served simply considering the four extreme values generated by plugging in the extremes of the two input ranges. The resulting range won't represent a uniform distribution, more triangular perhaps, but it will give you hard limits, and in engineering that may be more relevant.
Wrt your working so far, the formula you used for derivative of arctan assumes the angle is in radians. If it's in degrees you'll need to make a correction for that. It should be dimensionless - because angles are dimensionless.
 
richyw said:
the formula I am getting is

[tex]\sqrt{ \left(\frac{L}{x^2+L^2}\delta x\right)^2+\left(-\frac{x}{x^2+L^2}\delta L\right)^2}[/tex]

I'm getting a very small error though. less than the number of significant digits. this seems a bit off...

The expression you derived is only valid if θ is measured in radians. When using the standard calculus formulas for derivatives of trig functions, the angles are assumed to be in radians. For example d(tanθ)/dθ = sec2θ assumes θ is in radians. So, your expression gives the error in θ measured in radians. You can then convert to degrees if you wish.
 
Last edited: