www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第8章/li8_9.m

    clear;
M=32;  %设置M的数目
fc=12;  %载波频率
fd=2;  %信号采样速率
fs=56;  %采样频率
x=randint(100,1,M);  %随机信号
y=dmod(x,fc,fd,fs,'qask',M);  %使用M元QASK调制方式
y_noisy=y+0.01*randn(fs/fd*100,1);  %添加高斯白噪声
z=ddemod(y_noisy,fc,fd,fs,'qask',M); %解调信号
s=symerr(x,z);  %计算符号误差率
t=0.1:0.1:10;
subplot(211);
plot(t,x);  %绘制原信号
title('原信号');
subplot(212);
plot(t,z);  %绘制调制解调后的信号
title('经调制解调后的信号');
figure;
modmap('qask',M);  %绘制星座图