基于单片机的脑波治疗仪的软件设计研究

其中Ω为调制信号的频率，即脑波频率，频率范围为1Hz～30Hz；w0为载波的频率，范围为90Hz～800Hz，初始频率设定为330Hz。输出调幅波波形如图3所示。

(2)调幅波的载波频率调节
载波频率即脑波音乐的音调，将音调调整到人耳敏感的范围，会使治疗效果更加明显。人们可以听到的声振动频率范围为20Hz～20kHz，设计要求频率调节范围为90Hz～800Hz。载波频率计算如式(5)：

fw0=1/nT （5）
载波的频率由载波的采样点点数n和取点时间间隔T决定。通过改变采样点点数来调节需要大量的采样点，同时由于取点数量的改变必然带来较大的相位截断误差。本设计通过改变T来调节频率，这样保证了采样点不漏取，同时采样点点数不用取太多。调频过程中，首先计算所需频率对应的时间T，然后将其赋给取点时间间隔参数，这样载波就被跳到所需频率。
(3)调幅波的调制信号频率调节
在设计中，调幅波的调制信号即脑波音乐的包络，它确定了脑波音乐的特征(尤其是治疗的内容)。调制波的频率调节是在一定的脑波频率范围内(1Hz～30Hz)变化，其频率随时间变化的时间、频率参数已经固化在ROM中，不能任意调节。调制信号频率调节与载波频率调节的原理相同，可以参考公式5。本设计中调制信号随时间的频率变化参数已经换算成时间参数，程序运行时采用查表方式调用该时间参数来改变调制信号的频率。
(4)脑波音乐强度的调节
将脑波音乐信号vAM作为DA2的参考电压Vref，单片机的强度参数PW做为DA2的输入信号。DAl和DA2采用DAC0832，因为在设计中DAC都采用单极性反向电压输出方式，所以DA2的输出如式(6)：

这样脑波音乐强度在0~-vAM内变化。
2．2．2 双声道音乐的生成
如图4，模拟开关的输入口A2、B1接低电平，Al、B2同时接入脑波音乐信号，单片机通过对控制口l、2来控制两路开关的通断。当使左右两路的开关都始终打开就会产生左右两路同步输出；当单片机以与调制信号相同的频率控制模拟开关交替开、关就会产生左右两路信号交替输出。这两路输出信号经过放大滤波电路输出至耳机，这种交替／同步的立体声作用到人脑将引导大脑起到双脑同步的有效治疗效果。图5就是采集到的左右交替输出的信号。

2．3 节目时间、节目指示灯显示部分
节目时间显示部分是由单片机调用内部定时器进行分钟计时，输出时间信号到LED数码管进行显示。节目时间在开机时显示0，在运行过程中，当节目时间到了预定时间将调用报警子程序发出报警声，并将节目结束标志置位。节目指示灯显示是在节目选择子程序中由单片机控制节目指示发光二极管灯的亮灭，指示当前输出哪个节目。

3 结果分析
通过实验观察及示波器测量，设备能输出有效频率的脑波音乐，满足预计的参数设置，波形无明显失真。图5是采集到的的波形。

但是实际输出波形还是有一些误差，其主要来自：①由于内部波形存储器中存储的正弦幅度值是用二进制表示的，对于越过存储器字长的正弦幅度值必须进行量化处理，这样就引入了量化误差。幅度量化主要有两种方式，即舍入量化和截尾量化，这里采用截尾量化方式。②DAC非理想转换特性所引起的杂散水平，DAC的非理想特性包括：DAC的有限分辨位数；差分、积分的非线性；D／A转换过程中的瞬间毛刺；时钟泄露；数字噪声馈通；转换速率受限等。③调制信号频率调节过程中正弦波的各采样点会有漏取点的情况发生，这样就不可避免地产生相位截断误差。
针对以上问题，提出一些改进的方法：①应用更高位数的单片机和高性能的数模转换器。②增大ROM容量及采样点数，提高幅度分辨率。