|
|
Vektor mit zufälligem Offset elementweise erweitern |
|
divB |

Forum-Anfänger
|
 |
Beiträge: 38
|
 |
|
 |
Anmeldedatum: 23.10.09
|
 |
|
 |
Wohnort: ---
|
 |
|
 |
Version: ---
|
 |
|
|
 |
|
Verfasst am: 16.06.2011, 22:35
Titel: Vektor mit zufälligem Offset elementweise erweitern
|
 |
Hallo,
Ich hab eine ziemlich einfache Aufgabe: Ich habe einen Datenvektor wie z.B.:
Dieser Vektor soll nun elementweise um einen Faktor R erweitert werden (also nicht repmat), z.B.: R=4:
Das habe ich so implementiert:
Allerdings ist das ziemlich langsam.
Doch ist das leider noch nicht alles. Jetzt soll der Vektor mit einem normalverteiltem Offset erweitert werden "jitter". Zum Beispiel mit sigma²=0.1 könnten die Offsets sein: [ 0 -1 0 1]. Das obige Beispiel soll dann expandieren zu:
Derzeit verwende ich:
Doch das ist unbrauchbar: Mit einem x von 19000 Elementen und R=100 benötigt ein Aufruf auf einem modernen Dual Xeon mit 4 Kernen mit 3.4Ghz und 8GB RAM viele Minuten.
Wie könnte man das (brauchbar) beschleunigen?
LG
divB
|
|
|
|
|
_Peter_ |

Moderator
|
 |
Beiträge: 537
|
 |
|
 |
Anmeldedatum: 08.12.10
|
 |
|
 |
Wohnort: ---
|
 |
|
 |
Version: 7.10, 2010a
|
 |
|
|
 |
|
Verfasst am: 17.06.2011, 13:00
Titel:
|
 |
Hallo divB,
für den ersten teil würde ich folgende Alternatvie vorschlagen:
Vielleicht hilft dir das schon.
_________________
Gruß
Peter
_________________
goMatlab-Knigge - dran gehalten?!
Schon in den FAQ gesucht? Oder der MATLAB Hilfe?
Ist vielleicht bei den Skripten oder den Tutorials was für dich dabei?
|
|
|
Jan S |

Moderator
|
 |
Beiträge: 11.057
|
 |
|
 |
Anmeldedatum: 08.07.10
|
 |
|
 |
Wohnort: Heidelberg
|
 |
|
 |
Version: 2009a, 2016b
|
 |
|
|
 |
|
Verfasst am: 18.06.2011, 23:29
Titel: Re: Vektor mit zufälligem Offset elementweise erweitern
|
 |
Hallo divB,
Das zeitraubende an deinem Programm ist das wiederholte Anwachsen eines Vectors. In jeder Iteration muss ein neuer Vektor allociert werden und die Werte des alten werden kopiert. Insgesamt werden also etwa sum(1:100:19000*100) = 1.8*10^10 Bytes temporär angefordert und belegt. Da hilft auch ein Gigahertz-Multicore-Prozessor nicht viel.
Zunächst eine gleichmäßige Verteilung:
;
Der REPMAT-Ansatz von Peter ist aber effizienter. Vielleicht noch eine Vereinfachung:
Nun zum Jitter: Es ist viel effizienter den Ausgabe-Vektor zu Beginn zu allocieren und dann die Werte einzutragen:
Ich schätze es wäre noch schneller, wenn man einen Index-Vektor erzeugt:
Leider kann ich das gerade nicht testen, weil auf meinem Netz-Rechner kein Matlab läuft. Also: Viel Erfolg beim Debuggen!
Gruß, Jan
|
|
|
divB |
Themenstarter

Forum-Anfänger
|
 |
Beiträge: 38
|
 |
|
 |
Anmeldedatum: 23.10.09
|
 |
|
 |
Wohnort: ---
|
 |
|
 |
Version: ---
|
 |
|
|
 |
|
Verfasst am: 22.06.2011, 20:15
Titel:
|
 |
Super, vielen Dank, das hilft mir sehr weiter!
LG
divB
|
|
|
|
|
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.
|
|