一种微孔板加样指示装置的制作方法

本实用新型涉及微孔板加样设备技术领域，具体为一种微孔板加样指示装置。

背景技术：

在生物化学或是卫生防疫检验等实验室，常常需要用到微孔板做检测试验，然后通过酶标仪检测各孔的吸光度，该微孔板是一种经事先包理专用于放置待测样本的透明塑料板，板上有多排大小均匀一致的小孔，通常有40孔、48孔、55孔或96孔不同规格，每个孔可盛放微量的溶液，由于微孔板上的孔排列较为紧密，且邻孔之间没有明显的标志，因此手动加样或是大批量加样的过程中常常容易加错孔，尤其是孔内试剂无颜色变化的情况下，加样出错率大为增加，降低了试验效率，鉴于此，我们提出一种微孔板加样指示装置。

技术实现要素：

本实用新型的目的在于提供一种微孔板加样指示装置，以解决上述背景技术中提出的手动加样或是大批量加样的过程中常常容易加错孔，尤其是孔内试剂无颜色变化的情况下，加样出错率大为增加，降低了试验效率的问题。

为实现上述目的，本实用新型提供如下技术方案：

一种微孔板加样指示装置，包括底座，所述底座上表面等间距设有若干指示灯，所述底座上方设有置放盒，所述置放盒上等间距开设有若干置放槽，所述底座的后侧设有置放板，所述置放板上表面开设有枪槽，所述枪槽内部嵌设有加样枪。

优选的，所述指示灯与所述置放槽的数量一一对应。

优选的，所述加样枪包括枪筒，所述枪筒侧壁靠近端部的位置设有弧形板，所述枪筒内安装有注射杆，所述注射杆的端部设有活塞，所述注射杆外套设有弹簧，所述弹簧的两端分别安装有第一固定环和第二固定环，所述注射杆侧壁顶端安装有压力传感器。

优选的，所述弧形板与所述枪筒为一体成型结构。

优选的，所述活塞与所述枪筒的内壁活动连接。

优选的，所述底座内部安装有感应装置。

优选的，所述置放盒右侧壁靠近后侧的位置安装有显示面板。

优选的，所述显示面板内部安装有控制系统，所述控制系统包括控制模块、按键模块、显示模块、传感模块和电源模块，

所述传感模块用于对压力信号和手掌信号进行采集；

所述按键模块用于对信号数据计数和工作模式进行切换；

所述控制模块用于对采集的信号数据进行分析和处理；

所述显示模块用于对分析和处理后的信号数据进行显示；

所述电源模块用于对电路进行供电，确保电路正常运行。

优选的，所述传感模块采用fdc2214芯片设计。

优选的，所述显示面板内部还设有定义成功加样模块，用于监测加样过程。

与现有技术相比，本实用新型的有益效果是：

1、本实用新型通过设计的置放板，可将加样枪置放与枪槽内，便于收纳放置，为再次使用提供便利，其次通过设计的加样枪，便于加样，并通过弹簧增加需要推动注射杆压力，避免压力过小，易造成加样过量的问题，提高实用性；

2、本实用新型通过传感模块用于对压力信号和手掌信号进行采集，通过按键模块用于对信号数据计数和工作模式进行切换，通过控制模块用于对采集的信号数据进行分析和处理，有利于降低加样出错率，提高实验效率，通过显示模块用于对分析和处理后的信号数据进行显示，通过电源模块用于对电路进行供电，确保电路正常运行。

3、本实用新型通过定义成功加样模块，若加样过程中出错，可手动退回前一步，其次优势中突出无线悬浮感应，不需要额外操作，无需改变用户操作习惯。

附图说明

图1为本实用新型的整体结构示意图；

图2为本实用新型中底座的结构示意图；

图3为本实用新型中加样枪的连接爆炸结构示意图；

图4为本实用新型中控制系统的工作框图；

图5为本实用新型中感应装置的工作原理图；

图6为本实用新型中单片机的引脚图；

图7为本实用新型中单片机与fdc2214芯片的工作原理图；

图8为本实用新型中显示模块的工作原理图；

图9为本实用新型中按键模块的工作原理图；

图10为本实用新型的算法流程图。

图中：1、底座；10、指示灯；11、感应装置；2、置放盒；20、置放槽；3、显示面板；4、置放板；40、枪槽；5、加样枪；50、枪筒；51、弧形板；52、注射杆；53、活塞；54、第一固定环；55、弹簧；56、第二固定环；57、压力传感器。

