ইইজি স্লো প্রবাহের জন্য 0.01 Hz এর হাই পাস ফিল্টার

Question

আমি খুব ধীর চালিকা থেকে মুক্তি পেতে ইলেক্ট্রোয়েন্সফ্লোগ্রাফিক (ইইজি) ডেটাসেটের জন্য একটি উচ্চ পাস ফিল্টার তৈরি করার চেষ্টা করছি। তবে এই ডেটা পরিবেশন করা গবেষণার জন্য ০.০ হার্জেডের প্রায় ফ্রিকোয়েন্সিগুলি খুব গুরুত্বপূর্ণ important

আমি মতলব ব্যবহার করি এবং আমি বিভিন্ন জিনিস চেষ্টা করেছিলাম। দেখে মনে হচ্ছে একটি খুব কম ডিবি অ্যাটেন্যুয়েশন, 3 ডিবি এর এফআইআর ফিল্টারগুলির সাথে ঠিক আছে। উদাহরণস্বরূপ আমি এটি ব্যবহার করেছি:

h=fdesign.highpass('Fst,Fp,Ast,Ap',0.005,0.01,3,1,250); % i needed to cut 
d=design(h,'butter');
fvtool(d)

তবে এই ফিল্টারটি ডেটার শুরুতে একটি ডিসি ড্রিফ্ট তৈরি করে এবং আমি এই প্রথম ডেটাপয়েন্টগুলি আলগা করতে চাই না।

আমি আপনার পৃষ্ঠাগুলির মাধ্যমে পড়েছি যে একটি চালাক ধারণাটি হ'ল একটি কম পাস ফিল্টার ডিজাইন করা যা পরে আমি আসল তথ্য থেকে বিয়োগ করতে পারি। আমি এটি একটি ইক্যুরিপল ফিল্টার দিয়ে করেছি এবং এটি কাজ করেছে তবে একটি বেসলাইন লিফট রেখে গেছে। বিশ্লেষণ প্রোটোকলের এটি পরবর্তী ধাপে আমি এখনই গড় অপসারণ চালাতে চাই না।

কোন পরামর্শ?

Comments

তাহলে সমস্যাটি কেবল ফিল্টারটির প্রাথমিক প্রতিক্রিয়া? এই পূর্বনির্ধারিত ডেটা যা আপনি কেবল শূন্যগুলির সাথে প্রিপেন্ড করতে পারেন বা এটি রিয়েলটাইম?

— এন্ডোলিথ

লোপাস এবং বিয়োগের পদ্ধতিটি আপনি যে সমস্যাটি দেখছেন তাতেও ভুগবেন। সমস্যাটি হ'ল প্রতিটি ফিল্টারটির বিলম্ব হয়। আপনি চান তার মতো ধারালো কাট অফ সহ ফিল্টারগুলিতে খুব বেশি বিলম্ব হতে পারে। আপনি কেন পরের গড় অপসারণ করবেন তা আমি নিশ্চিত নই; আপনি যে ফিল্টারটি দেখিয়েছেন তার ডিসি লাভটি শূন্য, সুতরাং অপসারণের জন্য কোনও প্রশংসনীয় গড় থাকবে না (ফিল্টারটির বিলম্বের কারণে আউটপুট শুরুর সময় ক্ষণস্থায়ী সময় ব্যতীত)।

— জেসন আর

Answer 1

আপনি একটি টাইম মেশিনের জন্য জিজ্ঞাসা করছেন: তে ডিসি বাম্প এড়াতে আপনাকে ফিল্টারটির অবস্থাটি জানতে হবে যেমন রেকর্ডিং শুরুর আগে এটি চলছে। আপনি কোন ধরণের ফিল্টার ব্যবহার করেন না কেন তা ঘটবে না।$t=0$

