Skip to content
Snippets Groups Projects
Select Git revision
  • develop
  • master
  • v2.3.4
3 results

board_xml_led.c

Blame
    • Forked from PSP Fanclub / AVRSimV2
    Source project has a limited visibility.
    Code owners
    Assign users and groups as approvers for specific file changes. Learn more.
    verstaerker.py 4.15 KiB
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
    100
    101
    102
    103
    104
    105
    106
    107
    108
    109
    110
    111
    112
    113
    114
    115
    116
    117
    118
    119
    120
    121
    122
    123
    124
    125
    126
    127 
    

    

    

    import numpy as np
import matplotlib.pyplot as plt
import scipy.odr as sodr
import scipy.stats as scs
from my_methods import *
import scipy.optimize

delta_U = 0.00027#20/(2**12) /np.sqrt(12)
errorscale = 1

print('\n########## Verstaerker ohne C3 ##########')
verstaerker = parseCSV('verstaerker_ohne_c3')
t = verstaerker[:, 0]
U_a = verstaerker[:, 1]
U_b = verstaerker[:, 2]

sigma_t = (t[1] - t[0])/np.sqrt(12) * np.ones_like(t)
sigma_U_a = delta_U * np.ones_like(U_a)
sigma_U_b = delta_U * np.ones_like(U_b)

errorscale = 50
plt.errorbar(t, U_a, xerr = sigma_t, yerr = sigma_U_a * errorscale, label = 'U_A', linestyle = '')
plt.errorbar(t, U_b, xerr = sigma_t, yerr = sigma_U_b * errorscale, label = 'U_B', linestyle = '')
plt.ylabel(r'$U [V]$')
plt.xlabel(r'$t [ms]$')
plt.title(r'Messdaten Gleichrichter')
plt.grid()
plt.legend()
plt.tight_layout()
plt.savefig('output/verstaerker_oszi_ohne_c3.pdf')
plt.clf()

print('########## fitting sine ##########')
# odr does not converge here :(
def sine2(x, a,b,c,d):
    return sine([a,b,c,d], x)

popt_a, pcov_a = scipy.optimize.curve_fit(sine2, t, U_a, sigma = sigma_U_a, p0 = [2*np.pi, np.pi/2, 0.5, 0])
popt_b, pcov_b = scipy.optimize.curve_fit(sine2, t, U_b, sigma = sigma_U_b, p0 = [2*np.pi, -np.pi/2, 2, 0])
chi_a, p_a = chi_sq_dof_p(t, U_a, xerr = sigma_t, yerr = sigma_U_a, beta = popt_a, model = sine)
chi_b, p_b = chi_sq_dof_p(t, U_b, xerr = sigma_t, yerr = sigma_U_b, beta = popt_b, model = sine)
sigma_a = np.sqrt(pcov_a.diagonal())
sigma_b = np.sqrt(pcov_b.diagonal())

print('## U_a ##')
print(popt_a)
print(sigma_a)
print(chi_a, p_a)
print('## U_b ##')
print(popt_b)
print(sigma_b)
print(chi_b, p_b)

beta = popt_b[2] / popt_a[2]
sigma_beta = np.sqrt((beta*sigma_b[2]/popt_b[2])**2 + (beta*sigma_a[2]/popt_a[2])**2)

print(f'\n\t Verstärkung: {round(beta, sigma_beta)}')

errorscale = 1
plt.errorbar(t, U_a, xerr = sigma_t, yerr = sigma_U_a * errorscale, label = 'U_A', linestyle = '')
plt.errorbar(t, U_b, xerr = sigma_t, yerr = sigma_U_b * errorscale, label = 'U_B', linestyle = '')
plt.plot(t, sine(popt_a, t), label = 'a')
plt.plot(t, sine(popt_b, t), label = 'b')
plt.legend()
plt.grid()
plt.savefig('./output/verstaerker_ohne_c3_fit.pdf')
plt.cla()

print('\n\n########## Verstaerker mit C3 ##########')
verstaerker = parseCSV('verstaerker_mit_c3')
t = verstaerker[:, 0]
U_a = verstaerker[:, 1]
U_b = verstaerker[:, 2]

sigma_t = (t[1] - t[0])/np.sqrt(12) * np.ones_like(t)
sigma_U_a = delta_U * np.ones_like(U_a)
sigma_U_b = delta_U * np.ones_like(U_b)

errorscale = 50
plt.errorbar(t, U_a, xerr = sigma_t, yerr = sigma_U_a * errorscale, label = 'U_A', linestyle = '')
plt.errorbar(t, U_b, xerr = sigma_t, yerr = sigma_U_b * errorscale, label = 'U_B', linestyle = '')
plt.ylabel(r'$U [V]$')
plt.xlabel(r'$t [ms]$')
plt.title(r'Messdaten Gleichrichter')
plt.grid()
plt.legend()
plt.tight_layout()
plt.savefig('output/verstaerker_oszi_ohne_c3.pdf')
plt.clf()

print('########## fitting sine ##########')
# odr does not converge here :(
def sine2(x, a,b,c,d):
    return sine([a,b,c,d], x)

popt_a, pcov_a = scipy.optimize.curve_fit(sine2, t, U_a, sigma = sigma_U_a, p0 = [2*np.pi, np.pi/2, 0.05, 0])
popt_b, pcov_b = scipy.optimize.curve_fit(sine2, t, U_b, sigma = sigma_U_b, p0 = [2*np.pi, -np.pi/2, 3.55, -0.3])
chi_a, p_a = chi_sq_dof_p(t, U_a, xerr = sigma_t, yerr = sigma_U_a, beta = popt_a, model = sine)
chi_b, p_b = chi_sq_dof_p(t, U_b, xerr = sigma_t, yerr = sigma_U_b, beta = popt_b, model = sine)
sigma_a = np.sqrt(pcov_a.diagonal())
sigma_b = np.sqrt(pcov_b.diagonal())

print('## U_a ##')
print(popt_a)
print(sigma_a)
print(chi_a, p_a)
print('## U_b ##')
print(popt_b)
print(sigma_b)
print(chi_b, p_b)

beta = popt_b[2] / popt_a[2]
sigma_beta = np.sqrt((beta*sigma_b[2]/popt_b[2])**2 + (beta*sigma_a[2]/popt_a[2])**2)

print(f'\n\t Verstärkung: {round(beta, sigma_beta)}')

errorscale = 1
plt.errorbar(t, U_a, xerr = sigma_t, yerr = sigma_U_a * errorscale, label = 'U_A', linestyle = '')
plt.errorbar(t, U_b, xerr = sigma_t, yerr = sigma_U_b * errorscale, label = 'U_B', linestyle = '')
plt.plot(t, sine(popt_a, t), label = 'a')
plt.plot(t, sine(popt_b, t), label = 'b')
plt.legend()
plt.grid()
plt.savefig('./output/verstaerker_mit_c3_fit.pdf')
plt.show()
plt.cla()