matlab-(含教程)可见光通信的matlab仿真，仿真输出误码率空间分布,信噪比空间分布

clc; clear; close all; warning off; addpath(genpath(pwd)); [length,width] = meshgrid(0.5:0.05:5); height = 3; q=1.602176487E-19; Aroom = length.*width; ARX = 50; d = 3; psi = 30; %入射角 psic = 70; %接收机视角 no_of_LEDs = 4; phi_half = 60; pLED = 0.02; Ro = real(RO(phi_half)); % Lambertian radiant intensity, Ro ptx = pLED*Ro; reflectivity_of_floor = 0.15; reflectivity_of_ceiling = 0.8; reflectivity_of_wall = 0.7; floor_area = Aroom; ceiling_area = Aroom; wall_area = length*d; pTotalLEDs = no_of_LEDs*pLED; g = 6; TF = 1; noisebandwidth_factor = 0.562; dataRate = 115200; R=0.55; hlos = HLOS(ARX,d,psi,psic,Ro); % Line-of-sight transfer function, H_LOS pn = 3.3*5.84E-3; prxlos = PRXLOS(no_of_LEDs,ptx,hlos); %Total power of LEDs, PrxLOS amplifier_current = 0.01; Ba = 4.5E6; rhoF = floor_area .* reflectivity_of_floor; rhoC = ceiling_area .* reflectivity_of_ceiling; rhoW = wall_area .* reflectivity_of_wall; rho = (1./Aroom) .* (rhoF + rhoC + rhoW); Ii = rho .* (pTotalLEDs./Aroom); intensity = Ii./(1-rho); received_diffused_power = ARX*intensity; diffused_channel_loss = received_diffused_power/pTotalLEDs; prx = (prxlos+received_diffused_power)*TF*g; PDoutput_current = prx*R; Bn = noisebandwidth_factor*dataRate; omegashot = 2*q*R*(prx+pn)*Bn; omegaamp = amplifier_current^2*Ba; omegatotal = omegashot + omegaamp; %Signal-to-Noise Ratio SNR = (R*prx)^2./omegatotal; SNRdb = 10*log10(SNR); % Error Function Q = (1/2)*erfc((omegatotal./(R*prx.^2)./sqrt(2))); % Bit Error Rate BER = Q*sqrt(SNR); % Plot 3D meshc(length,width,SNRdb) figure(2) meshc(length,width,BER)