本腳本測(cè)試軟件：①M(fèi)ATLAB R2019b；②SleepSign2-OBI

下載地址：https://wwxr.lanzoum.com/b04dsjohg?密碼:blct

歷史版本：https://wwxr.lanzoum.com/b04dsp7sb?密碼:ayzl

目錄

• 運(yùn)行腳本

• 自定義文件信號(hào)采樣率和單元格時(shí)間

• 繪制信號(hào)典型圖

• 繪制時(shí)頻譜圖

• 參考資料

運(yùn)行腳本

使用MATLAB打開(kāi)UntitledProgram.m，或雙擊UntitledProgram.exe直接運(yùn)行。

運(yùn)行腳本打開(kāi)程序界面，點(diǎn)擊右下角@Instruction可跳轉(zhuǎn)使用說(shuō)明。

自定義文件信號(hào)采樣率和單元格時(shí)間

File Sampling Rate：

? ? ? ?文件采樣率，即信號(hào)幅值文件中每秒顯示的信號(hào)個(gè)數(shù)。在腦波記錄軟件VitalRecorder中可以更改，通常為128Hz。

Epoch Time：

? ? ? ?單元格時(shí)間，在使用腦波分析軟件SleepSign導(dǎo)出數(shù)據(jù)前可以更改，根據(jù)實(shí)際填寫(xiě)即可，默認(rèn)為4s。

繪制信號(hào)典型圖

點(diǎn)擊Typical Graph，選擇使用SleepSign的text output功能導(dǎo)出的幅值文件。導(dǎo)出一列幅值即可，導(dǎo)出多列也只能繪制第一列幅值。

其中頻域圖像較密集，放大可見(jiàn)柱體。

點(diǎn)擊“文件”→“另存為”，保存為svg矢量圖，可導(dǎo)入Adobe Illustrator進(jìn)一步修補(bǔ)美化。

繪制時(shí)頻譜圖

點(diǎn)擊Spectrogram，選擇使用SleepSign的text output和FFT功能導(dǎo)出幅值與能譜文件。

其中幅值文件需導(dǎo)出兩列信號(hào)，第一列為EEG，第二列為EMG。

依次導(dǎo)入所需文件，幅值圖和時(shí)頻圖繪制互不影響，如果只需要時(shí)頻圖，隨便打開(kāi)一份幅值信號(hào)即可。

時(shí)頻圖可選擇繪制原圖或擬合后的平滑圖，參考效果如下：

N越大耗時(shí)越長(zhǎng)，繪制時(shí)會(huì)顯示等待條。但由于刷新等待條進(jìn)度的動(dòng)畫(huà)需耗費(fèi)5s左右，在N數(shù)低時(shí)較浪費(fèi)時(shí)間，故這版代碼暫不予顯示。

點(diǎn)擊“文件”→“另存為”，保存為svg矢量圖，可導(dǎo)入Adobe Illustrator進(jìn)一步修補(bǔ)美化。

參考資料

[1]?MathWorks官網(wǎng)：https://ww2.mathworks.cn/help/matlab/referencelist.html

[2]?知乎《傅里葉分析之掐死教程（完整版）》：https://zhuanlan.zhihu.com/p/19763358

[3]?CSDN《【STM32F407的DSP教程】第27章 FFT的示波器應(yīng)用》：https://blog.csdn.net/Simon223/article/details/106402717

[4]?bilibili《matlab實(shí)現(xiàn)FFT算法》：https://www.bilibili.com/video/BV1i14y1K7Bq

[5]?知乎《[振動(dòng)與測(cè)試 2] 什么是PSD（功率譜密度）》：https://zhuanlan.zhihu.com/p/49328001

[6]?知乎《信號(hào)頻域分析方法的理解（頻譜、能量譜、功率譜）》：https://zhuanlan.zhihu.com/p/441559677

[7]?CSDN《Matlab打包可獨(dú)立運(yùn)行的exe可執(zhí)行文件》：https://blog.csdn.net/thlzzz/article/details/109697423

[8]?CSDN《Matlab中imagesc用法》：https://blog.csdn.net/qq_21449473/article/details/123083741

[9]?CSDN《matlab添加進(jìn)度條（waitbar）》：https://blog.csdn.net/weixin_43465015/article/details/89294079

由于從零學(xué)起參考內(nèi)容太過(guò)瑣碎，大多數(shù)函數(shù)在官網(wǎng)均可查，故此處僅列出部分參考資料。

我盡可能將每段代碼的含義在腳本文件中用中文備注出來(lái)，詳見(jiàn)每段代碼“%”后的內(nèi)容。