Εργασίες 2006. εργασία Γ

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

Για οποιοδήποτε πρόβλημα με register/login, στείλτε email στο contact@thmmy.gr.

[BOBoMASTORAS]

Και κάτι ακόμα. Σχεδίασα το φάσμα του input2.wav κάνοντας διαδοχικά fft και abs. Ο fft μας βγάζει όμως ένα φάσμα που είναι συμμετρικό ως προς το 0. Για να πιάσει το σήμα την περιοχή (0-22ΚΗz) θα χρειαστεί να κάνουμε ένα fftshift. Τότε όμως π.χ. η συχνότητα 11KHz-500Ηz θα έχει την ίδια ένταση με τη συχνότητα 11ΚΗz+500Ηz! Είναι αυτό λογικό; Γιατί σε όλους τους players ο spectrum analyzer βγάζει σαν εικόνα ένα σήμα που δεν είναι συμμετρικό! Και λογικά να το σκεφτεί κανείς, δεν είναι δυνατόν μια χαμηλή συχνότητα να έχει υποχρεωτικά την ίδια ένταση με μια υψηλή! Αφού από διαφορετικές πηγές παράγεται η καθεμία! Ελπίζω να καταλαβαίνετε την απορία μου ...

Όποιος ξέρει κάτι σχετικό θα τον παρακαλούσα να με διαφωτίσει...

όταν κάνεις fft ουσιαστικά βλέπεις τις θετικές αρμονικές. Αυτό που βλέπεις και σου φαίνεται συμμετρικό είναι απλά η επέκταση του σήματος για να είναι περιοδικό που κάνει η fft. οπότε με fft η συχνότητα 0 αντιστοιχεί στο δείγμα 1. Κάνοντας fftshift ουσιαστικά οι συχνοτητες που είναι πάνω από FS/2 πάνε πίσω από τη συχνότητα μηδέν κατοπτρικά και πλέον το πρώτο δείγμα του αντιστοιχεί στη συχνότητα -Fs/2 και το τελευταίο στην fs/2. δηλαδή όντως είναι κεντραρισμένο το φάσμα και αυτό γίνεται για καλύτερη εποπτία.

Αυτό που λές για τους spectrum analyzer μάλλον συμβαίνει δίοτι κάνουν fourier για ήχο που είναι πραγματικό σήμα όπότε ο fourier έιναι άρτιος και δεν έχει νόημα να δείχνουν τις αρνητικές συχνόητες αφού είναι ίδιες. Αντίθετα η fft του matlab κάνει fourier και σε μιγαδικά σήματα οπότε εμφανίζει όλο το φάσμα.

ελπίζω να καταλάβατε τι ήθελα να πω γιατι η ώρα δεν είναι η καλύτερη......

[marauber]

Μάλιστα. Δηλαδή αν κατάλαβα καλά, για να δούμε το σήμα όπως πραγματικά είναι (δηλαδή στο 0 να αντιστοιχεί όντως η συχνότητα 0 κ.ο.κ.) πρέπει να πάρουμε μόνο το πρώτο από τα δύο ίδια κομματάκια του fft. Σωστά; Τότε όμως πως γίνεται το σήμα μας να φαίνεται ότι έχει συχνότητες πολύ πάνω από τα 22KHz; (αυτό που έγραψα για την περιοχή 0-22KΗz στο προηγούμενο μήνυμα είναι άκυρο).
Σίγουρα η ένταση τους δεν είναι τόσο ισχυρή όσο των πρώτων αλλά και πάλι μου φαίνεται σημαντική. Στο αυτί δεν ακούγονται βέβαια αλλά η σπατάλη χώρου πρέπει να είναι σημαντική.

[BOBoMASTORAS]

βασικά εγώ έκανα

samples=wavread...

H=fft(samples)
και μετά plot(abs(fftshift(H)))

ενώ από κάνω είχα βάλει κλίμακα από -22.05 έως 22,05 KHz

[abc]

βασικά εγώ έκανα

samples=wavread...

H=fft(samples)
και μετά plot(abs(fftshift(H)))

ενώ από κάνω είχα βάλει κλίμακα από -22.05 έως 22,05 KHz

Για όσους θέλουν μόνο τα θετικά, κρατήστε μόνο τα στοιχεία του πίνακα που είναι από τη μέση και μετά. Έτσι θα λυθεί και η απορία του Μαραού Μπερέ

Θα βλέπεις μόνο τα θετικά και θα'σαι ήσυχος.

