前提・実現したいこと
MatlabでPAM変調のサンプルプログラムを実行したいのです.
本のコードをそのまま書いたものなのですが上手く実行することができません.
図の出力が出来ません.
実行手順など詳しくご教授願いたいです.
どなたかご教授よろしくお願いいたします。
■■な機能を実装中に以下のエラーメッセージが発生しました。
該当のソースコード
Matlab
clear all
close all
snr1=0:1:12;
for i=1:length(snr1)
mont_err(i)=mont_4PAM(snr1(i));
end
snr2=0:0.1:12;
for i=1:length(snr2)
snr_linear=10^(snr2(i)/10);
the_err(i)=3/4erfc(sqrt(4/5snr_linear)/sqrt(2));
end
semilogy(snr1,mont_err,'*');hold on
semilogy(snr2,the_err);
title('PAM信号の誤り率')
xlabel('SNR(dB)')
ylabel('BER')
legend('実験値','理論値');
function p=mont_4PAM(snr)
dis=1;
snr_linear=10^(snr/10);
sig=sqrt((5dis^2)/(4snr_linear));
N=10000;
for i=1:N
temp=rand;
if (temp<0.25)
d(i)=0;
elseif (temp<0.5)
d(i)=1;
elseif (temp<0.75)
d(i)=2;
else
d(i)=3;
end
end
num_err=0;
for i=1:N
if(d(i)==0)
r=-3dis+gn1(sig);
elseif(d(i)==1)
r=-dis+gn1(sig);
elseif(d(i)==2)
r=dis+gn1(sig);
else
r=3dis+gn1(sig);
end
if(r<-2dis)
dd=0;
elseif(r<0)
dd=1;
elseif(r<2dis)
dd=2;
else
dd=3;
end
if(dd~=d(i))
num_err=num_err+1;
end
end
p=num_err/N;
end
あなたの回答
tips
プレビュー