معلومة

كيفية تقليل الضوضاء والتخلص من اتجاه 1 / f لملف PSD في تجربة مزامنة EEG

كيفية تقليل الضوضاء والتخلص من اتجاه 1 / f لملف PSD في تجربة مزامنة EEG



We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

أنا أقدر بعض المساعدة منكم يا رفاق.

أقوم حاليًا بتحليل تجربة حيث قمت بتوصيل المشاركين إلى هذه الترددات:

الشرط 1: 4.16 هرتز -8.33 هرتز
الشرط 2: 2.77 هرتز -8.33 هرتز
الشرط 3: 2.4 هرتز

لكل حالة مررت بـ 8 تجارب مدتها 26 ثانية عند Fs = 500Hz. أقوم بتحليل البيانات باستخدام FieldTrip وبعض صناديق الأدوات في Matlab. لقد حسبت السلطة و PSD استخدام pwelch تعمل بعد معالجة البيانات في Fieldtrip.

أريد بياناتي بوحدات ديسيبل والتخلص من الاتجاه $ 1 / f ^ gamma $ حيث $ 0 < gamma <2 $ (Demanuele et al.، 2007)

لدي مشكلتان:

1) لا أعرف كيفية إزالة الاتجاه 1 / f بشكل صحيح لأنني لا أفهم كيف أن المقالة التي أذكرها تفعل ذلك (انظر أدناه) في مجال التردد (لقد حاولت تركيب $ 1 / f ^ gamma $ مع أفضل $ gamma $ ، باستخدام lsqcurvefit وظيفة. باستخدام ميدفيلت 1 مع الطلب 15 يعطيني أفضل النتائج. على الرغم من أنني أعتقد أنه يجب أن تكون هناك طريقة أفضل ، إلا أنه يمكن الدفاع عن ذلك في مقال).

"تم حساب المخطط الطيفي لقناة EEG واحدة. ثم تم العثور على الوسيط عبر جميع النوافذ الزمنية لكل نقطة تردد. وبالتالي ، تم الحصول على رسم بياني لمتوسط ​​قيمة PSD لكل تردد. تم تكرار نفس الإجراء لجميع قنوات EEG. بعد ذلك ، تم حساب الوسيط العام لمتوسط ​​منحنيات PSD لجميع القنوات. تم إجراء نفس الشيء لكل من حالة "المهمة" و "الراحة". ثم تم حساب المتوسط ​​عبر جميع الظروف وتم استخدام هذا كمنحنى التطبيع الخاص بنا . "


تعديل اعتقدت أنه قد يكون من المفيد تقديم الكود الذي صنعته. آسف على أوجه القصور التي قد تجدها لأنني لست مطورًا:

٪ اختر ملف matfiles / موضوع تحميل الموضوع 5٪ حيث (القناة ، العينات ، التجربة ، الشرط) عبارة عن عينات مصفوفة 4D = numel (subject5 (1،:، 1،1)) ؛ ٪ عدد العينات٪ عدد العينات لكل نافذة ، كلما قل عدد النوافذ كلما كان دقة التردد أفضل وأسوأ دقة زمنية. لقد جربت نوافذ مختلفة. nsc = أرضية (عينات / 4.5) ؛ nov = أرضية (nsc * 0.75) ؛ ٪ تداخل بنسبة 75٪ fr = 0.5: 0.05: 15 ؛ ٪ foi fs = 500 ؛ ٪ التكرار لأخذ العينات٪ تهيئة متغير متوسط ​​ u200b  u200b = أصفار (291،28،21،3) ؛ شروط الحلقة٪ للشرط = 1: الطول (الموضوع 5 (1،1،1، :))٪ محاولات الحلقة للتجربة = 1: الطول (الموضوع 5 (1،1 ،: ، الشرط))٪ قنوات الحلقة لـ chan = 1: length (subject5 (:، 1، trial، cond))٪ حساب المخطط الطيفي لكل قناة ، ويعطي التجربة والشرط٪ t (الوقت) و ps (تقدير PSD لكل تكرار ووقت بوحدات السعة ^ 2 / هرتز) [~، ~، t، ps] = مخطط طيفي (موضوع 5 (تشان:، تجربة، شرطية)، طرق (nsc)، nov، fr، fs) ؛ ٪ حساب الوسيط عبر النقاط الزمنية لكل قناة ، التجربة والشرط medps (:، chan، trial، cond) = الوسيط (ps، 2)؛ end end٪ حساب الوسيط عبر القنوات لكل تجربة والشرط medchan = الوسيط (medps ، 2) ؛ ٪ هذا ينتج مصفوفة freqx1channelxtrialsxcond medchan = permute (medchan، [1 3 4 2]) ؛ ٪ الحصول على 4D إلى 3D ٪ حساب الوسيط على التجارب لكل حالة medtrial = الوسيط (medchan ، 2) ؛ وسطي = تبديل (وسطي ، [1 3 2]) ؛ ٪ احصل على مصفوفة ثنائية الأبعاد freqxcond٪ المتوسط ​​فوق الظروف هو منحنى التطبيع الخاص بنا منحنى التطبيع = متوسط ​​(متوسط ​​، 2) ؛

2) بمجرد أن أقوم بتحويل البيانات إلى وحدات ديسيبل (أي 10 دولارات * لوغاريتم 10 دولارات) ، فإنها تصبح صاخبة حقًا ، مما يجعل القمم في foi (الترددات ذات الأهمية) كبيرة مثل الخلفية أو تقريبًا. لقد بحثت في SNR ولكن في الحقيقة لا أعرف كيفية تطبيقها هنا. الشيء الوحيد الذي فعلته وما زلت لا يزيل كل الضوضاء هو ميدفيلت 1 المذكورة في 1).

آمل أن يساعدني أي منكم لديه خبرة في هذا في حل هذه المشكلة!


من الواضح أن الورقة التي ذكرتها تقوم ببساطة بتقييم الوسيط في طيف مخطط كهربية الدماغ. لنفترض أنه بالنسبة لجميع قنواتك ، تحصل على مصفوفة ثلاثية الأبعاد PSD ، مع البعد 1 رقم القناة ، والبعد 2 النوافذ الزمنية ، والبعد 3 الترددات. في الشفرة الكاذبة ، تحسب المقالة M = الوسيط (الوسيط (PSD ، 2) ، 1) ، ثم تزيل الاتجاه باستخدام CorrectedPSD = PSD - M. وهو في الواقع بسيط إلى حد ما.


شاهد الفيديو: écrire en arabe sur photoshop. الكتابة باللغة العربية في الفوتوشوب (أغسطس 2022).