Δλδ:

samples=wavread...

H = fftshift(fft(samples));
H = H(length(H)/2:length(H)-1);
plot(abs(H));

Τι λέει καλέ;
Το'χω?

Αυτή η φωνή του των Helloween σφυρίζει μέσα μου.

[Turambar]

      EAGLE FLY FREE
      Music & Lyrics: Weikath

      People are in a big confusion
      They don't like their constitutions
      Everyday they draw conclusions
      And they're still prepared for war

      Some can say what's ineffective
      Some make up themselves atractive
      Build up things they call protective
      Well you life seems quite bizarre

      In the sky a migthy eagle
      Doesn't care 'bout what's illegal
      On its wings the rainbow's light
      It's flying to eternity

      Eagle fly free
      Let people see
      Just make it your own way
      Leave time behind
      Follow thw sign
      Together we'll fly someday

      Hey, we think so supersonic
      And we made our bombs atomic
      Or the better quite neutronic
      But the poor don't see a dime

      Nowadays the air's polluted
      Ancient people persecuted
      That's what mankind contributed
      To create a better time

      In the sky a migthy eagle
      Doesn't care 'bout what's illegal
      On its wings the rainbow's light
      It's flying to eternity

      Eagle fly free
      Let people see
      Just make it your own way
      Leave time behind
      Follow the sign
      Together we'll fly someday
      Together we'll fly someday
      Forever we'll fly
      Together we'll fly
      Forever we'll fly someday

[elisabeth]

ΝΑ ΚΑΝΩ ΜΙΑ ΑΛΛΗ ΕΡΩΤΗΣΗ;  ΕΓΩ ΕΧΩ ΑΛΛΟ ΠΡΟΒΛΗΜΑ..ΠΩΣ ΒΡΙΣΚΟΥΜΕ ΚΑΝΟΝΙΚΟΠΟΙΗΜΕΝΗ ΣΥΧΝΟΤΗΤΑ;ΑΝ ΓΙΑ ΠΑΡΑΔΕΙΓΜΑ ΧΡΗΣΙΜΟΠΟΙΗΣΩ ΠΑΡΑΘΥΡΟ HAMMING ΓΙΑ  Wn ΜΠΟΡΩ ΝΑ ΒΑΛΩ  2*pi*f/44100? ΜΙΑ ΔΙΑΙΡΕΣΗ ΜΕ ΤΗΝ ΣΥΧΝΟΤΗΤΑ ΔΕΙΓΜΑΤΟΛΗΨΙΑΣ ΑΡΚΕΙ?

Wn is a number between 0 and 1, where 1 corresponds to the Nyquist frequency.
Nyquist frequency = One-half the sampling frequency. Some toolbox functions normalize this value to 1.

Δηλαδή όλα είναι κανονικοποιημένα ως προς το μισό της συχνότητας δειγματοληψίας.

Για τυχαίο f:
2*pi*f/22050

All of the filter design functions operate with normalized frequencies, so they do not require the system sampling rate as an extra input argument. This toolbox uses the convention that unit frequency is the Nyquist frequency, defined as half the sampling frequency. The normalized frequency, therefore, is always in the interval 0  f  1. For a system with a 1000 Hz sampling frequency, 300 Hz is 300/500 = 0.6. To convert normalized frequency to angular frequency around the unit circle, multiply by . To convert normalized frequency back to hertz, multiply by half the sample frequency.

normalized cutoff frequency Wn

άρα σύμφωνα με αυτά δεν πρέπει να είναι Wn=fc/22050?
άσε που  το δοκίμασα με 2*pi*f/22050 και έχει συχνότητα αποκοπής γύρω στα 500 hz ενώ απο το φάσμα φαίνεται ότι θα έπρεπε να είναι 2000-2500 περίπου

[abc]

άρα σύμφωνα με αυτά δεν πρέπει να είναι Wn=fc/22050?
άσε που  το δοκίμασα με 2*pi*f/22050 και έχει συχνότητα αποκοπής γύρω στα 500 hz ενώ απο το φάσμα φαίνεται ότι θα έπρεπε να είναι 2000-2500 περίπου

Άσε Ελίζαμπεθ.

Μέγκα πατάτα, στον παρονομαστή ξέχασα να μετατρέψω το ω σε f.
Θα είναι δλδ:

w=2*pi*f/(pi*fs)=2*f/fs

Επομένως, για fc = 2000 Hz
ωc = .09

Συγγνώμη αν έκαψα κανέναν.