程序如下：
N=input('Enter number of symbols>');
snrdB_min=-3;snrdB_max=8;
snrdB=snrdB_min:1:snrdB_max;
snr=10.^(snrdB/10);
len_snr=length(snrdB);
for j=1:len_snr
sigma=sqrt(1/(2*snr(j)));
Ne=0;
for k=1:N
d=round(rand(1));
x_d=2*d-1;
n_d=sigma*randn(1);
y_d=x_d+n_d;
if y_d>0
d_est=1;
else
d_est=0;
end
if(d_est~=d)
Ne=Ne+1;
end
end
errors(j)=Ne;
ber_sim(j)=errors(j)/N;
end
ber_theor=qfunc(sqrt(2*snr));
semilogy(snrdB,ber_theor,snrdB,ber_sim,'o')
axis([snrdB_min snrdB_max 0.0001 1])
xlabel('SNR in dB')
ylabel('BER')
legend('Theoretical','Simulation')
其中有个定义函数 qfunc为：
function y=qfunc(x)
y=(1/2)*erfc(x/sqrt(2));
要回答的问题：
1、程序中哪行代码相当于BPSK?测试函数X=randn(N),试用plot(x),hist(x)
去体会高斯分布。帮助信息：helprandn,helpplot,helphist
2、程序中哪行代码表示信号经过了噪声信道？噪声强度由哪个参数来控制？
3改变参数N观察结果随N如何变化？并分析原因3、改变参数N,观察结果随N如何变化？并分析原因。
4、为什么曲线是单调降的？
5、当N太少时，比如50，仿真结果的数据点数有什么变化？为什么？真什什
6、所有代码和参数都保持不变，运行多次，每次结果是否相同？为什么？
7、若X~N(0,1)，Y=aX则服从什么分布？程序中哪几行代码体现了该思想？
谢谢啦，各位高手急救啊！！