Cvičení 3.0: Vygenerujte si signál
, jehož úhlová rychlost
3.1
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
.
W0 = pi/4; W1 = 3*pi/4; DW = W1-W0; N = 500; n = 1:N;
Wchp = W0*ones(1,N); % konstantni frekvence W0 x = sin(W0*n);
Wchp = DW/N*n + W0; % linearni chirp x = sin(DW/N*n.^2/2 + W0*n);
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);
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