observations every 8 or 9 weeks-interpolation
3 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
I have the cell matrix
A={
'28/12/08' [2.3187] [1.1517]
'01/03/09' [2.2314] [1.1540]
'03/05/09' [2.2094] [1.1702]
'28/06/09' [2.2666] [1.1762]
'30/08/09' [2.2240] [1.1860]
'01/11/09' [2.2387] [1.1741]
'03/01/10' [2.3229] [1.1729]
'07/03/10' [2.2253] [1.1559]
'09/05/10' [2.2504] [1.1678]
'04/07/10' [2.2673] [1.1810]
'05/09/10' [2.1828] [1.1816]
'07/11/10' [2.1271] [1.1660]
'02/01/11' [2.3143] [1.1674]
'06/03/11' [2.1860] [1.1544]
'08/05/11' [2.1770] [1.1704]
'03/07/11' [2.1832] [1.1758]
'04/09/11' [2.1578] [1.1766]
'06/11/11' [2.1612] [1.1710]}
the observations represent a 8 or 9 week average. And I want to convert these values to estimated monthly averages via interpolation.My thought is the following
%convert the mm/yy to number dates
xi=datenum(A(:,1),'mm/yy')
B=cell2mat(A(:,[2,3]));
for c = 2:3
B(:,c-1) = interp1(xi(1:2:end),cell2mat(A(:,c)),xi);
end
But I have the feeling that my approach is not correct. Any suggestions?
thank you in advance
2 Commenti
Risposta accettata
Andrei Bobrov
il 14 Ago 2012
Modificato: Andrei Bobrov
il 14 Ago 2012
One way:
d1 = datenum(A(:,1),'dd/mm/yy');
[y,m,d] = datevec(d1(1));
d2 = datenum(y,m,d+ (0:diff(datenum(A([1,end],1),'dd/mm/yy')))');
v = arrayfun(@(x)interp1(d1,cell2mat(A(:,x)),d2),2:size(A,2),'un',0);
[a,b] = datevec(d2);
[a1,b1,c] = unique([a,b],'rows');
[x,y] = ndgrid(c,1:size(A,2)-1);
out1 = accumarray([x(:),y(:)], cat(1,v{:}),[],@(x){mean(x)});
out = [cellstr(datestr(d2(b1),'mm/yy')),out1];
second way:
d1 = datenum(A(:,1),'dd/mm/yy');
[y,m,d] = datevec(d1([1,end]));
d3 = datenum(y(1),m(1)+(1:diff([y,m])*[12;1]+1)',1)-1;
out = [cellstr(datestr(d3,'mm/yy')), cell(numel(d3),size(A,2)-1)];
for j1 = 2:size(out,2)
out(:,j1) = num2cell(interp1(d1,cell2mat(A(:,j1)),d3,'linear','extrap'));
end
etc.
0 Commenti
Più risposte (1)
Walter Roberson
il 14 Ago 2012
Note that cell2mat(A(:,c)) is going to be the same as B(:,c-1)
You need to have another look at how many x and y values you are passing into interp1()
2 Commenti
Walter Roberson
il 14 Ago 2012
xi(1:2:end) is half the length of cell2mat(A(:,c)) but interp1() requires the first two arguments to be the same length.
Vedere anche
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!