具体实施方式

下面将结合本实用新型实施例中的附图，对本实用新型实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本实用新型一部分实施例，而不是全部的实施例。基于本实用新型中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本实用新型保护的范围。

在本实用新型的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本实用新型和简化描述，而不是指示或暗示所指的设备或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本实用新型的限制。

此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本实用新型的描述中，“若干”的含义是两个或两个以上，除非另有明确具体的限定。

实施例1：

一种微孔板加样指示装置，如图1至图3所示，包括底座1，底座1上方设有置放盒2，底座1的后侧设有置放板4，置放板4上表面开设有枪槽40，枪槽40内部嵌设有加样枪5，加样枪5包括枪筒50，枪筒50侧壁靠近端部的位置设有弧形板51，枪筒50内安装有注射杆52，注射杆52的端部设有活塞53，注射杆52外套设有弹簧55，弹簧55的两端分别安装有第一固定环54和第二固定环56。

本实施例中，弧形板51与枪筒50为一体成型结构，保证弧形板51与枪筒50之间连接的稳定性，便于后续加样，弧形板51的下表面呈圆弧状，有利于提高手指与弧形板51接触时的舒适性。

进一步的，活塞53与枪筒50的内壁活动连接，当加样时，可推动注射杆52，使注射杆52推动活塞53沿枪筒50的内壁移动，推注枪筒50内部的溶液。

具体的，第一固定环54通过胶水粘接至注射杆52侧壁上，第二固定环56通过胶水粘接至枪筒50的内壁上，当推动注射杆52时，注射杆52带动第一固定环54同步移动，使弹簧55压缩，弹簧55压缩产生的弹力产生反作用，作用第一固定环54和注射杆52，增加需要推动注射杆52的压力，避免压力过小，造成加样过量的问题，提高实用性。

除此之外，当不使用加样枪5，可将加样枪5放置于枪槽40内部，便于收纳放置，为再次使用提供便利。

本实施例的微孔板加样指示装置在使用时，在加样时，首先，推动注射杆52，使注射杆52推动活塞53沿枪筒50的内壁移动，此时，注射杆52带动第一固定环54同步移动，使弹簧55压缩，进行推注枪筒50内部的溶液，实现加样，而弹簧55压缩产生的弹力产生反作用力，作用第一固定环54和注射杆52，增加需要推动注射杆52的压力，避免压力过小，造成加样过量的问题，提高实用性。

实施例2：

作为本实用新型的第二种实施例，为了降低出错率，提高实验效率，本实用新型人设置控制系统，作为一种优选实施例，如图4至图10所示，底座1内部安装有感应装置11，底座1上表面等间距设有若干指示灯10，置放盒2上等间距开设有若干置放槽20，注射杆52侧壁顶端安装有压力传感器57，置放盒2右侧壁靠近后侧的位置安装有显示面板3，显示面板3内部安装有控制系统，控制系统包括控制模块、按键模块、显示模块、传感模块和电源模块，传感模块用于对压力信号和手掌信号进行采集；按键模块用于对信号数据计数和工作模式进行切换；控制模块用于对采集的信号数据进行分析和处理；显示模块用于对分析和处理后的信号数据进行显示；电源模块用于对电路进行供电，确保电路正常运行。

本实施例中，控制模块采用采用msp430f5529lp单片机芯片设计，该芯片晶振最高可达25mhz，处理速度快，内置硬件乘法器、比较器、定时器、a/d转换等，具备usb集线器功能，可以一个usb口实现多种usb应用，支持在线调试，在线下载程序，具备5v转3.3v功能，既可以满足最小系统的供电也可以满足5v需求的扩展应用，且msp430f5529lp配合ccs编程使用较为方便。

传感模块采用fdc2214芯片设计，fdc2214传感器是具有抗电磁干扰架构、高达28位分辨率、高速、多通道的电容数字转换器，支持的传感器频率为10khz至10mhz，支持宽激励频率范围，可为系统设计带来灵活性。

进一步的，指示灯10与置放槽20的数量一一对应，指示灯10采用led指示灯，当将试管置于置放槽20内，试管加样后，led指示灯点亮，便于后续观察加样情况，为后续实验提供便利。

压力传感器57上安装有2.4ghz无线通通讯模块，2.4ghz无线通通讯模块采用nrf24l01模块设计，通过nrf24l01模块实现控制系统与手机无线通信。

