डेटा के लिए चेबीशेव श्रृंखला के कम से कम वर्ग फिट करने के लिए, PythonNumpy में chebyshev.chebfit() का उपयोग करें। यह विधि चेबीशेव गुणांकों को निम्न से उच्च तक क्रमित करती है। यदि y 2-D था, y के कॉलम k में डेटा के लिए गुणांक कॉलम k में हैं। पैरामीटर, x, M नमूना (डेटा) बिंदुओं (x[i], y[i]) के x-निर्देशांक हैं।
पैरामीटर, y नमूना बिंदुओं के y-निर्देशांक हैं। एक ही एक्स-निर्देशांक साझा करने वाले नमूना बिंदुओं के कई सेट (स्वतंत्र रूप से) एक कॉल के साथ पॉलीफिट के लिए फिट हो सकते हैं, जिसमें y 2-डारे के लिए पास किया जाता है जिसमें प्रति कॉलम एक डेटा सेट होता है। पैरामीटर, डिग्री फिटिंगबहुपद की डिग्री है। यदि deg एक एकल पूर्णांक है, तो deg'th पद तक के सभी पदों को फिट में शामिल किया गया है।
पैरामीटर, rcond फिट की सापेक्ष स्थिति संख्या है। सबसे बड़े एकवचन मान के सापेक्ष rcond से छोटे एकवचन मान पर ध्यान नहीं दिया जाएगा। डिफ़ॉल्ट मान लेन (x) * ईपीएस है, जहां ईपीएस प्लेटफॉर्म के फ्लोट प्रकार की सापेक्ष सटीकता है, ज्यादातर मामलों में लगभग 2e-16। पैरामीटर, पूर्ण वापसी मूल्य की प्रकृति का निर्धारण करने वाला स्विच है। जब गलत (डिफ़ॉल्ट) केवल गुणांक लौटाए जाते हैं; जब सच है, एकवचन मूल्य अपघटन से नैदानिक जानकारी भी लौटा दी जाती है।
पैरामीटर, डब्ल्यू वजन हैं। यदि कोई नहीं, तो भार w[i] अवर्गीकृत अवशिष्ट y[i]- y_hat[i] पर x[i] पर लागू होता है। आदर्श रूप से वज़न को चुना जाता है ताकि उत्पादों की त्रुटियों w[i]*y[i] सभी में समान भिन्नता हो। व्युत्क्रम-विचरण भार का उपयोग करते समय, w[i] =1/sigma(y[i]) का उपयोग करें। डिफ़ॉल्ट मान कोई नहीं है।
कदम
सबसे पहले, आवश्यक पुस्तकालय आयात करें -
import numpy as np from numpy.polynomial import chebyshev as C
x-निर्देशांक -
x = np.linspace(-1,1,51)
x-निर्देशांक प्रदर्शित करें -
print("X Co-ordinate...\n",x) y-निर्देशांक -
y = x**3 - x + np.random.randn(len(x))
print("\nY Co-ordinate...\n",y) डेटा के लिए चेबीशेव श्रृंखला के कम से कम वर्ग फिट करने के लिए, PythonNumpy में chebyshev.chebfit() का उपयोग करें। यह विधि चेबीशेव गुणांकों को निम्न से उच्च तक क्रमित करती है। यदि y 2-D था, y के कॉलम k में डेटा के लिए गुणांक कॉलम k -
में हैंc, stats = C.chebfit(x,y,3,full=True)
print("\nResult...\n",c)
print("\nResult...\n",stats) उदाहरण
import numpy as np
from numpy.polynomial import chebyshev as C
# The x-coordinate
x = np.linspace(-1,1,51)
# Display the x-coordinate
print("X Co-ordinate...\n",x)
# The y-coordinate
y = x**3 - x + np.random.randn(len(x))
print("\nY Co-ordinate...\n",y)
# To get the least-squares fit of Chebyshev series to data, use the chebyshev.chebfit() in Python Numpy
c, stats = C.chebfit(x,y,3,full=True)
print("\nResult...\n",c)
print("\nResult...\n",stats) आउटपुट
X Co-ordinate... [-1. -0.96 -0.92 -0.88 -0.84 -0.8 -0.76 -0.72 -0.68 -0.64 -0.6 -0.56 -0.52 -0.48 -0.44 -0.4 -0.36 -0.32 -0.28 -0.24 -0.2 -0.16 -0.12 -0.08 -0.04 0. 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4 0.44 0.48 0.52 0.56 0.6 0.64 0.68 0.72 0.76 0.8 0.84 0.88 0.92 0.96 1. ] Y Co-ordinate... [ 0.04578661 -0.41009751 -0.59839355 -0.86942574 1.19418042 -0.53671972 -0.71247683 0.7118818 -0.09274183 1.46114141 -0.40189463 -0.84017206 -1.00618725 -0.7191427 -0.48005631 -0.28661328 0.58161734 2.62382626 -0.56256678 0.92925678 1.68074305 0.97381262 1.22568804 1.71884192 1.03080843 0.55990935 0.29117168 -0.63718482 0.49396313 -0.32920431 1.16682261 0.90746863 -1.0058597 0.54972961 -1.06040041 -0.11828954 -0.51446299 -1.97932024 -0.91902371 -0.31859977 -1.16124938 0.31809796 0.54940462 -1.11008331 1.04918751 -2.60742632 -1.07242746 0.54313779 -0.3440979 -0.28234564 0.46429998] Result... [-0.12730537 -0.08699379 -0.4211565 0.32959334] Result... [array([43.34485511]), 4, array([1.20144978, 1.19227163, 0.76058422, 0.74600162]), 1.1324274851176597e-14]