Esempio n. 1
0
def plot():
	
	x = 0.01
	x_values = []
	while x <= 10.0:
		x_values.append(x)
		x += 0.01


	# x*sin(x) and fourier series of it:
	y_fourier = []
	for x in x_values:
		f_x = fourier[0]/2

		for n in range (1,16):
			f_x += fourier[n*2 - 1]*cos.number(n*omega*x) + fourier[n*2]*sin.number(n*omega*x)
		
		y_fourier.append(f_x)

	xsinx = []
	for x in x_values:
		xsinx.append( x*sin.number(x) )

	plt.plot(x_values, xsinx, label="g(x) = xsinx" )
	plt.plot(x_values, y_fourier, label="Fourier of g(x)" )




	# song part:

	y_song = []
	for x in x_values:
		f_x = current[0]/2

		for n in range (1,16):
			f_x += current[n*2 - 1]*cos.number(n*omega*x) + current[n*2]*sin.number(n*omega*x)
		
		y_song.append(f_x)

	y_unknown = []
	for x in x_values:
		f_x = unknown[0]/2

		for n in range (1,16):
			f_x += unknown[n*2 - 1]*cos.number(n*omega*x) + unknown[n*2]*sin.number(n*omega*x)
		
		y_unknown.append(f_x)

	plt.plot(x_values, y_song, label="Current song" )
	plt.plot(x_values, y_unknown, label="Unknown song" )
	
	
	plt.legend(bbox_to_anchor=(1.05, 1), loc=1, borderaxespad=0.)
	plt.show()
Esempio n. 2
0
def anFunction(x, n):
    return unknown_f(x) * cos.number(omega * n * x)
Esempio n. 3
0
def anFunction(x, n):
	return x*sin.number(x)*cos.number(omega*n*x)
Esempio n. 4
0
def cosRootFunction(s):
    return s*cos.number(cube_root.cube_root(s))