CN103809198B - 一种测量粒子加速器束流位置的方法 - Google Patents

一种测量粒子加速器束流位置的方法 Download PDF

Info

Publication number
CN103809198B
CN103809198B CN201410080631.9A CN201410080631A CN103809198B CN 103809198 B CN103809198 B CN 103809198B CN 201410080631 A CN201410080631 A CN 201410080631A CN 103809198 B CN103809198 B CN 103809198B
Authority
CN
China
Prior art keywords
signal
electrode
matrix
passage
beam position
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410080631.9A
Other languages
English (en)
Other versions
CN103809198A (zh
Inventor
赖龙伟
冷用斌
阎映炳
陈之初
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Alliance Investment Ltd
Shanghai Institute of Applied Physics of CAS
Original Assignee
Shanghai Institute of Applied Physics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Institute of Applied Physics of CAS filed Critical Shanghai Institute of Applied Physics of CAS
Priority to CN201410080631.9A priority Critical patent/CN103809198B/zh
Publication of CN103809198A publication Critical patent/CN103809198A/zh
Application granted granted Critical
Publication of CN103809198B publication Critical patent/CN103809198B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

本发明涉及一种测量粒子加速器束流位置的方法，其包括：步骤S1，通过N个通道对条带束流位置探头的N个电极在束流通过所述探头时输出的N个电极信号进行采样；步骤S2，将所述N个电极信号组成一M*N的矩阵X，并对所述矩阵X进行主成分分析；步骤S3，计算得到所述N个电极分别对应的束流信号感应幅度，并通过所述N个通道分别输出N个所述束流信号感应幅度；步骤S4，对所述N个通道输出的N个所述束流信号感应幅度进行差比和运算，以获得束流的位置信号。本发明通过综合多个电极探头的电极信号组成矩阵，以利用主成分分析法对这些电极信号同时进行处理，从而不仅能够有效消除各电极信号中的噪声，而且还能有效提高束流位置的测量精度。

Description

一种测量粒子加速器束流位置的方法
技术领域
本发明涉及加速器物理束流诊断领域，尤其涉及一种测量粒子加速器束流位置的方法。
背景技术
对于多电极的束流位置探测器，以图1所示直线电子加速器上的带有四个电极A、B、C、D的条带束流位置探测器1（探头）（BPM）为例（图1中标记“1”表示沿穿过纸面方向的束流），其设计频率达到几百兆赫兹，为获得束团通过时的位置信号，需要对探头的四个电极输出的信号先分别进行处理之后再做如式（1）、（2）的差比和运算，从而获得水平和垂直方向的位置信息X、Y。
X=kX(B-D)/(B+D)（1）；
Y=ky(A-C)/(A+C)（2）；
其中，kX、ky为标定系数。
当前处理电极信号的处理方法是先利用由射频器件搭建的前端（即主要进行放大、衰减、混频、滤波等处理，使信号适合进行数字化采样处理的处理模块，其为探头输出之后的第一个处理模块）对各个探头的电极信号分别进行带通滤波处理，然后以一百多兆采样率进行带通欠采样，接着用正交解调提取信号幅度包络并计算积分，从而得到用于进行后续差比和运算的数值。
由于四通道不可避免存在信号网络的插损和对信号的采集处理增益的不一致性问题，而以上处理方法是对各通道分别进行信号处理，因此无法解决该问题，从而会导致位置计算结果存在误差。
发明内容
为了解决上述现有技术存在的问题，本发明旨在提供一种经济、简单且现实可行的测量粒子加速器束流位置的方法，以提高束流位置的测量质量及精度。
本发明所述的一种测量粒子加速器束流位置的方法，其包括以下步骤：
步骤S1，通过N个通道对条带束流位置探头的N个电极在束流通过所述探头时输出的N个电极信号进行采样，其中，每个电极信号包含M个不同的基模信号；
步骤S2，将所述N个电极信号组成一M*N的矩阵X，其中，以所述每个电极信号的M个不同的基模信号为列向量，并对所述矩阵X进行主成分分析，以获得由各个所述基模信号的强度组成的M*N的对角线矩阵S、以及由各个所述电极对各个基模信号的感应强度系数组成的N*N的矩阵V；
步骤S3，根据以下公式计算得到所述N个电极分别对应的束流信号感应幅度，并通过所述N个通道分别输出N个所述束流信号感应幅度，
( V i , 1 * S 1,1 ) 2 + ( V i , 2 * S 2,2 ) 2 ,
其中，Vi,1和Vi,2为所述矩阵V中第i列向量前两个数据，1≤i≤N，S1,1和S2,2为所述对角线矩阵S中对角线上最大的前两个数据；
步骤S4，对所述N个通道输出的N个所述束流信号感应幅度进行差比和运算，以获得束流的位置信号；
其中，M为自然数，N为偶数，M大于N。
在上述的测量粒子加速器束流位置的方法中，所述步骤S1包括：先通过射频前端对所述N个电极信号进行带通滤波，然后对经过带通滤波后的信号进行带通采样。
在上述的测量粒子加速器束流位置的方法中，所述步骤S2包括：对所述矩阵X进行矩阵奇异值分解，以获得所述对角线矩阵S和所述矩阵V。
由于采用了上述的技术解决方案，本发明并不是对每个电极探头采样的数字信号单独处理，而是充分利用了多探头带来的系统处理增益，通过综合多个电极探头的电极信号组成矩阵，以利用主成分分析法对这些电极信号同时进行处理，从而不仅能够有效消除各电极信号中的噪声，而且仅需要对经过主成分分析后获得的数据进行常规的幅度计算，即可直接获得束流通过时各电极探头所对应的束流信号感应幅度，进而有效了提高束流位置的测量精度。
附图说明
图1是现有技术中条带束流位置探测器的原理图。
具体实施方式
下面结合附图，给出本发明的较佳实施例，并予以详细描述。
本发明，即一种测量粒子加速器束流位置的方法，其包括以下步骤：
步骤S1，通过N个通道对条带束流位置探头的N个电极在束流通过探头时输出的N个电极信号进行带通采样，并且在采样前，先通过射频前端对N个电极信号进行带通滤波，其中，每个电极信号包含M个不同的基模信号；
步骤S2，将N个电极信号组成一M*N的矩阵X，其中，以每个电极信号的M个不同的基模信号为列向量，并对矩阵X进行主成分分析，即，对矩阵X进行矩阵奇异值分解（SVD分解），以获得由各个基模信号的强度组成的M*N的对角线矩阵S、以及由各个电极对各个基模信号的感应强度系数组成的N*N的矩阵V；
步骤S3，根据以下公式计算得到N个电极分别对应的束流信号感应幅度，并通过N个通道分别输出N个束流信号感应幅度，
( V i , 1 * S 1,1 ) 2 + ( V i , 2 * S 2,2 ) 2 - - - ( 3 )
其中，Vi,1和Vi,2为矩阵V中第i列向量前两个数据，1≤i≤N，S1,1和S2,2为对角线矩阵S中对角线上最大的前两个数据；
步骤S4，对N个通道输出的N个束流信号感应幅度进行差比和运算，以获得束流的位置信号；
其中，M为自然数，N为偶数，M大于N。
下面对本发明的原理进行详细说明。
由于理想的电极对同一个束团（或束团串）的响应是完全相同的，因此所有电极在时域中获得的信号波形应该完全相似。有束团通过时，各电极由于相对于束团位置的不同、信号网络插损或通道增益的不同，输出信号可以表示为一个基矢量与不同比例系数的乘积，而信号中与基矢量不一致的部分均可以认为是系统测量误差引入的噪声。
以带有N（例如N=4）个电极的条带束流位置探头（BPM）为例，首先，经过射频前端带通滤波处理和带通采样之后，输出的N个电极信号组成了矩阵X=[X1,X2,…,XN]，其中，X1，X2,…,XN分别是N个电极在束流（束团）通过探头时采样到的一组数据组成的行向量，假设每组数据包含有M（例如M=1024）点，即M个不同的基模信号，则矩阵X为M*N矩阵（例如1024*4矩阵）。
然后，对矩阵X进行矩阵奇异值分解（SVD）分解（SVD分解是数学上经常使用一种现有方法），即：
X=U*S*V’（4）
式（4）中，矩阵U为M*M矩阵（例如1024*1024矩阵），是由矩阵X分解出的时间向量，矩阵U中的M个列向量代表电极信号的M个不同的基模信号，其中前N个基模信号是从电极信号分解出的有效信号，而后M-N个（例如1020个）基模信号则是无意义的数学计算信号；矩阵S为M*N对角线矩阵（1024*4对角线矩阵），其具有N个按从大到小顺序排列的特征值，其他位都为0，特征值的大小代表了该基模信号的强度，根据特征值大小即可分析出电极信号中的主成分基模信号（即，特征值大的为主成分）；矩阵V’是矩阵V的转置矩阵，矩阵V为N*N矩阵（例如4*4矩阵），是由矩阵X分解出的空间向量，矩阵V中的N个列向量代表前N个基模信号在个N个电极上的感应强度系数。
在N个有效的基模信号中，信号强度最大的第一、二个基模信号是主成分信号，即分别为束流通过时探头感应信号的正弦和余弦部分，两者强度大小是由信号的相位决定的，而其他的基模信号（例如第三、第四个基模信号）则分别是噪声信号和反射信号。因此，第i个通道（1≤i≤N）感应的束流信号强度可通过上式（3）得出。
最后，采用常用的差比和运算处理方法即可获得束流的位置信号，此处不再赘述。
在本发明中，主成分分析法是一种已有的分析方法，被广泛应用在模式识别、数据挖掘等领域。但在加速器领域，现有技术中技术人员仅利用该方法进行全环的BPM数据相关性分析，而从未用于各BPM的信号处理。BPM各通道的信号处理最常用的还是软件无线电中的正交解调，该方法有些噪声无法消除，在算法研究过程中发现四通道之间可以利用主成分分析算法提取出主成分信号，有效消除噪声的干扰。由于各学术领域之间的交流并不是很及时，有些在其他领域广泛应用的算法在加速器领域很可能还没有发现它的用处。
综上所述，本发明针对经过射频前端处理后的数字信号在原有设备基础上利用现有的主成分分析法对多电极探头输出的信号组成的矩阵进行主成分分析处理，而没有采用通常的IQ解调方法计算信号幅度包络，因此不需要做正交解调、滤波、抽取、幅度计算等一系列数字信号处理，而只需实现SVD分解和幅度计算就可获得束流通BPM过时探头感应的信号幅度，从而有效消除了噪声和反射信号等对系统性能的影响，提高了信号处理质量，进而有效了提高束流位置的测量精度。
以上所述的，仅为本发明的较佳实施例，并非用以限定本发明的范围，本发明的上述实施例还可以做出各种变化。即凡是依据本发明申请的权利要求书及说明书内容所作的简单、等效变化与修饰，皆落入本发明专利的权利要求保护范围。本发明未详尽描述的均为常规技术内容。

Claims (2)

1.一种测量粒子加速器束流位置的方法，其特征在于，所述方法包括以下步骤：
步骤S1，通过N个通道对条带束流位置探头的N个电极在束流通过所述探头时输出的N个电极信号进行采样，其中，每个电极信号包含M个不同的基模信号；
步骤S2，将所述N个电极信号组成一M*N的矩阵X，其中，以所述每个电极信号的M个不同的基模信号为列向量，并对所述矩阵X进行矩阵奇异值分解，以获得由各个所述基模信号的强度组成的M*N的对角线矩阵S、以及由各个所述电极对各个基模信号的感应强度系数组成的N*N的矩阵V；
步骤S3，根据以下公式计算得到所述N个电极分别对应的束流信号感应幅度，并通过所述N个通道分别输出N个所述束流信号感应幅度，
( V i , 1 * S 1 , 1 ) 2 + ( V i , 2 * S 2 , 2 ) 2 ,
其中，Vi,1和Vi,2为所述矩阵V中第i列向量前两个数据，1≤i≤N，S1,1和S2,2为所述对角线矩阵S中对角线上最大的前两个数据；
步骤S4，对所述N个通道输出的N个所述束流信号感应幅度进行差比和运算，以获得束流的位置信号；
其中，M为自然数，N为偶数，M大于N。
2.根据权利要求1所述的测量粒子加速器束流位置的方法，其特征在于，所述步骤S1包括：先通过射频前端对所述N个电极信号进行带通滤波，然后对经过带通滤波后的信号进行带通采样。
CN201410080631.9A 2014-03-06 2014-03-06 一种测量粒子加速器束流位置的方法 Active CN103809198B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410080631.9A CN103809198B (zh) 2014-03-06 2014-03-06 一种测量粒子加速器束流位置的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410080631.9A CN103809198B (zh) 2014-03-06 2014-03-06 一种测量粒子加速器束流位置的方法

Publications (2)

Publication Number Publication Date
CN103809198A CN103809198A (zh) 2014-05-21
CN103809198B true CN103809198B (zh) 2016-01-20

Family

ID=50706247

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410080631.9A Active CN103809198B (zh) 2014-03-06 2014-03-06 一种测量粒子加速器束流位置的方法

Country Status (1)

Country Link
CN (1) CN103809198B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104180824A (zh) * 2014-08-18 2014-12-03 中国科学院上海应用物理研究所 一种基于主成分分析算法提高探头测量精度的方法
CN104166152B (zh) * 2014-08-18 2016-07-06 中国科学院上海应用物理研究所 一种检测粒子加速器束流异常的自触发方法
CN104506293B (zh) * 2014-12-30 2017-09-15 中国科学院上海应用物理研究所 载波抑制射频前端和方法、束流位置测量系统和方法
CN105676261B (zh) * 2016-01-13 2018-07-13 中国科学院上海应用物理研究所 一种测量粒子加速器束流流强的系统和方法
CN106330617A (zh) * 2016-09-18 2017-01-11 中国科学院上海应用物理研究所 基于主成分分析算法的信号传输质量诊断方法
CN106855634A (zh) * 2016-12-29 2017-06-16 中国科学技术大学 一种紧凑结构的钮扣型束流位置检测器组件
CN107456664B (zh) * 2017-08-22 2021-05-11 哈尔滨工程大学 一种半球形可伸缩3d光纤加速器晨检仪
CN108897033A (zh) * 2018-06-15 2018-11-27 东莞中子科学中心 用于加速器的束流位置探测器
CN110261886B (zh) * 2019-04-29 2020-12-25 中国科学院高能物理研究所 利用数字bpm测量直线加速器束流位置的方法及系统
CN112698096B (zh) * 2020-12-09 2023-07-25 中国科学院上海高等研究院 一种逐束团三维位置测量系统及测量方法
CN113433582B (zh) * 2021-05-13 2022-08-23 上海交通大学 一种x射线球管束流诊断方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1918484A (zh) * 2004-02-10 2007-02-21 松下电器产业株式会社 束测量装置和使用束测量装置的束测量方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3819402B2 (ja) * 2004-07-14 2006-09-06 三菱電機株式会社 ビーム位置測定システム
JP3860184B2 (ja) * 2004-07-14 2006-12-20 三菱電機株式会社 ビーム位置測定システム
JP5164110B2 (ja) * 2008-09-05 2013-03-13 独立行政法人放射線医学総合研究所 粒子線治療装置のビーム位置モニタ及びビーム位置測定方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1918484A (zh) * 2004-02-10 2007-02-21 松下电器产业株式会社 束测量装置和使用束测量装置的束测量方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Model-independent analysis of the Fermilab Tevatron turn-by-turn beam position monitor measurements;A.V.Petrenko;《PHYSICAL REVIEW SPECIAL TOPICS-ACCELERATORS AND BEAMS》;20110930;第14卷(第9期);第092801-1-12页 *
Principal Component Analysis of Cavity Beam Position Monitor Signals;Y.I.Kim;《HOURNAL OF INSTRUMENTATION》;20140228;第9卷;第1-12页 *
数字束流位置信号处理算法优化;赖龙伟 等;《强激光与粒子束》;20130131;第25卷(第1期);第189-192页 *
用束流位置监测器测量上海光源束流寿命;阎映炳 等;《强激光与粒子束》;20120131;第24卷(第1期);第109-113页 *

Also Published As

Publication number Publication date
CN103809198A (zh) 2014-05-21

Similar Documents

Publication Publication Date Title
CN103809198B (zh) 一种测量粒子加速器束流位置的方法
Azais et al. Spike detection from inaccurate samplings
CN101919695B (zh) 一种基于小波变换的心电信号qrs波检测方法
CN102749381B (zh) 线性扫描极谱曲线重叠峰分离方法
CN110197209A (zh) 一种基于多特征融合的辐射源识别方法
Fouda et al. Experimental chaos detection in the Duffing oscillator
CN105629290B (zh) 一种数字核脉冲信号墨西哥草帽小波成形方法
CN104062642A (zh) 一种对激光雷达波形数据进行高斯回波分解的方法
CN102707252B (zh) 一种去除高速取样示波器时基抖动的方法
CN104215833A (zh) 电力系统频率测量方法及装置
CN104655965B (zh) 一种电力系统中相量测量方法
CN103245830A (zh) 一种结合ar谱估计与非线性优化的间谐波检测方法
Zhou et al. Geomagnetic sensor noise reduction for improving calibration compensation accuracy based on improved HHT algorithm
CN105717490B (zh) 基于时频分析的lfm信号分离及参数估计方法
CN104181508B (zh) 基于压缩感知的威胁雷达信号检测方法
CN104021275A (zh) 一种新型心电图相似性度量方法
Washizawa et al. A flexible method for envelope estimation in empirical mode decomposition
CN106501582B (zh) 信号基线处理方法及信号基线处理设备
CN104180824A (zh) 一种基于主成分分析算法提高探头测量精度的方法
CN107341519A (zh) 一种基于多分辨率分析的支持向量机识别优化方法
CN105022917B (zh) 一种信号精确提取与处理方法
CN111008356B (zh) 一种基于WTSVD算法扣除背景的γ能谱集分析方法
CN106598910A (zh) 一种emd端点效应抑制方法及系统
CN104268838A (zh) 一种面向超光谱数据库的小波去噪算法
CN105405444B (zh) 一种在Odd-DFT域对含噪正弦信号进行参数估计的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20140521

Assignee: SHANGHAI AIPUQIANG PARTICLE EQUIPMENT Co.,Ltd.

Assignor: SHANGHAI INSTITUTE OF APPLIED PHYSICS, CHINESE ACADEMY OF SCIENCES

Contract record no.: X2019310000019

Denomination of invention: Method for measuring particle accelerator beam position

Granted publication date: 20160120

License type: Exclusive License

Record date: 20191113

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221212

Address after: 201800 Shanghai city Jiading District Road No. 2019

Patentee after: SHANGHAI INSTITUTE OF APPLIED PHYSICS, CHINESE ACADEMY OF SCIENCES

Patentee after: Shanghai Alliance Investment Ltd.

Address before: 201800 Shanghai city Jiading District Road No. 2019

Patentee before: SHANGHAI INSTITUTE OF APPLIED PHYSICS, CHINESE ACADEMY OF SCIENCES