~Nožka

20111004-Tue-stridavy-vykon.py

   1 #!/usr/bin/env python
   2 # -*- coding: utf8 -*-
   3 # Soubor:  20111004-Tue-stridavy-vykon.py
   4 # Datum:   04.10.2011 08:21
   5 # Autor:   Marek Nožka, nozka <@t> spseol <d.t> cz
   6 # Licence: GNU/GPL 
   7 
   8 import pylab as pl
   9 from pylab import pi, e, sin
  10 
  11 print pi
  12 print e    
  13 
  14 f=77      # frekvence
  15 fi = 50   # fázový posun
  16 fi = input("zadej fázový posun: ")
  17 ##############################################
  18 # vytvořím si časovou osu
  19 t = []
  20 krok = 0.2
  21 n = 0 # číslo vzorku
  22 vzorek = 0
  23 while vzorek<=20:
  24     t.append(vzorek)
  25     n = n+1
  26     vzorek = krok * n
  27 ###############################################
  28 # vytvořím vzorky namětí
  29 u = []
  30 n = 0
  31 while n < len(t):
  32     u.append( sin( 2*pi*f*t[n]/1000) )
  33     n = n + 1
  34 
  35 ###############################################
  36 # vytvořím vzorky proudu
  37 i = []
  38 n = 0
  39 while n < len(t):
  40     i.append( sin( 2*pi*f*t[n]/1000 + fi*pi/180) )
  41     n = n + 1
  42     
  43 ###############################################
  44 # vypočet výkonu
  45 p = []
  46 n = 0
  47 while n < len(t):
  48     p.append( u[n] * i[n] )
  49     n = n+1
  50 
  51 pl.plot(t,u)
  52 pl.plot(t,i)
  53 pl.plot(t,p)
  54 pl.grid()
  55 pl.show()

`--> stáhnout