এখানে একটি কৌতুক যা এর পারে আপনার জন্য উপকারে আসা। ধরা যাক আপনার একটি ট্যাপ এফআইআর ফিল্টার রয়েছে (বা একটি আইআইআর ফিল্টার যার অনুপ্রেরণাগুলি নমুনার পরে 0 পর্যায়ে পর্যাপ্ত হয়ে যায়)। আপনার সিগন্যালের প্রথম নমুনাগুলি নিন , সেগুলি বিপরীত করুন এবং এটি আপনার সিগন্যালের শুরুতে প্রিপেন্ড করুন। কার্যত, আমরা ফিল্টারিংয়ের উদ্দেশ্যে একটি নতুন সংকেত সংজ্ঞায়িত করছি । আপনি যখন এটি ফিল্টার করেন, ফিল্টার রাষ্ট্রটি ইতিমধ্যে প্রাথমিক সিগন্যালের সাথে "প্রাইমড" হয়ে উঠবে: আউটপুটটির প্রথম নমুনাগুলি মুছুন , এবং ডিসি বাম্প আশাবাদী চলে যেতে হবে।$N$$N$$N$$g(t)=g(-t)=f(t)$$N$

বলা বাহুল্য, সেই স্কিমটি নিয়ে প্রচুর গোছা রয়েছে - প্রারম্ভিকদের জন্য, রিয়েলটাইম প্রসেসিংয়ের জন্য এটি করা আরও শক্ত এবং যদি সিগন্যাল ডেরিভেটিভ চেয়ে বেশি হয় তবে আপনার প্রথম কয়েকটি নমুনা নিয়ে বিশ্বাসের সমস্যা থাকতে পারে। তবে এটিকে গলদা থেকে মুক্তি দেওয়া উচিত।$t=0$

Answer 2

ফিল্টফিল্ট ফাংশনটি একবার দেখুন। এটি জিরো-ফেজ প্রতিক্রিয়া এবং একটি নিখুঁত পদক্ষেপ প্রতিক্রিয়া দেয়। বিশেষত, ফিল্টফিল্ট ফিল্টারিংয়ের পদক্ষেপের প্রতিক্রিয়া সম্পত্তি আপনার সমস্যার সমাধান করতে পারে।

Answer 3

আমি মিঃ rtollert এর কৌশলটি ব্যাখ্যা করেছিলাম এবং ভেবেছিলাম এটি আমার পক্ষে সবচেয়ে ভাল।

আপনি যদি অবিচ্ছিন্ন নমুনা / হার্ডওয়্যার ক্ষতিপূরণের পথে চলে যান তবে দক্ষতা আপনাকে ডিসি ব্লকারদের প্রতি আগ্রহী করে তুলতে পারে র্যান্ডি ইয়েটস এবং রিচার্ড লায়ন্স দ্বারা ডিসিএস টিপস এন্ড ট্রিকস মার্চ ২০০৮

Answer 4

আমি এমন কিছু চেষ্টা করেছি যা নির্দিষ্ট এম্প্লিফায়ারের জন্য ভালভাবে কাজ করেছিল। ম্যাটল্যাব কোডটি এখানে: https://sites.google.com/site/marialstavrinou/home/dc-offset-removal-filter-in-matlab ।

Answer 5

আপনি নিখুঁত ডিসি মান সম্পর্কে যত্নশীল না?

আপনার ডেটাতে কেবল স্ট্যাটিক অফসেটটি যুক্ত করবেন না কেন যাতে প্রথম ডেটা পয়েন্টটি শূন্যে থাকে?

অবশ্যই, আপনাকে প্রত্যেকটি ডেটাপয়েন্টে অফসেট যুক্ত করতে হবে তবে আপনি ফিল্টারটির সাথে যে বড় পদক্ষেপ-প্রতিক্রিয়া করছেন তা সম্পূর্ণরূপে এড়ানো হবে।

কার্যকরভাবে ফিল্টার শূন্য থেকে শুরু হয়। অতএব, যখন আপনার ডেটা শুরু হয়, ফিল্টার শূন্য থেকে আপনার যে কোনও ডি সি লেভেল পর্যন্ত একটি বড় সিঁড়ি ধাপ দেখে।

সিঁড়ি ধাপটি সরাতে কেবল অফসেট যুক্ত করুন।

Answer 6

গড় অপসারণ করবেন না কেন? পরবর্তী পদক্ষেপ যা এটি "আবার" করে কেবল কার্যকর হয় না (এটি শূন্যকে বাদ দিয়ে দেবে)।

Comments

এটি আরও একটি প্রশ্নের, এবং সত্যই প্রশ্নের উত্তর নয়। এটি একটি পর্যায়ে মন্তব্য হিসাবে উপযুক্ত হতে পারে, যদিও প্রশ্নটি এত পুরানো তবে এই মুহুর্তে এটি কেবল অর্থহীন বলে মনে হচ্ছে।

— স্যাম মালুনি