चरित्र संबंधी कार्य संगणना और यादृच्छिक चर के अंतर को वास्तव में आसान बना सकते हैं। Mathematica में सांख्यिकीय वितरण के साथ काम करने के लिए बहुत सारे कार्य हैं, जिसमें एक बिलियन को इसके विशेषता फ़ंक्शन में बदलना शामिल है।
मैं इसे दो ठोस उदाहरणों के साथ चित्रित करना चाहता हूं: (1) मान लीजिए कि आप पक्षों की भिन्न संख्याओं के साथ पासा के संग्रह का परिणाम निर्धारित करना चाहते थे , उदाहरण के लिए, दो छह-पक्षीय पासा और एक आठ-पक्षीय मर (यानी) रोल करें , 2d6 + d8 )? या (2) मान लें कि आप दो पासा रोल (जैसे, d6-d6 ) के अंतर को खोजना चाहते थे ?
X fφX(t)fφX(t)=F{f}(t)=E[eitX]
XYfghX+Yh(n)=(f∗g)(n)=∑∞m=−∞f(m)g(n−m)
हम फूरियर ट्रांसफ़र के कन्वेन्शन प्रॉपर्टी का उपयोग कर सकते हैं ताकि इसे और अधिक आसानी से फ़ीचर कार्यों के संदर्भ में देखा जा सके:
φX+Y(t)XYφX(t)φY(t)
यह गणितज्ञ समारोह एक तरफा मरने के लिए विशेषता कार्य करेगा:
MakeCf [s_]: =
मॉड्यूल [{सीएफ़},
Cf: = CharacteristicFunction [डिसक्रीट यूनीफॉर्मडिस्ट्रिएशन [{1, s}],
टी];
सीएफ़]
वितरण की पीएमएफ को इसकी विशेषता फ़ंक्शन से पुनर्प्राप्त किया जा सकता है, क्योंकि फूरियर ट्रांसफॉर्म उल्टे हैं। यहाँ यह करने के लिए गणित कोड है:
RecoverPmf [Cf_]: =
मॉड्यूल [{एफ},
F [y_]: = SeriesCoeffic [Cf /। t -> -I * लॉग [x], {x, 0, y}];
एफ]
हमारे उदाहरण को जारी रखते हुए, F को pmf होने दें जो 2d6 + d8 के परिणाम देता है।
F := RecoverPmf[MakeCf[6]^2 MakeCf[8]]
62⋅8=288S={3,…,20}20=2⋅6+8
में: = एफ / @ रेंज [3, 20]
आउट = {1/288, 1/96, 1/48, 5/144, 5/96, 7/96, 13/144, 5/48, 1/9, 1/9, \
5/48, 13/144, 7/96, 5/96, 5/144, 1/48, 1/96, 1/284}
यदि आप परिणामों की संख्या जानना चाहते हैं जो कि 10 की राशि है, तो गणना करें
में: = 6 ^ 2 8 एफ [10]
आउट = 30
XYfghX−Yh(n)=(f⋆g)(n)=∑∞m=−∞f(m)g(n+m)
हम फूरियर ट्रांसफॉर्म की क्रॉस-सहसंबंधी संपत्ति का उपयोग विशेषता कार्यों के संदर्भ में इसे और अधिक सरल बनाने के लिए कर सकते हैं:
φX−Y(t)X,YφX(t)φY(−t)
इसलिए, D6-d6 के pmf G को खोजने के लिए गणितज्ञ का उपयोग कर:
G := RecoverPmf[MakeCf[6] (MakeCf[6] /. t -> -t)]
62=36S={−5,…,5}−5=1−66−1=5
में: = जी / @ रेंज [-5, 5]
आउट = {1/36, 1/18, 1/12, 1/9, 5/36, 1/6, 5/36, 1/9, 1/12, 1/18, 1/36}।