clear all nsymbol=100000; T=1; fs=100; ts=1/fs; t=0:ts:T-ts; %时间向量 fc=10; %载波频率 c=sqrt(2/T)*cos(2*pi*fc*t); %载波信号 M=4; %8-PAM graycode=[0 1 3 2]; %Gray编码规则 EsN0=0:15; %信噪比,Es/N0 snr1=10.^(EsN0/10); %信噪比转换为线性值 msg=randi(2,1,nsymbol); %消息数据 msg1=graycode(msg+1); %Gray映射 msgmod=pammod(msg1,M).'; %基带4-PAM调制 tx=msgmod*c; %载波调制 tx1=reshape(tx.',1,length(msgmod)*length(c)); spow=norm(tx1).^2/nsymbol; %求每个符号的平均功率 for indx=1:length(EsN0) sigma=sqrt(spow/(2*snr1(indx))); %根据符号功率求噪声功率 ll=length(tx1); rx=tx1+sigma*randn(1,ll); %加入高斯白噪声 rx1=reshape(rx,length(c),length(msgmod)); y=(c*rx1)/length(c); %相关运算 y1=pamdemod(y,M); %PAM解调 decmsg=graycode(y1+1); [err,ber(indx)]=biterr(msg,decmsg,log2(M)); %误比特率 [err,ser(indx)]=symerr(msg,decmsg); %误符号率 end semilogy(EsN0,ber,'-ko',EsN0,ser,'-k*',EsN0,1.5*qfunc(sqrt(0.4*snr1))); title('4-PAM载波调制信号在AWGN信道下的性能') xlabel('Es/N0');ylabel('误比特率和误符号率') legend('误比特率','误符号率','理论误符号率')上述程序怎么改为8PAM,请列出程序和理论误符号率公式
时间: 2024-01-04 21:04:06 浏览: 221
以下是修改后的程序和理论误符号率公式:
```matlab
clear all
nsymbol=100000;
T=1;
fs=100;
ts=1/fs;
t=0:ts:T-ts; %时间向量
fc=10; %载波频率
c=sqrt(2/T)*cos(2*pi*fc*t); %载波信号
M=8; %8-PAM
graycode=[0 1 3 2 6 7 5 4]; %Gray编码规则
EsN0=0:15; %信噪比,Es/N0
snr1=10.^(EsN0/10); %信噪比转换为线性值
msg=randi(2,1,nsymbol); %消息数据
msg1=graycode(msg+1); %Gray映射
msgmod=pammod(msg1,M).'; %基带8-PAM调制
tx=msgmod*c; %载波调制
tx1=reshape(tx.',1,length(msgmod)*length(c));
spow=norm(tx1).^2/nsymbol; %求每个符号的平均功率
for indx=1:length(EsN0)
sigma=sqrt(spow/(3*snr1(indx))); %根据符号功率求噪声功率
ll=length(tx1);
rx=tx1+sigma*randn(1,ll); %加入高斯白噪声
rx1=reshape(rx,length(c),length(msgmod));
y=(c*rx1)/length(c); %相关运算
y1=pamdemod(y,M); %PAM解调
decmsg=graycode(y1+1);
[err,ber(indx)]=biterr(msg,decmsg,log2(M));
%误比特率
[err,ser(indx)]=symerr(msg,decmsg); %误符号率
end
semilogy(EsN0,ber,'-ko',EsN0,ser,'-k*',EsN0,3/7*qfunc(sqrt(0.4/3*snr1)));
title('8-PAM载波调制信号在AWGN信道下的性能')
xlabel('Es/N0');
ylabel('误比特率和误符号率')
legend('误比特率','误符号率','理论误符号率')
```
其中,理论误符号率公式为:
$$P_{s}=\frac{1}{log_2(M)}\sum_{i=0}^{M-1}\sum_{j=0,j\neq i}^{M-1}Q\left(\sqrt{\frac{3}{M^2-1}\frac{d_{ij}^2}{N_0}}\right)$$
其中,$M$为调制阶数,$d_{ij}$为第 $i$ 个符号与第 $j$ 个符号之间的欧几里得距离,$N_0$为单侧带宽为 $B$ 的白噪声功率谱密度。在本程序中,$B=2/T$,因此 $N_0=0.4/3$。$Q(x)$ 为高斯 Q 函数。
阅读全文