x=[1 1 1 0 0 0 1 0];
stem(x);  %繪制基帶信號的波形圖
tone=8;
fc=20;fd=10;fs=500;
y=dmod(x,fc,fd,fs,'psk',2,tone);
figure
plot(y) %繪制調制后的波形圖 
a=fft(x,1024);
f=(0:length(a)-1)*fs/length(a)-fs/2;
figure
plot(f,abs(a));  %調制前的頻譜圖
b=fft(y,1024);
f=(0:length(b)-1)*fs/length(b)-fs/2;
figure
plot(f,abs(b)); %調制后的頻譜圖
y=dmod(x,fc,fd,fs,'psk',2,tone);
ynoisy=y+10*awgn(y,20,0);  % 加入高斯噪聲
figure
plot(ynoisy) %  繪制加入高斯噪聲后調制前的波形圖
f=(0:length(ynoisy)-1)*fs/length(ynoisy)-fs/2;
figure
plot(f,abs(fft(ynoisy))); %繪制加入高斯噪聲調制波形圖