krzysiek414
Tłuste zero
Dołączył: 24 Paź 2006
Posty: 52
Przeczytał: 0 tematów
|
 Wysłany: Wto 21:16, 25 Lis 2008 Temat postu: |
|
|
Nie jestem tego pewny oto kod:
| Cytat: | fc=2000;
Tp=0.0001;
fp=1/Tp;
t=[-0.01:Tp:0.0099];
N=length(t);
m=sinc(400*t).^2;
c=cos(2*pi*fc*t);
s=c.*m;
n=0:N-1;
f=(n/N)*fp;
SC=(fft(s)./length(s));
mod=abs(SC);
faz=angle(SC);
M=(fft(m)./length(m));
mod1=abs(M);
%demo dodatkowo demodulacja
c1=cos(2*pi*fc*t);
z=s.*c1;
wn=2*fc/fp;
[B,A]=butter(4,wn);
y=filter(B,A,z);
subplot(2,2,1);
plot(t,m);
title('informacyjny');
subplot(2,2,2);
plot(t,c);
title('nosna');
subplot(2,2,3);
plot(t,s);
title('zmodulowany sygnal czasowa postac');
subplot(2,2,4);
plot(t,y);
title('sygnal po demodulacji');
figure;
subplot(2,1,1);
stem(f(1:100),mod(1:100));
title('widmo amplitudowe zmodulowanego');
subplot(2,1,2);
stem(f(1:100),faz(1:100));
title('faza zmodulowanego');
figure;
stem(f,mod1);
title('widmo informacyjnego'); |
|
|
