t_start=clock
for snr=0:5:30
    runtime=10000*50/480000;
    sim('ofdm1');
    Pe(snr+5)=bee(1);
end
figure;                                                     %繪制誤碼性能曲線
semilogy(snr,Pe,'r-');
title('OFDM導引信道估計系統碼性能曲線');
xlabel('Eb/No(dB)');
ylabel('BER');
grid on;