Jump to content
The mkiv Supra Owners Club

Help with differentiation


martini

Recommended Posts

It's been a good five years since I did any differentiation, so I wondered if someone could help me with a relatively easy equation...

 

I need to differentiate this with respect to RHm:

 

Gc=(1+(gm^3-1)((1-RHm)RHt)/(1-RHt)RHm)^1/3

 

It's probably much easier to read here:

 

http://img.skitch.com/20100208-gr6mrqj9uygh81bitbc2xn7dqe.jpg

 

Help on this would be greatly appreciated :)

Link to comment
Share on other sites

You should change the equation really and change some of the variables with single letters to help you engineer that.

 

if y = Gcorr, a=Gmeas, x=RHmeas, b=RHtarget then you've got

 

y=(1+(a^3-1) (b/1-b)(1-x/x))^1/3

 

if you're holding a and b as constant you're looking to diff dy/dx then so it's simple differentiation

 

u=1+(a^3-1)(b/1-b)(1-x/x) to get du/dx = -(a^3-1)(b/1-b)(1/x^2) and so y=u^(1/3) you can work out dy/dx and then mutiply them together.

 

Sorry doing this without being able to draw it out... and haven't done diff since my a-level maths which is criky... about 10 years ago?

Link to comment
Share on other sites

I know we have some clever people on here but only a scientist would be able to work that one out.

 

Oh... You are a scientist :D

 

Cheeky! We need a mathematician, not an atmospheric physicist! :)

I'm rusty too but...

 

the first question is, are RHt and gm constants w.r.t. RHm?

 

After a little thinking about this, I managed to do it. I must have been between 5 and 8 years since I did anything like this... it's been a good reminder!

 

You should change the equation really and change some of the variables with single letters to help you engineer that.

 

if y = Gcorr, a=Gmeas, x=RHmeas, b=RHtarget then you've got

 

y=(1+(a^3-1) (b/1-b)(1-x/x))^1/3

 

if you're holding a and b as constant you're looking to diff dy/dx then so it's simple differentiation

 

u=1+(a^3-1)(b/1-b)(1-x/x) to get du/dx = -(a^3-1)(b/1-b)(1/x^2) and so y=u^(1/3) you can work out dy/dx and then mutiply them together.

 

Sorry doing this without being able to draw it out... and haven't done diff since my a-level maths which is criky... about 10 years ago?

 

I went away and thought about it for a while, and got the same result as you :)

 

Similar to Newton and Leibniz? ha!

 

OK, I've managed it and put it into code. Thanks for the input guys (especially caseys :thumbs:)

 

a=gfavg//GFmeas

b=awtarget//Target RH

h=awmeas//Measured RH

c=((a^3)-1)*(b/(1-b))

dGdRH=-(1/3)*((1+c*((1/h)-1))^(-2/3))*(c/(h^2))

Link to comment
Share on other sites

Cheeky! We need a mathematician, not an atmospheric physicist! :)

 

 

After a little thinking about this, I managed to do it. I must have been between 5 and 8 years since I did anything like this... it's been a good reminder!

 

 

 

I went away and thought about it for a while, and got the same result as you :)

 

Similar to Newton and Leibniz? ha!

 

OK, I've managed it and put it into code. Thanks for the input guys (especially caseys :thumbs:)

 

a=gfavg//GFmeas

b=awtarget//Target RH

h=awmeas//Measured RH

c=((a^3)-1)*(b/(1-b))

dGdRH=-(1/3)*((1+c*((1/h)-1))^(-2/3))*(c/(h^2))

 

Phew, I wondered if I was fully accurate though, as jet lagged and calculus is not my strong point (I love mechanical/physics maths a lot more).

 

I mean ask me about orbital mechanics or such and you've got me interested :)

 

Just out of curiosity what's this for? I assume something academic? What you studying?

 

Also I have failed to find the significance of 18571... bar it being a largish prime. It's not a 'weird number' is it??

Link to comment
Share on other sites

Phew, I wondered if I was fully accurate though, as jet lagged and calculus is not my strong point (I love mechanical/physics maths a lot more).

 

I mean ask me about orbital mechanics or such and you've got me interested :)

 

Just out of curiosity what's this for? I assume something academic? What you studying?

 

Also I have failed to find the significance of 18571... bar it being a largish prime. It's not a 'weird number' is it??

 

I too prefer orbital mechanics. Especially trying to calculate gravitational forces between 3 bodies ;)

 

This is in an attempt to try and attribute error bars to some measurements I've made in the field.

 

I'm doing a Ph.D in atmospheric physics, and I have taken measurements of particle water uptake from multiple field sites. At this time, there is a lot of published material on how this may change cloud behaviour and so on, but nobody has actually tried propagating measurement errors through into global models (for these instruments).

 

Thus, for my first paper (that I'm currently writing - I'll show you all when it's submitted! :yawn:) I want to try to propagate these errors in a kosher manner.

 

This is the code I've written this afternoon for a snippet of what I've been up to:

 

http://img.skitch.com/20100208-kjcm2wya5reag6s4wsw9nqm3q6.jpg

 

As I said, when I finally get my first paper published, I'll show you lot :D

 

quick refresher and it's coming back.

 

x^n => nx^(n-1)

 

Chain Rule, product rule - I'm remembering stuff (:blink:)

 

Mike

 

That's it mate!

 

Also I have failed to find the significance of 18571... bar it being a largish prime. It's not a 'weird number' is it??

 

Nothing fun unfortunately, I'll make the next sig more interesting!

Link to comment
Share on other sites

V. impressive martini. I do hope they pay you for all this malark :)

 

Funny you should say that; Government's cut funding for all research departments and I'll be out on money starting next month... :(

 

will write up the thesis and then I'll be hunting for a real job! :)

Link to comment
Share on other sites

Phew, I wondered if I was fully accurate though, as jet lagged and calculus is not my strong point (I love mechanical/physics maths a lot more).

 

 

errr... just for the record kids it's actually Calculus ;) :D

 

Yep, I'm in the short term memory loss club too... I did Maths or the first year of my Psychology degree... and then promptly dumped it in Year 2 :D :blink:

 

Ahem :)

 

Notice you took part in the Psychometrics, but is it in fact any good or a load of old tosh?

 

http://www.mkivsupra.net/vbb/showthread.php?t=158016

Link to comment
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Restore formatting

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

×
×
  • Create New...

Important Information

We have placed cookies on your device to help make this website better. You can adjust your cookie settings, otherwise we'll assume you're okay to continue. You might also be interested in our Guidelines, Privacy Policy and Terms of Use.