मैं इस पेपर में वर्णित मॉडल को पुन: पेश करने का प्रयास कर रहा हूं https://hal.archives-ouvertes.fr/file/index/docid/683477/filename/clarinette-logique-8.pdf

यहां एक विधि है जो त्रिज्या a और लंबाई L के सिलेंडर के लिए स्थानांतरण मैट्रिक्स लौटाती है।

# density of air
rho = 1.2250

# viscosity of air
eta = 18.5

# transfer matrix for a cylinder
def Tcyl(a, L):
    rv = a * math.sqrt(rho*omega/eta)
    Z0 = (rho * constants.c) / (constants.pi * a**2)
    Zc = Z0 * complex(1 + 0.369/rv, -(0.369/rv + 1.149/rv**2))
    gamma = k * complex(1.045/rv + 1.080/rv**2, 1 + 1.045/rv)
    return
       np.matrix([[cmath.cosh(gamma*L),Zc*cmath.sinh(gamma*L)],
       [(1.0/Zc)*cmath.sinh(gamma*L),cmath.cosh(gamma*L)]], dtype=complex)

इनपुट प्रतिबाधा की गणना विभिन्न आवृत्तियों के लिए निम्नानुसार की जाती है।

for f in range(1,4000):
    
    # angular frequency
    omega = 2*constants.pi*f

    # wave number
    k = omega/constants.c

    # transfer matrix for a cylinder of radius 7.5mm and length 100mm
    T = Tcyl(7.5, 100.0)

    # radiation impedance
    Zr = complex(0,0)

    # [P,U] vector (pression and velocity)
    T = np.dot(T,np.matrix([[Zr],[complex(1,0)]], dtype=complex))

    # input impedance (Z = P/U)
    Z = T.item(0)/T.item(1)

खेलने की आवृत्ति समीकरण Im[Z]=0 को संतुष्ट करती है। Z के काल्पनिक भाग को प्लॉट करते समय मुझे निम्न आकृति मिलती है: गलत प्रतिबाधा

यह स्पष्ट रूप से गलत है क्योंकि अपेक्षित आउटपुट कुछ इस तरह होना चाहिए: सही प्रतिबाधा

मैं क्या गलत कर रहा हूं? धन्यवाद।

0
pamplemousse 21 मार्च 2021, 05:53
(1) प्रश्न में कोड दिखाएं कि आप इसे कैसे प्लॉट करते हैं। (2) यदि वेरिएबल्स "ओमेगा" और "के" का उपयोग "टीसीएल" में किया जाता है, तो वे फ़ंक्शन के पैरामीटर भी होने चाहिए।
 – 
Michael Butscher
21 मार्च 2021, 06:05
आप सही हैं ओमेगा और के फ़ंक्शन के पैरामीटर होने चाहिए
 – 
pamplemousse
21 मार्च 2021, 21:27

1 उत्तर

सबसे बढ़िया उत्तर

आपके पास

    Z0 = (rho * constants.c) / (constants.pi * a**2)

शहनाई की प्रतिबाधा ध्वनि की गति पर निर्भर करती है, न कि प्रकाश की गति पर। constants.c को 343 से बदलें और आपके परिणाम करीब होंगे। मुझे अभी भी यकीन नहीं है कि यह बिल्कुल सही है, लेकिन करीब है।

एक शहनाई वादक के रूप में, मैं लोगों को यह सोचने की कोशिश करता हूं कि मेरी उंगलियां प्रकाश की तरह तेज चलती हैं, लेकिन ऐसा नहीं है।

0
Tim Roberts 21 मार्च 2021, 06:10
हा हा यह बहुत शर्मनाक है। यह अब पूरी तरह से काम करता है! शुक्रिया !!!
 – 
pamplemousse
21 मार्च 2021, 21:28
कुछ अन्य सुधार: हवा की चिपचिपाहट 0.0000181 होनी चाहिए, त्रिज्या और लंबाई मीटर में व्यक्त की जानी चाहिए।
 – 
pamplemousse
21 मार्च 2021, 23:04