Witam. Chciałbym zrobić w matlabie Przesunięcie częstotliwości i kilka herców tak aby nie pwstawało sprzężenie zwrotne pomiędzy mikrofonem a głośnikiem. Sam algorytm przesuwania znalazłem w helpie ale mam proplem z odczytem próbek z karty dźwiękowej tak aby nie tracić żadnej. Próbuje robić na eventach ale po pierwszym evencie wszystko się mi zatrzymuje....
Witam, potrzebuję napisać skrypt(m-file) matlaba który moduluje sinusoidalnie częstotliwość. Mam to zrobić "na piechotę". Oto mój skrypt: fpr=50000; dt=1/fpr; N=10000; n=0:N-1; t=dt*n; fs=100; fm=10; a=(2*pi*fs*t); b=(2*pi*fm*t); x1=230*sin(a); plot(t,x1); pause xm=sin(b); ff=abs(xm); plot(t,ff); pause f=fs*ff; x2=230*sin(2*pi.*f.*t); plot(t,x2); pause...
Witam. Zakładam ten wątek pomimo tego, że istnieje dużo podobnych tematów poruszanych w całej sieci ale niestety nie udało mi się z dowiedzieć w jaki sposób mogę rozwiązać mój problem choć spędziłem nad tym już kilka długich wieczorów. Przedstawie problem od początku. A więc mam za zadanie utworzyć w matlabie funkcję, która będzie rozpoznawać wysokość...
Witam. Mam gotowe próbki sygnału, które ładuję funkcją load. Nie znam częstotliwości próbkowania tego sygnału, jedyne co posiadam, to te wartości. Jak dokonać jego analizy częstotliwościowej?
Witam. Chodzi o pewny algorytm w Matlabie. Dotyczy on badania powierzchni za pomocą urządzenia o pewnej częstotliwości próbkowania, przez którą częstotliwość powstają błędy odczytu (praktycznie niemożliwe do uniknięcia). Posiadamy profil powierzchni powiększonej do tego stopnia że nie jesteśmy w stanie odczytać jego kształtu ani falowania materiału,...
Witam, mam problem z rysowaniem wykresu FFT w Matlabie, otóż mam funkcje: x=(2-400*t).*cos(2000*pi*t) Dla tej funkcji dobrałem częstotliwość próbkowanie f=16kHz, i przyjąłem N=16 próbek na okres, czyli analizując 2 i 4 okresy kolejno otrzymam 32 i 64 próbki. Kod mojego programu w matlabie jest następujący: >> t=linspace(0,2*0.001,10000); >> A=(2-400*t)...
Próbuję w Matlabie zasymulować demodulator amplitudy oparty na transformacie Hilberta.W skrócie metodę można opisać jako sqrt (HT(n)^2+s(n)^2), gdzie HT jest transformatą Hilberta, a s(n) sygnałem odebranym. Trzeba wziąć pod uwagę opóźnienie wprowadzane przez HT. W przypadku filtru przyczynowego o długości N=2M+1 próbek opóźnienie wynosi M próbek. N-ta...
przykładowy programik z fft clear all close all clc tp=.001; t=0:tp:.112; f=50; y=sin(2*pi*f*t)+5*sin(2*pi*3*f*t+8); %próbki z sumy sinusoid p=fft(y); %szybka transformata furiera g=abs(p); figure plot(t,y) % tak wyglądają próbki N=length(t); figure stem((0:1:N-1)/tp/N,g/N*2) %tu mogłem coś spieprzyć xlabel('czestotliwosc') ylabel('amplituda') jeden...
Mam za zadanie zrobić filtr IIR dla zadanych wartości (w matlabie). Dane: f=1000 fs= Za mało informacji. Co to jest f=1000? Prawdopodobnie częstotliwość próbkowania, którą przyjęło się oznaczać fs (od sampling frequency). Co to zatem jest fs? Może wektory sygnałów testowych o częstotliwościach 100, 300 i 400, którymi masz sprawdzić filtr? Masz zrobić...
Witam serdecznie. Mam zadanie, z którym utknąłem. Otóż napisałem kod, który odtwarza przebieg na podstawie danych zaimportowantch w formacie .dat z oscyloskopu Siglent SDS1104X-E, następnie robię FFT tego sygnalu i wyświetlam w oknie oba rysunki. Jednak apetyt rośnie w miarę jedzenia i chcialbym zrobić rozkład widma w 3D jak na załączonym obrazku i...
Potrzebuje zaprojektowac filtr antyaliasingowy w Matlabie. Tzn mam zrobic skrypt, w którym wykorzystam filtr dolnoperzpustowy do tłumienia wyższych częstotliwości powyżej częstotliwości Nyquista. Czy ma ktoś może materiały na ten temat ,podobny skrypt,albo chociaż wie jak trzeba napisać ten skrypt,żeby mnie nakierować w tym temacie. Z góry dziękuję...
Witam. Moja wiedza na temat matlaba jest słaba, jednak potrzebuje wykreślić w nim charakterystykę filtru. To co posiadam to równanie różnicowe ( lub równoznaczne w domenie Z). Aby wykreślić charakterystyki filtru należy za " Z " podstawić "e^jw" ( "e" do potęgi j*omega). I tyle w teorii. Wychodzi skomplikowane równanie które przy dowolnej zmianie trzeba...
Dzień dobry, ciekawy pomysł, raczkuję jeszcze w MATLABie ale dodam od siebie jak ja to widzę aby rozwiązać to w miarę dynamicznie. Od strony MATLAB-a dostajemy funkcje: wavread/wavwrite do odczytu/zapisu danych z plików na dysku. Statycznie można po prostu wczytać pliki, odjąć 2 funkcje, zapisać. Aby było to dynamiczne można by wykorzystać potoki, zaś...
Witam wszystkich forumowiczów. Mam problem z filtracją w Matlabie. Mam napisany kod filtru pasmowoprzepustowego metodą eliptyczną. Program działa jednakże wydaje mi się, że wyniki nie są do końca poprawne. Prosiłbym o sprawdzenie kodu i określenie w czym leży problem. Widmo sygnału po filtracji jest w zasadzie takie samo jak przed obróbką, co wskazywałoby...
- kazde wymuszenie ma swoje parametry. Ty ustawiles V1 jako sinus 0 10 5k. Czy parametr 5k to rezystancja wewnetrzna zrodla wymuszajacego czy moze maksymalna wartosc amplitudy 5kV ? Tak jak zauważył (at)elektryku5 10 to wartość amplitudy sinusoidalnej, a 5k to częstotliwość. Jak dałem tylko jeden warunek początkowy to już wartość prądu oscyluje się...
Witam Do kolegi "kedar50". Co masz na myśli pisząc "usunąć przeciek"? Piku przy 100Hz nie otrzymasz bezpośrednio, bo nie masz danych dla tej częstotliwości (w siatce częstotliwości nie ma tej wartości). A)Jeżeli chcesz zostawić tylko pik dominujący, to możesz wyszukać maksimum lokalne (tutaj około 97Hz) i "wyzerować otoczenie". ""czyścisz tak długo,...
Ja tam nie wiem, jak to w Matlabie wygląda, ale w MathCADzie by to jakoś tak wyglądało (zakładamy, że wszystkie indeksy tablic zaczynają się od zera) 0. Próbka dźwięku to wektor (V) 1. wybieramy ilość próbek do przetworzenia przez jedną transformatę (N) 2. Tworzymy macierz M o odpowiedniej wielkości: 2.1 liczba wierszy - dla transformaty liczonej z...
Witam, chcę odfiltrować z sygnału pewną częstotliwość przy pomocy fft. Wykonuję filtrowanie w matlabie. Problem pojawia się przy sprawdzeniu. Jeśli analizuję cały filtrowany sygnał wszystko jest w porządku ale jeśli poddam analizie fragment sygnału odfiltrowana częstotliwość pojawia się ponownie. Myślałem że jest to kwestia nieprecycyjnego wyznaczenia...
Do wykonania sa dwa projekty - bez Waszej pomocy moze sie to nie udac... Temat projektu: Cyfrowe przetwarzanie sygnałów Opis: Twoim zadaniem jest zapoznanie się z tematyką szybkiego przekształcenia Fouriera oraz projektowania filtrów cyfrowych przy użyciu oprogramowania Matlab. Zadania: 1.Zaprojektuj filtr o nieskończonej odpowiedzi impulsowej metodą...
Do tego się mierzy kontrast pomiędzy mało oddalonymi punktami - jakaś transformata Fouriera, obcięcie niskich częstotliwości, i policzenie normy - przy ostrym nastawieniu ma maksimum.
To jakie są dokładne wytyczne co do tego zadania, które kolega chce rozwiązać? To mże po zrobieniu transformaty fft po prostu wyzerować pierwszy, ostatni lub co drugi element tablicy częstotliwości.
Coś tu nie gra. Skoro analizujesz przebieg sygnału słyszalnego dla ucha ludzkiego, gdzie częstotliwość graniczna górna to 20 kHz to powinieneś próbkować z częstotliwością przynajmniej 2x większą czyli 40 kHz zgodnie z twierdzeniem Nyquista. Co to ma wspólnego z zadanym pytaniem?
Witam wszystkich. Mam do napisania program w Matlabie generujący różne dźwięki. Założenie jest takie, że użytkownik naciskając dany klawisz wywołuje konkretny dźwięk. Czyli ma to być coś w stylu symulatora pianina, czy czegoś w tym rodzaju. Jako, że nigdy nie miałem styczności z Matlabem nie wiem, jak w tym środowisku wygenerować dźwięki o określonej...
Witam, W ramach realizacji zadania laboratoryjnego, mam zaprojektować filtr FIR. Ponieważ moja wiedza w tym temacie jest uboga, a instrukcja do wykonania niewiele bogatsza, proszę o wytłumaczenie, jak to zrobić, może ktoś zaproponuje kod. Poniżej polecenie z instrukcji: "Wywołanie tej funkcji ma postać : b = remez (N, f, mag) Gdzie argumenty wejściowe...
Witam mam następujący problem: wygenerować impuls jednostkowy o długości 200 próbek; narysować wykres sygnału zaznaczając próbki kropkami nie bardzo wiem jak za to się wziąć, z poleceniem: wygenerować impuls jednostkowy o czasie trwania(tj. długości) 0.137 s przy częstotliwości próbkowania fs=1750fs=1750 Hz ktoś wysłał mi taki kod - nwm co to znaczy,...
Witam, zaczynam się uczyć sygnałów i to jest moje pierwsze polecenie. Czy mógłby ktoś krok po kroku napisac i wytłumaczyć mi ten program ? Napisać funkcję która wykona kwantyzację dowolnego sygnału wejściowego, przy podanej (dowolnej) liczbie bitów przetwornika. Jako wynik funkcja ma zwraciać sygnał skwantowany i sygnał błędu kwantyzacji. Zakładamy...
Witam otóż mam mały problem, gdyż muszę zaprojektować filtr FIR typu pierwszego low pass (LP) w programie MATLAB. Mam już praktycznie cały kod programu, lecz gdy wpisuje dane np. rząd filtru 31, częstotliwość próbkowania 1000, częstotliwość 200 to za każdym razem wyskakuje mi ten sam błąd: Error using + Matrix dimensions must agree. Error in fir1lp...
Z okresu próbkowania i ilości próbek transformaty. Na szybko z manuala do matlaba. df = Fs/NFFT; df- kwant częstotliwości Fs- częstotliwość próbkowania NFFT, ilość próbek transformaty. Sprawdź to, mogłem się pomylić. Ps. :) nie zapomnij modułów znormalizować. Dodano po 13 aha, dla jasności NFFT to nie jest ilość próbek sygnału, tylko ilu próbowa jest...
Witam. Mam pewien przebieg dyskretny i analizuje 5 okresów składowej o częstotliwości f1. oto kod w Matlabie: [syntax=matlab] clc clear all; close all; A1=1; A2=-0.5; A3=0.3; A4=0.2; A5=0.1; f1=100; f2=2*f1; f3=3*f1; f4=4*f1; f5=5*f1; fi1=0; fi2=(pi*30)/180; fi3=(pi*45)/180; fi4=(pi*60)/180; fi5=0; N=128; t=0: (1/f1)/128:(5*(1/f1)); x1=A1*sin(2*pi*f1*t)+fi1;...
Wszystko zależy od tego czy jest potrzebny wobulator, czy charakterystyka filtru. Do tej pory nie używałem Matlaba w bezpośrednim połączeniu ze światem zewnętrznym (bo uzależniłem się od LabViewa). Jeżeli jednak masz opanowane wystawianie sygnału i pobieranie próbek, lub całych rekordów danych to nie widzę przeszkód. Do realizacji "wobulatora" wystarczy...
Jak w temacie - mam problem z za modelowaniem w matlabie siusoidy prostowanej dwu-połówkowo. Jest może jakoś możliwa do wywołania - lub jakiś zbliżony sygnał ? Wykonałem to w następujący sposób: y = (A)*abs(sin(2*pi*x*F)); gdzie : A - amplituda x=0:1/128:1; F - częstotliwość Problem polega na tym, ze po wywołaniu odpowiedniej funkcji do liczenia widma...
Witam, mam problem z wczytaniem pliku audio i zmianie wartości w wektor binarny. Mój kod to funkcja o następującej formie: [syntax=matlab] function [ y, fp ] = audiot( x ) [sound fp]=audioread(x); dod=abs(sound) dzwiek=(uint16(dec2bin(dod))) [/syntax] Problem polega na tym, że wynik takiego zabiegu to macierz o wartości 48, co chyba nie jest prawidłowym...
(at)wesolyyyy zagapiłem się, w zadaniu jest podane, że to rezystancja wew. transformatora. Pewnie mogę sobie dobrać jakąś. Tylko teraz mam problem z tego w matlaba. Jakiś pomysł jak mógłbym wyliczyć prąd na kondensatorze. Z napięciem jest dosyć łatwo. Jeśli < 0 to zmieniam znak. Ale za nic w świecie nie wiem jak obliczyć prąd. Próbuje coś z tym:...
Witam mam do zrobienia projekt w matlabie "Rozpoznawanie tonów sygnalizacji w telefonii analogowej " NIe mam pojęcia jak to zacząć wiem tylko że ma to być na zasadzie filtrów dla każdej częstotliwości. I problem rodzi się w tym że całkowicie nie potrefię pisać w matlabie. więc poszukuje chętnej osoby która mi w tym pomoże, a ja się chętnie z nią jakoś...
Witam, jestem nowy na tym forum i za bardzo nie wiedziałem w jakim dziale umieścić mój problem - jeżeli istnieje "bardziej odpowiedni" dział to przepraszam. Mam do rozwiązania następujący problem: Posługując się metodą okien czasowych zaprojektować filtr pasmowo-przepustowy o danych częstotliwościach odcięcia oraz częstotliwości próbkowania. Rząd filtru...
Witam, mam do wykonania projekt z cyfrowego przetwarzania sygnalow. Musze napisac program, ktory zmierzy czestotliwosc podanego sygnalu. Oto zalozenia projektowe: 21)Pomiar częstotliwości Do zadań należy: -zapoznanie się z wybranymi algorytmami pomiaru częstotliwości -wybranie zaprogramowanie wybranego algorytmu -wykonanie symulacji dla wybranych sygnałów...
Mam narysować FFT do sygnału, czytam i próbuję od kilku dni, nic mi nie wychodzi. Proszę o pomoc i jakiekolwiek wskazówki. W Matlabie robię to pierwszy raz, proszę o wyrozumiałość co do moich może (banalnych pewnie) pytań. Chcę uzyskać wykres gdzie na osi x będzie częstotliwość , na osi y będzie moc. 1) Mam problem z wczytaniem danych z pliku. Jak odwołać...
Muszę w matlabie zasymulować taki oto schemat: Jednak na samym początku mam jakiś błąd którego nie mogę wyłapać, od strefy nieczułości (Z->A), przynajmniej tak mi się wydaje. Nie wiem czy dobrze interpretuję tę strefę nieczułości. Jeśli mam dwa sygnały wejściowe X i Y, więc moje Z będzie 2, to po przejściu przez tę strefę nieczułości co powinnam...
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...
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. Męczę się z tym od kilku dni. Może zna ktoś jakiś prosty sposób, bo ja chyba zbyt komplikuje sobie sporawe, a zapewne są już opracowane metody. Chciałbym aby w simulinku w czasie trwania symulacji wyświetlała mi się aktualna częstotliwość sygnału. Obecnie sam napisałem algorytm w m.pliku, działający mniej więcej tak, że wykrywam sobie przejście...
Jak w temacie - mam problem z za modelowaniem w matlabie siusoidy tyrystorowej. Jest może jakoś możliwa do wywołania - lub jakiś zbliżony sygnał ? Najlepiej w podobny sposób gdyby się dało : y = (A)*square(2*pi*x*F); gdzie : A - amplituda x=0:1/128:1; F - częstotliwość Może, nie wiadomo o jaki sygnał mi chodzi - dlatego dodaję rysunek :
Hej. Czy ktoś byłby w stanie pomóc mi napisać kod w Matlabie do stworzenia takiej animacji: ? Mam definicje tych zmiennych, umiem wyświetlić te wykresy oddzielnie, ale mam problem z napisaniem kodu do takiej animacji. Zmienne są zależne od częstotliwości (lub czasu). Korzystać z funkcji refresh czy comet czy jeszcze jakiejś innej? Gdyby ktoś otrzymał...
Witam. Interesuje mnie następująca kwestia. Zaprojektowałem w dziedzinie czasu pewien układ, który wykonuje pewien algorytm. Algorytm ten wyznacza mi wartość liczby, a przebieg działania algorytmu w postaci wartości liczby w kolejnych informacjach jest tą moją " odpowiedzią układu". Chciałbym teraz zamodelować taki układ w dziedzinie częstotliwości....
Pierwsza czynność ---> wstęp teoretyczny nt. filtrów (np. wikipedia) Druga czynność --> oprogramowanie np. Matlab ale też programy freeware: free filter 2007, elesie, i wiele innych na sieci Jak juz wybierzesz matlaba, to HELP FDATOOL (filter design) i krok po kroku, metodą "obrazkową" zaprojektujesz FDP i wyznaczysz jego parametry oraz charakterystykę....
Witam Serdecznie Próbuje zasymulować celowe zakłócenia( nie zakłócenia środowiskowe i losowe) sygnału radiowego dla transmisji 802.11g Niestety jak to w życiu trafiłem na kilka problemów których nie potrafię rozwiązać i liczę na sugestie/pomysły bardziej doświadczonych osób. Model warstwy fizycznej 802.11g który wykorzystałem pochodzi z tond: Elementy...
Hmm dzięki dzięki za tą teorię ale ja mam zaprojektować ten filtr teoretycznie. Napisać skrypt w Matlabie, oparty na funkcji butter i to właśnie ten skrypt ma mi odfiltrowywać odpowiednie częstotliwości.
Od wygenerowania tonu podstawowego i zmodulowania jego częstotliwości sinusoidą ok. 7Hz (czyli zwykła modulacja FM wąskopasmowa) Dewiacja częstotliwości np. w zakresie +/- 0.25 tonu Kolejne tony powstają przez mnożenie lub dzielenia tonu "a" razkreślne (440 Hz) przez pierwiastek 12 stopnia z dwóch (1.059463094) Głębokość modulacji (jako wielkość odchyłki)...
Witam. Mój problem sprowadza się do tego jak zapiać sinusoidę koherentną za pomocą matlaba w formacie *.wav. W chciałem uzyskać częstotliwość 500 clear all; close all; % zerowanie N=44.1*10000; fp=44100*10000; Tp=1/fp; %okreslenie ilosci probek T1=N*Tp; % okres dla sinusoidy koherentnej f1=1/T1; %f koherentna t=(0:1:(N-1))*Tp...
Witam zna sie ktoś na funkcji fft w matlabie? znalazlem programik t = 0:0.001:0.6; x = exp(-100.*t); y = x ; plot(y(1:50)) Y = fft(y,512); Pyy = Y.* conj(Y) / 512; f = 1000*(0:256)/512; plot(f,Pyy(1:257)) jednak nie do końca go rozumiem ja chciałbym przeprowadzić analize na zbiorze danych dyskretnych (1000) ale nie umiem dopasować częstotliwości tak...
Oto kod programu napisany przez mojego znajomoego function [mixedVector, fadeOut, fadeIn] = crossfade(leftSig, rightSig, crossLength, nCase) % Przejście jednego sygnału w drugi za pomocą zmiany amplitud. % autor: Bernard Latanowicz latanb(at)wp.pl % nargin % leftSig: lewy wektor % rightSig: prawy wektor % crossLength:...
Analiza sygnału: x(t)=4×sin^3 8000πt×cos3000πt Zadania: Dobór parametrów próbkowania (częstotliwości oraz ilości próbek) Przeprowadzić szybką transformatę Fouriera Przekształcenie wzoru sygnału: x(t)=4×sin^3 8000πt×cos3000πt x(t)=4×sin^2 8000πt×(sin8000πt×cos3000π... x(t)=4×sin^2 8000πt×(1/2[sin(11000πt...
Witam. Ogólnie sprawa wygląda tak, że po wygenerowaniu sygnału FM w celu symulacji kanału dodaję do niego zakłócenia. Następnie przepuszczam tak zmodyfikowany sygnał przez filtr pasmowy(szerokość pasma ustalam na podstawie reguły Carsona - Bnd=2*(df+fm), df -dewiacja częstotliwości, fm częstotliwość fali modulującej). Cały problem polega na tym, że...
Witam napisałem program który ma za zadanie generować 128 próbek sinusoidy o częstotliwości 2000Hz z różnymi częstotliwościami próbkowania. Moim problem jest to, że owszem wychodzi 128 próbek ale łącznie na całym przebiegu co jest trochę bez sensu. Bo te próbki powinny być na okres. Niech ktoś mądrzejszy da jakieś wskazówki... Wrzucam skrypt. N=128;...
Witam, mam do napisania skrypt w matlabie, a kompletnie nie wiem jak to zrobić: Mam do wygenerowanie następujące sygnały w czasie 1. sinus o częstotliwości 30 Hz 2. impuls prostokątny 80 ms z okresem powtarzania 0.4 s 3. impuls sinusoidalny jak w punkcie 1 kluczowany impulsem prostokątnym z punktu 2 Za jakąkolwiek pomoc będę wdzięczny
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, pewnie to zadanie jest już przedawnione, ale pozwoliłem sobie zapuścić na nie 'scilaba' (=matlab free) kod: //--------------------------- clear; //czyszcze wszystkie zmienne //ilość punktów do przybliżenia generowania modułu transmitancji N=128; //próbkowanie fpr=480000; //zadany wektor współczynników filtra h=[0.0717 0.0451 -0.0956 -0.2705 0.5587...
Założyłem ze częstotliwości sie zgadzają. Na początek wywaliłem w edytorze tekstowym pierwsze linie z opisami zmiennych. load gl_logchirp_ff.txt % otwieranie zbioru load mb_logchirp_ff.txt g=gl_logchirp_ff; % zmiana nazwy wektora na krótszą m=mb_logchirp_ff; plot(g(:,1),g(:,2)); % podstawowy wykres, skale liniowe, 1 kolumna x , 2 kolumna y plot(g(:,1),g(:,3));...
Zasadniczo należy filtrować. W tym przypadku dolnoprzepustowo z odpowiednim dobraniem częstotliwości odcięcia. Filtracja (o ile pamiętam) y=filter(A,B,x) albo filtfilt. A i B to wielomiany licznika i mianownika transmitancji H(z). Otrzymuje się je projektując filtr wg. któregoś z mędrców np. butter, cheby1, cheby2, ellip. Pozostaje jeszcze sprawa używania...
Witam Mam do odfiltrowania z sygnału PWM trzy pierwsze harmoniczne (50, 100 i 150 Hz). Najlepiej użyć filtru środkowo-przepustowego, tak by przepuszczał wszystko w pasmie od 50 do 150Hz. Do tego celu chciałem użyć filtr z bloczku "Analog Filter Design" - nastawiłem częstotliwość dolną (50*2*pi rad/s) i górna (150*2*pi rad/s). Niestety to co uzyskałem...
Sinusa o zadanej częstotliwości robi się tak (sorry, że w C, ale ten język jest dla mnie bardziej naturalny): for (i = 0; i < ilosc_probek; i++) { tablica[i] = sin(2 * PI * i / ilosc_probek); } To daje próbki zmiennoprzecinkowe. Jeżeli mamy 10-bitowy przetwornik, to zakres jest (-512; +511), zatem mnożymy próbki...
Rozumiem teraz mniej więcej, plik txt pomógł. Kolejne pytanie - kolumna buffer zawiera kolejne próbki z bufora, a 10kHz to częstotliwość próbkowania? To znaczy kolumna ta zawiera kolejne odczyty z urządzenia w chwilach 0,1,2,... itd. Jeśli 10kHz jest częstotliwością próbkowania, to każdy kolejny odczyt następuje co 1/10000 sekundy? Czekam na potwierdzenie...
już wiem jeżeli chodzi o aliasing to jeden sygnał ma częstotliwość 5Hz drugi 9Hz sygnały miały zostać pomnożone x3=x1.*x2 czyli częstotliwości dodajemy 5+9=14Hz aby nie było aliasingu musimy przyjąć częstotliwość 2 x większą 14 x 2 = 28Hz ≈ 30Hz
Witam! Borykam się z takim oto problemem: "korzystając z narzędzia GUIDE zbudować wirtualny generator sygnałów, który pozwoli na wysłanie generowanego sygnału poprzez łącze RS-232 komputera. Generator powinien mieć możliwość generowania sygnału sinusoidalnego o zmiennej amplitudzie i częstotliwości lub sygnału prostokątnego o zmiennej częstotliwości,...
% --- Executes on button press in wykres. function wykres_Callback(hObject, eventdata, handles) % hObject handle to wykres (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) axes(handles.axes1) f= str2num( get( handles.edit1...
witam Szanownych Użytkowników, chciałbym prosić Was o pomoc, mianowicie, napisałem program, który nie jest do końca dobry, gdyż występują w nim 2 zasadnicze problemy: 1. Wykres sygnału przefiltrowanego i zdecymowanego jest "przesunięty" nieznacznie względem wykresu sygnału podstawowego. 2.(ważniejszy) nie mam pojęcia dlaczego program nie chce mi wyrysować...
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 Jestem w trakcie pisania symulacji na pracę inżynierską. Wykonuję symulację transmisji STBC-VBLAST dla symboli OFDM w kanale z zanikami. Nie będę się rozpisywał o samej pracy. Mój problem polega na tym, że w celu wyznaczenia macierzy kanału dziele sygnał wyjściowy przez wejściowy w dziedzinie częstotliwości. Teraz potrzebuję przekonwertować macierz...
Kolego, pytasz o podstawy z przetwarzania sygnałów. Pozwól, że przypomnę Ci transformaty sinusa i cosinusa: Myślę, że po przeczytaniu zrozumiesz skąd biorą się zespolone wartości DFT. Aby usłyszeć w miarę dokładnie potrzebujesz n=fs próbek, wtedy sygnał trwa 1 sekundę. Częstotliwość próbkowania naturalnie potrzebna jest podczas odtwarzania w matlabie.
a jaka częstotliwość próbkowania? N = 128; Tp = 0.5; t = 0:Tp:(N-1)*Tp; w = pi/20; y = sin(w*t); plot(t,y)
Nie, pytam, bo nie wiem jak się do tego zabrać. Może ktoś ma jakieś pomysły? Chcę przedstawić te dane (ruch ręki), żeby zaobserwować jakieś charakterystyczne zmiany w częstotliwościach. Przyszedł mi do głowy FFT. Ale może zamiast 3D raczej lepiej zrobić 2D: x i z oraz y i z? Czy to ma sens? Albo jakie inne funkcje matlaba mogę użyć, żeby zaobserwować...
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: NFFT = 2^nextpow2(L); % Next power of 2 from length of y Y = fft(y,NFFT)/L; f = Fs/2*linspace(0,1,NFFT/2+1);...
PWM: czy chodzi o to, że wartości 1001 - 1003 nie zgłaszają błędu? Sprzętowo PWM jest ośmiobitowy, zakres 0..1000 jest przeliczany na 0..255 (we * 100 / 392) i tak wychodzi. Czy możesz podać szczegóły odnośnie zależności czasowych? Częstotliwość PWM, czas cyklu sterownika (bez przerwań), częstotliwość i czas trwania impulsów z wentylatora, w jakim czasie...
Z zadanie mam wykreślić charakterystyki częstotliwościowe jeśli mam podane bieguny,zera i wzmocnienie. Bez używania funkcji zp2tf i freqz. Przy zamianie z dziedziny Z na częstotliwość korzystałem z wzorów z załącznika. Problem w tym że mój wynik nie pokrywa się z tym co ma być. Ktoś mi podpowie gdzie robię błąd? Mój kod:
Może mi ktoś wytłumaczyć, jak mam zaprojektować, by "Matlab" wyświetlił sygnał trwający 1 sekundę i 100000 punktów wykresu? Jak to się robi? Czytam różne skrypty i nie umiem znaleźć zależności. Albo wygenerować sygnał sinusoidalny o częstotliwości próbkowania 1000 i czasie trwania 75 ms. Fajnie by było, jakby ktoś to łatwo wytłumaczył :)
Jest tutaj dysonans nazewnictwa. Widmem sygnału zwykło się nazywać sygnał w dziedzinie częstotliwości, więc określenia "widmo zespolone" czy "widmo amplitudowe" jest czymś dla mnie dziwnym. Cóż jak zrobisz fft z sygnału to otrzymasz w wyniku wektor liczb zespolonych które możesz przedstawić na wykresie, ale to nie będzie widmo. Jeśli z takiego wektora...
Witam, W jaki sposób mo¿na prosto policzyć częstotliwość główną i harmoniczne, a także THD sygnału, dostępnego wyłacznie jako kolejne próbki - zapisane do pliku. Do wykorzystania możliwe są środowiska Matlaba i C, ale najpierw interesuje mnie teoria do tego tematu. Czy ktoś mógłby coś pomóc? Może ktoś robił coś podobnego i móglby sie podzielić...
Wrzuć dane do Matlaba i zobacz co wychodzi.
hm. Chwila przypomnienia.... Najłatwiej jest dobrać regulator typu PID na podstawie porównania transmitancji Twojego układu objętego jednostkową pętlą sprzężenia zwrotnego do równania oscylatora. W literaturze opisy występują często, może nie zwórciłeś na to uwagi, lub Twój prowadzący nie nakierował Cię na odpowienią. Mniejsza o większość, bo z chęcią...
hm, trochę to strzelanie do wróbla z armaty ale może to rozwiąże problem. Spróbuj wygenerować wektor modułu widma i przefiltrować go filtrem dolnoprzepustowym (eksperymentalnie dobierając częstotliwość odcięcia).
No w sumie wartość omegi nie powinna mieć tutaj nic do znaczenia bo i tak bilans prądów powinien się zgadzać dla każdej częstotliwości (prawo Kirchoffa (prądowe) dla węzła V2 się nie zgadzało). Z tego co wiem to są problemy w tej metodzie jeśli w gałęzi występuje idealne źródło napięcia. Wydaje mi się że trzeba zastosować metodę eliminacji idealnego...
Jeżeli ma to wyglądać analogicznie jak w tym linku to nie rozumiem problemu, zamieniasz tylko dane i przecież jest po problemie. Cały kod jest więc problemu brak. Częstotliwość 56000 to chyba nie jest potrzebna do wykresu. Pozdrawiam
Witam, mam małe pytanko do kolegów, orientujecie się może jak wycina się wybrane prążki. Próbuje w Matlabie wyciąc poleceniem filter ale nic mi to nie daje, pokazują mi się tylko dwie kreski nie moge nazwać tego prążkami gdyż zajmują one miejsce na początku i końcu osi częstotliwości. Głowie się już od miesiąca i jakoś mi to nie za bardzo wychodzi....
Witam. Prosił bym o pomoc z następującym zadaniem. Z sygnału: x=sin(2*pi*1600*t)+sin(2*pi*4000*t)+sin(... próbkowanego 44100, wyciągnąć pojedynczą częstotliwość za pomocą filtrów: butter, cheby1, cheby2 10-tego rzędu To akurat zrobiłem, ale mam porównać odpowiedzi impulsowe położenia zer i biegunów, charakterystyki częstotliwościowe, skuteczność...
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.
We współczesnychsystemach HVAC(HEATING,Ventilation,AirCondition)st... się regulatory dwupołożeniowe(termostaty)z ujemnym sprzężeniem zwrotnym wokół elementu przekaźnikowego.Sprzężenie zwrotne zmniejsza oscylacje wartości regulowanej kosztem zwiększenia częstotliwości przełączeń w układzie .Zaproponuj schemat Simulinka ilustrujący działanie zmodyfikowanego...
Dzięki za odpowiedź ale nie do końca mi o to chodziło. Mam wygenerować sygnał o stałej amplitudzie i zmiennej częstotliwości (najlepiej jakąś sinusoidę), który będzie tą częstotliwość zmieniał w sposób ciągły od jednej wartości do drugiej(np od 60 Hz do 70 Hz). Potrzebny mi taki sygnał po to by poddać go analizie za pomocą krótko-czasowej transformaty...
Witam Piszę skrypt, który ma na celu analizę częstotliwościową w pętli, w krótkim przesuwnym oknie czasowym(krótko-czasowa transformata fouriera). Celem analizy jest uzyskanie zależności częstotliwości od czasu. Wynik analizy należy porównać ze znaną zależnością częst. od czasu(znam ją gdyż sam generuje sygnał, który później analizuje). Na razie napisałem...
Cześć! Mam do zrobienia projekt filtrów. Filtry już zaprojektowałem i działają, problem tylko mam z wykreśleniem modułów ich funkcji transmitancji w skali logarytmicznej. Możecie powiedzieć jak to trzeba zrobić? Filtry są tak zrobione: % Tworzymy wektor częstotliwości filtru f4=[0 2950 3000 4000 4050 fNq]; ff4=f4/fNq; % Wektor pożądanych amplitud...
No ja właśnie też sobie przerobiłem program - już nie podaję ilości próbek - tylko częstotliwość próbkowania, oraz nie podaje długości okresu tylko częstotliwość i w ten sposób tworze zadany sygnał i potem widmo i też mam ok - wydaje mi się , ze przy poprzednim zapisie (jak podawałem ilość próbek, długość okresu itd i tak zapisany sygnał przechodził,...
Witam. Mam nadzieję, ze nie trafilem źle, jeśli temat jest w złym dziale, to przepraszam za problemy. Usiłuję przeprowadzić analizę częstotliwościową sygnału dźwiękowego w formacie .wav. Na własną rękę, korzystając z własnej wiedzy nie za bardzo udało mi się to zrobić. Zacząłem więc korzystać z poradników i tutoriali. Jednym z nich jest ten tutaj: On...
Z tego co widzę, to kolega nie rozpisał się zbyt dokładnie na podany temat, ale najpierw wypadałoby wywołać plik kk.m, a następne wywołania pliku kkk.m będą powodować kolejne wywołania generatora. Przyglądając się teamu schematowi odnoszę wrażenie, że sygnał wyjściowy z bloku Transfer Fcn jest nieprawidłowo połączony z sygnałem filtra. Jeden z nich...
Użyj znaczników syntax - popraw post. To co narysowałeś to jest Twój sygnał, z którego chcesz obliczyć widmo ? FFT daje informację o sygnale w dziedzinie częstotliwości. Rozumiem, że chcesz wykreślić PSD. Żeby je wykreślić trzeba podnieść wartość bezwzględną FFT do kwadratu i podzielić przez 2 pi. [syntax=bash]widmo = fft(y); PSD = abs(widmo).*abs(widmo)/2/pi;...
Mam sygnał prostokątny, który generuje za pomocą poniższego skryptu. function [t,x] = rect(A,f,fp,T,w) t=zeros(1,T*fp); x=t; tmpA=floor(w*fp/f); tmpB=floor((1-w)*fp/f); k=1; i=0; while k<=T*fp for j=1 : tmpA x(k)=A; t(k)=i; i=i+1/fp; k=k+1; end; for j=1 : tmpB x(k)=0;...
Przetestuję i pochwalę się wynikami bo doszedłem już do tego jak zrealizować tę funkcję na bloku Signum. Niestety nie jestem w stanie tego zrobić dzisiaj ale na pewno zrealizuję obie opcje przez weekend. Dziękuję za zainteresowanie tematem. Pozdrawiam Dodano po 1 Jednak znalazłem chwilkę czasu żeby przetestować blok o którym mówiłeś. Niestety nie ma...
Obiektem regulacji jest rzeczywisty elektryczny silnik liniowy, zaimplementowany w środowisku MATLAB & Simulink. Ktoś mi kiedyś powiedział, że częstotliwość dla regulatora nadrzędnego powinna w tym przypadku wynosić 1kHz a dla podrzędnego od 20 do 50 kHZ. A te 1kHz i 50kHz to częstotliwość próbkowania czy tak zwane crossover frequency dla regulatorów?
Po wywołaniu polecenia Nyquist otrzymujemy charakterystykę, na której dzięki odpowiedniej funkcji możemy zaznaczyć zapasy fazy i modułu. Problemem dla mnie jest to, że na tym wykresie zapas modułu podany jest w decybelach, a ja potrzebowałbym możliwości wypisania go w wartościach bezwzględnych. Czy i gdzie jest opcja, która pozwoli mi na taka zmianę...
Spróbuj tego: [syntax=cpp]%DANE: fp = 250;%częstotliwość próbkowania Tp = 1/fp; fs = 50;%częstotliwość sygnału do odfiltrowania [Hz] ws = 2*pi*fs; %transmitancja filtru: %H = 1 - 2*cos(ws*Tp)*z^(-1) + z^(-2); %współczynniki filtru: L = [1 -2*cos(ws*Tp) 1]; M = 1; sig2 = filter(L, M, sig2);[/syntax] To jest filtr zerujący. Ja używałem go do odfiltrowania...
Witam, mam z pozoru łatwy problem ale nie potrafię sobie z nim poradzić. Mam dwa pliki z pomiarami, napięcie w funkcji częstotliwości. Obydwa w tym samym zakresie częstotliwości tylko jeden z większym próbkowaniem, w konsekwencji mam pliki o różnej ilości komórek ale w tym samym zakresie częstotliwości. Czy da się w matlabie tak przeskalować ten krótszy...
Zgoda, ale od 0.5 na osi częstotliwości charakterystyka się nie pojawia (właściwie to pojawia, ale została ukryta zmianą zakresu fazy). Dlaczego?
