سحاسبه و رسم خصوصیات زمان و فرکانس

[hd_uni_pro]

سلام دوستان راستش من یک سوال دارم راجع به نرم افزار متلب ، من می خوام با استفاده ا ز سیگنال ECG در متلب کارهای پردازشی انجام بدم مثلا می خوام خصوصیاتی مانند :energyوaverage magnitude,zero crossing ,autocorrelation , signal segmentation و signal filtering و denoisingو نسبت سیگنال به نویز یا همون0 SNR رو حساب کنم با استفاده از سیگنال ECG .حال توی متلب چه دستورات هستش که این کارو بکنه البته من به هلپش رفتم ولی خیلی می پیچونه من فقط دستوراتش و نحوه کارشو رو می خوام البته من تونستم دادهای سایت فیزیونت رو وارد متلب کنم ولی برای محاسبه مواردی که گفتم سردرگمم لطفا کمک کنید باتشکر

[matrix12]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

سلام...

کد متلب برای انرژی:
e=c.*c;
E=sum(e);

کد متلب اتوکریلیشن:
autocorr(c);

و اینکه zero crossing رو اگه بخواین نمونه کد رو میتونم واستون بذارم؟

[hd_uni_pro]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

سلام
ممنون
نیکی و پرسش

[hd_uni_pro]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

ممنون .اما این طوری انرژی یک عدد میشه !یعنی سیگنال انرژی یک ECG همین طوری محاسبه می شود !؟

[matrix12]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

سلام

قاعدتا شما اول سیگنالتون رو فریم بندی میکنید و بعد انرژی هر کدوم از فریم ها رو حساب میکنید که میشه یه وکتور انرژی.چون هر محاسبه ای که بخواین انجام بدین مثل فوریه و ....باید روی فریم ها انجام بشه و محاسبه روی کل سیگنال بدرد نمیخوره.

[matrix12]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

این یه نمونه کد واسه نرخ قطع صفر هست.

p=signal % no DC signal
for i=1:47
if p(i).*p(i+1)<0
zero=1
stem(i,p(i))
end
end

که signal قبلا تعریف شده تو برنامه.

[hd_uni_pro]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

منظور از فریم بندی همون سگمنتیشن است دستورش چیه ؟ر ضمن برای محاسبه خواص بالا لازم سیگنال برای هر کدام مقدار دی سی اون حذف بشه !

[matrix12]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

راستش سگمنتیشن رو نمیدونم!شاید....اما فریم بندی همون پنجره بندی سیگنال هست که سه فریم متداول مورد استفاده مستطیلی،همینگ و هنینگ هستن.

فریم بندی رو میشه هم بصورت کد نویسی ارائه داد و هم با دستور window در متلب میشه پنجره تنظیماتش رو اورد که برای حالت دو باید بعد از فراخوانی سیگنال دستور window رو بزنیم .یه پنجره باز میشه که همه ی ویژگی هایی که میخوایم رو در اون تعیین میکنیم که نوع فریم و طول فریم هست و اون رو توی command window ذخیره میکنیم و سپس اون رو به سیگنال اعمال میکنیم.

کد فریم بندی رو هم سعی میکنم بین کدهام پیدا کنم بذارم براتون.

موفق باشید.

[hd_uni_pro]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

این یه نمونه کد واسه نرخ قطع صفر هست.

p=signal % no DC signal
for i=1:47
if p(i).*p(i+1)<0
zero=1
stem(i,p(i))
end
end

که signal قبلا تعریف شده تو برنامه

چرا i =1:47 است برای یک سیگنال ECG که ارایه ای 2*3600 داره فکر کنم با i =3600 و یا 1800 باشه .درسته در ضمن خط p=signal این سیگنال !منظورتون چه جور سیگنالی است مثلا این جورب:

a=val(1,

p=a
val یم ارایه ecg 1*3600 است

[hd_uni_pro]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

سلام دوستان من تمام مواردی که گفتم محاسبه کردم بات وجه به اینکه بیشتر دوستان فقط موضوع تاپیک رو نگاه می کنن و زورشون میاد کمک کنم تصمیم گرفتم اول مقاله رو کلا تو سایت اپلود کنم ولی نظرم عوض شد حالا اگر کسی می خواد یک پیغام خصوصی بفرسته تا براش اپلود کنم ((البته از زحمات کاربر محترم (فرشته) هم ممنون هستم

[matrix12]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

سلام

الان دیدم سوالاتتون رو تو پست یکی مونده به اخر...اومدم جواب بدم که دیدم کارتون به نتیجه رسیده.خداروشکر.

ولی خب شاید واقعا نمیدونسته کسی که کمکتون کنه.

اگر دوست داشتین مقاله رو بذارین همین جا که بقیه هم استفاده کنن.

موفق باشید

[HamEd276]

پاسخ : سحاسبه و رسم خصوصیات زمان و فرکانس

:sad: :sad: :sad: :sad: :sad: :sad: :sad: :sad: :sad: :sad: :sad: درود :sad: :sad: :sad: :sad: :sad: :sad:
دوستان دستورِ فریم بندیِ سیگنال و محاسبه FFT و انرزی و اندازه هر فریم در متلب دستورِ اش چیه ؟! ؟؟؟؟؟؟؟؟؟؟؟؟؟؟؟؟؟؟؟
گیرم بخدا :sad: :sad: :sad: :sad:
ای خداااااااااااااااااااااااااااااا