ostatnio na wykladzie mialam cos takiego..nie bardzo rozumiem o co tu chodzi..wykladowca pokazal nam program w matalbie na wygenerowanie syganalu y(t)=Acos(wt)..mamy podane parametry generowania syganlu..po uruchomieniu pojawiaja sie 4 wykresy..moglby ktos wytlumaczyc mi o co dokladnie tam biega..???jak rozumiec te wykresy?? dziekuje bardzo za pomoc..
clear all
clc
Gr=8;
A1=Gr*50
Usk1=A1/(sqrt(2))
Fi1=2*pi/3;
n=128;
Ts=1/128;
k= 0:n-1;
fo=8; % czestotliwosc podstawowa
D=50; % skladowa stala
fprintf('skladowa podstawowa %g Hz\n', fo)
y=D+A1*cos(2*pi*fo*k*Ts+Fi1)
mocty=y*y'/n
uskt=sqrt(mocty)
figure (1)
plot(y(1:n));
xlabel('czas [s]'),ylabel('wartosc chwilowa')
title('wykres sygnalu y(t)')
% widmo=fft(y)/n; %algorytm FFT z Signal Proccessing
for r=1:n
for p=1:n;
wsp(r,p)=cos(2*pi*(r-1)*(p-1)/n)-j*sin(2*pi*(r-1)*(p-1)/n);
end
end
widmo=wsp*(y')/n;
modul=abs(widmo);
mocfy=modul'*modul
uskfy=sqrt(mocfy)
freq=k/(n*Ts);
figure (2)
stem(freq(1:n), modul(1:n))
title('modul widma')
xlabel('czestotliwosc [Hx]'), ylabel('modul')
ph=angle(widmo);
figure (3)
stem(freq(1:n),ph(1:n))
title ('kat fazowy,wykres zawiera tez bledne wartosci')
xlabel('czestotliwosc [Hz]'), ylabel('kat fazowy')
for kk=1:n;
if modul(kk)>0.001
ph(kk)=(180/pi)*angle(widmo(kk));
elseif modul(kk) < 0.001
end
end
figure (4)
stem(freq(1:n),ph(1:n))
title('poprawiony wykres kata fazowego w stopniach')
xlabel('nr probki widma'), ylabel('kat fazowy')