पायथन

Question

के साथ एक रैखिक प्रणाली मॉडलिंग मैं पाइथन के साथ बंद-लूप, रैखिक, समय-परिवर्तनीय प्रणाली (विशेष रूप से लॉक पीएलएल सन्निकटन) अनुकरण/मॉडल करना चाहता हूं।

मॉडल के भीतर प्रत्येक उप-ब्लॉक में ज्ञात स्थानांतरण फ़ंक्शन है जो जटिल आवृत्ति H(s) = K/(s * tau + 1) के संदर्भ में दिया जाता है। मॉडल का उपयोग करके, मैं देखना चाहता हूं कि सिस्टम प्रतिक्रिया के साथ-साथ शोर प्रतिक्रिया प्रभावित होती है क्योंकि पैरामीटर (उदा। वीसीओ लाभ) बदल जाते हैं। इसमें बोड प्लॉट्स और रूट-लोकस प्लॉट्स का उपयोग शामिल होगा।

नौकरी पाने के लिए मुझे क्या पाइथन मॉड्यूल की तलाश करनी चाहिए?

Answer 1

मुझे पता है कि यह थोड़ा पुराना है, लेकिन एक खोज ने मुझे इस प्रश्न पर लाया। मैंने इसे एक साथ रखा जब मुझे इसके लिए एक अच्छा मॉड्यूल नहीं मिला। यह ज्यादा नहीं है, लेकिन अगर कोई और खुद को यहां पाता है तो यह अच्छी शुरुआत है।

import matplotlib.pylab as plt 
import numpy as np 
import scipy.signal 

def bode(G,f=np.arange(.01,100,.01)): 
    plt.figure() 
    jw = 2*np.pi*f*1j 
    y = np.polyval(G.num, jw)/np.polyval(G.den, jw) 
    mag = 20.0*np.log10(abs(y)) 
    phase = np.arctan2(y.imag, y.real)*180.0/np.pi % 360 

    plt.subplot(211) 
    #plt.semilogx(jw.imag, mag) 
    plt.semilogx(f,mag) 
    plt.grid() 
    plt.gca().xaxis.grid(True, which='minor') 

    plt.ylabel(r'Magnitude (db)') 

    plt.subplot(212) 
    #plt.semilogx(jw.imag, phase) 
    plt.semilogx(f,phase) 
    plt.grid() 
    plt.gca().xaxis.grid(True, which='minor') 
    plt.ylabel(r'Phase (deg)') 
    plt.yticks(np.arange(0, phase.min()-30, -30)) 

    return mag, phase 

f=scipy.signal.lti([1],[1,1]) 
bode(f) 

संपादित करें: मैं यहाँ वापस हूँ क्योंकि किसी को इस उत्तर upvoted, आप Control Systems Library कोशिश करनी चाहिए। उन्होंने मिलान वाक्यविन्यास और सबकुछ के साथ मैटलैब नियंत्रण प्रणाली टूलबॉक्स का बड़ा हिस्सा लागू किया है।

Answer 2

scipy और numpy मॉड्यूल आपके आवेदन के लिए उपयुक्त हैं।

Answer 3

जैसा कि @ मैट ने कहा, मुझे पता है कि यह पुराना है। लेकिन यह मेरी पहली Google हिट के रूप में आया, इसलिए मैं इसे संपादित करना चाहता था।

आप मॉडल रैखिक, समय परिवर्तक प्रणाली के लिए scipy.signal.lti का उपयोग कर सकते हैं। यह आपको lti.bode देता है।

एच (एस) = (जैसा कि^2 + बीएस + सी)/(डीएस^2 + ईएस + एफ) के रूप में एक आवेग प्रतिक्रिया के लिए, आप h = scipy.signal.lti([A,B,C],[D,E,F]) दर्ज करेंगे। बोड प्लॉट प्राप्त करने के लिए, आप plot(*h.bode()[:2]) करेंगे।

Answer 4

http://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.bode.html के अनुसार अब आप इस का उपयोग कर सकते हैं:

from scipy import signal 
import matplotlib.pyplot as plt 

s1 = signal.lti([1], [1, 1]) 
w, mag, phase = signal.bode(s1) 

plt.figure() 
plt.semilogx(w, mag) # bode magnitude plot 
plt.figure() 
plt.semilogx(w, phase) # bode phase plot 
plt.show() 

Answer 5

मैं इस तरह से बाहर काम कर रहे, अजगर-नियंत्रण का उपयोग बोडे भूखंडों मिला है।

from matplotlib.pyplot import * # Grab MATLAB plotting functions 
from control.matlab import * # MATLAB-like functions 


# Transfer functions for dynamics 
G_modele = tf([1], [13500, 345, 1]); 

# Use state space versions 
G_modele = tf2ss(G_modele); 

figure(1); 
bode(G_modele, dB=1); 
show(); 

कोड मुख्य रूप से इस उदाहरण जो बहुत व्यापक

http://www.cds.caltech.edu/~murray/wiki/index.php/Python-control/Example:_Vertical_takeoff_and_landing_aircraft