Gemeinsamen Weg zur Generierung von endlichen geometrischen Reihe in MATLAB

Angenommen ich habe eine Zahl a, und ich will Vektor [ 1 , a , a^2 , ... , a^N ]. Ich benutze [ 1 , cumprod( a * ones( 1 , N - 1 ) ) ] code. Was ist die beste (und vermutlich effizienter) Weg, es zu tun?

InformationsquelleAutor UncleAli | 2011-06-10
  1. 13

    Was a.^[0:N] ?

    • Ich fürchte, es müssen O(N^2) Produkte, während cumprod braucht nur O(N) Produkte
    • Definitiv die kompakte Art
    • Versuchen profiling beide Methoden. Aber ja, es gebunden ist, langsamer zu sein.
    • Mein Weg ist ca. 2x schneller. Dennoch, es ist akzeptabel, weil der große Lesbarkeit zu verbessern. Danke.
    • Ein paar Dinge zu beachten... 1) ich würde nicht erwarten, dass a.^(0:N) O(N^2) Abhängigkeit, wahrscheinlich mehr wie O(NlogN), da Stell ich mir die Potenzierung erfolgt mit ein effizienter Algorithmus. 2) seien Sie nicht zu überrascht, wenn die Ergebnisse der beiden unterschiedlichen Methoden nicht *genau gleich bei der Verwendung von floating-point-Werte für a. Da Sie mit zwei unterschiedlichen zugrunde liegenden algorithmen, mit verschiedenen Aufträgen von mathematischen Operationen, sehen Sie möglicherweise kleine Unterschiede auf die Bestellung der floating-point-Präzision.
    InformationsquelleAutor ThibThib
  2. 2

    ThibThib s Antwort ist absolut korrekt, aber es nicht verallgemeinern sehr leicht, wenn a passiert mit einem Vektor. So als Startpunkt:

    > a= 2
a =  2
> n= 3
n =  3
> a.^[0: n]
ans =
   1   2   4   8

    Konnte man nun nutzen auch die integrierte Funktion vander (obwohl die Reihenfolge ist anders, aber das ist leicht behoben, falls erforderlich), zu produzieren:

    > vander(a, n+ 1)
ans =
   8   4   2   1

    Sowie mit Vektor geschätzt a:

    > a= [2; 3; 4];
> vander(a, n+ 1)
ans =
   8    4    2    1
  27    9    3    1
  64   16    4    1
    InformationsquelleAutor eat
Schreibe einen Kommentar