No i bardzo ładnie wąsy, tylko pamiętaj, że jeśli sygnał ma N próbek, to widmo ma N/2 + 1 próbek, to co jest od połowy to lustrzane odbicie widma. Nie wiem czemu liczysz np długość y skoro y definiujesz potem. Daj na początku programu clear all; W parametrze fft nie musisz nic podawać (ani tym bardziej nic przeliczać na ten cel). Abs liczysz trzy razy...
A próbowałeś sobie odsłuchać potem efektu swojego filtru z użyciem polecenia soundsc(z)? Bo niestety, ale jest tylko jakiś szum, czyli filtr działa źle. Chyba, że ja coś niepoprawnie robię. Po ustaleniu parametrów filtra w oknie filterDesigner brałem edit>convert to single section, a potem file>export(coefficients). Dzięki temu miałem w workspace Num...
Udało mi się to zrobić fft dla fragmentu jednego z twoich przebiegów (wziąłem 1000 próbek z pierwszej kolumny - czyli A1, wprowadziłem do matlaba za pomocą xlsread oraz zrobiłem fft w wyniku czego otrzymałem nasępujące widmo http://obrazki.elektroda.net/33_12494691... daje się zauważyć piki co 50 Hz - być może w innych przebiegach będą wyraźniejsze....
ad 1) wygrzebalem wzor na szereg Fouriera sygnalu x=A|sin(ωt)| gdzie ω=(2*pi)/T czyli ustalona pulsacja (no i pi=3.1415... :-) ), wyglada on tak: X= (2A)/pi - ((4A)/pi) * ∑ Dodano po 2 2) Obliczyc i narysowac widmo ampl. i faz. sygnalu s(t)=sin(t), gdzie T/3 <= t <= (3T)/8 (1/24 okresu), Dla t z poza przedzialu to 0 oczywiscie...
Genario , jeżeli chodzi o matlaba, to nie za wiele pomogę, bo nie używam. :( A z tym skalowaniem to chodzi o coś takiego, że nie znasz częstotliwości próbkowania, ale mimo wszystko możesz przedstawić na skali częstotliwości poszczególnych prążków, właśnie w postaci znormalizowanej. mając widmo sygnału (rzeczywistego) o N próbkach składające się z N/2+1...
Witam, mam problem z widmem sygnału z matlabie. Środek wykresu jest pusty. Ma ktoś pomysł co może być powodem?
Możesz wykorzystać do tego celu szybką transformatę Fouriera w celu wyznaczenia częstotliwości wchodzących w skład badanego sygnału, a następnie porównujesz widmo częstotliwości występujących w sygnale z częstotliwościami składowymi danego kodu DTMF ( Dodano po 17 Znalazłem swój stary program z przetwarzania sygnałów: http://obrazki.elektroda.pl/2043492400_1...
stem(fplot(1:length(fplot)/10), widmo(1:length(fplot)/10)); Cokolwiek autor ma na myśli to działająca wersja będzie taka: stem(t, widmo); Tylko, że zamiast t należy stworzyć macierz częstotliwości. Niech help będzie Twoim przyjacielem: http://www.mathworks.com/help/matlab/ref... NFFT = 2^nextpow2(L); % Next power of 2 from length of y...
No ja mam taki większy program - ale idea jest ta sama - podaje się okres sygnału, amplitudę, typ (sinus np), liczbe okresów, od do ile pokazać : A= 3; %amplituda sygnału lo= 3; %liczba okresów P= 0; %zacznij na czas=0 TE = 6; %zakończ na czas = 6 lp=112; %liczba próbek o=3; %czas trwania okresu sygnalu w = pi*lo/TE; plp=TE/(lp-1); %obliczenie co jaki...
Zagladam, zagladam, tylko... brak natchnienia przez ostatni tydzien. Co do FFT w konsoli... własnie w tym jest problem, że trzeba przeskalowac os czestotliwosci... Matlab wyrysuje prążki, ale liczy sobie od 0 i ma sie to nijak do czestotliwosci. Dlatego wlasnie lepiej robic w PowerGUI, bo tam sie podaje Fundamental freq i matlab sam 'umiejscowi' prążki...
Bardzo nie precyzyjne pytanie, gdyż podejrzewam że chcesz mieć wykres amplitudy z widma częstotliwościowego? Matematyka pozwala na rozwiązania ciągłe, matlab już na to nie pozwala i trzeba stosować metody dyskretne, więc nie za bardzo da się coś policzyć w -inf i +inf, trzeba określić rzeczywisty zakres i ustalić krok obliczenia. Nie do końca wiem dlaczego...
Widmo rand nie zawiera tylko jeden prążek. Ta częstotliwość 0Hz po prostu dominuje nad pozostałymi. Funkcja rand generuje wartości od 0 do 1.0. FFT rozkłada ten sygnał na sinusoidy + składowa stała, która wynosi około 0.5 i widać ją wyraźnie jako prążek w zerze i jest to normalne zachowanie. Reszta to inne częstotliwości, które FFT wydedukował z podanego...
Witam.Jestem początkującym użytkownikiem języka Matlab i mam pewien problem z którym nie mogę sobie poradzić.Na projekt mam zadanie aby wygenerować sygnał sonaru o obwiedni prostokątnej tak by można było zmieniać jego parametry,dodać szum, obliczyć sygnał odebrany (oczywiście wszystko ma być zobrazowane i wykresy mają wyglądać no wiadomo "książkowo")...
witam mam problem musze napisać program w języku Matlab filtrujący zadanym filtrem sygnał dźwiękowy (plik wawa) .Program powinien wyświetlać widmo sygnału przed i po filtracji. Niemam zielonego pojęcia jak się do tego zabrać , może ktoś zgromadził jakieś ciekawe materiały którymi chciałby się podzielić zemną a może znajdzie się kogoś kto by mi napisał...
Jeśli zastosuję powyższy sposób, to nie otrzymam widma takiego jakie oczekuję. Załączam rysunek, na którym jest przebieg czasowy sygnału i jego widmo (jest to sygnał po modulacji QPSK). Widmo jest względem unormowanej częstotliwości próbkowania. Szukam instrukcji w Matlabie, która wykreśli takie widmo. http://obrazki.elektroda.net/16_12519939...
witam mam podobny problem musze napisać program w języku Matlab filtrujący zadanym filtrem sygnał dźwiękowy (plik wawa) .Program powinien wyświetlać widmo sygnału przed i po filtracji. Niemam zielonego pojęcia jak się do tego zabrać , może zgromadziłeś jakieś ciekawe materiały którymi chciałbyś się podzielić zemną a może znasz kogoś kto by mi napisał...
2 rzeczy które mi się rzuciły w oko to: - różnica w rozkładzie energii widma podczas uderzenia. Zwróć uwagę, że niektóre uderzenia mają widmo kończące się na 150-250 Hz z szybkimi "strzałami" powyżej tej granicy (rys 1, 2 i 4), poza tym każde uderzenie, ma po 2 takie strzały - koncentracja energii w czasie. W pozostałych rysunkach (3 i 5) te "strzały"...
cos wyszlo ale tez nie tak jak chcialem:(..jak zastosowac okno hamminga do wektora x??.. bo tych szumow widmo mi nieuwydatnia:( clear all clf x=FFT_x0_1(:,2); x=x-mean(x); % to brakowalo Fs=100; T=1/Fs; L=6009; t=(0:L-1)*T; NFFT=2^nextpow2(L); Y=fft(x,NFFT)/L; f = Fs/2*linspace(0,1,NFFT/2); plot(f,2*abs(Y(1:NFFT/2)))
Powinno wystarczyc przesuniecie calosci w wyzsze pasmo, przy czym widmo jest symetryczne. Jesli nie korzystales z funkcji fftshift to dolna polowe widma przesun w prawo a gorna w lewo newfft = zeros(1,n-1); newfft(1,n/2-1) = fft(1+shift, n/2-1+shift); newfft(n/2,n-1) = fft(n/2-shift, n-1-shift); i odwrotna fft. Wydaje mi sie ze powinno wystarczyc. oczywiscie...
Tak są to dwie różne rzeczy. Dlaczego to takie dziwne? Mając wynik w postaci sygnału wyjściowego mogę zmierzyć zawartość harmonicznych - to chyba nie ulega wątpliwości. Natomiast na podstawie sygnału wejściowego oraz wyjściowego matlab wykreśla charakterystykę przejściową - każdego stopnia z osobna, oczywiście później jeszcze dla pełniejszego obrazu...
Witam Możesz uściślić pytanie? Bo na tak postawione jak powyżej, to na przyklad: widmo_usrednione=mean(widmo) . Pozdrawiam Sanwa
Proponuje Matlab, ma gotowe funkcje do generowania i analizy sygnału. Do analizy użyj transformaty Fouriera - funkcja FFT.
Niestety, ale nie mam namiaru na nich :( Powiem Ci, że rozpoznawanie osób mówiących ten sam wyraz nie jest problemem. Oni to zrobili bez problemu, ale to była tylko tycia część całej magisterki... Głównym tematem ich pracy miałobyć rozpoznawanie słowa wymawianego przez różne osoby, np. bez względu na to czy wyraz "Cześć" powie dziecko/dorosły, kobieta/mężczyzna...
No to kolega się musi jeszcze zastanowić jaką metodą chciałby obliczać FFT i poszukać w sieci wzorów. FFT ma to do siebie że istnieje wiele metod przy czym im prostrza (mnie operacji) tym wiecej zakłóceń i błędów w widmie.
Wie ktoś, jak znaleźć w Matlabie widmo amplitudowe i fazowe sygnału okresowego, jak znam jego wzór?
Już mi się udało. Pokazuje, gdyby ktoś miał podobny problem. Wystarczy : y =abs(A*(sin(2*pi*x*F))); W sumie, niewielka różnica. Mam natomiast problem z jedno połówkowym - jeśli robię przy pomocy if-ów to w widmie fazowym są błędy... .
relacja z postepow: znalazlem funkcje fftshift ktora mi ladnie poukladala wykres [url=http://www.freewebtown.com/fisz/mir... sinc widmo tego sygnalu jednak nie jest widmem jaki powinien miec sinc [url=http://www.freewebtown.com/fisz/mir... mojego sygnalu wyglada tak a powinno byc [url=http://www.freewebtown.com/fisz/mir...
Cześć! Potrzebuję wygładzić widmo sygnału, którego już trasformate Fouriera sobie zrobiłam. Zastanawiam się czy w matlabie jest jakaś fukcja w Signal Processing Toolbox gotowa, czy bede musiała napisać to sama? Chodzi o metodę uśredniania harmonicznych widma względem sąsiednich harmonicznych. Z góry dzięki za jakieś odp
Mam nastepujacy skrypt dotyczący widma amplitudowego i fazowego : clc Fs = 16;% Częstotliwość próbkowania F=2; T=1; %czas końca A=2; %amplituda t = 0:1/Fs:T-(1/Fs); x = A*sin(2*pi*t*F); fpx=fft(x); mf=abs(fpx); p=unwrap(angle(fpx)); %widmo fazowe; ful=(0:length(x)-1)'*100/length(x); %oś częstotliwości stem(ful,mf) figure(2) plot(ful,p*180/pi); %faza...
Witam. Mam problem odnosnie filtrowania w matlabie. Oto zadanie: a) wygenerowac widmo transformatą fft z abs i bez abs, i porownac te widma b) przeprowadzic filtracje, tj. wyciac polowe skladowych gornych lub dolnych(wysokie lub niskie czestotliwosci) c) przeprowadzic filtracje, tj. wyciac parzyste lub nieparzyste skladowe sygnalu. Proszę o pomoc z...
No tak, w końcu widmo jest dodatnie, więc musisz wziąć moduł.
Czy macie pomysł, skąd mogło wyniknąć takie cudo jak obrazku? Już sporo czasu nad tym siedzę i nie mogę do niczego sensownego dojść. https://obrazki.elektroda.pl/5450278900_... [syntax=matlab][/syntax]
Powiem więcej (prawdą jest, że nie powinno Cię to dziwić jeśli chodziłeś na wykłady :P ), że jeśli sygnał jest parzysty w dziedzinie czasu to widmo jest czysto rzeczywiste, jeśli natomiast nieparzysty, to widmo jest czysto urojone. Dokładnie, musisz wykreślić próbki w dziedzinie czasu i voila, co prawda nie wiem co chce Twój ćwiczeniowiec. Co do obserwacji...
Cały czas uzyskiwałem rezultaty podobne do twoich - okazało się, że trzeba po prostu zmienić skalę i oglądać tylko początek wykresu, wtedy od biedy przypomina on ten uzyskany w pracy przedstawionej w linku! Tylko przeliczyłem jeszcze jednostkę na osi x clc; clear; fs = 10000; % probkowanie dane=importdata('pomiar5.dat'); y=dane.data(:,2);...
A gdzie tam wymiękam :) biorę się ostro za lekturę, pobrałem wykłady studentów z politechniki na temat modulacji cyfrowych i postaram się je przestudiować, matlab też się przyda do prezentacji widm sygnałów. Pozdrawiam i dziękuję za pomoc.!
close all clear all % częstotliwość próbkowania fp=44100; Tp=1/fp; % częstotliwość sinusoidy %f0=501; %gdy T0 jest całkowitą wielokrotnością Tp, to widmo nie będzie miało przecieków f0=1/(88*Tp); N=fp*10; t=(0:(N-1))*Tp; x=0.99*sin(2*pi*f0*t); wavwrite(x, 44100, 'Sintone.wav'); y=wavread('Sintone.wav');...
Witam. Czy pomoże mi ktoś z rozwiązaniem tego zadania? ) wyznaczyć i narysować widmo X1(f), (funkcja fft(..)),) dla znormalizowanej wartości częstotliwości oraz amplitudy x(t)=5cos(100πt)+3cos(300πt)+c... Ja próbowałem w ten sposób clear all; close all; fp=40; t=(0:fp-1)*1/fp; x=5*cos(100*pi*t)+3*cos(300*pi*t)+cos(50...
Cześć! Chcesz poprostu zobaczyć widmo sygnału czy je analizować? Jeśli to pierwsze to wizualizacje w Winampie, jeśli drugie to może Matlab. Pozdrawiam! Alfrednovi
Dany jest sygnał cyfrowy x. Jak poprawnie usunąć z widma sygnału X=fft(x) składową o danej częstoliwość - X(m) ,gdzie m-indeks? Środowisko MATLAB. Jedyne co mi do głowy przychodzi to wykonanie okienkowania przed fft, w celu wyeliminowania przecieku, a następnie wyzerowanie prążka widma X(m)=0; Czy jest to poprawnie wykonane?
Witam Mam zadanie, muszę wykonać program w dowolnym środowisko, jezyku,(najlepiej Matlab) który będzie obliczał współczynniki szeregu Fouriera(szereg trygonometryczny) dla jakiś przykładowych sygnałów, dodatkowo może być dodane widmo amplitudowe i fazowe sygnału oraz przykładowa aproksymacja na podstawie wyliczonych wspolczynnikow. Pomoże ktoś ?
"Na piechotę" możesz to zrobić przez wymnożenie widma przygotowaną maską o długości równej rozmiarowi wektora widma. Wartości elementów maski dla częstotliwości które mają zostać wycięte równe są 0, w przeciwnym wypadku 1. Trzeba tylko wiedzieć w jakim miejscu wektora zaczynają się konkretne częstotliwości. Może coś takiego będzie pomocne: Fs=25000;...
Mam problem i to duży:( Mam dokonać FFT dowolnej próbki w formacie .wav i wykreślić widmo amplitudowe do 16 harmonicznej a następnie przekonwertować .wav na wektor zrozumiały dla MATLABA!Zupełnie nie mam pojęcia jak się za to zabrać więc będę bardzo wdzięczny za każdą pomoc!!!!!!!
[syntax=matlab][/syntax] Witam mam do zrobienia na zajęcia zadanie z MATLABA, ale nigdy nie miałem do czynienia z tym programem, a prowadzący nie bardzo chce z tym pomóc. Macie jakieś wskazówki albo gotowy kod który rozwiązałby to zadanie? Wygenerować przebiegi czasowe i policzyć widma dla przebiegów: sinusoidalnego, złożenia kilku sinusoid, prostokątnego,...
Witam Podaję przykładowe zadanie : Impuls prostokątny - Policzyć tw. Furiera - Wykreślić widmo częstotliwościowe i amplitudowe abs(fft(x)) rysuje widmo amplitudowe angle(ffft(x) rysuje widmo fazowe >> t=0:0.01:1; >> w(1:25)=0; >> w(26:75)=1; >> w(76:101)=0; >> figure(1)...
Uśrednianie niekoherentne jest po krótce opisane w ksiące "Wprowadzenie do cyfrowego przetwarzania sygnałów". Generalnie chodzi o to, że należy zebrać n widm amplitudowych sygnału a następnie je uśrednić. Zysk SNR wynosi 10*log(sqr(n)).
Czy przedział czasowy dla x ma być na pewno taki jaki wpisałeś? Przedział dość szeroki a częstotliwość próbkowania niska. Zobacz wyniki dla np. x=[10E-5:0.1:1000]; i normalnej skali bo nie wiem czy to dobry pomysł żeby stosować logarytmiczną. Nie bardzo wiem co chciałeś zrobić w tym kodzie, jest część zmiennych z których w ogóle nie...
Witam wszystkich. Do tej pory moja styczność z Matlabem była praktycznie zerowa, więc pojawił się problem. Muszę wykreślić widmo harmonicznych prądu poprzez FFT. Zadany sygnał prądu mam zapisany jako plik w postaci tekstowej csv. Otworzony w excelu ma 10000 próbek. Na razie tylko podstawową funkcja fft w matlabie wykreśliłęm widmo aczkolwiek widzę,...
Witam, Potrzebuję wyznaczyć wzór na widmo następującego sygnału (falka z gaussowską obwiednią, najlepiej z jakiegoś źródła podać, bo ja nie mogłem znaleźć): http://obrazki.elektroda.net/86_12738757... gdzie: AF - amplituda (stała), t - czas (zmienna), t0 - przesunięcie (stała), tau - czas trwania (stała), f0 - częstotliwość drgań falki (stała)....
W jaki sposób w Matlabie można w najprostszy sposób uzyskac widmo analizowanego obrazu? Bo używam następującego zestawu instrukcji: A=imread('lena.bmp'); B=fft2(A); B=fftshift(B); imshow(abs(B)); i niestety dostaję białe pole. Na czym polega mój błąd. Pomocy!
Odpowiadam na goraco, i to po nocy, sorry jesli niezbyt precyzyjnie. Zgadzam sie ze zastosowanie cepstrum powinno dawac rezultaty przy detekcji dzwiecznosci, ale wiem to jedynie z teorii. proponuje najpierw sprawdzic pare rzeczy prostszymi metodami. 1. zobacz jak wyglada estymator mocy widmowej sygnalu konkretnie stosunek mocy sygnalu w piku do calkowitej...
widmo sygnału matlab analiza widmo matlab widmo sygnał matlab
junkers gasnąć zapalać klapa bagażnika laguna corsa naprawa
Lumel N14 - Przywracanie ustawień fabrycznych krok po kroku Poprawa zasięgu pilota napędu Genius 300 BAT