具体的，感应装置11采用一块10cm*15cm和一块15cm*15cm的铜板作为两路被测电容传感端，其中一块主要测手指信息，一块主要测手掌信息，与一块fdc2214连接，并将手指信息和手掌信息整合。

电容传感器识别手势原理

当手触摸在金属层上时，由于人体电场，用户和触摸屏表面形成以一个耦合电容，对于高频电流来说，电容是直接导体，会影响电路整体电容特性。

耦合电容的计算方式可以近似地用理想平行板公式

当d＝1cm,当d<1cm时，传感器比较敏感且收噪声影响大，当d>1cm时，传感器变化较小，不易判决。

当距离一定，介电常数一定，耦合电容的大小仅与接触面积有关，即可以有效地比较出接触手指的数量。

耦合电容与fdc2214中lc振荡电路的固定电容进行并联。

利用电容并联公式

没有接触形成电容板时，电容为0，接触时，电容大小上升，lc振荡频率下降，fdc2214通过基准频率计算lc振荡频率，并转化为28位的二进制数字值，

利用电容计算公式

得到最后的电容计算值，将非线性变化的lc振荡频率转换成与面积成线性关系的电容值。

手势数据滤波算法原理

使用iir加权滑动平均算法

ft＝(w1at-1+w2at-2+w3at-3+...+wnat-n)/n

wn为权重，n为窗口宽度，这里选择6，at为t时刻电容数据值。

加权滑动平均给固定周期内的每个变量值附以不同的权重，不同时间段的数据影响力不同。对获取的信号进行滤波，有效预测得到其发展方向和趋势，消除部分随机和随机波动的影响，又能高效地利用数据。

利用多距离投票算法

多距离投票算法将训练时所有的数据保存下来。采用欧式距离，曼哈顿距离和切比雪夫距离计算测试数据和训练数据的距离。

利用欧式距离公式：

利用曼哈顿距离公式：

利用切比雪夫距离公式：

得到三种距离后，通过权重为1：1：1的投票决定最后的手势判决结果。

fdc原始数据转化为电容程序

channeldata＝fdc2214_channelread(0,time)；

fsensor＝channeldata/50；

k＝(37513180.0/fsensor)；

k＝k*k；

capacitance＝k*10000；

fdc原始数据读取滤波程序

投票算法

除此之外，控制模块通过i2c与传感模块实现数据交互，i2c是一种较低速度的串行通讯总线i2c协议v2.1规定了100k,400k和3.4m三种速率(bps)。

控制模块与传感模块通过i2c的调通步骤：

步骤1：单片机(master)写fdc芯片(slave)的i2caddress，该地址有两个，addr引脚为高时为0x2b，addr引脚为低时为0x2a,选用后者，8位地址0x2a的二进制虽然是00101010，但需在a6-a0再加一位读写位(读是1，写是0)，所以这里应该写01010100(最后一位是写)，简单操作可以将0x2a先左移一位，再加0；

步骤2：单片机写从机寄存器地址，device_id的寄存器地址为0x7f，我们读出的fdc2214的device_id应该是0x3055，如果读到这个值了就说明通讯成功了；

步骤3：单片机(master)写fdc芯片(slave)的i2caddress，也是将0x2a先左移一位，不过这次是准备接收device_id了，所以是读，最后一位应为1，所以应该写01010101(最后一位是读)；

步骤4：单片机(master)将sda设为输入，读取fdc芯片(slave)返回值的高8位，对于fdc2214应该是0x30，读完后将sda设为输出，发送ack信号，即将sda拉低；

步骤5：单片机(master)将sda设为输入，读取fdc芯片(slave)返回值的低8位，对于fdc2214应该是0x55，读完后将sda设为输出，发送nack信号，即将sda置高。

控制模块与传感模块采用gpio口模拟的方式进行编程，实现后续单片机与fdc2214芯片之间的通讯。

为了实现上述时序的要求，需进行如下的编程操作：

步骤1确定时钟和延迟时间：

为了较为准确的实现上述时序，首先要写delay函数，开发环境中有自带的延迟函数，void_delay_cycles(unsignedlongcycles)；可以看出该函数是基于系统时钟的，所以第一步是设置系统时钟，编写如下所示的clock_init函数。

一般来说对于msp430f5529比较常用的系统时钟频率为8mhz和25mhz，所以我们可以在某头文件里宏定义cpu_f为8192000或25001984。

#definecpu_f((double)8192000)

#definecpu_f((double)25001984)

