16.04.2015
Identyfikacja układów
mechatronicznych
Byczek Marta
Projekt #3
Chrupczalski Karol
Grabowska Patrycja
1. Kod programu
clear all
close all
clc
[x,fs] =
wavread('C:\Users\Marcia\Desktop\Byczek_Chrupczalski_Grabowska\Byczek_Ch
rupczalski_Grabowska\whitenoise2012.wav');
y =
wavread('C:\Users\Marcia\Desktop\Byczek_Chrupczalski_Grabowska\Byczek_Ch
rupczalski_Grabowska\whitenoise_rec2012.wav');
X = fft(x);
Y = fft(y);
H = Y./X;
h= ifft(H);
Ha = angle(H);
Ha = unwrap(Ha);
f = linspace(0,fs/2,length(x)/2);
figure(4)
subplot(2,1,1), %ylabel('Amplituda'
semilogx(20*log(abs(H(1:length(H)/2))))
xlabel('Czestotliwosc [Hz]')
ylabel('Amplituda [dB]')
subplot(2,1,2)
semilogx(Ha(1:length(h)/2))
xlabel('Czestotliwosc [Hz]')
ylabel('Faza [rad]')
h=ifft(H);
t=0:(1/fs):((length(h)-1)/fs);
figure(2)
plot(t,h)
xlim([0.01, 0.015])
xlabel('Czas [s]')
ylabel('Amplituda [dB]')
grid
s2=555*340/44100
2. Wykresy Bodego
100
0
-100
-200
-300
0 1 2 3 4 5 6
10 10 10 10 10 10 10
Czestotliwosc [Hz]
1000
0
-1000
-2000
-3000
0 1 2 3 4 5 6
10 10 10 10 10 10 10
Czestotliwosc [Hz]
3. Odpowiedz
impulsowa
0.3
0.2
0.1
0
-0.1
-0.2
-0.3
-0.4
-0.5
0.01 0.0105 0.011 0.0115 0.012 0.0125 0.013 0.0135 0.014 0.0145 0.015
Czas [s]
Amplituda [dB]
Faza [rad]
Amplituda [dB]
4. Odległość mikrofonu od głośnika
Czas dotarcia sygnału do mikrofonu 12.46ms
Prędkość dz
więku 320 m/s
Odległość mikrofonu od głośnika 4.2789m
5. Odpowiedz
częstotliwościowa ma charakter filtra pasmowo-przepustowego. Pasmo
przenoszenia zawiera się w przedziale od 2.5kHz do 16.98kHz.
6. Odpowiedz
impulsowa ma kształt impulsu, po którym następują gasnące oscylacje. Impuls
pojawia się z opóz
nieniem równym 12.47ms, co wynika z odległości jaką ma przebyć fala
dz
więkowa.