Metoda przenoszenia średnio metryczna

Używając MATLAB, w jaki sposób można znaleźć 3-dniową średnią ruchu konkretnej kolumny macierzy i dodać średnią ruchoma do tej matrycy staram się obliczyć 3-dniową średnią ruchu od dołu do góry matrycy. Mam podany mój kod: Biorąc pod uwagę następujące macierzy a i maski: próbowałem wykonania komendy conv, ale otrzymuję błąd. Oto komenda conv, którą próbowałem użyć w drugiej kolumnie macierzy a: Wyjście I pragnienie jest podane w następującej matrycy: Jeśli masz jakieś sugestie, bardzo by to doceniłem. Dziękuję W kolumnie 2 matrycy a, obliczam średnią ruchu trzydniowego w następujący sposób i umieszczając wynik w kolumnie 4 matrycy a (zmieniłam nazwę matrycy jako 39desiredOutput39 tylko dla ilustracji). Średnia 3-dniowa z 17, 14, 11 wynosi 14, 3-dniowa średnia z 14, 11, 8 jest 11, 3-dniowa średnia z 11, 8, 5 jest równa 8, a średnia 3-dniowa z 8, 2 to 5. W czwartej kolumnie nie ma wartości w dolnych dwóch wierszach, ponieważ obliczenia dla 3-dniowej średniej ruchomej zaczynają się od dołu. Wyjście 39valid39 nie będzie wyświetlane do co najmniej 17, 14 i 11. Mamy nadzieję, że to ma sens ndash Aaron 12 czerwca 13 w 1:28 Ogólnie pomogłoby, gdybyś wykazał błąd. W tym przypadku robisz dwie rzeczy niewłaściwie: najpierw trzeba podzielić rozdzielczość na trzy (lub długość średniej ruchomej). Zauważ rozmiar c. Nie możesz po prostu zmieścić się w c. Typowym sposobem na uzyskanie średniej ruchomej byłoby użycie tego samego: ale to nie wygląda tak, jak chcesz. Zamiast tego musisz zmusić się do użycia kilku linii: Utworzono w środę, 08 października 2008 20:04 Zaktualizowano w czwartek, 14 marca 2017 01:29 Wpisany przez Batuhan Osmanoglu Odsłon: 41492 Przeprowadzka Średnia W Matlabu Często znajduję się w potrzebie uśredniając dane muszę trochę zmniejszyć hałas. Napisałem kilka funkcji, aby dokładnie to, czego chcę, ale matlab wbudowane w funkcję filtru działa całkiem nieźle. Tu napisać o uśrednieniu danych 1D i 2D. Filtr 1D może być realizowany przy użyciu funkcji filtra. Funkcja filtrująca wymaga co najmniej trzech parametrów wejściowych: współczynnika licznika dla filtra (b), współczynnika mianownika dla filtra (a) oraz oczywiście oczywiście danych (X). Filtr średniej prędkości można zdefiniować po prostu: W przypadku danych 2D możemy użyć funkcji filtru Matlabs2. Aby uzyskać więcej informacji na temat działania filtru, wpisz: Oto szybka i brudna implementacja filtru średniego 16 na 16. Najpierw musimy zdefiniować filtr. Ponieważ chcemy tylko równego udziału wszystkich sąsiadów, możemy użyć tych funkcji. Dzielimy wszystko 256 (1616), ponieważ nie chcemy zmieniać ogólnego poziomu (amplitudy) sygnału. Aby zastosować filtr, po prostu powiedzmy, że poniżej Poniżej przedstawiono wyniki fazy interferogramu SAR. W tym przypadku zakres znajduje się w osi Y a azymut jest odwzorowywany na osi X. Filtr miał szerokość 4 pikseli w zakresie i szerokość 16 pikseli w Azimuth. Należy obliczać średnią ruchomej w serii danych w pętli for. Muszę uzyskać średnią ruchomej przez dziewięć dni. Tablica Im w komputerach to 4 serie 365 wartości (M), które same są wartościami średnimi innego zestawu danych. Chcę wykreślić średnie wartości moich danych ze średnią ruchoma w jednym wykresie. I googled nieco o ruchu średnich i conv polecenia i znalazłem coś, co próbowałem wdrożyć w moim kodzie .:. Więc w zasadzie obliczyć moje średnie i spróbować to z (zła) średnia ruchoma. Wybrałem wartość wts tuż przy stronie matematyki, więc jest to błędne. (źródło: mathworks. nlhelpeconmoving-average-trend-estimation. html) Mój problem jest jednak taki, że nie rozumiem, co to jest. Czy ktoś mógłby wyjaśnić, jeśli ma coś wspólnego z wagami wartości: jest to nieważne w tym przypadku. Wszystkie wartości są ważone tak samo. A jeśli robię to całkowicie złe, mogę uzyskać pomoc z nim moje najszczersze podziękowania. zapytał 23 września o godzinie 19:05 Korzystając z conv jest doskonałym sposobem na zaimplementowanie średniej ruchomej. W kodzie, którego używasz, ws jest tym, ile ważysz każdą wartość (jak się domyślasz). suma tego wektora powinna zawsze być równa. Jeśli chcesz równomiernie wyważyć każdą wartość i zrobić filtr rozmiaru N, to chcesz to zrobić Używając prawidłowego argumentu w conv będzie powodować mniejsze wartości Ms niż masz w M. Użyj tego samego, jeśli nie masz nic przeciwko efektom zero wypełnienia. Jeśli masz skrzynkę narzędziową do przetwarzania sygnałów, możesz użyć cconv, jeśli chcesz spróbować średniej ruchomej. Coś, co trzeba przeczytać w dokumentacji conv i cconv, aby uzyskać więcej informacji, jeśli już nie masz. Pobierz plik movAv. m (zobacz także movAv2 - zaktualizowana wersja umożliwiająca ważenie) Opis Matlab zawiera funkcje movavg i tsmovavg (średnia ruchoma serii czasowej) w Financial Toolbox, movAv ma na celu powielenie podstawowych funkcji tych. Kod tutaj zapewnia dobry przykład zarządzania indeksami wewnątrz pętli, które mogą być mylące na początek. Ive celowo przechowywać kod krótkie i proste, aby zachować ten proces jasne. movAv wykonuje prostą średnią ruchową, która może być wykorzystywana do odzyskiwania hałaśliwych danych w niektórych sytuacjach. Działa, przyjmując średnią wejścia (y) w oknie czasu przesuwu, którego wielkość jest określona przez n. Im większy n, tym większa jest ilość efektu wygładzania efektu n względem długości wektora wejściowego y. i skutecznie (dobrze, sortuje) tworzy filtr częstotliwości dolnoprzepustowej - patrz sekcja przykłady i rozważania. Ponieważ ilość wygładzania zapewniona przez każdą wartość n jest względna względem długości wektora wejściowego, zawsze warto testować różne wartości, aby zobaczyć, co jest właściwe. Pamiętaj również, że n punktów są tracone po każdej średniej, jeśli n wynosi 100, pierwsze 99 punktów wektora wejściowego nie zawiera wystarczających danych dla średniej wartości 100pt. Można to uniknąć poprzez ułożenie średnich, na przykład poniższy kod i wykres porównują różne średnie okna długości. Zauważ, że płynność 1010pt jest porównywana z pojedynczą średnią 20pt. W obu przypadkach 20 punktów danych są tracone w całości. Utwórz xaxis x1: 0.01: 5 Wygeneruj hałas hałasuReps 4 reprezacja szumu (randn (1, stłum. (Numel (x) noiseReps)), noiseReps, 1) zmiana szumu (hałas, hałas 1, hałas) Odtworzenie hałasu ydata yexp x) 10noise (1: długość (x)) Średnie perfrom: y2 movAv (y, 10) 10 pkt y3 movAv (y2, 10) 1010 pkt y4 movAv (y, 20) 20 pkt y5 movAv (y, 40) 40 pkt y6 movAv (y, 100) 100 pkt Opis wykresu wykresu (x, y, y2, y3, y4, y5, y6) (dane surowe, średnia ruchoma 10pt, 1010pt, 20pt, 40pt, 100pt) xlabel (x) ylabel y) tytuł (porównanie średnich ruchomej) movAv. m kod biegu wyjściowego wyjście movAv (y, n) Pierwsza linia definiuje nazwy funkcji, wejścia i wyjścia. Wejście x powinno być wektorem danych do wykonywania średniej na, n powinno być liczbą punktów do wykonywania średniej nad wyjściem będzie zawierać uśrednione dane zwracane przez funkcję. Prealokacja wyjścia wyjściowegoNaN (1, numel (y)) Znajdź punkt środkowy n midPoint round (n2) Główna praca funkcji jest wykonywana w pętli for, ale zanim rozpoczniesz dwie rzeczy. Po pierwsze, wyjście jest wstępnie alokowane jako NaN, to służyło dwóm celom. Po pierwsze, prealokacja jest ogólnie dobrą praktyką, ponieważ zmniejsza żonglowanie pamięci, co Matlab musi zrobić, po drugie, bardzo łatwo umieścić uśrednione dane jako wyjściowe w takim samym rozmiarze, jak wektora wejściowego. Oznacza to, że ten sam xaxis może być używany później dla obu, co jest wygodne do spisu, alternatywnie NaN można usunąć później w jednej linii kodu (wyjście wyjściowe (Zmienna midPoint będzie używana do wyrównywania danych wektora wyjściowego. n 10, 10 punktów zostanie utracone, ponieważ w pierwszych 9 punktach wektora wejściowego nie ma wystarczająco dużo danych, aby uzyskać średnią z 10 punktów. Jak wyjście będzie krótsze niż wejście, to musi być wyrównane prawidłowo. należy użyć tak, aby na początku i na końcu była zgubna ilość danych na początku i na końcu, a dane wejściowe są ustawione w linii z wyjściami buforów NaN utworzonych podczas wstępnego przypisywania danych wyjściowych dla długości 1: y (y) - n over (a: b) ban Oblicz średnie wyjście (amidPoint) mean (y (a: b)) end W samej pętli for, średnia jest przejęta przez każdy kolejny segment wejścia. Pętla będzie działać na a. zdefiniowane jako 1 do długości wejścia (y), minus dane, które zostaną utracone (n).Jeżeli wejście jest równe 100 punktom ng i n wynosi 10, pętla będzie działać z (a) od 1 do 90. Oznacza to, że pierwszy indeks segmentu ma być uśredniony. Drugi indeks (b) jest po prostu-1. Więc na pierwszej iteracji, a1. n10. więc b 11-1 10. Pierwsza średnia przejęła y (a: b). lub x (1:10). Średnia tego segmentu, która jest pojedynczą wartością, jest przechowywana na wyjściu w indeksie amidPoint. lub 156. Podczas drugiej iteracji, a2. b 210-1 11. więc średnia jest przejęta przez x (2:11) i zapisana na wyjściu (7). Przy ostatniej iteracji pętli dla wejścia o długości 100, a91. b 9010-1 100, więc średnia jest przejęta przez x (91: 100) i zapisana na wyjściu (95). Pozostawia to wyjście z wartością n (10) NaN w indeksie (1: 5) i (96: 100). Przykłady i rozważania Poruszanie się średnimi jest przydatne w niektórych sytuacjach, ale nie zawsze jest to najlepszy wybór. Oto dwa przykłady, w których niekoniecznie są one optymalne. Kalibracja mikrofonu Ten zestaw danych reprezentuje poziomy każdej częstotliwości generowanej przez głośnik i zarejestrowany przez mikrofon o znanej liniowej odpowiedzi. Wyjście głośnika zmienia się z częstotliwością, ale możemy poprawić tę zmianę za pomocą danych kalibracji - wyjście może być regulowane na poziomie w celu uwzględnienia fluktuacji kalibracji. Zwróć uwagę, że surowe dane są hałaśliwe - oznacza to, że niewielka zmiana częstotliwości wydaje się wymagać dużego, niepoprawnego, zmian poziomu w celu uwzględnienia. Czy jest to realistyczne Lub czy jest to produkt środowiska zapisu W tym przypadku rozsądne jest zastosowanie średniej ruchomej, która wygładzi krzywą poziomu, aby uzyskać krzywą kalibracji, która jest nieco mniej nieregularna. Ale dlaczego nie jest to optymalne w tym przykładzie? Więcej danych byłoby lepiej - wielokrotne kalibracje uśrednione razem zniszczą hałas w systemie (tak długo, jak jego losowo) i zapewniają krzywiznę z mniej subtelnym szczegółem utracone. Średnia ruchoma może tylko przybliżać tę sytuację i może usunąć niektóre z wyższych częstotliwości i szczytów z krzywej, która naprawdę istnieje. Sine waves Wykorzystanie średniej ruchomej na falach sinusoidalnych podkreśla dwa punkty: Ogólny problem polegający na wybraniu rozsądnej liczby punktów do osiągnięcia średniej. To proste, ale istnieją skuteczniejsze metody analizy sygnału niż uśrednione sygnały oscylacyjne w dziedzinie czasowej. Na tym wykresie pierwsza fala sinusoidy jest wykreślana na niebiesko. Hałas jest dodawany i wykreślany jako krzywa pomarańczowa. Średnia ruchoma jest wykonywana w różnych punktach, aby sprawdzić, czy oryginalna fala może zostać odzyskana. 5 i 10 punktów daje rozsądne wyniki, ale nie usuwaj całkowicie szumu, gdzie coraz większe punkty zaczynają tracić szczegóły amplitudy, gdy średnia rozciąga się na różnych fazach (pamiętaj, że fala oscyluje wokół zera, a średnia -1) Alternatywnym podejściem byłoby zbudowanie filtra dolnoprzepustowego niż można go zastosować do sygnału w dziedzinie częstotliwości. Nie będę wchodzić w szczegóły, ponieważ wykraczało poza zakres tego artykułu, ale ponieważ szum jest znacznie wyższy niż częstotliwość podstawowa fal, byłoby w tym przypadku dość łatwo utworzyć filtr dolnoprzepustowy niż usunie wysoką częstotliwość hałas.