W zeszłym tygodniu pożegnaliśmy Steve’a Jobbs’a, natomiast dzisiaj w świat obiegła wiadomość o śmierci Dennis’a Ritchie.
Mówi się, że nieszczęścia chodzą trójkami więc strach pomyśleć jaką jeszcze wiadomość możemy usłyszeć….
To nie jest dobry rok dla IT
Wysłane przez itguy
w 13 października 2011
Komentarze są wyłączone
Zmiany.
Wysłane przez itguy
w 11 października 2011
3 komentarze
Krótki wpis z kategorii prywatnych 
- Możliwa zmiana pracy – czeka mnie jeszcze tylko jedna rozmowa (finałowa) i może w końcu zacznę zarabiać prawdziwe pieniądze.
- Miałem swój debiut jako operator kamery i drugi reżyser teledysku dla pewnego heavy metalowego zespołu. Po skończeniu prac nad montażem postaram się zaprezentować wyniki.
- Do mojej IT nory dołączyła nowa postać:
Piractwo też zaskakuje
Wysłane przez itguy
w 22 września 2011
Brak komentarzy
Bo HelpDesk to nie ma prawa się nudzić
Wysłane przez itguy
w 21 września 2011
Brak komentarzy
Hałas wentylatora czyli Szybka Transformata Fouriera (FFT) w domowych warunkach
Wysłane przez itguy
w 4 września 2011
Brak komentarzy
Wstęp.
„back to basics” czyli mały powrót do tematyki moich studiów połączony z „za” dużą ilością wolnego czasu na urlopie.
Chłodząc się ostatnio za pomocą małego wiatraczka przypomniały mi się zajęcia na katedrze wibroakustyki poświęcone analizie hałasu oraz drgań od urządzeń mechanicznych – a zwłaszcza wentylatorów. Postanowiłem przeprowadzić mały pomiar a następnie przeprowadzić jego analizę. Oczywiście mówimy tu o totalnej amatorce – przecież nie dysponuję laboratorium – ale idea pozostaje ta sama.
Tor pomiarowy.
Jak już wspominałem nie dysponuję laboratorium (komora bezechowa, czułe mikrofony, oscyloskop) więc wszystko zostało sprowadzone do absolutnego minimum.
Tor pomiarowy wyglądał następująco:
Jak już widać, że za mikrofon posłużyła funkcja dyktafonu z telefonu komórkowego Nokia E5. Nagrałem 5s materiału w najlepszej jakości (jaką tylko oferował telefon).
Następnie rozpakowałem plik mp4 do wav i przeszedłem do analizy.
Obliczenia.
Mając już plik wav pozostała już szybka obróbka w Matlabie.
Pierwszy plik to tylko pomocnicza funkcja wyliczająca kilka podstawowych parametrów sygnału wejściowego.
function [suma,srednia,sr_geom, wariancja] = obliczenia(x) %funkcja licząca sume elementow, srednia i wariancje suma = 0; suma_kw = 0; sr_geom = 1; for i=1:(length(x)) suma = suma + x(i); sr_geom = sr_geom*((abs(x(i)))^(1/length(x))); suma_kw = suma_kw + (x(i))^2; end srednia = suma/(length(x)); wariancja = suma_kw/(length(x) - 1);
Poniższy kod jest odpowiedzialny za całą resztę.
%Analiza danych z pliku wave
%Wczytanie pliku wave
[dane_wav,fs,format]=wavread('Wentylator.wav');
%wyrysowanie danych z pliku
wykresy_sygnalow=figure;
set(gcf,'Name', 'Wykresy wave');
plot(dane_wav);
%Wyliczanie podstawowych parametrow sygnalu
amplituda = (max(dane_wav) - min(dane_wav))/2;
maxim=max(dane_wav);
minim=min(dane_wav);
[suma srednia srednia_geom wariancja]=obliczenia(dane_wav);
wyniki(1) = maxim;
wyniki(2) = minim;
wyniki(3) = suma;
wyniki(4) = srednia;
wyniki(5) = srednia_geom;
wyniki(6) = wariancja;
wyniki(7) = amplituda;
%Wypisanie wynikow na ekran
wyniki'
%Wylicznie Szybkiej Transformaty Fouriera dla zadanego sygnalu
y1 = fft(dane_wav);
fs=2*fs;
for k=1:length(y1),
y1(k) = sqrt(real(y1(k))*real(y1(k))+imag(y1(k))*imag(y1(k)));
end;
zm = round(length(y1)/2);
y1 = y1(1:zm);
clear zm;
elementy = length(y1);
krok = (fs/2)/(elementy-1);
wfs1 = 0:krok:(fs/2);
%Wyrysowanie wykresu w dziedzinie czestotliwosci
wykresy_fft=figure;
set(gcf,'Name', 'Wykresy fft wave');
plot(wfs1,y1);
Pierwszy krok to wczytanie pliku wave do przestrzeni roboczej matlaba. Wykonujemy to za pomocą funkcjiwavread.
Następnie mamy wyświetlenie badanego sygnału na ekranie oraz wypisanie parametrów sygnału. Od linijki 28 następuje wyliczenie transformaty Fouriera-czyli reprezentacja sygnału w dziedzinie częstotliwości (linijka 29 to nic innego tylko wynik twierdzenia o próbkowaniu [tw. Kotielnikowa-Shannona / warunek Nyquista] ).
Na samym końcu następuje wyrysowanie wyników na ekran.
Wyniki.
Sygnał wejściowy z pliku wave:
Analizowany sygnał w dziedzinie częstotliwości:
Oczywiście należy pamiętać jak prymitywny był tor pomiarowy i jak wiele”śmieci” zarejestrowaliśmy (turbulencje związane z obudową wentylatora, drgania poprzez luzy na konstrukcji, zakłócenia z zewnątrz) jednakże na podstawie uzyskanych wyników widać w jakim zakresie wentylator generuje hałas.
P.S.
Plik z nagraniem pracy wentylatora znajduje się tutaj: Wentylator.