Summe einer Zeitreihe berechnen

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

eve

Forum-Anfänger


	Beiträge: 41

	Anmeldedatum: 18.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 02.09.2009, 13:12 Titel: Summe einer Zeitreihe berechnen

Hallo zusammen!
Ich würd gern einzelne Werte eines Jahres aufsummieren und anschließend eine Graphik erstellen, auf der die summierten Werte pro Jahr dargestellt sind.
Leider plottet er mir einen riesigen Wert nur für 1982...Wo ist der Fehler? Könnt ihr mit evtl. weiterhelfen?

Code:

%Einlesen
[istnr, YYYYMMDD, ns] = ...
textread('XXX.txt','%f %f %f',-1,'headerlines',1);

%Datumsumwandlung

YYYY = fix(YYYYMMDD*1E-4);
MM = fix((YYYYMMDD-YYYY*1E4)*1E-2);
DD = fix(YYYYMMDD-YYYY*1E4-MM*1E2);

sdate = datenum(YYYY,MM,DD);
date = datevec(sdate);
year = (1982:2008); %Vektor der Jahre 1982-2008 erzeugen

%Übersichtsplot
if ns == -1;
ns = NaN;
elseif ns == 9999;
ns = NaN;
else
for i = 1:length(year); %von 1982-2008
ind_year = find(date(:,1) == 1982); %anfangen bei 1982
ns_sum(i) = sum(ns(ind_year)); %die Werte des Jahres aufzusummieren
ns_sum_mm(i) = ns_sum(i)./10; %Division durch 10
end
figure ('PaperPositionMode','auto');axes; %graphik erstellen
bar(date(i),ns_sum(i));
%set(h,'BarWidth',0.4);
set(gca,'Layer','top');
box on;grid on;xlim(date);datetick('x','keeplimits');
end

Funktion ohne Link?

Die Daten sehen folgendermaßen aus:
14800 19821107 -1
14800 19821108 66
14800 19821109 148

Grüße
eve

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 04.09.2009, 10:51 Titel:

Hallo eve,

versuch diesen leicht umgeschriebenen Code mal.

Code:

%Einlesen
[istnr, YYYYMMDD, ns] = ...
textread('XXX.txt','%f %f %f',-1,'headerlines',1);

%Datumsumwandlung

YYYY = fix(YYYYMMDD*1E-4);
MM = fix((YYYYMMDD-YYYY*1E4)*1E-2);
DD = fix(YYYYMMDD-YYYY*1E4-MM*1E2);

sdate = datenum(YYYY,MM,DD);
date = datevec(sdate);
year = (1982:2008); %Vektor der Jahre 1982-2008 erzeugen

%Übersichtsplot
ii=1; % Indexvariable
if ns == -1;
ns = NaN;
elseif ns == 9999;
ns = NaN;
else
for j=1982:2008;
ind_year = find(date(:,1) == j); %anfangen bei 1982
a = sum(ns(ind_year)); % Zwischenvariable
ns_sum(ii,:) = [j a]; % Matrix der Jahressummen
ns_sum_mm=ns_sum(:,2)./10;
ii=ii+1;
end;
figure ('PaperPositionMode','auto');axes; %graphik erstellen
bar(ns_sum(:,1),ns_sum(:,2));
%set(h,'BarWidth',0.4);
set(gca,'Layer','top');
box on;grid on;xlim(date);datetick('x','keeplimits');
end

Funktion ohne Link?

Mit meinen Beipieldaten hat es problemlos funktioniert und sollte im Ergebnis hoffentlich deinen Vorstellungen entsprechen.

Ciao


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.