Jest tutaj jakiś specjalista od Matlaba?

14.01.2009

13:56

[1]

QrKo_ [ ]

Jest tutaj jakiś specjalista od Matlaba?

Mam problem i potrzebowalbym malej pomocy :)

Otóż mam wygenerowany sygnał składajacy się z kilku harmonicznych + dodanego na koncu szumu, o takich parametrach:
A - dowolne, dla kazdej harmonicznej inne, zakres 10-100
fi - dowolne, dla kazdej harmonicznej inne, zakres 0-90
w - poczatkowe 1,21kHz = 1210, z kazda nastepna mnozone *2, *3, *4 etc...

Po przepuszczeniu calosci przez filtr FFT, otrzymuje cos takiego jak na obrazku (9 slupkow, po jednym dla kazdej harmonicznej + odbicie lustrzane) ->
Moim zadaniem jest potraktowac sygnal 3 filtracjami:
1. Tak aby odfiltrowac harmoniczna ktora reprezentuje najnizszy slupek.
2. Tak aby odfiltrowac harmoniczna ktora reprezentuje najwyzszy slupek.
3. Tak aby odfiltrowac harmoniczne ktore reprezentują co drugie słupki (parzyste albo nieparzyste, obojetne).

Niestety nie mam pojecia jak sie za to zabrac, jezeli ktos moglby w jakikolwiek pomoc bylbym dozgonnie wdzieczny.

Kod calego sygnalu:

fs=25000;
t=0:1/fs:5;
w=1210;
A0=10;
fi0=pi/2;
y0=A0*sin(2*pi*w*t+fi0);

w1=w*2;
A1=90;
fi1=pi/6;
y1=A1*sin(2*pi*w1*t+fi1);

w2=w*3;
A2=55;
fi2=pi/3;
y2=A1*sin(2*pi*w2*t+fi2);

w3=w*4;
A3=25;
fi3=pi/7;
y3=A3*sin(2*pi*w3*t+fi3);

w4=w*5;
A4=75;
fi4=pi/4;
y4=A4*sin(2*pi*w4*t+fi4);

w5=w*6;
A5=85;
fi5=pi/10;
y5=A5*sin(2*pi*w5*t+fi5);

w6=w*7;
A6=35;
fi6=pi/6;
y6=A6*sin(2*pi*w6*t+fi6);

w7=w*8;
A7=60;
fi7=pi/5;
y7=A7*sin(2*pi*w7*t+fi7);

w8=w*9;
A8=15;
fi8=pi/9;
y8=A8*sin(2*pi*w8*t+fi8);

power=0.5;
ran=sqrt(power)*randn(size(t));

h=y0+y1+y2+y3+y4+y5+y6+y7+y8+ran;

plot(t,h);

Hakim [ Generaďż˝ ]

Nie mam pojęcia o czym piszesz, ale w międzyczasie czekając na odpowiedź tutaj, polecam zapytać też na elektroda.pl. Tam raczej są tęgie głowy od Matlaba.