logo elektroda
logo elektroda
X
logo elektroda
REKLAMA
REKLAMA
Adblock/uBlockOrigin/AdGuard mogą powodować znikanie niektórych postów z powodu nowej reguły.

Odpowiedź Impulsowa ???

Sebastian. 06 Maj 2005 16:12 7505 7
REKLAMA
  • #1 1466896
    Sebastian.
    Poziom 11  
    Witam
    Może na początku napisze dokladnie w calosci co chce zrobic, mianowicie
    pisze programik ktory generuje z karty muz sygnał MLS(maximum length sequence), czyli przebieg pseudo losowy przyjmujacy tylko wartosci ±Umax. Nastepnie sygnal z wyjscia karty jest puszczany na jakies tam uzadzenie poczym z uzadzenia kierowany jest spowrotem do wejscia karty muzycznej. I tu zaczyna sie problem bowiem sygnal otrzymany z uządzenia musze tak przeksztalcic wzgledem orginalnego sygnalu generowanego aby otrzymac pojedynczą odpowiedz impulsowa tego uzadzenia i nie wiem jak tego dokonac jakie dzialania, algorytmy tego dokonają. Dodam jeszcze bo moze to bardziej zobrazuje o co mi chodzi iz juz widzialem taki program w działaniu i jezeli sygnal nie jest przepuszczany przez uzadzenie czyli takjakby sygnal orginalny jest przekrztalcany wzgledem orginalnego to otrzymujemy wykres sygnalu w dziedzinie czasu ktory wyglada nastepujaco pierwszy punkt ma jakas wartoasc a nastepne maja wartosc zerową i wlasnie taki efekt chce miedzyinnymi uzyskac. Nastepnie na podstawie tej odpowiedzi impulsowej mam obliczyc widmo sygnalu przy pomocy DFT ale z tym juz rade sobie dam.

    Jezeli cos w moim pytaniu jest niezrozumiale prosze napisac a postaram sie jeszcze bardziej doglebnie przedstawic problem.

    Bardzo prosze o jakies konkretne wskazowki jak taką odpowiedż impulsową otrzymać.
  • REKLAMA
  • Pomocny post
    #2 1473762
    kafar
    Poziom 15  
    Witam
    z twojego postu wynika, ze chcesz uzyskac transmitancje urzadzenia. Metoda MLS pozwala w prosty sposob usyskac odpowiedz impulsowa, mialem gdzies pdfa, jutro jak bede w pracy to go poszukam na kompie. Jezeli jednak masz zamiar na koncu obliczac tranformate Fouriera, to zastanwo sie czy nie zaczac od niej,a pozniej po prostu podziel widma zespolone sygnaly pomiary i otrzymany. Otrzymasz wtedy transmitancje, z ktorej pozniej mozesz otrzymac odpowiedz impulsowa. Najtrudniejsza rzecza w tym wszystkim jest jednal wtedy szybka transformata w obie strony. Z pozostalych algorytmow potrzebowal bys juz tylko dzielenie. Jezeli jednak musi byc stricte metoda MLS to poszukam tego pdfa na pewno gdzies go mam. W jakim srodowisku chcesz zrealizowac ten program?
  • REKLAMA
  • #3 1474947
    Sebastian.
    Poziom 11  
    na samym początku dziekuje za odpowiedz

    niestety musze na samym poczatku policzyc odpowiedz impulsową ze wzgledu iz programik bedzie tez sluzyl do pomiarow za pomoca mikrofonu tzn np sygnal z karty idzie na glosnik nastepnie zostaje on nagrany przez mikrofon i wraca do wejscia karty muz nastepnie sygnal powinien byc poddany dft w celu uzyskania widma jednakze przed obliczeniem dft w sygnale nie moze byc zadnych zanieczyszczen zewnecznych typu odbicia od scian pomieszczenia w ktorym dokonywany jest pomiar a dzieki odpowiedzi impulsowej bede mial czarno na bialym w postaci pików te odbicia od scian i mozna bedzie odciąc niepotrzebna(nieczysta) czesc sygnalu a nastepnie czysty glowny sygnal(impuls) poddac dft.

    Programik pisze w delphim w sumie to tak troche od tylu zaczolem bo napisane juz mam wyliczanie czesci urojonej, zeczywistej, fazy itp itd. no i stanolem na tejze wlasnie odpowiedzi impulsowej pozniej dojdzi jeszcze pare innych zeczy(problemow) do wyliczenia ale narazie staram sie skupic na jednej sprawie takze bardzo bym prosil o tego pdf-a
  • #4 1475278
    jarek_krakow
    Poziom 17  
    1)korelacyjna metoda identyfikacji
    możliwe ,że będzie o tym w książce S.Brandt "analiza danych"
    albo szukaj równanie wienera-hopfa
    2)i tak musisz wiedzieć jakie jest wymuszenie. po co ograniczasz się tylko do wartości +/- umax?
  • REKLAMA
  • Pomocny post
    #5 1475863
    kafar
    Poziom 15  
    1. Bardzo trudno bedzie ci odfiltrowac jakiekolwiek odbicia z odpowiedzi impulsowej, ja bym raczej na to nie liczyl. Wystarcza jezeli mikrofon bedzie znajdowal sie teoretycznie w polu bliskim zrodla. Praktyka pokazuje jednak, ze odleglosc graniczna w pomieszczeniu tak dokladnie realizowana. Dlatego mikrofon nalezy ustawiac jak najblizej. Domyslam sie, ze chodzi ci o pomiar miedzy innymi glosnikow lub zestawow, nie spotkalem sie aby profesjonalne systemu pomiarowe mialy korekcje odbic. Zreszta skad wiesz jaka wysokosc ma miec ten pik no i czy to bedzie pik, odbicie pojawi sie jako nowe zrodlo dzwieku posiadajace rowniez jakis konkretny czas zaniku, czas ten oczywiscie jest czasem poglosu pomieszczenia. No i dochodzi ostatni fakt czasu po jakim mialo by nastapic to odbicie w bardzo malych pomieszczeniach bedzie on bardzo maly.
    2. Tak jak obiecywalem znalazlem pdf dotyczace metody MLS. Wynika z nich ze potrzebna jest transformata Hadawarda. Ja nie wglebialem sie w to dokladnie poniewaz wykorzystywalem szum MLS jako sygnal o odpowiednim widmie.
    3.Tez pisalem w Delphiaczu.Masz juz cos do transformaty Fouriera?Ile masz zamiar uzyc probek? Ja w moim programie pracowalem nawet na 512 tys. ale moj komp P4 2,4 HT troche sie przy tym pocil o i duzo ramu.wtedy.
  • REKLAMA
  • #6 1476436
    Sebastian.
    Poziom 11  
    Dzieki za pomoc

    jezeli chodzi o odfiltrowanie to mnie zle zrozumiales ja nie bede tego filtrowal tylko poprostu dft poddany zostanie sygnal od początku glównego piku do momentu początku pierwszego odbicia a reszta wogule niebedzie brana pod uwage potocznie mowi sie na to bramkowanie.Zalanczm printscr takiej odpowiedzi z programu na ktorym sie wzoruje
    czerwone markery zaznaczaja wlasnie ten sygnal ktory ma byc poddany dft i takze bardzo ladnie widac pierwsze odbicie. I tez dlatego stosuje dft zamiast fft bo ilosc prubek sygnalu bedzie zmienna zalezna od ustawien markerow a przy tym takze nie bedzie tych prubek do obliczania duzo takze dft wystarczy bo przewaznie taki sygna zbramkowany nie trwa dluzej niz 10ms a przy prubkowaniu karty 48kHz bedzie to okolo 480 prubek.


    jedyna wada takiego pomiaru jest mala ilosc informacji dotyczacych niskich czestotliwosci gdzies do okolo 500Hz jednakze do zestrojenia zwrotnicy wzupelnosci wystarcza bo podzialow przewaznie dokonuje sie w granicach 3kHz . a do pomiaru niskich czestotliwosci w celu zestrojenia obudowy wykonuje sie pomiar w polu bliskim glosnika o ktorym mówiles, przewaznie 1cm od centrum membrany i nie podlega on wtedy zadnym czynnoscia bramkujacym. A dlatego jest on wykozystywany tylko do niskich czestotliwosci bo przy wyzszych, przewaznie powyzej 1,5kHz niestety juz nie jest prawdziwy ze wzgledu iz mebrana takiego 16cm glosnika w ruznych swoich punktach jest odpowiedzialna za przetwazanie ruznych czestotliwosci


    w programie mam narazie napisane obliczanie czesci zeczywistej i urojone za pomoca dft. Nastepnie co za tym idzie widmo fazowe i amplitudowe w dziedzinie czestotliwosci. Jest tez juz napisana idft
    oraz na obcych modulach nie pisanych przezemnie ubsluge generacji sygnalu i jego nagrywania puzniej moze jeszcze dopisze ETC(enrgy time curve) ale do tego musze miec najpierw odpowiedz impulsowa


    idft napisalem wczoraj na szybko bo chcialem sprawdzic jaka mi wyjdzie ta odpowiedz impulsowa wedle twoich wskazowek ale chyba cos zle robie lub zle zrozumialem bo nie wychodzi mi to co niby powinno a robie to tak

    obliczam czesc urojona i czesc rzeczywista generowanego sygnalu nastepnie to samo z sygnalem nagranym i przeprowadzam tarnsmitancje czyli??

    RE_transmitacji = RE_generowany / RE_nagrany
    IM_transmitacji = IM_generowany / IM_nagrany

    Nastepnie IDFT z re_transmitacji i im_transmitacji

    ale ta odpowiedz impulsowa nie za bardzo jest taka jak ma byc pewnie cos zle zrozumialem i prosze o wyrozumialosc bo jestem tylko po zawodówie i nie za bardzo mialem mozliwosc taka wiedze przyswoic tak wiec poprawcie mnie.
  • Pomocny post
    #7 1479494
    kafar
    Poziom 15  
    Na poczatku tu respekcic ze po zawodowie zajmujesz sie takimi rzeczami, ja sie nimi zajmowalem na studiach.
    Dzielenie liczb zespolonych to bardziej skomplikowana, postac a+bi, jest niestety nierozerwalana i dzielenie jest błędna, ale bez dodatkowych wywodów, które możesz znalżc w ksiązkach wzory powinny wygladac nastepujaca:
    przy zalozeniu ze dzielenie wyglada nstepujaco:(a+bi)/(c+di), wynik jest nastepujacy Re+Imi, gdzie
    Re=(ac-bd)/(c^2+d^2)
    Im=(cb-ad)/(c^2+d^2)

    Zastanawiales sie czy przy okienkoiwaniu 10 ms wystarczy ci rozdzielczosci?
  • #8 1482643
    Sebastian.
    Poziom 11  
    witam jeszcze raz bardzo Ci dziekuje
    zrobilem tak jak mowisz i na poczatku mialem problemy bo dobrze wychodzily tylko wyniki na podstawie sygnalów generowanych przes funkcje sinusa a inne byly lekko zpaprane powodem byl omylkowo przez cibie zamieniony + z - w jednym miejscu i brak - przed jednym z nawiasów ale poszukalem na sieci o dzieleniu liczb zespolonych i po wprowadzeniu zmian jest wkoncu tak jak ma byc i jestem bardzo zadowolony a bez twoich podpowiedzi bylo by mi bardzo ciezko do tego wszystkiego dojsc

    teraz jeszcze sprubuje zkumac ta transformate hadamarda no i zaimplementowac ja w delphim

    Tak zastanawialem sie nad rozdzielcoscia przy tak malej ilosci probek w sumie to siedzialem nad tym chyba z 2 miesiace i sa na to sposoby ja osobiscie uzyskalem berdzo zadawalajacy mnie efekt identyczny jak w programie na ktorym sie wzoruje
    i przedstawie moze w zalaczniku widmo ktore zostalo obliczone z odpowiedzi impulsowej ktora zalonczylem wczesniej o czasie bramkowania 6,7ms
    efekt jak dla mnie jest dosc dobry (sygnal wygląda byc moze jakby byl wygladzany ale ja jednak zapewniam ze nie jest)


    Aha jeszcze chcialem sie zapytac
    1. w pewnym programie zauwazylem ze wylicza on jakies 3 rodzaje fazy pierwsz nosila nazwe normal zapewne obliczana w sposob standardowy czyli
    atan(im/re) jesli re>0
    atan(im/re)+pi jesli re<0

    ale byla tez faza o nazwie "minimum" i "excess" co one soba reprezentuja i jak sie je wylicza.

    2.w jaki sposob mozna obliczyc czasowe opuznienie odtwazania sygnalu dla danej czestotliwosci podane w ms.
REKLAMA