यह sinc मेरे LP-turn-HP फ़िल्टर आउटपुट में क्या कर रहा है?


10

मैं एलपी-इसे छानकर और मूल फिल्टर से आउटपुट घटाकर एक सिग्नल को एचपी-फ़िल्टर करने की कोशिश कर रहा हूं। एलपी-फ़िल्टर किए गए सिग्नल को प्राप्त करने के लिए, मैं 16,000-नमूना सिग्नल के लिए तीसरे क्रम के बटरवर्थ फ़िल्टर का उपयोग करके आगे-पीछे फ़िल्टर किया गया। कटऑफ आवृत्ति 0.5 हर्ट्ज @ 250 हर्ट्ज थी। कोई खिड़की का उपयोग नहीं किया गया था।

पहले 85% या इससे अधिक नमूनों में, वास्तव में एचपी-फ़िल्टर किए गए सिग्नल के साथ इस पद्धति का आउटपुट पूरी तरह से ओवरलैप होता है। हालांकि, आउटपुट के अंत में, एक सुपरसुनीक सिनकी जैसी घटना होती है।

  • इसका क्या कारण होता है?
  • इस विधि को संभव बनाने के लिए मैं क्या बदल सकता हूं?

यहाँ दो संकेत ओवरले हैं। नीला एक HP- फ़िल्टर्ड है, और लाल एक घटाव द्वारा प्राप्त किया गया है।

एचपी और मूल-एलपी सिग्नल

यहां प्रति नमूना पूर्ण अंतर है। वे कभी बराबर नहीं होते; न्यूनतम अंतर 8e-12 है। 87% मान 0.01 से नीचे हैं।

दो तरीकों के बीच अंतर

(अगर एलपी फिल्टर के से एचपी फिल्टर प्राप्त करने का एक बेहतर तरीका है , तो मुझे यह जानने में खुशी होगी)।एच(z)


क्या आप MATLAB कोड को शामिल कर सकते हैं जिसका उपयोग आप प्लॉट बनाने के लिए कर रहे हैं? डेटा महत्वपूर्ण नहीं है; यदि आवश्यक हो तो आप यादृच्छिक नमूने के साथ इसे अनुकरण कर सकते हैं।
जेसन आर

जवाबों:


0

मिथिंक्स @ hotpaw2 सही रास्ते पर है। मेरा अनुमान है कि filtfilt रिवर्स फिल्टर के लिए प्रारंभिक शर्तों को काफी सही नहीं संभाल रहा है। यह एक बहुत कठिन समस्या की तरह लगता है - मुझे नहीं लगता कि मैं गलत तरीके से इसे प्राप्त करने के लिए मैथवर्क्स को दोष दे सकता हूं ...

यह फिक्स मेरे लिए काफी आसान लगता है: जब तक कि साइनक दूर नहीं हो जाता, तब तक अपने सिग्नल को ज़ीरो-पैड करें और फ़िल्टर के बाद (अब नॉनज़रो!) पैडिंग को हटा दें।


0

आप एक के लिए एक पीछे की ओर-आगे फिल्टर कर रहे हैं, और दूसरे के लिए एक आगे-केवल फिल्टर। इस प्रकार दो फिल्टर एक या दोनों किनारों पर अलग-अलग व्यवहार करेंगे।

यदि आपका सिग्नल डेटा विंडो की लंबाई में दृढ़ता से स्थिर और बिल्कुल आवधिक है, तो आप डेटा को परिपत्र रूप से चारों ओर लपेटने का प्रयास कर सकते हैं और अपने आवेग प्रतिक्रिया लंबाई द्वारा अपने फ़िल्टर को शुरू कर सकते हैं।


क्षमा करें यदि मैं पर्याप्त स्पष्ट नहीं था। जिस सिग्नल की मैं तुलना कर रहा हूं (नीला) वह भी आगे-पीछे फ़िल्टर ( filtfiltऔर butter(3,0.5/250,high)) है। संकेत ईसीजी से है कि मैं बिल्कुल आवधिक होने की उम्मीद नहीं कर सकता।
एंड्रियास

क्या मैं कुछ गलत समझ रहा हूं? मैं upvote को नहीं समझ सकता क्योंकि पहला पैराग्राफ स्पष्ट रूप से गलत है।
एंड्रियास
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.