CODEWITHSUNDEEP
numerical-methodsnumerical-integrationnumerical-computing
pinky
pinky

Reputation: 1

multivariate Fox-H numerical computation

I have some troubles in numerically calculating PDF below: SNR of particular system model

This PDF is the expression of SNR of practical RIS-aided system: https://ieeexplore.ieee.org/document/9387559 (expression A.8)

Authors mentioned some sources of computation of multivariate Fox-H type functions, and one of most popular is presented below: https://github.com/melayadi/multivarFoxH/blob/master/multiFoxH.py

Below there is a definition of multivariate Fox-H function:
Multivariate Fox-H Function

I adjusted test_case according to parameters in the expression, but having errors with executing ( a = (1-kmin)/kmin ).

import numpy as np
from scipy.special import gamma
import math

if __name__ == '__main__':
    xi = 1.5
    kmin = 0.5
    mn = [(0, 2)] + [(0, 2)] + [(0, 1)] + [(0, 1)]    
    pq = [(2, 1)] + [(3, 0)] + [(1, 1)] + [(1, 0)]
    a = [[(0, 1), (0, 1), (0, 2)], [(0, 1/2)], [(0, 1/(2*xi))]]
    b = [[], [(1/2, 1/2)], []]
    c = [(1, -2, 0, -1/(2*xi)), (0, 0, -1, -1)]
    d = [(-1, 0, -1, -1)]

    for i in np.arange(0.1, 1, 0.1):
        z1 = [(kmin/i)**2, 1, (kmin/(1-kmin))**(1/(2*xi))]
        params1 = z1, mn, pq, c, d, a, b
        H1 = np.real(compMultiFoxH(params1, nsubdivisions=20))
        H1 = 1/(math.sqrt(math.pi)*xi*i) * H1

Hope you'll help me to find out why it is not working.

Upvotes: 0

Views: 15

Answers (0)

Related Questions