Interpolate percentages - maintain sum of rows =1, maintain all values >0

8 次查看（过去 30 天）
Edward Byers
Edward Byers 2013-4-20
Hi
I would like to interpolate data that are percentages, that change at various points in time.
For each row (year) in the data, the sum of the columns must =100%
i.e.
A B C D Sum
1990 10 40 40 10 100
1995 5 35 32 18 100
1998 0 32 44 24 100
2005 0 37 55 8 100
The data above is just an example.
I have been using interp1, but understand that there is no linkage between the columns
When I use linear interpolation, this works fine - all points are above 0, sum of each row = 100.
Using other methods:
  • cubic - all points are above 0, but the rows do not sum 100.
  • spline - some points go negative, but rows sum 100.
The problem is that the approximation oscillates around 0
What other interpolation method can I use (that isn't linear) where I can set conditions to prevent the approximation going below 0?
Thanks
  3 个评论
显示 1更早的评论
Edward Byers
Edward Byers 2013-4-21
Yes, just looking to find values in intermediate years.
It's not actually past data, am actually making assumptions about the percentage splits of fuel vectors for future projections of energy, i.e. 2010 2025 2050 2075.... i.e. what will total energy demand be in 2050 for transport if the fuel split is 25% diesel, 45% electric, 15% hydrogen.. etc and overall passenger demand is x trillion passenger kms.
I could do it linearly, but we know that these things change slowly over time, so a big change of direction between 2049 and 2050 doesn't make the graphs look so good!
Mel
Mel 2013-4-21
If you are trying to predict future fuel usage, you may want to include other factors. For example, unless they can make non-diesel generators and ship them to the Arctic fairly inexpensively, you will always have a baseline diesel amount (I know it is not a lot of people, but Canada's Arctic pretty much runs off of diesel at the moment). If you are looking at a particular area (instead of the world), you may want to consider other factors in that area.
If your rows do not sum to 100, you could always multiply them by a scalar quantity so that they do (sort of like normalizing the sum).
Maybe look at what shape your initial data has and see what form fits it best (semilinear, linear, quasilinear, quadratic, exponential, etc.) To keep data from following an exact curve (so it fluctuates from year to year), you can add in random normal noise with an appropriate std.
You could also set up a projection scheme using differential equations, such that if the usage of a particular fuel went up a lot recently it will continue to or not continue to, depending on how you set up your equations.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Matt J
Matt J 2013-4-21
编辑：Matt J 2013-4-21
Raised cosine interpolation would be a possibility. For data xi spaced apart by 1, this would be
y(x)=sum_i y_i h(x-xi)
where
h(x)= 0.5*cos(x/pi)+0.5, abs(x)<=1
= 0, otherwise
but AFAIK there is no pre-packaged MATLAB function that offers it. You would have to code it from scratch. What's wrong with linear interpolation? What do you mean by "oscillates around zero"?
  1 个评论
Matt J
Matt J 2013-4-21
For non-uniformly spaced data the extension is fairly simple. E.g., interpolating at x with x1<=x<=x2, one would do
delta= x2-x1;
weight= 0.5*cos(x/pi/delta)+0.5;
y=weight*y1+(1-weight)*y2

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Interpolation 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by