数字基带调制

Fc=10;   %载频
Fs=40;   %系统采样频率
Fd=1;    %码速率
N=Fs/Fd;
df=10;
numSymb=25;%进行仿真的信息代码个数
M=2;       %进制数
SNRpBit=60;%信噪比
SNR=SNRpBit/log2(M);
seed=[12345 54321];
numPlot=25;
%%===============产生25个二进制随机码
x=randsrc(numSymb,1,[0:M-1]);%产生25个二进制随机码
figure(1)
stem([0:numPlot-1],x(1:numPlot),'bx');
title('二进制随机序列')
xlabel('Time');
ylabel('Amplitude');
%%==============ASK调制
y1=dmod(x,Fc,Fd,Fs,'ask',M,df);
numModPlot=numPlot*Fs;
t=[0:numModPlot-1]./Fs;
figure(2)
plot(t,y1(1:length(t)),'b-');
axis([min(t) max(t) -1.5 1.5]);
title('ASK调制后的信号')
xlabel('Time');
ylabel('Amplitude');
%%==============FSK调制
y2=dmod(x,Fc,Fd,Fs,'fsk',M,df);
numModPlot=numPlot*Fs;
t=[0:numModPlot-1]./Fs;
figure(3)
plot(t,y2(1:length(t)),'b-');
axis([min(t) max(t) -1.5 1.5]);
title('FSK调制后的信号')
xlabel('Time');
ylabel('Amplitude');
%%===============PSK调制
y3=dmod(x,Fc,Fd,Fs,'psk',M,df);
numModPlot=numPlot*Fs;
t=[0:numModPlot-1]./Fs;
figure(4)
plot(t,y3(1:length(t)),'b-');
axis([min(t) max(t) -1.5 1.5]);
title('PSK调制后的信号')
xlabel('Time');
ylabel('Amplitude');
%%===============MSK调制
y4=dmod(x,Fc,Fd,Fs,'msk',M,df);
numModPlot=numPlot*Fs;
t=[0:numModPlot-1]./Fs;
figure(5)
plot(t,y4(1:length(t)),'b-');
axis([min(t) max(t) -1.5 1.5]);
title('MSK调制后的信号')
xlabel('Time');
ylabel('Amplitude');