চিত্র: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টি পাতা এই ফাইল ব্যবহার করে:

অধি-উপাত্ত