中国生物医学工程学报杂志
Chinese Journal of Biomedical Engineering 중국생물의학공정학보
- 主管单位: 中国科学技术协会
- 主办单位: 中国生物医学工程学会
- 影响因子: 0.61
- 审稿时间: 1-3个月
- 国际刊号: 0258-8021
- 国内刊号: 11-2057/R
- 论文标题 期刊级别 审稿状态
-
利用心率的去趋势波动分析提取睡眠呼吸暂停低通气综合症患者的睡眠结构信息
以经基于心率变异性的睡眠分期大都是针对健康人进行的,拟利用心率变异性分析睡眠呼吸暂停低通气综合症(SAHS)病人的睡眠结构.根据微动敏感式床垫式睡眠监测系统(MSMSMS)所得10名健康者的心率和30名不同程度SAHS患者的心率,利用去趋势波动分析(DFA)方法,计算他们在各睡眠阶段的标度指数.实验结果显示,SAHS患者各个睡眠阶段的标度指数具有明显差异,与健康者也有不同,为基于心率变异性对SAHS患者进行睡眠分期以及SAHS患者的判定提供了依据.
-
组织工程化纳米细菌纤维素管体内生物相容性研究
通过肌肉植入试验,初步探讨组织工程化纳米细菌纤维素管的体内生物相容性.利用茶水发酵法培养红茶菌,制备纳米细菌纤维素并使之成管,扫描电镜观察其表面结构.将消毒后的纤维素管植入大鼠大腿肌肉内贴近坐骨神经处,分别于第1、2、6周各随机抽取6只大鼠灌流固定、取材,做大体观察、组织学观察及血液学检查.扫描电镜观察显示,细菌纤维素管壁由细菌纤维素条相互交织,形成不规则的、疏松的网状或絮状多孔结构.材料植入后第1周,可见材料周围有较厚的结缔组织包裹,与周围组织粘连较明显,光镜下有较多的中性粒细胞和少量的单核细胞、淋巴细胞浸润,同时周围肌肉组织炎症细胞浸润明显,以中性粒细胞为主;植入后第2周结缔组织包裹明显缩小,管壁周围的炎性反应明显减轻,以中性粒细胞、淋巴细胞为主,周围肌肉组织的炎性反应亦明显减轻,同时偶可见嗜酸性粒细胞和巨噬细胞;植入后第6周,材料与周围组织分界清楚,仍然有结缔组织包裹,炎细胞数量明显减少,以淋巴细胞为主,周围肌肉组织未见明显炎症细胞浸润,同时材料开始轻度降解.血常规和血生化指标检测显示,各时间点与正常大鼠数据间的差异无统计学意义.组织工程化纳米细菌纤维素管在短期观察中,具有良好的体内生物相容性.
-
基于改进支持向量机方法的蛋白质相互作用预测
蛋白质与蛋白质相互作用研究是蛋白质组学的重要研究内容之一.本研究采用支持向量机学习方法,将氨基酸物理化学特性和序列信息方法相结合构建支持向量,选取DIP数据库中的酵母表达蛋白序列进行蛋白质相互作用预测.在34 000对酵母表达蛋白实验数据中,预测准确率达到83.72%,而单独运用基于氨基酸物理化学特性的方法和基于序列信息的方法预测准确率分别为75.86%和79.63%.在提高预测准确率的同时通过引入离散信息度量函数(FDOD)减少支持向量的维数,使支持向量学习时间缩短,提高相互作用预测的速度.
-
应用搅拌式生物反应器大量扩增人胎盘间充质干细胞的实验研究
探讨应用搅拌式生物反应器培养技术体外扩增人胎盘间充质干细胞(hPDMCs)的方法.以 hPDMCs为种子细胞,利用微载体悬浮法在搅拌式生物反应器内进行体外培养,同时设置静态培养的培养瓶为对照组,观察细胞的扩增速率、细胞代谢(乳酸、葡萄糖、培养液的pH值)的变化,检测生物反应器培养前、后干细胞的表面标记物.hPDMCs 在搅拌式生物反应器内每代可以扩增 (10.55±1.62) 倍,明显高于对照组 (6.10±0.11) 倍的扩增值(P<0.05),且细胞生长代谢指标优于对照组;流式细胞仪检测应用搅拌式生物反应器,培养后细胞表面标记物的表达率无明显改变(P>0.05). 搅拌式生物反应器能够提供良好的细胞生长环境,建立用于 hPDMCs 的大量扩增的体外三维培养体系.
-
颜色辨识状态下ERP时频特征与认知能力的相关性研究
研究不同程度认知能力人群的事件相关电位(ERP)的特征,从中提取出能有效识别血管性痴呆(VD)的电生理指标.选择正常青年组、VD组以及与VD同一年龄段的正常老年组作为被试人群,采用视觉通道诱发的oddball实验模式,对ERP信号运用平滑伪魏格纳分布(SPWVD)进行时频分析.通过不同被试人群的群延迟(GD)和瞬时频率(IF)的对比发现,大脑额、中央和顶区ERP信号的群延迟可反映3组不同被试之间认知能力的差异,为VD患者的早期诊断识别提供一个新的途径.
-
一种基于Gene Ontology注释信息的基因选择算法
基因选择算法是辅助生物学分析重要的方法之一,但这类统计学算法受样本量相对基因数目过少的困扰.提出一种结合Gene Ontology(GO)注释信息的基因选择算法,用GO注释接近基因的方差的加权平均进行修正,增强小样本量下对总体的估计,进而寻找差异表达基因.将该算法与其他5种常见算法对比,以选择出的基因为特征构建分类器,以分类器的可靠性作为衡量算法的标准.3组芯片实验的结果表明,该算法在小样本情况下具有一定优势.亦有Pubmed文献证明,该算法可以鉴别出其他算法未曾发现的致病基因.该方法所建立起来的框架,是把生物学注释信息引入算法改进的一种有效尝试.
-
从肿瘤基因表达数据挖掘分类规则的研究
基于肿瘤基因表达数据,利用信息科学的方法和技术建立肿瘤预测分类模型,对肿瘤基因表达模式研究和肿瘤的诊断识别具有重要意义.本研究提出一种从肿瘤基因表达数据中直接挖掘分类规则建立肿瘤预测分类器的方法.该方法首先抽取实验样本集,分别找出标记肿瘤和正常组织样本的分类特征,由此生成可预测样本类别的分类规则,对每个未知类别样本,按照置信度高原则,选择一个分类规则作为预测结构.本研究的实验数据来自Broad Institute的前列腺癌基因表达数据,实验结果显示该方法的预测精度在90%以上,且同时获得了大量结构透明的分类预测规则,表明本研究的方法是可行的和有效的.
-
多尺度交叉熵及其在血压与肾交感神经活动耦合模式中的应用
多尺度熵反映了单个时间序列在多个尺度上的复杂性,本研究提出了多尺度交叉熵的概念和算法,并探索两个时间序列在多个尺度上的耦合行为.首先,将其应用到一个混沌动力学系统的分析中,即Henon-Henon映射.然后,应用多尺度交叉熵分析清醒和麻醉两种状态下血压与肾交感神经活动的非线性耦合模式.结果表明,肾神经活动和血压之间耦合可能具有分形动力学的特征,其复杂度在麻醉后显著降低,而且两种状态下的耦合模式具有明显的差异.
-
手机辐射与人体头颅相互作用的仿真分析
手机与人体相互作用中涉及两方面问题,一是人体吸收的电磁剂量是否超过限值标准从而对人体造成伤害;二是人体对手机的辐射特性会产生什么影响.采用时域有限差分(FDTD)法经数值仿真进行探讨.结果表明,当手机内置PIFA天线、辐射频率为1 800 MHz、辐射功率为1 W时,人体头颅的总体平均比吸收率未超过ICNIRP颁布的公众基本限值标准,但头颅内任意10 g组织中平均比吸收率的大值却超过了公众基本限值标准.同时,头颅的影响使手机的辐射特性发生变化,当手机贴于左耳时,手机天线的方向图较之于处在自由空间时的方向图而言发生了向上、向左的偏斜,并且辐射到空间的功率约为辐射总功率的60%,被头颅吸收的功率约为辐射总功率的40%.
-
旋转自动插入大肠内窥镜研究
现有的大肠内窥镜由于插入技术复杂,操作困难,检查过程中常给病人造成痛苦,有时甚至出现肠穿孔等医疗事故.针对这一难题,提出了一种具有螺旋结构的旋转自动插入大肠内窥镜新技术方案,阐述了螺旋结构与大肠壁摩擦产生推动力的机理,并建立了摩擦推进动力理论模型.研究设计了内窥镜的主体-镜体和螺旋结构,针对实际要求选择了多种材料制作镜体和螺旋结构,利用人体大肠模型进行了插入试验,试验结果为:内窥镜在直筒状大肠模型中以8 r/min的转速旋转插入时可实现推进速度约为30~45 cm/min,表明所设计的螺旋结构可在大肠模型中,通过旋转获得一定的推进力带动镜体行进.此项研究为今后进一步开发这种新型大肠内窥镜打下了基础.
-
牙周膜细胞在网格型与无纺型聚乳酸纳米纤维支架材料上体外培养的研究比较
研究对比牙周膜细胞在无纺型和网格型聚乳酸纳米纤维膜上的生长行为,探讨支架结构对细胞生长的影响.采用静电纺丝技术,用金属平板或金属网分别接收,得到无纺型和网格型聚乳酸纳米纤维膜;通过SEM观察两种支架形貌差异,并测试比较它们的力学性能.通过MTT测试和SEM观察,比较无纺型和网格型纳米纤维膜对细胞生长的影响.实验结果:网格型膜的纤维直径平均为500~600 nm;无纺型膜的纤维直径大于网格型膜,平均直径约为700 nm,但网格型膜的拉伸断裂应变略大.牙周细胞与支架联合培养的MTT结果显示,与在聚苯乙烯(TCPS) 培养板上的培养比较,两种纳米纤维膜都显示出促进细胞增殖的效果,其中网格膜的促进效果比无纺膜更加明显.SEM观察的结果显示,细胞无法进入无纺型膜内部生长,而网格型膜中由疏松纤维堆积形成的大孔结构则非常有利于细胞进入支架内部,细胞在后者上生长良好.因此,网格型纳米纤维支架是一种优于纤维为完全无纺排布的支架,更适用于组织工程研究.
-
基于双基因分析的结肠癌标志基因选择
针对结肠癌标志基因选择问题,引入一种高得分对(TSP)方法,处理一组包含40个肿瘤和22个正常样本的结肠癌微阵列数据,得到标志基因对(VIP,DARS)并构建双基因分类器.结果显示,该分类器对62例样本的分类准确率可达93.55%.进一步,利用荧光实时定量PCR在10个独立样本上检验所选基因对的分类效果,正确率为100%,表明该方法能有效地选择结肠癌标志基因,对临床诊断有积极的意义.
-
胶囊内窥镜无线能量传输系统的人体安全性研究
采用电磁感应方式为体内胶囊内窥镜供能时,人体必然暴露于发射线圈产生的强电磁场中.本研究采用数值计算的方法,在包含62种人体组织、分辨率为0.33 mm×0.33 mm×2 mm的三维真实人体仿真模型上求解人体不同组织在该电磁场中的感应电流密度和比吸收率,并根据国际非电离辐射防护委员会安全性标准对其进行安全性评估.结果表明,发射功率为25 W,谐振频率为266.5 kHz的能量传输系统产生的电磁场对人体是安全的.本研究为无线能量传输系统奠定了生物安全性基础.
-
利用超声回波特征无创监测微波热疗的实验研究
寻找一种合适的方法区别正常组织与坏死组织,是癌微波热疗应用中所面临的重要难题.本研究探讨微波热疗过程中组织回波时域的变化,以寻找合适的参数,无创监测热疗组织状态的变化.在原理分析的基础上设计了实验系统,采用微波介入的加热方式,采集超声回波信号,对变性区域做回波幅值时域积分、幅值平方时域积分(能量)和损伤数字减影积分,并分析它们的变化趋势.实验结果表明:生物组织回波的不同方式积分值与加热过程中的组织变性具有相关性,微波热疗前后组织回波的幅值、能量、不均匀性均有增大的趋势.利用超声背向散射积分和反应不均匀度的参数进行监测,可得到较理想的结果.
-
基于自相关约束的信号波形估计及其在诱发电位提取中的应用
从观测数据中估计源信号的波形一直是信号处理领域中的一个重要问题.利用源信号的时间自相关函数构造一个非线性方程组,借助大规模数值计算方法,把原来直接估计源信号这样一个较困难的问题,转化为估计迭代初始值与源信号时间自相关函数的问题.使用小波模极大值法,估计迭代初始值与源信号时间自相关函数,再对方程组迭代求解,收敛后的结果即为对源信号的估计.将该方法应用于诱发电位的单通道、单次提取实验,在信噪比为0 dB时,提取结果与源信号的相关系数约为0.93,对幅度与潜伏期的估计都比较准确.实验结果表明,算法性能受迭代初始值估计精度与源信号时间自相关函数估计精度的影响,受前者影响较小,而受后者的影响相对较大.
-
基于先验信息的肺部电阻抗成像算法
电阻抗层析成像 (EIT) 系统逆问题求解的欠定性是导致重构图像分辨率低的主要原因之一.结合先验信息,改善逆问题的不适定性,提高算法鲁棒性是提高成像质量可行的办法.本研究基于Tikhonov正则化算法,结合肺部组织结构、器官电导率分布参数以及肺部呼吸动态结构变化等先验信息,构建正则化矩阵,进行EIT图像重构.研究结果表明,结合先验信息的Tikhonov正则化图像重构算法减小了成像相对误差,改进了图像质量,EIT用于区域性肺部通气变化的检测和监护是可行的.
-
一种循环频率域内的心音包络合成方法
在循环频率域内分析心音的特性,提出一种新的心音包络合成方法.采用循环统计量,研究心音的循环平稳特性,把心音包络表示为循环频率的线性和,线性系数由循环谱估计.以心音分裂为例,进行包络提取仿真.在加性高斯白噪声和随机干扰环境下,对各种包络提取方法进行对比分析.仿真结果表明,心音包络可以在循环频率域内合成.心音具有显著的循环平稳特征,心音包络可以表示成循环频率的线性和.所提出的方法不仅对心音包络的合成与分解给出物理解释,而且具有一定的抗噪声和抗干扰能力.
-
采用多尺度多级组合分类器快速定位乳腺X片中的感兴趣区域
乳腺癌是妇女常见恶性肿瘤之一,早期诊断和治疗是降低乳腺癌患者死亡率的关键. 微钙化是乳腺癌早期的一个重要标志,因此快速准确地找出乳腺X光片中含有微钙化簇的感兴趣区域(ROI)是成功诊断的第一步.乳腺X光片中含有大量无病变区域和少量微钙化区域,形成了一种典型的不对称分类问题.本研究结合大量无病变区域的信息训练多级组合分类器,并借助多尺度方法加快筛选速度,以定位ROI.在真实的数字化X线乳腺照片上的实验表明,该方法在无漏检的情况下,可以排除92.64%的正常区域,而且基于Matlab处理,对于每幅图片的平均处理时间仅为7 s.
-
基于贝叶斯分类器的脉象自动识别方法
传统脉诊依靠医生按压腕部挠动脉脉搏进行脉象识别,具有很强的主观性和模糊性,其准确性与可靠性依赖于医生个人的主观感觉与经验积累,缺乏客观和量化的诊断指标.针对脉象信号复杂性以及脉象特征与脉象类别之间非线性等特点,提出一种基于贝叶斯分类器的脉象自动识别方法,并据此建立脉象定量诊断模型.首先,提取脉象信号的特征参数,创建脉象特征参数-脉象类别数据库,采用少数类合成过采样技术SMOTE结合Tomek links的方法,对数据库进行均衡,使不同的脉象类别具有大致相同的样本;然后基于均衡后的数据库学习贝叶斯网络结构,将得到的马尔可夫毯选择为特征集合并作为贝叶斯分类器的输入,创建脉象信号与类别之间的映射关系模型.通过创建的脉象样本数据库和交叉验证方法,对所提出的方法进行验证.结果表明:所提出的方法可有效识别脉象类型,对于脉位、脉率和脉律的预测准确率都超过90%,是一种有效的脉象定量诊断方法.
-
陡脉冲电场对体外人肝癌细胞SMMC-7721的诱导凋亡作用
为研究陡脉冲电场的诱导凋亡作用,以体外人肝癌细胞SMMC-7721为实验对象,采用Annexin V-FITC/PI双染色并利用流式细胞术检测磷脂酰丝氨酸(PS)外翻,采用琼脂糖凝胶电泳检测DNA降解片段化.实验发现,陡脉冲电场能有效(P<0.05)地使PS外翻,并且微秒级陡脉冲(μsSPEF)比纳秒级陡脉冲(nsSPEF)更有效(P<0.05);陡脉冲电场能使DNA降解片段化,且nsSPEF比μsSPEF更有效.实验结果表明,陡脉冲电场能有效诱导体外人肝癌细胞凋亡,其机理与陡脉冲电场非热效应及其频谱特性有关.μsSPEF主要作用于细胞膜,而可能通过外源性通路诱导细胞凋亡;nsSPEF主要作用于细胞核和线粒体等胞内细胞器,可能通过内源性通路诱导细胞凋亡.实验结果为陡脉冲电场杀伤肿瘤细胞的机制和参数选择提供了依据.
-
大鼠腓肠肌细胞介电谱测量和Cole-Cole数学模型分析
测量大鼠腓肠肌细胞介电谱,利用Cole-Cole方程的曲线拟合建立Cole-Cole数学模型参数.在100 Hz~100 MHz频率范围,使用Agilent 4294A阻抗分析仪对电场与肌纤维方向呈平行和垂直的大鼠腓肠肌进行了阻抗测量,并利用Cole-Cole方程的数值计算,对细胞介电谱、Cole-Cole图、介电损耗因子频谱、电导率虚部频谱和损耗角正切频谱进行曲线拟合分析和残差分析,比较电场与肌纤维方向呈平行或垂直两种情况的介电行为区别.结果表明:骨骼肌细胞的介电弛豫具有两个特征频率,其平行方向的特征频率低于垂直方向;平行方向的介电增量高于垂直方向的介电增量;平行方向的低频极限量高于垂直方向的低频极限量;平行方向和垂直方向的介电常数高频极限量相等;平行方向和垂直方向的电导率高频极限量存在一定差距.平行方向的电导率虚部频谱和损耗角正切分别高于垂直方向的值,但是,平行方向的电导率虚部大值低于其垂直方向的值.腓肠肌细胞的介电常数和电导率存在频率依存性和各向异性,其介电行为存在α和β两个介电散射,并满足Cole-Cole数学模型.
-
基于弹性小球动态平衡的血管中心线和直径的测量
血管中心线及直径的测量是血管类疾病诊断和微创介入治疗的重要依据.当前研究主要依据图像细化算法,但是细化算法缺乏计算结果的一致性.提出一种基于弹性小球动态平衡的血管中心线和直径的测量方法.首先在血管内部以交互的方式输入一条多段线(polyline)来确定感兴趣的血管部位,然后对多段线离散化为一组点,通过对离散点逐步膨胀为小球以及在小球与血管壁碰撞的两个力影响下,对小球和多段线同时变形,后在平衡状态下得到的多段线即是血管中心线,小球直径就是该处血管直径.实验中,通过对模拟血管图像和临床DSA图像的量化测量及与常用方法对比,表明所提出算法能够精确地得到血管中心线与血管直径.
-
后碰撞中人体颈椎骨的有限元分析
对建立并验证的全颈椎有限元模型进行仿真以分析后碰撞中颈部的响应及其损伤机制.模型包括C1-T1的8块椎骨以及椎间盘、韧带和肌肉等软组织,共35 042个节点和22 618个单元.运用后碰撞志愿者实验数据对模型进行仿真,得到各椎骨间相对转角、各椎骨应力、应变曲线及应力云图.颈部S形曲线、C形曲线得到体现;大应力、应变主要发生在颈部过度伸展过程中,上位颈椎中C2应力、应变较大,而下位颈椎中T1数值大,受力较大值发生在椎体上下表面及其关节突关节面上,易产生多部位骨折.分析结果可为人体颈部损伤的预防、诊断和治疗提供理论依据.
-
一种乳腺钼靶X线图像的软拷贝显示优化方法
数字乳腺钼靶X线图像比胶片乳腺X线图像能更有效地表现和检测乳腺病变,本研究提出了针对DICOM 格式乳腺钼靶X线图像的高效软拷贝显示方法.通过提取图像中的乳腺区域获得对初始窗口参数的准确估计,使乳腺组织在图像初始显示时就有比较好的对比度.利用乳腺区域分割结果,交互调窗操作时只更新乳腺区域内部像素的显示灰度,背景区域像素灰度保持不变;用整数运算和移位运算实现了多项式插值函数,优化了频繁调用的图像缩放操作.算法能显著提高乳腺钼靶X线图像软拷贝显示效率.
-
基于光力学实验方法的生物膜材料性能的研究
为研究细胞培养膜在外力作用下的应变分布,将影栅云纹法应用于培养膜的力学性能测试中,实现对培养膜进行光力学实验研究.实验可得到生物培养膜的离面位移,与相移技术结合能够快速得到培养膜的三维形貌,并利用干涉云纹图得到培养膜的应变分布.应用影栅云纹法,可测量低弹模材料的大变形,较好地解决了大变形力学量的测试问题,不仅对观察附着其上的细胞生长和发育规律有一定的指导意义,而且促进了光力学实验方法在生物工程领域中的应用.
-
电阻抗成像交叉测量协议的灵敏度性能
电阻抗重构图像技术对噪声敏感.测量数据中比较小的误差可能在终图像上导致很大的误差.测量技术应提高对被测场域的灵敏度.在二维圆形均匀场域和相邻电流激励条件下,研究相邻和交叉测量协议的灵敏度性能.相比相邻测量协议,交叉测量协议提高整体灵敏度,同时也提高了中心区域的灵敏度比例值.交叉测量的合成灵敏度峰值是相邻测量的4.54倍.中心处相邻测量的灵敏度比例值是10.7%,交叉测量时为27.2%.实验结果表明交叉测量协议在灵敏度性能方面优于相邻测量协议.
-
基于纤维增强超弹性基体模型的骨骼肌有限元模拟
为了进一步研究骨骼肌复杂的机械动力行为,采用有限元方法对骨骼肌行为进行力学分析.将骨骼肌视为肌纤维和超弹性基体构成的复合材料,纤维应力分为主动和被动两部分,主动应力依赖于疲劳状态、激励波形、收缩速度和现时长度,而被动应力与普通材料相同,仅依赖于现时应变.利用商用软件ABAQUS,对蛙大腿肌进行数值模拟,比较在周期激励作用下,不同时刻的单轴应变和代表性节点的位移分布规律.结果显示,骨骼肌表现出大变形特性,单轴大应变达到69%.而且由于引入了疲劳函数,所采用的数值分析技术可模拟骨骼肌的主动与被动机械行为,且容易实现,结果稳定可靠.
-
运动在模拟失重尾吊大鼠模型中的应用及对肌骨系统改善的研究进展
大鼠作为模拟人类骨骼状态的常用实验动物模型,被广泛用来探寻对抗骨质疏松的有效措施,其中运动疗法既安全又有效.宇航员长期处于微重力环境会发生失重性骨质疏松,相应的大鼠尾吊模型作为地面模拟微重力效应的实验模型已得到广泛应用.文中综述了应用在尾吊和非尾吊大鼠实验中的运动训练,力图涵盖运动训练对抗大鼠骨质疏松的国内外研究进展,并对今后模拟微重力下大鼠运动方式的研究方向做了展望.
-
深部脑刺激对神经精神疾病的治疗与未来展望
随着社会老龄化和日渐加深的心理压力等因素的影响,世界范围内包括帕金森症、肌张力失常、强迫症、抑郁症、癫痫等与神经回路异常相关的神经精神系统疾病患者的数量急剧增加,目前全球患者已逾五亿六千万人.深部脑刺激技术(DBS)已广泛成为治疗这些疾病的重要方法.尽管DBS的作用机制尚不完全明了,但由于其具有微损伤、可回复和可调节的优点,已成为一项可靠、有效的神经外科手术方法.文中对DBS装置、作用机制、疾病治疗和发展方向等方面进行阐述,同时简述了其他治疗方法(如经颅磁刺激),并对目前广泛应用于在体神经回路研究的光感基因神经调控技术(光遗传技术)对神经精神系统疾病的干预进行了展望.
年 | 期数 |
2018 | 01 02 03 04 05 06 |
2017 | 01 02 03 04 05 06 |
2016 | 01 02 03 04 05 06 |
2015 | 01 02 03 04 05 06 |
2014 | 01 02 03 04 05 06 |
2013 | 01 02 03 04 05 06 |
2012 | 01 02 03 04 05 06 |
2011 | 01 02 03 04 05 06 |
2010 | 01 02 03 04 05 06 |
2009 | 01 02 03 04 05 06 |
2008 | 01 02 03 04 05 06 |
2007 | 01 02 03 04 05 06 |
2006 | 01 02 03 04 05 06 |
2005 | 01 02 03 04 05 06 |
2004 | 01 02 03 04 05 06 |
2003 | 01 02 03 04 05 06 |
2002 | 01 02 03 04 05 06 |
2001 | 01 02 03 04 05 06 |
2000 | 01 02 03 04 |
1999 | 01 02 |
1998 | 01 02 |
1994 | 04 |