% % 仿真主程序
% clear all;
% error = [0 0 0 0 0 0;
% 0 0 0 0 0 0;
% 0 0 0 0 0 0;
% 0 0 0 0 0 0;
% 0 0 0 0 0 0;
% 0 0 0 0 0 0];
% for j=1:6;
% for k=1:6;
% for i=1:1000;
% pr_part; % 产生原始信号
% interp_part; % 插值
% channel_part; % 加信道H(t),考虑可见光特性
% noise_part; % 加噪声
% final_part; % 判决
% error_bit = sum(nonzeros(abs(final-pr_st))); % 误码个数
% error(j,k) = error(j,k) + error_bit;
% if mod(i,10)==0;
% disp(['(',num2str(j),',',num2str(k),')',' ',num2str(i/10),'%']);
% end
% end
% ber(j,k) = error(j,k)./100000 %误码率
% end
% end
%
% % 注:实际仿真时,需要增加序列长度(1000以上),改变不同信噪比,求出误码率曲线。
clear all;
pr_part; % 产生原始信号
interp_part; % 插值
channel_part; % 加信道H(t),考虑可见光特性
noise_part; % 加噪声
final_part; % 判决
error_bit = sum(nonzeros(abs(final-pr_st))) % 误码个数
ber = error_bit./length(final)