कुशल दृढ़ संकल्प (आर में)

Question

मैं कनवल्शन की गणना / मूल्यांकन करना चाहता हूं

$$g(x)=\int_D f(x-t) \phi(t) dt,$$

जहां एक घनत्व है और कॉम्पैक्ट समर्थन साथ एक चिकनी कार्य है । समापन बंद-रूप में उपलब्ध नहीं है और मुझे इसे संख्यात्मक रूप से एकीकृत करने की आवश्यकता है। मेरा प्रश्न है: क्या ऐसा करने का एक कुशल तरीका है? मैं इसे आर में लागू करना चाहता हूं, इसलिए, मैं यह देखना चाहूंगा कि क्या कमांड का उपयोग करने से बेहतर तरीका है ।$f$$\phi$$D$ integrate()

Comments

परिस्थितियों के आधार पर मैं आमतौर पर 2 डिब्बे की एक बड़ी शक्ति का विवेक करता हूं, और तेजी से फूरियर रूपांतरण ( ?fft) का उपयोग करता हूं convolve। Fft दृष्टिकोण सेट करने के लिए थोड़ा अधिक काम करता है, लेकिन बेहतर है यदि आपको कई बार कुछ के साथ मनाने की आवश्यकता है। कभी-कभी सही तर्क सेटिंग्स का पता लगाने में कुछ समय लगता है।

— Glen_b -Reinstate मोनिका

@Glen_b धन्यवाद। एक univariate फ़ंक्शन , मुझे लगता है कि प्रत्यक्ष एकीकरण तेज हो सकता है, फिर। $f$

— कुक

आपने एक कुशल तरीका पूछा - एफएफटी वास्तव में तेज है ; इसके लिए बस थोड़े से सेटअप की आवश्यकता होती है (शून्यकरण, शून्य के साथ पैडिंग)।

— Glen_b -Reinstate मोनिका

@Glen_b हाँ, मैं मानता हूँ कि एफएफटी वास्तव में तेज़ है, लेकिन पिछला चरण प्रक्रिया को धीमा कर सकता है। मैं वैसे भी दोनों तरीकों की तुलना करूंगा। धन्यवाद।

— कुक

मुझे convolveइस उद्देश्य के लिए कई बार उपयोग करना याद है । सरल कार्यशील उदाहरण आंकड़े .stackexchange.com / a / 41263 , आँकड़े . stackexchange.com/a/41255 और आँकड़े . stackexchange.com/a/49444 पर दिखाई देते हैं ।

— व्हीबर

Answer 1

क्या आपने इसके लिए समर्पित R पैकेजों पर एक नज़र डाली है? जैसे convolve, https://stat.ethz.ch/R-manual/R-devel/library/stats/html/.net/vv.html

Comments

धन्यवाद। अगर मैं ठीक से समझूं, तो यह पैकेज दो कार्यों के बजाय संख्यात्मक अनुक्रमों को समझाने के लिए है। क्या मैं कुछ भूल रहा हूँ?

— कुक 13

ठीक है, मैंने उस पैकेज का उपयोग नहीं किया है, लेकिन यदि आप आर में कन्वेंशन करने की कोशिश कर रहे हैं, तो मूल रूप से आपके पास दो संख्यात्मक अनुक्रम होंगे जो आपके कार्यों से मान होंगे, क्या वे नहीं? मुझे नहीं लगता कि R आपको फ़ंक्शन की परिभाषा के साथ काम करने देता है।

— jnnavarro 13