clc
clear all
%---------扩频通信过程---------------
%---------只仿真一个频点的情况---------------
%---------初始化---------------
Ts=0.00001; fs=1/Ts;
EndTime=2-Ts; %2s
%---------产生信息序列(双极性不归零码)---------------
Tm=0.25; fm=1/Tm; %码率
[u,time] = gensig('square',2*Tm,EndTime,Ts);
y = 2*(u-0.5);
figure(1);
plot(time,y);
title('信息序列');
xlabel('time (seconds)');
axis([0 2 -2 2]);
%---------FSK调制---------------
T0=0.1; f0=1/T0;
T1=0.2; f1=1/T1;
[u0,time]=gensig('sin',T0,EndTime,Ts);
[u1,time]=gensig('sin',T1,EndTime,Ts);
y0=u0.*sign(-y+1);
y1=u1.*sign(y+1);
SignalFSK=y0+y1; % FSK信号
%---------FSK调制的频谱---------------
nfft=fs+1;
Y = fft(SignalFSK,nfft);
PSignalFSK = Y.* conj(Y)/nfft;
f = fs*(0:nfft/2)/nfft;
figure(2);
plot(f,PSignalFSK(1:nfft/2+1));
title('FSK调制后的频谱');
xlabel('frequency (Hz)');
axis([0 100 -inf inf]);
%---------FSK调制后,低通滤波---------------
cof_low=fir1(64,25/fs);
SignalFSK_l=filter(cof_low,1,SignalFSK);
figure(3)