Question on manually modelling multipath & coherent/correlated signals for DOA estimation

I am trying to model a DOA algorithm in detecting multiple coherent signal sources such as in a multipath environment. I have the following lines of codes which are 3 signals coming in at different DOAs for a 4 element ULA.

clear all
close all
clc

format long
rng('default')

fc = 5.500e9; %center frequency
M = 4; %number of array elements
c = physconst('Lightspeed');
lambda = c/fc; %wavelength
d = 0.5*lambda; %element spacing

%DOA & SNR Model
true_ang = [10 17 45];

ang = true_ang*(pi/180); %DOA in radians
N = length(ang); %number of signals

snrdb = 1; %in dB
snr = db2pow(snrdb) %convert db to pow
sig_pow = 1; %watt
noise_pow = sig_pow/snr %watt

K = 1024; %snapshots

%Array Steering Vector
elementPos = linspace(0,(d*(M-1)),M);
A = steervec(elementPos/lambda,true_ang) %MxN matrix

%Signal and Noise Derivations
s = sqrt(sig_pow)*sign(randn(N,K))
n = sqrt(noise_pow)*sign(randn(M,K))

I need help in defining a multipath scenario where the signals impinging on an array is correlated due to multipath. Can anyone assist? I believe it has something to do with manipulating the 's' and 'n' term of my code?