Mam wygenerowany przebieg modulacji FM i potrzebuję wyznaczyć jego widmo. Jak to poprawnie zrobić?
clear all
close all
Fc = 10; %syg nosnej
fs = 1; %f syg modulujacy
t = 0:0.0001:3; %okres probkwoania
beta=5; %indeks modulacji
Fp = 4*(Fc); %czestotliwosc probkowania
y = sin(2*pi*Fc*t);
subplot(4,1,1)
plot (t, y,'k');
title('Sygnal nosnej');
m=sin(2*pi*fs*t);
hold on
subplot(4,1,2)
plot(t,m,'k')
title('Sygnal modulujacy');
subplot(4,1,3)
x=sin(2*pi*Fc*t + beta*(sin( 2*pi*fs*t)));
hold on
plot(t,x,'k');
title('Sygnal FM');
widmo=abs(fft(x));
%widmo=widmo/500;
subplot(4,1,4);
widmo=abs(fft(y));
widmo=widmo/500;
subplot(2,2,4);
stem(fplot(1:length(fplot)/10), widmo(1:length(fplot)/10));
xlabel('czas');
ylabel('Amplituda');
title('widmo');
grid on;
