मैं सर्किट विश्लेषण के लिए MATLAB का थोड़ा बहुत उपयोग करता हूं। कभी-कभी मैं इसे मसाले के लिए पसंद करता हूं, अन्य बार मैं मसाला पसंद करता हूं, मेरे मूड और आवश्यकताओं पर निर्भर करता है।
ये निम्नलिखित चरण हैं:
- 1: सर्किट के लैपल्स ट्रांसफॉर्म को लें
- 2: हस्तांतरण समारोह प्राप्त करते हैं
- 3: MATLAB फ़ंक्शन का उपयोग करके प्लॉट / विश्लेषण करें। बोडे, आवेग, freqresp और इतने पर।
मेरे द्वारा पाया जाने वाला सबसे कठिन हिस्सा लाप्लास ट्रांसफ़ॉर्म को लेना और आपके ट्रांसफ़र फ़ंक्शन समीकरण को प्राप्त करना है।
इंटरनेट पर लाप्लास लेने के कई उदाहरण और पाठ्य पुस्तकें हैं। संक्षेप में यहाँ उद्देश्य के रूप में समीकरण प्राप्त करना है
एच( रों ) = एक रों2+ b s + cघरों2+ ई s + च
एसीघच
ऐसा करने के लिए आप सभी निष्क्रिय तत्वों को जटिल बाधाओं में परिवर्तित करते हैं। यानी कि
- सी = 1 / एससी
- आर = आर
- एल = एसएल
अगला Vout / Vin के रूप में आपके सर्किट के लिए एक समीकरण प्राप्त करता है।
के रूप में एक सरल कम पास फिल्टर के लिए:
Vin -------R-------------- Vout
|
C
|
------------------------------
इससे उपज होगी:
वीओ यू टीवीमैं एन= एस सीआर + एस सी
MATLAB के लिए संख्या और मांद के रूप में उपरोक्त समीकरण लिखें:
num = [C 0];
den = [C R];
फिर किसी भी मैटलैब फ़ंक्शन का उपयोग करने पर अनुसरण करें जिसे आप ट्रांसफर फ़ंक्शन (बोड), पोल जीरो आरेख और इतने पर विश्लेषण करना चाहते हैं।
नीचे दिए गए फिल्टर का एक उदाहरण है जो मैं हाल ही में खेल रहा था और मूल्यों को ट्यून करने की कोशिश कर रहा था:
R1 = 20e3;
C1 = 235e-9;
R2 = 2e3;
C2 = 22e-9;
num = [2*R2*C1 0];
den = [C1*R1*C2*R2*2 (2*C1*R1 + C2*2*R2) 2];
g = tf(num,den);
P = bodeoptions; % Set phase visiblity to off and frequency units to Hz in options
P.FreqUnits = 'Hz'; % Create plot with the options specified by P
bode(g,P);
%[num,den] = eqtflength(num,den); % Make lengths equal
%[z,p,k] = tf2zp(num,den) % Obtain zero-pole-gain form