clearclcEs_No = 0:15;    % 信噪比  dBBER = Es_No;   % 誤比特率，長度與Es_No相同%   循環執行仿真程序for ii = 1:length(Es_No)    EsNo = Es_No(ii);    sim('OFDMModulate');    sim('MC_CDMASim');    load rxerror;    rxerror(1,:) = 0;    error=sum(sum(rxerror));    BER(ii) = error/(1.2e5);endhold ongrid onsemilogy(Es_No,BER,'r+:');save data_for_plot BER Es_No% %%%計算MQAM的理論誤碼率曲線% %%%MQAM的理論誤碼率公式——P=[2*（L-1）/L]*Q[sqrt（[6*log2(L)/(L^2-1)]*(Eb/n0)）]  % M = 16;% L = sqrt(M);% for x = 1:16%     y(x) = x-1;%     EbNo = 10^(y(x)/10);%     P(x) = 2*(L-1)/L*q(sqrt(6*log2(L)/(L^2-1)*EbNo));% end% semilogy(y,P,'k:');%%%計算BPSK的理論誤碼率曲線%%%P=Q（sqrt(2*Eb/no)）for ii = 1:11    Ebno_indB(ii) = ii-1;    Ebno = 10^(Ebno_indB(ii)/10);    Pe(ii) = Q(sqrt(2*Ebno));endsemilogy(Ebno_indB,Pe,'k-');