有了较为准确的系统时钟，可以继续宏定义延时毫秒和延时微秒的子程序，计算方法也很简单，一个时钟周期是1/cpu_fs，所以延时xμs的话，需要(x/(1/cpu_f))/1000000个时钟周期，即cpu_f*(double)x/1000000.0。

#definedelay_us(x)__delay_cycles((long)(cpu_f*(double)x/1000000.0))

#definedelay_ms(x)__delay_cycles((long)(cpu_f*(double)x/1000.0))

步骤2规定i2c的引脚和主要流程函数：

在i2c.h中定义相关引脚和函数原型声明：

(一)引脚定义：

#definei2cdirp1dir

#definei2coutp1out

#definei2cinp1in

#definesclbit2

#definesdabit3

p1口，p1.2对应scl，p1.3对应sda

(二)将常用的操作用宏定义为更容易理解的符号，包括scl/sda的输入输出的切换，scl/sda的置高置低，sda输入的比特位，如下所示：

#definescl_out()i2cdir|＝scl//scl脚输出

#definescl_h()i2cout|＝scl//scl拉高

#definescl_l()i2cout&＝～scl//scl拉低

#definesda_out()i2cdir|＝sda//sda脚输出

#definesda_h()i2cout|＝sda//sda拉高

#definesda_l()i2cout&＝～sda//sda拉低

#definesda_in()i2cdir&＝～sda//sda脚输入

#definesda_dat()(i2cin&sda)//sda输入数据

(三)i2c流程函数声明：

这里所谓的流程函数指，我们上文写到的一个步骤(frame)里可能出现的基本操作，包括起始(start)操作、停止(stop)操作、接收ack信号、发送ack信号、写一个字节的数据、接收一个字节的数据。

voidi2c_start(void)；//起始操作

voidi2c_stop(void)；//结束操作

voidi2c_sendack(ucharack)；//发送ack

uchari2c_recvack(void)；//接收ack

voidi2c_sendbyte(uchardat)；//写字节

uchari2c_recvbyte(void)；//收字节

在i2c.c中实现这几个基本操作的函数

a.起始操作：

b.停止操作:

c.发送ack:

注意这里的ack有两种，在master读单字节时，收到数据后直接发送nack信号，即将sda置高，如果是读双字节，高8位收到后，master先发送ack信号，即将sda置低，收到低八位后再发送nack信号将sda置高，因此函数中需要加入一个参数ack，ack＝0时发送ack,ack＝1时发送nack。

d.接收ack:

e.发送字节

d.接收字节

步骤3在基本的流程函数基础上进一步定义完整的i2c读写操作：

主要包括单字节写，双字节写，单字节读，双字节读

voidbytewrite(ucharslave_address,ucharreg_address,ucharreg_data)；//单字节写

voidwordwrite(ucharslave_address,ucharreg_address,ucharreg_datam,ucharreg_datal)；//双字节写

ucharbyteread(ucharslave_address,ucharreg_address)；//单字节读

intwordread(ucharslave_address,ucharreg_address)；//双字节读

device_id的步骤明显是一个双字节读的过程，会写这个，其他操作应该是很简单的，双字节读的程序如下所示：

步骤4“刷脸”函数及其他fdc2214的操作函数

在fdc2214.c和fdc2214.h中定义更高层次的操作函数，包括刷脸函数和其他操作函数。

(一)刷脸函数：

刷脸函数fdc2214_begin()可以直接调用上述的wordread()函数，将i2c地址和寄存器地址(0x7f)作为参数传入，就可以读到device_id了，下述程序还兼顾了device_id等于0x3054的fdc2114和fdc2112。

i2c地址应该把0x2a先左移一位，下面是宏定义。

#definefdc2214_i2c_addr_0(0x2a<<1)

#definefdc2214_device_id0x7f

(二)定义其他fdc2214操作函数，原型函数如下所示。

voidfdc2214_init(void)；//初始化

ulongfdc2214_readdata(ucharslave_address,ucharreg_addressf,ucharreg_addresss)；//读数据

ulongfdc2214_channelread(ucharchannel)；//通道读取

voidfdc2214_channelinit(fdcstru_channel*fdcchannel,ulongintethreshold,ulongdthreshold,ucharleak,ucharchannel)；//通道初始化

charfdc2214_objectdetect(fdcstru_channel*fdcchannel,ucharleak,ucharchannel)；//物体探测

