急求QAM在多径信道下的性能分析。怎么把下面的程序的高斯噪声改成瑞利噪声，且能够显示误码率曲线图？谢谢

%QAM调制误码率计算及仿真%第一步：设置相关参数Fd=2400; %信号数据率Fc=4800; %载波频率Fs=12000; %调制信号的采样频率N=Fs/Fd;modulation='qask';M=4;while M<16; M=M*4;L=sqrt(M); k=log2(M); %正交振幅调制参数 SNRpBit=0:2:30; %设置仿真信噪比范围 SNR=SNRpBit+10*log10(k); symbPerIter=2048; %设置每次迭代符号的次数 iters=3; %迭代次数 expSymErrs=30; %设置预期误差符号数目 numSymbTot=symbPerIter*iters; rand('state',56789*10^10); %设置均匀随机数产生的“随机种子” randn('state',98765*10^5); %设置高斯噪声产生的“随机种子” %计算理论的误码率PE(M),并绘制误码率PE(M)曲线 warning('off','comm : bsolete：mod'); warning('off','comm : bsolete : demod'); errSym=zeros(1,1); ratSym=zeros(1,1); ratio=zeros(1,1); for idx2=1:length(SNR) %仿真不同的SNR值时的误码率PE(M) idx=1; while((idx<=iters)||(sum(errSym)<=expSymErrs)) %第二步：创建信号 msg_orig=randsrc(symbPerIter,1,(0:M-1)); %产生(0,M-1)之间的信息序列 %第三步：数字带通调制 msg_tx=dmod(msg_orig,Fc,Fd,Fs,modulation,M); %第四步：对调制信号添加高斯噪声，由于是带通信号，所以噪声功率只要一半 msg_rx=awgn(msg_tx,SNR(idx2)-10*log10(0.5.*N),'measured',[],'dB'); %第五步：解调带通调制信号 msg_demod=ddemod(msg_rx,Fc,Fd,Fs,modulation,M); %计算本次迭代的误码率PE(M) [errSym(idx),ratSym(idx)]=symerr(msg_orig,msg_demod); idx=idx+1; pause(0.01); end %第六步：计算各次迭代后的平均误码率PE(M) ratio(idx2, : =mean(ratSym); %绘制仿真的误码率PE(M)曲线 semilogy(SNRpBit(1:size(ratio(:,1),1)),ratio(:,1),'r : ' ); legend('仿真的P_{E}(M)',0); pause(0.01); endendhold off;