#2.2
f = 1000
fp = 80000
Tp = 1/fp
dur = 0.005
A = 1
phase = 0
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)

plt.plot(1000*t,y)
plt.xlabel("Czas[ms]")
plt.ylabel('Amplituda')
plt.title('Ampliduta=1')
plt.grid()
plt.ylim(-5,5)
plt.axis([0,3,-5,5])
plt.show()

#
f = 1000
fp = 80000
Tp = 1/fp
dur = 0.005
A = 2
phase = 0
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)
plt.title('Ampliduta=2')
plt.ylabel('Amplituda')
plt.ylim(-5,5)
plt.plot(1000*t,y)
plt.xlabel("Czas[s]")
plt.grid()
plt.axis([0,3,-5,5])
plt.show()

#
f = 1000
fp = 80000
Tp = 1/fp
dur = 0.005
A = 5
phase = 0
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)
plt.title('Ampliduta=5')
plt.ylabel('Amplituda')
plt.ylim(-5,5)
plt.plot(1000*t,y)
plt.xlabel("Czas[s]")
plt.grid()
plt.axis([0,3,-5,5])
plt.show()

#
#
#
#
#2.2
f = 1000
fp = 2500000
Tp = 1/fp
dur = 0.005
A = 1
phase = 0
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)

plt.plot(1000*t,y)
plt.ylabel('Amplituda')
plt.title('f=1000Hz')
plt.xlabel("Czas[s]")
plt.grid()
plt.ylim(-1,1)
plt.axis([0,3,-1,1])
plt.show()

#
f = 2000
fp = 2500000
Tp = 1/fp
dur = 0.005
A = 1
phase = 0
t = np.arange(0,dur,Tp)
plt.ylabel('Amplituda')
y = A*np.sin(2*np.pi*f*t + phase)
plt.ylim(-1,1)
plt.title('f=2000Hz')
plt.plot(1000*t,y)
plt.xlabel("Czas[s]")
plt.grid()
plt.axis([0,3,-1,1])
plt.show()

#
f = 5000
fp = 2500000
Tp = 1/fp
dur = 0.005
A = 1
phase = 0
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)
plt.title('f=5000Hz')
plt.ylabel('Amplituda')
plt.ylim(-1,1)
plt.plot(1000*t,y)
plt.xlabel("Czas[s]")
plt.axis([0,3,-1,1])
plt.grid()
plt.show()

#
#
#
#
f = 1000
fp = 80000
Tp = 1/fp
dur = 0.005
A = 1
phase = 0
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)
plt.ylabel('Amplituda')
plt.title('ϕ=0 rad')
plt.plot(1000*t,y)
plt.xlabel("Czas[s]")
plt.ylim(-1,1)
plt.axis([0,3,-1,1])
plt.grid()
plt.show()

#
f = 1000
fp = 80000
Tp = 1/fp
dur = 0.005
A = 1
phase = -1*pi/2
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)
plt.ylabel('Amplituda')
plt.title('ϕ=-π/2 rad')
plt.ylim(-1,1)
plt.plot(1000*t,y)
plt.xlabel("Czas[s]")
plt.grid()
plt.axis([0,3,-1,1])
plt.show()

#
f = 1000
fp = 80000
Tp = 1/fp
dur = 0.005
A = 1
phase = 1*pi/2
t = np.arange(0,dur,Tp)
y = A*np.sin(2*np.pi*f*t + phase)
plt.ylabel('Amplituda')
plt.title('ϕ=π/2 rad')
plt.ylim(-1,1)
plt.plot(1000*t,y)
plt.xlabel("Czas[s]")
plt.grid()
plt.axis([0,3,-1,1])
plt.show()
ipd.Audio(y, rate=fp)