在main.c的主循环while(1)中反复调用刷脸函数即可,可以加入一个延时，这样如果通讯正常的话，拿逻辑分析仪看波形的时候，就可以每10ms观察到一个如说明书所示的双字节读取时序。小使用nivirtualbench测试仪中的逻辑分析仪功能，观察时序，结果如下所示。

while(1)

{

fdc2214_begin()；

delay_ms(10)；

}

时序没问题，程序会在测试仪的液晶屏上显示“通讯成功”的字样，如果时序有问题，仔细核对时序，尤其是要观察frame1中fdc2214有没有返回ack信号，即在master写完i2c地址后，释放sda控制权后，fdc2214有没有将sda拉低。

值得注意的是，按键模块采用独立式按键。此类键盘采用独立键盘扫描方式，电路设计简单，编程相对容易，占用cpu时间较少，显示模块采用oled显示屏，i2c通信设计，该显示屏具有自发光，不需背光源、对比度高、厚度薄、视角广、刷新速度快、可用于挠曲性面板、使用温度范围广、构造及制程较简单等特点，电源模块采用采用锂电池供电，充电效率高，比较轻便，使用寿命时间长，同样容量的锂电池，有效使用的能量和续行里程比蓄电池多。

值得说明的是，本实施例中涉及到的压力传感器57可以采用型号为hh8204的微型测力传感器，其配套电源和控制开关也可由厂家提供，除此之外，本实用新型中涉及到电路和电子元器件以及模块均为现有技术，本领域技术人员完全可以实现，无需赘言，本实用新型保护的内容也不涉及对于内部结构和方法的改进。

本实施例的微孔板加样指示装置在使用时，通过传感模块用于对压力信号和手掌信号进行采集，通过按键模块用于对信号数据计数和工作模式进行切换，通过控制模块用于对采集的信号数据进行分析和处理，有利于降低加样出错率，提高实验效率，通过显示模块用于对分析和处理后的信号数据进行显示，通过电源模块用于对电路进行供电，确保电路正常运行。

实施例3：

作为本实用新型的第三种实施例，为了便于监测加样过程，提高加样的准确性，本实用新型人设置定义成功加样模块，作为一种优选实施例，显示面板3内部还设有定义成功加样模块，用于监测加样过程。

本实施例中，定义成功加样模块采用自回归模型对监测的数据信号进行处理计算。

进一步，底座1上表面设有铜板，铜板将置放槽20全部覆盖，便于后续监测加样过程，预防加样出错。

具体的，ar模型(自回归模型)的引入：

由单摆系统，设xt为第t次摆动过程中的摆幅，根据物理原理，第t次的摆幅xt由前一次的摆幅xt-1决定，即有xt＝a1xt-1，考虑到空气振动的影响，我们往往假设

xt＝a1xt-1+εt，t≥1(1)

其中，随机干扰εt～n(0，σ²)。

设初始时刻x0＝1,现在取不同的a1和σ值进行实验，参数a1对序列的稳定性起到决定性的作用，而噪声强度σ²决定了序列的波动程度，在这里，我们称模型(1)为一阶自回归模型，更一般地，可以考虑序列值xt可由前p个时刻的序列值及当前的噪声表出，即

xt＝a1xt-1+a2xt-t+…+apxt-p+εt(2)

其中，aj为参数，{εt}为白噪声，为了显示序列值为随机变量，这里使用xt而不是xt。

ar模型的定义：

如果{εt}白噪声，服从n(0，σ²)，a0,a1,...,ap(ap≠0)为实数，就称p阶差分方程

是一个p阶自回归模型，简称ar(p)模型，称a＝(a0，a1，...，ap)^t是ar(p)模型中的自回归系数，满足ar(p)模型(3)的时间序列{xt}称为ar(p)序列，当a0＝0时，称为零均值ar(p)序列，即

需要指出的是，对于a0≠0的情况，我们可以通过零均值化的手段把一般的ar(p)序列变为零均值ar(p)序列。

ar序列的建模：

对于给定的时间序列{xt}，对其进行建模

步骤①对序列作白噪声检验，若经检验判定序列为白噪声，建模结束，否则转步骤②；

步骤②对序列作平稳性检验，若经检验判定为非平稳，则进行序列的平稳化处理，转步骤①，否则转步骤③；

步骤③对模型进行识别，估计其参数，转步骤④；

步骤④检验模型的适用性，若检验通过，则得到拟合模型并可对序列做预测，否则转步骤③。

ar模型的判定：

