哪位高人能幫我改壹下這個MATLAB程序，需要運行出結果來，謝謝了！！！程序已給出（但是有錯誤）。

function [ray]=ray_sig(zhongzi)貼進去，保存；

2 新建m文件，把函數

function [p]=wubitlv(snr_in_dB,Lc,W0)貼進去，保存；

3 在command window裏輸入

echo on

Lc=20;

W0=1; %窄帶幹擾信號的中心頻率

SNRindB=0:2:30;

for i=1:length(SNRindB)

smld_err_prb(i)=wubitlv(SNRindB(i),Lc,W0);

echo off

end;

% 繪制仿真結果曲線，誤碼率用對數表示

echo on

semilogy(SNRindB,smld_err_prb);

axis([0,10,0,10^-1]);

xlabel('信噪比Eb/N0');

ylabel('誤碼率Pb');

title('DS系統在窄帶幹擾下的誤碼率曲線')

%%%即可，可以出圖！