চিত্র:Polynomial fitting.png
ডিপথট থেকে
Polynomial_fitting.png (৭০৩ × ৫৬৪ পিক্সেল, ফাইলের আকার: ৭১ কিলোবাইট, এমআইএমই ধরন: image/png)
This plot has been created using python. The codes are as follows:
import numpy as np, matplotlib.pyplot as pl
from sys import argv
def main(N,sig,order):
print "Fitting %sth order polynomial to data with N=%s and sigma=%s"% (str(order),str(N),str(sig))
n=np.arange(-1,1,1./(N/2))
x=1./(1.+n**2.)
w=np.random.randn(len(x))*sig**2
f=x+w
pl.plot(n,x,'r',label='$(1+t^2)^{-1}$')
pl.plot(n,f,'b',label='$(1+t^2)^{-1}+w(t)$')
H = np.zeros((N,order+1))
p = []
for i in range(N):
for j in range(order+1):
H[i,j] = n[i]**j
H=np.matrix(H)
f=np.transpose(np.matrix(f))
p = np.linalg.inv(np.transpose(H)*H)*np.transpose(H)*f
f_est=np.zeros((N))
for i in range(N):
for j in range(len(p)): f_est[i] += p[j]*n[i]**j
pl.plot(n,np.transpose(f_est),'g',label='Polynomial fit')
pl.ylabel('Actual function, noisy data, and polynomial fit')
pl.xlabel('$t$')
pl.title("Fitting a "+str(order)+"th order polynomial to a data with $\sigma$="+str(sig))
leg = pl.legend(loc='upper right', fancybox=True)
leg.get_frame().set_alpha(0.5)
pl.savefig('polynomial_fitting.png', bbox_inches='tight', transparent=True)
if __name__=='__main__':
main(int(argv[1]),float(argv[2]),int(argv[3]))
ফাইলের ইতিহাস
যেকোন তারিখ/সময়ে ক্লিক করে দেখুন ফাইলটি তখন কী অবস্থায় ছিল।
তারিখ/সময় | সংক্ষেপচিত্র | মাত্রা | ব্যবহারকারী | মন্তব্য | |
---|---|---|---|---|---|
বর্তমান | ২০:১৪, ১৫ ডিসেম্বর ২০১৩ | ৭০৩ × ৫৬৪ (৭১ কিলোবাইট) | Muhammad2017 (আলোচনা | অবদান) | {{pyplots}} |
আপনি এই ফাইলটি প্রতিস্থাপন করতে পারবেন না।
ফাইলের ব্যবহার
নিম্নলিখিত 2টি পাতা এই ফাইল ব্যবহার করে: