differentiating the Veloclinic model for nonlinear least-squares fitting

Two posts ago I described fitting the linear CP model with an iteratively weighted least-square fit to approximate an envelope fit. The weights were either 1/t2 for points falling under the CP curve or 10k/t2 for points falling above the CP curve. This did a decent job in the example I showed of matching the fit using a 2-point method used in Golden Cheetah.

But of greater interest is using this method for my modified Veloclinic model. That model should do much better at fitting the full time spectrum of the power-duration data. The challenge there, however, is the model is nonlinear.

Nonlinear least squares fitting is described by Wolfram. In my last posts, I added weighting to the method described on the Wolfram page. Nonlinear least square fitting is a matter of navigating a hyperdimensional space, looking for the point where all error terms go to zero (where the model perfectly fits the data). There is no such point, so you search for it until improvement slows below some threshold, then declare victory. But the search involves finding the direction to zero, and that requires tracking the slope, and that requires derivatives. No problem -- it's easy to differentiate the Veloclinic model.

The model as I modified it consists of two components. First, the anaerobic component:

P1 (τ1 / t) [1 − exp(−t / τ1) ] ,

and then an aerobic component, which I modified as follows:

P2 (τ2 / t)1/2 [1 − exp(−[t/τ2]1/2) ] .

I'll call the total modeled power f, a function of time.

There's four parameters: P1, τ1, P2, and τ2. The model is clearly linear on P1 and P2. It's the τ1 and τ2 which are nonlinear.

First, the derivatives of f on P1 and P2 are trivial since the dependence is linear:

∂ f / ∂ P1 = [τ1 / t] [1 − exp(−t / τ1)]

∂ f / ∂ P2 = [τ2 / t]1/2 [1 − exp(−[t / τ2]1/2)]

The derivatives of the τ terms is a bit more complicated:

∂ f / ∂ τ1 = P1 [ (τ1 / t) (1 - exp[-t / τ1]) - exp(-t / τ1) ] / τ1

∂ f / ∂ τ2 = 1/2 P2 [ (τ2 / t)1/2 (1 - exp[-(t / τ2)1/2]) - exp(-[t / τ2]1/2) ] / τ2

I checked these numerically and made some corrections versus the originally posted derivatives.

Comments

Post a Comment

Popular posts from this blog

Marin Avenue (Berkeley)

Image
On Saturday, I got an email from Tim telling letting me know that Murphy's Spring Classics was going to start up Marin Ave. I'd missed that. Murphy announces his courses the day before, and indeed it was already well into the night before when the course announcement went out. I'd skipped over the details, jumping straight to middle and end games. Marin? I laughed at the irony, because not long before when he'd told me he was choosing to repeat the Nifty Twn Fifty which we'd done together last year versus selecting Murphy Mack's Stage Mullett two-day, I responded I didn't think I had it in me to ride Marin Ave two years in succession. Well, so much for that. I was committed, so I went downstairs to swap my 11-26 cassette for my Recon 12-27. The name is easily overlooked: "Marin Ave" fails to strike the same impression as, for example, "Redwood Gulch" or "China Grade". But among those in the know, Marin Ave is infamous i
Read more

hummingbird feeder physics

Image
I couldn't figure out how the hummingbird feeder worked. Why didn't it overflow? Forces need to balance, of course. Neglecting surface tension, there liquid level is higher in the inner reservoir than in the feeding chamber, so there must be a corresponding pressure difference. Suppose the pressure in the inner chamber were zero. Then the column height difference would need to be atmospheric pressure / (density of liquid × gravity). But this is over 9 meters! Obviously the height difference is only approximately 1% of this. So the pressure difference inside versus outside is only approximately 1%. The inside is only slightly below atmosphere. So air is getting in. How? Does it diffuse through the liquid? If this were the dominant mechanism, it wouldn't take long for the pressure inside to go from 99% to 99.3%, for example, which should be plenty to push the column of liquid down in the inside chamber and thus push liquid out through the holes. It would over
Read more

Gallium Pro geometry

Image
On Caltrain yesterday, I was surprised to see among the usual set of commuting bikes a SRAM Red-equipped Argon Gallium Pro. The bike immediately stole your attention. The design was one which has been trending: clear-coated carbon with plenty of text detail so the viewer is sure to realize he's gazing upon the result of advanced, proprietary engineering and not at yet another paint job on the same old OEM frames being pumped out of the same old Taiwanese factories. But it worked: the bottom bracket area was huge, the downtube a large-diameter "inverted Kamm tail" design seemingly designed to maximize wind resistance, the top tube a broad, eccentric shape which screamed "vertically compliant yet torsionally stiff". All it lacked were the pencil-think seat stays Cervelo popularized, but this was a machine designed for stiffness over anything else. It seemed dramatic overkill, since the bike was small (the Argon "XS" I suspect). There wasn't muc
Read more