You can refer to the following code for the implementation.
% Generate the transmitted signal using the communication beamformer
s_t = wc*(HBR'*ome'*h_RU + h_BU)*(h_RU'*ome*HBR + h_BU')*wc';
% Add noise to the transmitted signal
r_target = awgn(s_t, sigma2);
% Apply the matched filter using the sensing beamformer
template = conj(ws*(HBR'*ome'*h_RU + h_BU));
filtered_signal = conv(r_target, template, 'same');
% Compute the numerator and denominator of the SINR equation
numerator = wc*(HBR'*ome'*h_RU + h_BU)*(h_RU'*ome*HBR + h_BU')*wc';
denominator = ws*(HBR'*ome'*h_RU + h_BU)*(h_RU'*ome*HBR + h_BU')*ws';
% Compute the SINR
SINR = abs(filtered_signal).^2 ./ (denominator + var(filtered_signal) - abs(numerator).^2 ./ (denominator + var(r_target)));
% End of code
You may need to change some parameters to values that are relevant to your model.
