3

I have these curves,

life curve

From this curve I can determine the life of a prop shaft due to gyroscopic forces at different yaw angles and certain speeds. I performed curve fitting on data points to get accurate high order polynomials for this interval of yaw angles. The polynomials are as follows,

y_150 = @(x) 22*((x-23)/4.9)^4 - 48*((x-23)/4.9)^3 + 27*((x-23)/4.9)^2 - 37*((x-23)/4.9) + 40;

y_200 = @(x) 11*((x-19)/4.8)^4 - 48*((x-19)/4.8)^3 + 73*((x-19)/4.8)^2 - 72*((x-19)/4.8) + 48;

y_212 = @(x) 23*((x-19)/4.8)^4 - 43*((x-19)/4.8)^3 + 22*((x-19)/4.8)^2 - 40*((x-19)/4.8) + 41;

But what about at 180 knots? Or 205 knots? Can I do some sort of 3 dimensional interpolation to account for different speeds? Since it is not considered good enough to use the closest speed value.

I would appreciate ANY ideas or comments on this problem.

hbaderts
  • 14,136
  • 4
  • 41
  • 48
Dan_space
  • 31
  • 1
  • 2
    I think this question is better suited for [stats.stackexchange.com](http://stats.stackexchange.com/), because it is not about programming but rather about how to do multivariate interpolation. However, a [quick search](http://stats.stackexchange.com/questions/11/multivariate-interpolation-approaches) shows that you might not get the answer you want. – dasdingonesin Feb 25 '16 at 08:46
  • 1
    have you tried anything yet? I am thinking of something like [griddatan](http://ch.mathworks.com/help/matlab/ref/griddatan.html) or [Scattered Data Interpolation](http://ch.mathworks.com/help/matlab/scattered-data-interpolation.html) – marco wassmer Feb 25 '16 at 11:19
  • Thank you. That is a great suggestion. I'm not sure how to use it yet. But it looks like it can interpolate higher dimensions. – Dan_space Feb 25 '16 at 11:57

1 Answers1

0

Part of the problem is you don't have estimates for the quality of the individual fits. All your curves go asymptotic, so the constant terms are likely to have large uncertainties.

The first thing I'd do is re-run the fits but cut off the data at, say, 16 or 18 degrees to avoid the asymptotic region.

Next, the zero-point values for 200 and 212 are almost certainly not identical -- your fitting coefficients look suspiciously to have been rounded to integers, which is not a good idea. If you can re-calcuate, then at least you will have three points on a curve for each Yaw angle, which makes sensible fitting along the Lifetime axis more likely to occur.

Carl Witthoft
  • 20,573
  • 9
  • 43
  • 73