哪位大神可以那帮忙看下matlab程序吗？就是最好给每行注释下？我是小白不太懂这个，再此感谢好心人

Question

哪位大神可以那帮忙看下matlab程序吗？就是最好给每行注释下？我是小白不太懂这个，再此感谢好心人function [p]= cm_sm41(snr_in_dB)
N=10000; d=1; %设定步长
Eav=10*d.^2;
snr=10.^(snr_in_dB/10); %设定信噪比
sgma=sqrt(Eav./(8*snr));
M=16;
for i=1:N
temp=rand; %设定随机数
dsource(i)=1+floor(M*temp);
end;
mapping=[-3*d 3*d; -d 3*d; d 3*d;
3*d 3*d; -3*d d; -d d;
d d ; 3*d d; -3*d -d;
-d -d; d -d; 3*d -d;
-3*d -3*d; -d -3*d; d -3*d;
3*d -3*d]; %设定矩阵
for i=1:N
qam_sig(i,:)=mapping(dsource(i),:);
end;
scatterplot(qam_sig);
for i=1:N
[n(1) n(2)]=gauss(sgma);
r(i,:)=qam_sig(i,:)+n;
end
scatterplot(r); %绘制散点图
numoferr=0;
for i=1:N
for j=1:M
metrics(j)=(r(i,1)-mapping(j,1))^2+(r(i,2)-mapping(j,2))^2;
end
[min_metric decis]=min(metrics);
if (decis~=dsource(i))
numoferr=numoferr+1;
end
end
p=numoferr/(N);

主函数 clear %test1.m
SNRindB1=0:1:10;
SNRindB2=0:0.1:10;
M=16;
k=log2(M);
for i=1:length(SNRindB1)
smld_err_prb(i)= cm_sm41(SNRindB1(i));
echo off;
end;
echo on;
for i=1:length(SNRindB2)
SNR=exp(SNRindB2(i)*log(10)/10); theo_err_prb(i)=3*qfunc(sqrt(4/5*SNR))-9/4*(qfunc(sqrt(4/5*SNR)).^2); echo off;
end;
echo on;
semilogy(SNRindB1,smld_err_prb,'ro-'); %绘图
hold
semilogy(SNRindB2,theo_err_prb);

自定义gauss函数
func

Answer 1

应该生成完整文件路径的算法，.一般表示当前路径下，..一般是表示上一级路径，是构成相对路径的表达方式