|
|
|
Polyfit in MonteCarlo Simulation |
|
| Sonk |
Forum-Anfänger
|
 |
Beiträge: 10
|
|
|
|
Anmeldedatum: 19.05.10
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 19.05.2010, 08:48
Titel: Polyfit in MonteCarlo Simulation
|
 |
| |
 |
 |
Hallo,
Ich versuche gerade anhand von MonteCarlo Simulationen ein statistisches Verfahren zu testen. Die erzeugte Matrix umfasst in Spaltenlänge x Stichproben mit jewels 300 Werten in 300 Zeilen. Zuerst habe ich meine Routine mit 100 MonteCarlo Simulationen, also einer 300*100 Matrix getestet, was auch reibungslos funktioniert hat. Erhöhe ich jedoch die Zahl der Simulationen auf über 100, bricht Matlab die Routine mit der folgenden Meldung ab:
??? Error using ==> polyfit at 48
X and Y vectors must be the same size.
Error in ==> MMWUTMC_Analyse at 21
yp=polyfit(time(:,1),MC_Vektor(:,i),1); %Ermittlung der
Regressionsgeraden der gesamten Stichprobe
Kann mir eventuell jemand helfen?
Die Routine sieht folgendermaßen aus:
[[/code]
[Anz_y Anz_x]=size(MC_Vektor);
Ergebnis_p=NaN(Anz_y,Anz_x);
Ergebnis_h=NaN(Anz_y,Anz_x);
for i=1:Anz_x
time=[1:length(MC_Vektor)]'; %Erzeugung eines Zeitvektors zur Elimination der Trendkomponente
yp=polyfit(time(:,1),MC_Vektor(:,i),1); %Ermittlung der Regressionsgeraden der gesamten Stichprobe
bt=yp(1,1).*time(:,1);
Time_w1=[time,[MC_Vektor(:,i)-bt]];
for j=8:Anz_y-8
x=Time_w1(1:j,2); %Bestimmung des Auswertefenste x
y=Time_w1(j+1:Anz_y,2); %Bestimmung des Auswertefenste y
[p(j,1),h(j,1)]=ranksum(x,y,0.05); %Durchführung des Wilcoxon_Rangsummentest (Verteilungsunabhängig)
end
Ergebnis_p(:,i)=[p;NaN(8,1)];
Ergebnis_h(:,i)=[h;NaN(8,1)];
end
|
|
|
|
|
|
| Thomas84 |
Forum-Meister
|
|
Beiträge: 546
|
|
|
|
Anmeldedatum: 10.02.10
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 19.05.2010, 10:37
Titel:
|
|
statt
folgendes versuchen:
|
|
|
|
| Sonk |
Themenstarter
Forum-Anfänger
|
|
Beiträge: 10
|
|
|
|
Anmeldedatum: 19.05.10
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
|
Verfasst am: 19.05.2010, 11:01
Titel:
|
|
Dankeschön,hab den Fehler eben auch gefunden. 
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
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 - 2025
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.
|
|
