本腳本測試軟件：①MATLAB R2019b；②SleepSign2-OBI

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

目錄

• 運行腳本

• 自定義文件信號采樣率和單元格時間

• 繪制信號典型圖

• 繪制時頻譜圖

• 參考資料

運行腳本

雙擊UntitledProgram.exe，或使用MATLAB打開UntitledProgram.m。

運行腳本打開程序界面，點擊右下角@Instruction可跳轉(zhuǎn)使用說明。

自定義文件信號采樣率和單元格時間

File Sampling Rate：

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

Epoch Time：

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

繪制信號典型圖

點擊Typical Graph，選擇使用SleepSign的text output功能導(dǎo)出的幅值文件。

導(dǎo)出一列幅值即可，導(dǎo)出多列也只能繪制第一列幅值。

其中頻域圖像較密集，放大可見柱體。

點擊“文件”→“另存為”，保存為svg矢量圖。

可以先在MATLAB里畫好草圖，再根據(jù)需求將svg圖片導(dǎo)入Adobe Illustrator進一步修補美化。

繪制時頻譜圖

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

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

依次導(dǎo)入所需文件。

在Input處輸入時頻圖需要平滑的倍數(shù)，不宜太大以免卡死。

可根據(jù)實際情況選擇繪制原圖或擬合后的平滑圖，參考效果如下：

原圖（n=1，左）；平滑圖（n=10，右）

幅值圖和時頻圖繪制互不影響，如果只需要時頻圖，隨便打開一份幅值信號即可。時頻圖顏色深度可在坐標(biāo)區(qū)屬性中自行修改。

點擊“文件”→“另存為”，保存為svg矢量圖。

?

可以先在MATLAB里畫好草圖，再根據(jù)需求將svg圖片導(dǎo)入Adobe Illustrator進一步修補美化。

參考資料

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

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

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

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

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

[6]?知乎《[振動與測試 2] 什么是PSD（功率譜密度）》：https://zhuanlan.zhihu.com/p/49328001

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

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

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