m06_qpsk.m

matlab的经典讲义

% 第一步：设置参数
M = 8; % Number of symbols in alphabet
len = 10000; % Number of symbols in the original message
Fd = 1; % 频率分辨力
Fs = 3; % 调制信号的采样频率

% 第二步：创建输入信号
signal = randint(len,1,M); % 这里输入一个随机信号

% 第三步：调制信号
modsignal(:,1) = dmodce(signal,Fd,Fs,'qask',M);
inphase = [-3:2:3 -3:2:3];
quad = [ones(1,4), -1*ones(1,4)];
modsignal(:,2) = dmodce(signal,Fd,Fs,'qask/arb',inphase,quad);

%第四步：添加噪声
noisy = modsignal+.5*randn(len*Fs/Fd,2)+j*.5*randn(len*Fs/Fd,2);

% 第五步：调制恢复的信号
newsignal(:,1) = ddemodce(noisy(:,1),Fd,Fs,'qask',M);
newsignal(:,2) = ddemodce(noisy(:,2),Fd,Fs,'qask/arb',inphase,quad);

%第六步：计算和显示比特误差率
[num,rate] = biterr(newsignal,signal); 
disp('Bit error rates for the two constellations used here')
disp('----------------------------------------------------')
disp(['Gray code constellation:     ', num2str(rate(1))])
disp(['Non-Gray code constellation: ', num2str(rate(2))])

% 第七步：绘制图形
modmap('qask',M);