CN110673523B - 一种平衡三进制与索引算法结合的约瑟夫森结阵驱动方法 - Google Patents
一种平衡三进制与索引算法结合的约瑟夫森结阵驱动方法 Download PDFInfo
- Publication number
- CN110673523B CN110673523B CN201910903099.9A CN201910903099A CN110673523B CN 110673523 B CN110673523 B CN 110673523B CN 201910903099 A CN201910903099 A CN 201910903099A CN 110673523 B CN110673523 B CN 110673523B
- Authority
- CN
- China
- Prior art keywords
- bias
- matrix
- junction
- bias state
- junction array
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/04—Programme control other than numerical control, i.e. in sequence controllers or logic controllers
- G05B19/042—Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
- G05B19/0423—Input/output
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/25—Pc structure of the system
- G05B2219/25257—Microcontroller
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Superconductor Devices And Manufacturing Methods Thereof (AREA)
Abstract
本发明提供了一种基于SNS型约瑟夫森结阵的驱动方法,属于计量领域。根据SNS型约瑟夫森结阵的驱动原理以及结阵分段的特殊性,介绍以平衡三进制与索引算法相结合的驱动算法。首先,将对应的电压值转换成所需的结阵数;其次,根据约瑟夫森结阵的分段特点,使用平衡三进制算法高效并快速地将所需结阵数转换成偏置状态,将有序的偏置状态组合成偏置状态矩阵。再次,当发生捕获磁通的现象时,平衡三进制算法不再适用,索引算法可以替代平衡三进制驱动算法,在捕获磁通的状态下依旧的得到有效的计算结果。最后,输出所需的偏置状态矩阵。
Description
技术领域
本发明属于计量领域,具体涉及一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法。
背景技术
可编程约瑟夫森(PJVS)交流量子电压的合成是通过软件编程的方式,采用模数变换输出偏置电流,控制不同约瑟夫森结阵的偏置状态,实现随时间变化的量子电压台阶。目前,国际上PJVS结阵的驱动方法以二进制分段驱动为主,PJVS结阵为一系列二进制排列的约瑟夫森结阵的分段组合构成,即各段结阵的单结数量按1,2,4,8......排列。新型SNS型2V结阵共分20段:2916、108、36、972、12、4、324、8742、8744、8744、8744、8744、8742、324、4、12、972、36、108、2916,共61204个结。该型结阵分段采用非二进制方式,且每段结所含结的个数都是偶数,对前7段或者后7段结由小到大排序并除以4后,可以得到一个三进制的序列,传统二进制分段的驱动方法已经无法计算输出PJVS结阵的偏置状态组合方式;同时,PJVS结阵在工作过程中,可能受到外部的干扰信号的影响,造成结阵捕获磁通现象,失去超导量子台阶,此时基于三进制的算法同样不适用。
发明内容
本发明的目的在于解决上述现有方法中存在的难题,通过研究结阵的驱动原理以及结阵分段的特点,提供一种平衡三进制与索引算法相结合的可编程约瑟夫森结阵驱动方法,从而能够随时间有序地输出DAC控制结阵的偏置状态,实现交流量子电压的精确合成。在捕获磁通的状态下,需要在计算时关闭失去超导特性的结阵分段时,通过索引算法可以输出偏置状态,作为平衡三进制驱动算法的补充。
本发明是通过以下技术方案实现的:
一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,包括:
步骤1,初始化偏置状态矩阵B1、B2,初始化两路正弦波信号,分别为信号1和信号2,设置台阶数、微波频率与信号参数。
步骤2,采用N个台阶数对两路正弦波信号进行离散,将每个离散点的电压值转换成所需结阵数输入至向量n1(i)、n2(i)。
步骤3,判断是否发生捕获磁通现象,若未发生捕获磁通现象则进入步骤4,否则,进入步骤5。
步骤4,采用平衡三进制算法计算量子电压序列对应的偏置状态矩阵。
步骤5,计算时关闭异常结段,采用索引算法计算偏置状态矩阵。
步骤6,输出偏置状态矩阵B1、B2。
所述步骤1是这样实现的:
初始化偏置状态矩阵,将两个偏置状态矩阵B1、B2分别置零,设置台阶数N=20,矩阵维度分别为N×9与N×11;设置微波的频率以及信号1与信号2的幅值与相位。
所述偏置状态矩阵每个元素的取值范围是{-1,0,1},在初始化时矩阵所有元素置零,所述N为离散台阶个数。
所述步骤2是这样实现的:
根据已知的正弦信号ysin1、ysin2分别从零点开始等间隔取N个点的电压值,并将这些电压值通过舍入和取整得到与之最接近的偶数值,这个偶数值即为所述的所需结阵数,对于两路信号第i个点的所需结阵数分别为n1(i)与n2(i)
所述步骤3是这样实现的:
检查所述SNS型约瑟夫森结阵是否发生捕获磁通现象,若出现,则进入步骤4,若未发生,则进入步骤5。
所述捕获磁通现象包括:①某一段或者几段结的正负台阶或者0台阶宽度变窄;②某一段或者几段结无正负台阶或者0台阶;③输出的某一量子电压台阶具有一定的斜率,量子电压台阶不稳定;
所述步骤4是这样实现的:
(41)判断偶数n1、n2的正负,分别存储在sign_1与sign_2中,所述sign_1与sign_2取值范围为{1-1},并将n1、n2取绝对值;
(42)以n2为例,n1同理,判断n2能否被4整除,若能进入(43),若不能进入(44);
(43)N2为n2整除4得到,通过循环判断N2的绝对值是否大于1093得到需要加入几段分段数为8744的结阵分段,所述分段数为8744的结阵分段有第9段、第10段以及第11段。这样就能够确定所述第9到11段结的偏置状态。
(44)N2通过N2=(n2-8742)/4得到,此时将第八段结的偏置状态设置为正偏,即B2(8)=1,通过循环判断所述第9到11段结的偏置状态。
(45)确定第9到第11段结偏置状态后,确定前7段结的偏置状态,采用平衡三进制计算,令A=N2,循环七次,Rem为A与3的余数,判断余数绝对值是否等于2。若等于2,则A=(A-Rem)/3+sign(Rem),B(i)=-sign(Rem);否则A=(A-Rem)/3,B(i)=Rem。所述N2是经过之前确定9至11段结运算后的N2,所述Rem为A与3的余数,所述sign(Rem)为Rem的正负性,取值是1或-1,所述B(i)为第i段结的偏置状态。
所述步骤5是这样实现的:
(51)分别将两个结阵的所有情况遍历,分别有177147(3^11)和19683(3^9)种情况。以177147为例,将177146对应的十进制数转换成三进制的序列后,序列每一个元素减一存储至B1_all,将B1_all的每一行点乘JJS_1得到所有情况的偏置结总数B1M,将结果B1M与B1_all存储至索引矩阵search_1中。
所述B1_all为[-1-1-1-1-1-1-1-1-1-1-1]到[1 1 1 1 1 1 1 1 1 1 1]的所有向量组成的偏置矩阵情况,B1M为组合偏置结总数的所有情况矩阵,JJS_1为结段向量,JJS_1=[2916 108 36 972 12 4 324 8742 8744 8744 8744];
(52)分别输入所需的结阵数n1、n2,在索引矩阵中找到与所需结阵数差值最小的结阵数,并且找到对应的偏置矩阵。
所述步骤6是这样实现的:
将步骤4或步骤5运行的结果输出,输出两个偏置矩阵B1、B2。
与现有技术相比,本发明的有益效果是:
根据SNS型约瑟夫森结阵分段数的特点,提出了平衡三进制算法,快速、有效地合成了控制约瑟夫森结阵输出量子电压的偏置矩阵;同时,提出了索引算法,在约瑟夫森结阵发生捕获磁通的现象时能够计算关闭异常结阵并输出量子电压偏置矩阵。为实现两路交流量子电压的输出提供了高效完备的驱动方法。
附图说明
图1为平衡三进制算法框图
图2为索引算法框图
图3为驱动方法总框图
具体实施方式
下面结合附图对本发明作进一步详细描述:
利用所述基于SNS型约瑟夫森结阵的驱动方法,包括:
(1)初始化和设置参数:设置微波频率microwave_freq为17.8GHz;约瑟夫森常数Kj_90=483597.9;台阶数N=20;信号1与信号2的幅值分别为Amp_1=1V和Amp_2=0.5V、相位分别为phase_sin_1=0°和phase_sin_2=60°。
(2)离散正弦信号与计算所需结阵数:以信号1为例,离散成二十个点,即:j=1:1:N;ysin1=Amp_1*sin((j-1)*2*π/N+phase_sin_1/360*2*π);
计算所需结的个数:n1=round(ysin1*Kj_90/microwave_freq);即可计算20个点分别所需的结个数。
(3)判断是否发生捕获磁通现象:扫描结阵的I-V特性,观察超导量子台阶是否正常。
(4)采用平衡三进制计算:首先判断第九到第十一段结的偏置状态,再对前七段结采用平衡三进制计算。具体算法流程如图1所示。即包括如下步骤:
(41)判断偶数n1、n2的正负,分别存储在sign_1与sign_2中,所述sign_1与sign_2取值范围为{1-1},并将n1、n2取绝对值;
(42)以n2为例,n1同理,判断n2能否被4整除,若能进入(43),若不能进入(44);
(43)N2为n2整除4得到,通过循环判断N2的绝对值是否大于1093得到需要加入几段分段数为8744的结阵分段,所述分段数为8744的结阵分段有第9段、第10段以及第11段。这样就能够确定所述第9到11段结的偏置状态。
(44)N2通过N2=(n2-8742)/4得到,此时将第八段结的偏置状态设置为正偏,即B2(8)=1,通过循环判断所述第9到11段结的偏置状态。
(45)确定第9到第11段结偏置状态后,确定前7段结的偏置状态,采用平衡三进制计算,令A=N2,循环七次,Rem为A与3的余数,判断余数绝对值是否等于2。若等于2,则A=(A-Rem)/3+sign(Rem),B(i)=-sign(Rem);否则A=(A-Rem)/3,B(i)=Rem。所述N2是经过之前确定9至11段结运算后的N2,所述Rem为A与3的余数,所述sign(Rem)为Rem的正负性,取值是1或-1,所述B(i)为第i段结的偏置状态。
(5)采用索引算法计算:分别将电压和电流的所有情况的偏置矩阵穷举出来,电压结阵与电流矩阵分别有177147与19683种偏置矩阵。以电压结阵为例,将0到177146中所有的数从十进制转换成三进制并且将矩阵每个元素减1,能得到从[-1-1-1-1-1-1-1-1-1-1-1]到[1 1 1 1 1 1 1 1 1 1 1]所有的偏置矩阵,构成索引矩阵search_1,当输入所需要的结阵数后,通过在索引矩阵中寻找最佳偏置状态找到偏置矩阵。具体算法流程如图2所,即包括如下步骤:
(51)分别将两个结阵的所有情况遍历,分别有177147(3^11)和19683(3^9)种情况。以177147为例,将177146对应的十进制数转换成三进制的序列后,序列每一个元素减一存储至B1_all,将B1_all的每一行点乘JJS_1得到所有情况的偏置结总数B1M,将结果B1M与B1_all存储至索引矩阵search_1中;
所述B1_all为[-1-1-1-1-1-1-1-1-1-1-1]到[1 1 1 1 1 1 1 1 1 1 1]的所有向量组成的偏置矩阵情况,B1M为组合偏置结总数的所有情况矩阵,JJS_1为结段向量,JJS_1=[2916 108 36 972 12 4 324 8742 8744 8744 8744];
(52)分别输入所需的结阵数n1、n2,在索引矩阵中找到与所需结阵数差值最小的结阵数,并且找到对应的偏置矩阵。
(6)输出偏置状态矩阵:分别调用两个算法的matlab的m文件,输出偏置状态矩阵。离散二十个点的用平衡三进制和索引算法输出的偏置状态矩阵如表1、表2所示。所述表中数据-1、0、1分别表示第i段结阵在第j个离散点处的偏置状态。
表1平衡三进制输出偏置状态矩阵
表2索引算法输出偏置状态矩阵
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (7)
1.一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,包括:
步骤1,初始化偏置状态矩阵B1、B2及初始化两路正弦波信号,两路正弦波信号分别为信号1和信号2;设置台阶数N、微波频率与信号参数,
步骤2,采用N个台阶数对两路正弦波信号进行离散,将每个离散点的电压值转换成所需结阵数输入至向量n1(i)、n2(i);
步骤3,判断是否发生捕获磁通现象,若未发生捕获磁通现象则进入步骤4,否则,进入步骤5;
步骤4,采用平衡三进制算法计算量子电压序列对应的偏置状态矩阵;
步骤5,计算时关闭异常结段,采用索引算法计算偏置状态矩阵;
步骤6,输出偏置状态矩阵B1、B2。
2.根据权利要求1所述的一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,其中所述步骤1是这样实现的:
初始化偏置状态矩阵,将两个偏置状态矩阵B1、B2分别置零,矩阵维度分别为N×9与N×11,设置微波的频率以及信号1与信号2的幅值与相位;所述偏置状态矩阵每个元素的取值范围是{-1,0,1},在初始化时矩阵所有元素置零,所述N为离散点数。
3.根据权利要求1所述的一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,其中所述步骤2是这样实现的:
根据已知的正弦信号ysin1、ysin2分别从零点开始等间隔取N个点的电压值,并将这些电压值通过舍入和取整得到与之最接近的偶数值,这个偶数值即为所述的所需结阵数,对于两路信号第i个点的所需结阵数分别为n1(i)与n2(i)。
4.根据权利要求1所述的一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,其中所述步骤3是这样实现的:
检查瑟夫森结阵是否发生捕获磁通现象,若出现,则进入步骤4,若未发生,则进入步骤5;
所述捕获磁通现象包括:①某一段或者几段结的正负台阶或者0台阶宽度变窄;②某一段或者几段结无正负台阶或者0台阶;③输出的某一量子电压台阶具有一定的斜率,量子电压台阶不稳定。
5.根据权利要求4所述的一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,所述步骤4是这样实现的:
(41)判断偶数n1、n2的正负,分别存储在sign_1与sign_2中,所述sign_1与sign_2取值范围为{1-1},并将n1、n2取绝对值;
(42)以n2为例,n1同理,判断n2能否被4整除,若能进入(43),若不能进入(44);
(43)N2为n2整除4得到,通过循环判断N2的绝对值是否大于1093得到需要加入几段分段数为8744的结阵分段,所述分段数为8744的结阵分段有第9段、第10段以及第11段,这样就能够确定所述第9到11段结阵的偏置状态;
(44)N2通过N2=(n2-8742)/4得到,此时将第八段结的偏置状态设置为正偏,即B2(8)=1,通过循环判断所述第9到11段结的偏置状态;
(45)确定第9到第11段结偏置状态后,确定前7段结的偏置状态,采用平衡三进制计算,令A=N2,循环七次,Rem为A与3的余数,判断余数绝对值是否等于2;若等于2,则A=(A-Rem)/3+sign(Rem),B(i)=-sign(Rem);否则A=(A-Rem)/3,B(i)=Rem,所述N2是经过之前确定9至11段结运算后的N2,所述Rem为A与3的余数,所述sign(Rem)为Rem的正负性,取值是1或-1,所述B(i)为第i段结的偏置状态。
6.根据权利要求4所述的一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,所述步骤5是这样实现的:
(51)分别将两个结阵的所有情况遍历,分别有177147(3^11)和19683(3^9)种情况,以177147为例,将177146对应的十进制数转换成三进制的序列后,序列每一个元素减一存储至B1_all,将B1_all的每一行点乘JJS_1得到所有情况的偏置结总数B1M,将结果B1M与B1_all存储至搜索矩阵search_1中;
所述B1_all为[-1-1-1-1-1-1-1-1-1-1-1]到[1 1 1 1 1 1 1 1 1 1 1]的所有向量组成的偏置矩阵情况,B1M为组合偏置结总数的所有情况矩阵,JJS_1为结段向量,JJS_1=[2916 108 36 972 12 4 324 8742 8744 8744 8744];
(52)分别输入所需的结阵数n1、n2,在搜索矩阵中找到与所需结阵数差值最小的结阵数,并且找到对应的偏置矩阵。
7.根据权利要求4所述的一种平衡三进制与索引算法相结合的约瑟夫森结阵驱动方法,所述步骤6是这样实现的:
将步骤4或步骤5运行的结果输出,输出两个偏置矩阵B1、B2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910903099.9A CN110673523B (zh) | 2019-09-24 | 2019-09-24 | 一种平衡三进制与索引算法结合的约瑟夫森结阵驱动方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910903099.9A CN110673523B (zh) | 2019-09-24 | 2019-09-24 | 一种平衡三进制与索引算法结合的约瑟夫森结阵驱动方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110673523A CN110673523A (zh) | 2020-01-10 |
CN110673523B true CN110673523B (zh) | 2020-11-13 |
Family
ID=69077449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910903099.9A Active CN110673523B (zh) | 2019-09-24 | 2019-09-24 | 一种平衡三进制与索引算法结合的约瑟夫森结阵驱动方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110673523B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111628072B (zh) * | 2020-05-22 | 2020-11-27 | 中国计量科学研究院 | 约瑟夫森结阵偏置组合计算方法、电子设备及介质 |
CN112632894B (zh) * | 2020-12-30 | 2021-11-23 | 中国计量科学研究院 | 约瑟夫森结阵电压输出方法、装置、设备和存储介质 |
CN114153257B (zh) * | 2021-11-18 | 2023-04-11 | 国网陕西省电力公司营销服务中心(计量中心) | 一种双路可编程约瑟夫森结阵驱动方法及其装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095639A (zh) * | 2019-04-04 | 2019-08-06 | 中国计量科学研究院 | 一种产生交流量子电压的装置及方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3834635B2 (ja) * | 2002-11-18 | 2006-10-18 | 独立行政法人産業技術総合研究所 | 電圧標準装置 |
JP5414031B2 (ja) * | 2009-02-26 | 2014-02-12 | 独立行政法人産業技術総合研究所 | ジョセフソン素子を用いた交流電圧発生回路及び交流電圧標準回路 |
CN203313163U (zh) * | 2013-05-28 | 2013-11-27 | 中国计量科学研究院 | 基于可编程约瑟夫森结阵的动态信号合成装置 |
CN106771556B (zh) * | 2016-12-23 | 2019-02-22 | 中国计量科学研究院 | 一种基于量子技术的交流功率差分测量系统及方法 |
-
2019
- 2019-09-24 CN CN201910903099.9A patent/CN110673523B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110095639A (zh) * | 2019-04-04 | 2019-08-06 | 中国计量科学研究院 | 一种产生交流量子电压的装置及方法 |
Non-Patent Citations (1)
Title |
---|
基于可编程SNS 型约瑟夫森结阵的驱动系统研究;刘志尧;《电测与仪表》;20180425;第55卷(第8期);第125-130页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110673523A (zh) | 2020-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110673523B (zh) | 一种平衡三进制与索引算法结合的约瑟夫森结阵驱动方法 | |
Bertsekas et al. | Relaxation methods for minimum cost ordinary and generalized network flow problems | |
Zhou et al. | A multi-objective tabu search algorithm based on decomposition for multi-objective unconstrained binary quadratic programming problem | |
CN106787695A (zh) | 一种动态响应优化的开关电源控制方法 | |
Soeken et al. | Optimizing the mapping of reversible circuits to four-valued quantum gate circuits | |
CN108228536A (zh) | 使用FPGA实现Hermitian矩阵分解的方法 | |
Muthavarapu et al. | An efficient sorting algorithm for capacitor voltage balance of modular multilevel converter with space vector pulsewidth modulation | |
Maji et al. | FMACA: A fuzzy cellular automata based pattern classifier | |
Gal et al. | Improvements to the bacterial memetic algorithm used for fuzzy rule base extraction | |
CN105337522A (zh) | 一种模块化多电平换流器的双载波调制方法 | |
US20210294304A1 (en) | Control of matrix converters using machine learning | |
CN106230293B (zh) | 集中式多电平调制策略的状态机型脉冲分配方法 | |
Elsayed et al. | A three-strategy based differential evolution algorithm for constrained optimization | |
Weintraub | The shortest and the K‐shortest routes as assignment problems | |
Zholbaryssov et al. | On the use of entropy-like functions for solving the network flow problem with application to electric power systems | |
CN106849955A (zh) | 一种脉冲密度调制方法及脉冲密度值信号转换电路 | |
Liu | An exact algorithm for the biobjective 0-1 linear knapsack problem with a single continuous variable | |
Lee et al. | A distributed proportional load sharing controller for dc microgrids | |
Mizobe et al. | Strawberry optimization for multi-objective knapsack problem | |
Mahmoud | New Approach for Controlling PV-PV Differential Power Processing Converters | |
CN106295150A (zh) | 一种基因调控网络的优化控制方法 | |
US11656988B2 (en) | Memory device and operation method thereof | |
CN117134639B (zh) | 换流桥臂电路控制方法、装置、电子设备和计算机设备 | |
Antczak | EXACTNESS OF THE ABSOLUTE VALUE EXACT PENALTY FUNCTION METHOD FOR OPTIMIZATION PROBLEMS VIA UPPER CONVEXIFICATORS | |
Sankur et al. | Extremum seeking over a discrete action space |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |