Next: 3.2 LMS prediktor, a
Up: 3. Cvičení 3: Odhad
Previous: 3. Cvičení 3: Odhad
Subsections
Cvičení 3.0: Vygenerujte si signál
, jehož úhlová rychlost
3.1
- A
- je konstantní
. Výsledný signál pak můžeme
popsat jako průmět rovnoměrného pohybu po kružnici.
Pro úhel dostáváme
- B
- se mění lineárně
Výsledný signál pak můžeme
popsat jako průmět rovnoměrně zrychleného pohybu po kružnici.
Pro úhel dostáváme
- C
- se mění harmonicky
kde
Výsledný signál pak můžeme
popsat jako průmět nerovnoměrně zrychleného pohybu po kružnici.
Pro úhel dostáváme
Uvedené parametry zvolte např.
,
,
,
,
Jelikož jsme omezeni na systémy diskrétní v čase, musíme uvažovat
vzorkování. Předpokládejte, že vzorkovací perioda je 1 s.
Dobu vzorování uvažujte shodnou s .
Nejprve zadání parametrů:
W0 = pi/4;
W1 = 3*pi/4;
DW = W1-W0;
N = 500;
n = 1:N;
- A
- Generování harmonického signálu, úhlovou frekvenci si
taky někam uložte (budeme ji potřebovat pro porovnání s odhadem)
Wchp = W0*ones(1,N); % konstantni frekvence W0
x = sin(W0*n);
- B
- Lineární chirp
Wchp = DW/N*n + W0; % linearni chirp
x = sin(DW/N*n.^2/2 + W0*n);
- C
- Harmonický chirp
Wchp = DW/2*sin(2*pi/N*n) + (W0+W1)/2; % harmonicky chirp
x = sin(-DW/2/(2*pi/N)*cos(2*pi/N*n)+(W0+W1)/2*n);
Vykreslení průběhů (u specgram již opravíme časovou osu):
figure(1);
subplot(3,1,1);
plot(1:N,x);
subplot(3,1,2);
wsize = 50; % delka okna pro specgram
[B,F,T] = specgram(x,wsize,[],[],0);
imagesc(T*2+wsize/2,F,20*log10(abs(B))); % napraveni casove osy u specgram
axis('xy');
subplot(3,1,3);
plot(1:N,Wchp/pi); % normalizovana frekvence
Figure 3.1:
Lineární chirp: časový průběh , spektrogram, průběh úhlové frekvence
|
Next: 3.2 LMS prediktor, a
Up: 3. Cvičení 3: Odhad
Previous: 3. Cvičení 3: Odhad
Mirek
2006-12-12