对于观测到的时间序列，若通过白噪声检验确定为非白噪声，且经平稳性检验确定为平稳后，我们常根据相关系数和偏相关系数来识别模型，判断该问题是否适用ar模型建模，以及大致确定阶数p。

可通过下面的代码，计算自相关系数(autocorrelationfunction,saf)和偏自相关系数(partialautocorrelationfunction,pacf)

fromstatsmodels.tsa.stattoolsimportacf，pacf

#pacf计算偏自相关系数

#acf计算自相关系数

如果一个时间序列满足以下两个条件

acf具有拖尾性，即acf(k)不会在k大于某个常数之后就恒等0；

pacf具有截尾性，即pacf(k)在k>p时变为0。

第二个条件还可以用来确定阶数p，考虑到存在随机误差的存在，因此pacf在p阶延迟后未必严格为0，而是在0附近的小范围内波动，具体来说

设k阶偏自相关系数为ak，若阶数大于p大部分的偏自相关系数满足下式，则ar模型的阶数取p，其中n表示样本序列长度

例如，对于模型xt＝0.9xt-1-0.3xt-2+εt，它的acf和pacf自相关系数呈现一定的周期性，故判定为拖尾，偏自相关系数2步后截尾，因此，我们可以尝试使用ar(2)模型来建模。

ar模型的参数估计，ar模型的参数估计采用最小二乘估计。

对于样本序列{xt}，当j≥p+1时，记白噪声εj的估计为

通常称为残差，我们的优化目标是使得残差平方和

达到最小，我们称使上式达到最小的为ar(p)模型中自回归系数a1，a2，……，ap的估计。

记

得到如下线性方程组

y＝xa+ε(7)

于是式(6)的目标函数可表示为

s(a)＝(y-xa)^t(y-xa)＝y^ty-2y^txa+a^tx^txa(8)

上式对参数a求导并令其为0，可得

因此，参数a的最小二乘估计为

此时，误差方差的最小二乘估计

ar模型的定阶

在对ar模型识别时，根据其样本偏自相关系数的截尾步数，可初步得到ar模型的阶数p，然而，此时建立的ar(p)未必是最优的，一个好的模型通常要求残差序列方差较小，同时模型页相对简单，即要求阶数较低，因此我们需要一些准则来比较不同阶数的模型之间的优劣，从而确定最合适的阶数，下面给出两种常用的定阶准则。

fpe准则

最终预报误差准则，简称为fpe准则，其判据就是最终预报误差最小，设ar(p)为拟合模型，是序列的各阶样本自协方差函数，其最终预报误差可表示为

在具体应用时，通常是分别建立从低阶到高阶的ar模型，并计算出相应的fpe的值，由此确定使fpe达到最小的p值。

贝叶斯信息准则

定义

使得bic达到最小值的p即为该准则下的最优ar模型的阶数。

ar模型的检验

在模型拟合之后需要进行模型的检验，主要分为两部分

有效性检验：检验拟合模型对序列中信息的提取是否充分

显著性检验：检验模型中的个参数是否显著为0，从而判断拟合魔心是否可以进一步简化。

模型的有效性检验

一个好的拟合模型应该能够提取观测值序列中几乎所有的样本相关信息，即残差序列应该为白噪声序列，因此，模型的有效性检验即为残差序列的白噪声检验，如果残差序列是白噪声，那么理论赏其延迟任意阶的自相关系数为0，考虑数据的偏差，那么绝大多数应该在0附近的范围内，通常在95％的置信水平(2倍标准差以内。

参数的显著性检验

这一部分的目标是，删除那些不显著参数使模型结构最为精简，对于模型参数aj(j＝1,...,p)的检验，其原假设和备择假设分别为

检验统计量为t统计量

在给定的显著水平α下，当检验统计量t大部分位于分点t1-α/2，或该统计量的p值小于α时，则可以1-α的置信水平拒绝原假设，认为模型参数显著，反之，则不能显著拒绝参数为0的假设。

在实施过程中，为了方便加样，通过定义成功加样模块，若加样过程中出错，可手动退回前一步。

以上显示和描述了本实用新型的基本原理、主要特征和本实用新型的优点。本行业的技术人员应该了解，本实用新型不受上述实施例的限制，上述实施例和说明书中描述的仅为本实用新型的优选例，并不用来限制本实用新型，在不脱离本实用新型精神和范围的前提下，本实用新型还会有各种变化和改进，这些变化和改进都落入要求保护的本实用新型范围内。本实用新型要求保护范围由所附的权利要求书及其等效物界定。