İntegral dönüşümleri tersine çevirmek için sayısal yöntemler?

Sayısal olarak aşağıdaki integral dönüşümü ters çevirmeye çalışıyorum:

F (y) = \int_{0}^{\infty} y \exp [- \frac{1}{2} (y^{2} + x^{2})] I_{0} (x y) f (x) d x

$F(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x$

Bu nedenle, belirli bir için yaklaşık değerine ihtiyacım var: $F(y)$ $f(x)$

ve gerçek ve pozitiftir $f(x)$ $F(y)$ (sürekli olasılık dağılımlarıdır)
gerçek ve pozitiftir $x,y$ (büyüklüklerdir)

Şu anda bunu yapmak için çok dağınık ve kaba bir güç yöntemi var:

Bir dizi nokta üzerinde ve eğriyi tanımlarım , kesikli noktaların değerleri rastgele örnekleme ile 'tahmin edilir', bu da tahmin edilen bir . Yazdığım temel genetik algoritma, tahmin edilen ve ölçülen dizisi arasındaki farkı en aza indirir . Daha sonra algoritmanın tersine çevirmeye cevabım olduğu değerini alıyorum. $f(x)$ $F(y)$ $F(y)$ $f(x)$

Bu yaklaşım bazı basit durumlar için oldukça iyi çalışıyor, ancak bana dağınık geliyor ve özellikle sağlam değil.

Birisi bana bu sorunu çözmenin daha iyi yolları hakkında rehberlik verebilir mi?

Zaman ayırdığınız için teşekkürler!

[bilgisayar bilimlerinde yayınlanmıştır]

— CBowman
kaynak

Oldukça basit bir yöntem, fonksiyon uzayında bir temel seçmek ve integral dönüşümü bir matrise dönüştürmek olacaktır. Sonra matrisi ters çevirebilirsiniz.

Eğer ortonormal baz fonksiyonları bazı kümesi gerekir: Matematiksel olarak burada nasıl çalıştığı aşağıda . (Onların da normalize olmadan kurtulmak, ancak bu şekilde açıklamak için daha kolay.) Birimdik araçlarının iç çarpım , nerede $T_i(x)$ $\langle T_i,T_j\rangle = \delta_{ij}$

\begin{matrix} (1) & ⟨ T_{i}, T_{j} ⟩ \equiv \int_{a}^{b} W (x) T_{i} (x) T_{j} (x) d x = δ_{i j} \end{matrix}

$\langle T_i,T_j\rangle \equiv \int_a^b W(x) T_i(x)T_j(x)\,\mathrm{d}x = \delta_{ij}\tag{1}$

Burada bir ağırlık fonksiyonudur. O ve limitler ve seçiminize bağlıdır . Hangi temel fonksiyon setinin kullanılacağını seçtikten sonra, limitleri ve ağırlık fonksiyonunu programınıza sabit olarak kodlayabilirsiniz. $W(x)$ $a$ $b$ $T_i$

Ortonormallik kullanarak, ve gibi herhangi bir işlevi, bu temel işlevlerin doğrusal kombinasyonları olarak ifade edebilirsiniz: $f(x)$ $F(y)$

\begin{aligned} (2) & f (x) & = \sum_{i} c_{i} T_{i} (x) & F (y) & = \sum_{j} C_{j} T_{j} (y) \end{aligned}

$\begin{align} f(x) &= \sum_{i} c_i T_i(x) & F(y) &= \sum_{j} C_j T_j(y) \tag{2} \end{align}$

burada katsayılar şu şekilde hesaplanır:

\begin{aligned} (3) & c_{i} & = ⟨ f, T_{i} ⟩ = \int_{a}^{b} W (x) f (x) T_{i} (x) d x \\ (4) & C_{j} & = ⟨ F, T_{j} ⟩ = \int_{a}^{b} W (y) F (y) T_{j} (y) d y \end{aligned}

$\begin{align} c_i &= \langle f,T_i\rangle = \int_a^b W(x) f(x)T_i(x)\,\mathrm{d}x \tag{3}\\ C_j &= \langle F,T_j\rangle = \int_a^b W(y) F(y)T_j(y)\,\mathrm{d}y \tag{4} \end{align}$

Bu ifadelerin katsayıların tanımlarıyla tutarlı olduğunu doğrulayabilirsiniz, denk. (2) ve ortonormallik, ek. (1).

Şimdi, temel fonksiyonların her birinin dönüşümünü hesaplayın; diyelim . $\tilde{T}_i(y)$

{\tilde{T}}_{i} (y) \equiv \int_{0}^{\infty} y \exp [- \frac{1}{2} (y^{2} + x^{2})] I_{0} (x y) T_{i} (x) d x

$\tilde{T}_i(y) \equiv \int_0^\infty y \exp\biggl[-\frac{1}{2}(y^2 + x^2)\biggr]I_0(xy)T_i(x)\,\mathrm{d}x$

bir işlevdir ve bu nedenle onuveile yaptığımız gibi temel işlevlerin doğrusal bir kombinasyonu olarak ifade edebilirsiniz: $\tilde{T}_i(y)$ $f(x)$ $F(y)$

{\tilde{T}}_{i} (y) = \sum_{k} A_{i k} T_{k} (y)

$\tilde{T}_i(y) = \sum_k A_{ik}T_k(y)$

burada matris elemanlarının bulduğumuz aynı şekilde belirlenir ve , yukarıda: $A_{ik}$ $c_i$ $C_j$

\begin{matrix} (5) & A_{i k} = ⟨ {\tilde{T}}_{i}, T_{k} ⟩ = \int_{a}^{b} W (y) {\tilde{T}}_{i} (y) T_{k} (y) d y \end{matrix}

$A_{ik} = \langle \tilde{T}_i, T_k\rangle = \int_a^b W(y) \tilde{T}_i(y)T_k(y)\,\mathrm{d}y\tag{5}$

$i$ $k$ $T_i(x)$ $W(x)$

$A_{ik}$ $c_i$ $C_j$ $f(x)$ $F(y)$

\begin{aligned} \overset{F (y)}{\overset{⏞}{\sum_{j} C_{j} T_{j} (y)}} & = \int_{0}^{\infty} y \exp [- \frac{1}{2} (y^{2} + x^{2})] I_{0} (x y) \overset{f (x)}{\overset{⏞}{\sum_{i} c_{i} T_{i} (x)}} d x \\ = \sum_{i} c_{i} \int_{0}^{\infty} y \exp [- \frac{1}{2} (y^{2} + x^{2})] I_{0} (x y) T_{i} (x) d x \\ = \sum_{i} c_{i} \sum_{k} A_{i k} T_{k} (y) \end{aligned}

$\begin{align} \overbrace{\sum_{j} C_j T_j(y)}^{F(y)} &= \int_0^\infty y \exp\biggl[-\frac{1}{2}(y^2 + x^2)\biggr]I_0(xy)\overbrace{\sum_{i} c_i T_i(x)}^{f(x)}\,\mathrm{d}x\\ &= \sum_{i} c_i \int_0^\infty y \exp\biggl[-\frac{1}{2}(y^2 + x^2)\biggr]I_0(xy)T_i(x)\,\mathrm{d}x\\ &= \sum_{i} c_i \sum_{k} A_{ik}T_k(y) \end{align}$

$C_{\ell}$ $T_{\ell}$

\begin{aligned} ⟨ (\sum_{j} C_{j} T_{j}), T_{ℓ} ⟩ & = ⟨ (\sum_{i} c_{i} \sum_{k} A_{i k} T_{k}), T_{ℓ} ⟩ \\ \int_{a}^{b} W (y) \sum_{j} C_{j} T_{j} (y) T_{ℓ} (y) d y & = \int_{a}^{b} W (y) \sum_{i} c_{i} \sum_{j} A_{i k} T_{k} (y) T_{ℓ} (y) d y \\ \sum_{j} C_{j} \int_{a}^{b} W (y) T_{j} (y) T_{ℓ} (y) d y & = \sum_{i} c_{i} \sum_{k} A_{i k} \int_{a}^{b} W (y) T_{k} (y) T_{ℓ} (y) d y \\ \sum_{j} C_{j} δ_{j ℓ} & = \sum_{i} c_{i} \sum_{k} A_{i k} δ_{k ℓ} \\ C_{ℓ} & = \sum_{i} c_{i} A_{i ℓ} \end{aligned}

$\begin{align} \left\langle \biggl(\sum_{j} C_j T_j\biggr), T_\ell\right\rangle &= \left\langle \biggl(\sum_{i} c_i \sum_{k} A_{ik}T_k\biggr), T_\ell\right\rangle \\ \int_a^b W(y) \sum_{j} C_j T_j(y)T_\ell(y)\,\mathrm{d}y &= \int_a^b W(y) \sum_{i} c_i \sum_{j} A_{ik}T_k(y)T_\ell(y)\,\mathrm{d}y \\ \sum_{j} C_j \int_a^b W(y) T_j(y)T_\ell(y)\,\mathrm{d}y &= \sum_{i} c_i \sum_{k} A_{ik}\int_a^b W(y) T_k(y)T_\ell(y)\,\mathrm{d}y \\ \sum_{j} C_j \delta_{j\ell} &= \sum_{i} c_i \sum_{k} A_{ik}\delta_{k\ell} \\ C_\ell &= \sum_{i} c_i A_{i\ell} \end{align}$

$\ell$ $C_j$

$C_j$ $c_i$ $A_{ij}$ $c_i$ $A_{ij}$ $C_j$ $F(y)$

$F(y)$ $C_j$

C_{j} = \sum_{i} c_{i} A_{i j}

$C_j = \sum_{i} c_i A_{ij}$

$A$

$i$ $j$ $1$ $N$ $N$ $f(x)$ $T_1(x),\ldots,T_N(x)$ $1$ $M$ $F(y)$ $T_1(y),\ldots,T_M(y)$ $M = N$ $M$ $N$ $N$ $c_i$ $A$ $M\times N$ $A_{11}$ $A_{NM}$

$[-1,1]$ $T_i$ $W(x) = \frac{1}{\sqrt{1 - x^2}}$ $a = -1$ $b = 1$ $\langle T_i, T_j\rangle = \delta_{ij}\pi/2$ $i = j\neq 0$ $\langle T_0, T_0\rangle = \pi$

— David Z
kaynak