Tyto stránky již nejsou udržovány. Obsah je postupně přesouván/aktualizován na adrese chytrosti.marrek.cz.
1575/14246
1 #!/usr/bin/env python
2 # -*- coding: utf8 -*-
3 # Soubor: 20120514-Mon-wien.py
4 # Datum: 14.05.2012 08:58
5 # Autor: Marek Nožka, nozka <@t> spseol <d.t> cz
6 # Licence: GNU/GPL
7 # Úloha: Přenosová ch. Wienova článku
8 ###################################################3
9 import matplotlib as mpl
10 mpl.use('Qt4Agg')
11 #mpl.rc('font',**{'family':'sans-serif','sans-serif':['Helvetica']})
12 ## for Palatino and other serif fonts use:
13 #mpl.rc('font',**{'family':'serif','serif':['Palatino']})
14 #mpl.rc('lines', linewidth=2, color='r')
15
16 mpl.rcParams['text.usetex']=True
17 mpl.rcParams['text.latex.unicode']=True
18
19 import numpy as num
20 import scipy as sci
21 import pylab as lab
22 from pylab import pi
23
24 R = 1e3
25 C = 100e-9
26
27 # frekvenční osa
28 f = lab.linspace(10,100e3,7e3)
29 w = 2*pi*f
30
31 # Impedance
32 Z1 = R + 1./(1j*w*C)
33
34 Y2 = 1./R + 1j*w*C
35 Z2 = 1./Y2
36
37 Au = Z2/(Z1+Z2)
38
39 # amplitudový přenos
40 lab.figure()
41 lab.semilogx(f,20*num.log10(abs(Au)))
42 lab.grid()
43 # fázový přenos
44 lab.figure()
45 lab.semilogx(f,lab.angle(Au))
46 lab.grid()
47
48 # dohromady
49 lab.figure()
50 lab.subplot(211)
51 lab.semilogx(f,20*num.log10(abs(Au)))
52 lab.grid()
53 lab.subplot(212)
54 lab.semilogx(f,lab.angle(Au))
55 lab.grid()
56
57
58 lab.show()