~Nožka

20120514-Mon-wien.py

   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()

`--> stáhnout