python
Dashboard
My Repos
Compilers
Python Online
Node JS Online
Golang Online
codepy
Login
My Repos
Sign Out
Online Python Interpreter
Stop
Run
from sympy import * from numpy import arange import matplotlib.pyplot as plt T1,T2,w =symbols('T1 T2 w',real=True) z=factor((T1*w*I+1)*(T2*w*I+1)*w*I+1) print("Характеристический многочлен замкнутой системы -\n%s"%z) zr=re(z) zm=im(z) print("Действительная часть Re= %s"%zr) print("Мнимая часть Im= %s"%zm) x=[zr.subs({T1:0.01,T2:0.1,w:q}) for q in arange(0,100,0.1)] y=[zm.subs({T1:0.01,T2:0.1,w:q}) for q in arange(0,100,0.1)] plt.axis([-150.0, 10.0, -15.0, 15.0]) plt.plot(x, y) plt.grid(True) plt.show()
Share this code with others
Public
Clear
My Repos
Repo
Lang
Login
Register
Login
Create a free account. No Credit card info required.
I agree with the Codepy
Term of Service
Sign Up