1.本发明实施例涉及滤波技术,尤其涉及一种狭窄带通滤波方法、装置、设备及存储介质。
背景技术:
2.在电力系统发生故障后,继电保护主要是通过对故障后电流和电压信号的基频分量进行提取分析从而能够快速定位发生故障的具体位置。
3.常用的基频分量提取方法是采用全波傅里叶变换算法和半波傅里叶变换算法。全波傅里叶变换算法的数据窗需要一个周期,导致了响应速度比减慢,无法满足有些继电保护的速动性要求。而半波傅里叶变换算法只需要半个周期的数据窗,虽然响应速度快,但是无法有效滤除直流分量和偶次谐波分量,同时对偶次谐波具有放大影响,对低频分量抑制效果较差,容易导致继电保护中故障的误判断,造成故障扩大化。
技术实现要素:
4.本发明实施例提供了一种狭窄带通滤波方法、装置、设备及存储介质,能够有效滤除直流分量和偶次谐波分量,满足继电保护的速动性要求,从而避免故障的误判断。
5.第一方面,本发明实施例提供了一种狭窄带通滤波方法,该方法包括:
6.在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,其中,所述两个极点的值是共轭关系,所述传递函数是通过零极点配置法得到的;
7.将所述新的传递函数通过z反变换得到差分方程;
8.通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值,所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
9.第二方面,本发明实施例提供了一种狭窄带通滤波装置,该装置包括:
10.传递函数确定模块,用于在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,其中,所述两个极点的值是共轭关系,所述传递函数是通过零极点配置法得到的;
11.差分方程确定模块,用于将所述新的传递函数通过z反变换得到差分方程;
12.提取值确定模块,用于通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值,所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
13.第三方面,本发明实施例提供了一种计算机设备,该计算机设备包括:
14.一个或多个处理器;
15.存储装置,用于存储一个或多个程序;
16.当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现本发明任意实施例所述的狭窄带通滤波方法。
17.第四方面,本发明实施例提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现本发明任意实施例所述的狭窄带通滤波方法。
18.本发明实施例提供了一种狭窄带通滤波方法、装置、设备及存储介质,首先在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,然后将新的传递函数通过z反变换得到差分方程,最后通过差分方程对目标信号进行滤波,得到目标信号的基频分量提取值,能够有效滤除目标信号的直流分量和偶次谐波分量,满足继电保护的速动性要求,从而避免故障的误判断。
附图说明
19.图1为本发明实施例一提供的一种狭窄带通滤波方法的流程图;
20.图2为本发明实施例二提供的一种狭窄带通滤波方法的流程图;
21.图3为本发明实施例三提供的一种狭窄带通滤波装置的结构示意图;
22.图4为本发明实施例四提供的一种计算机设备的结构示意图。
具体实施方式
23.下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
24.实施例一
25.图1为本发明实施例一提供的一种狭窄带通滤波方法的流程图,本实施例可适用于对继电保护过程中电力系统故障后的电流信号或者电压信号进行滤波的情况。本实施例提供的狭窄带通滤波方法可以由本发明实施例提供的狭窄带通滤波装置来执行,该装置可以通过软件和/或硬件的方式实现,并集成在执行本方法的计算机设备中。
26.参见图1,本实施例的方法包括但不限于如下步骤:
27.s110,在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数。
28.其中,所述两个极点的值是共轭关系,所述传递函数是通过零极点配置法得到的。
29.考虑到要将某次谐波分量完全滤除,需要将该谐波分量设置为零点,要使某次谐波分量输出最大,需要将该谐波分量设置为极点,由此在相位角为0时,通过零点公式可以确定出传递函数中第一个零点的值,在相位角为π时通过零点公式可以确定出传递函数中第二个零点的值,然后通过极点公式确定出传递函数中两个极点的值,得到新的传递函数。
30.具体的,传递函数的推导过程如下所示:
31.(1)、一个n阶线性时不变系统的常系数线性差分方程可用下式表示:
[0032][0033]
其中,x表示自变量,y表示因变量,n表示大于或者等于1的整数,a0,a1…
,a
n
为常数,b1,b2,
…
,b
m
为常数。
[0034]
(2)、对上述差分方程进行z变换可以得到下式:
[0035]
y(z)=a0x(z) a1x(z)z-1
…
a
n
x(z)z-n
b1y(z)z-1
[0036]
b2y(z)z-2
…
b
m
y(z)z-m
[0037]
其中,z表示z变换的算子。
[0038]
(3)、对上式进行合并同类项可以得到下式:
[0039]
y(z)(1-b1z-1-b2z-2
-…-
b
m
z-m
)=x(z)(a0 a1z-1
a2z-2
…
a
n
z-n
)
[0040]
(4)获取传递函数:
[0041][0042]
其中,c
i
表示传递函数的零点,d
k
表示传递函数的极点,a表示增益。
[0043]
可选的,所述零点公式和极点公式为:
[0044]
c
i
=e
jωts
[0045]
d
k
=re
jωts
[0046]
其中,c
i
表示零点,i表示零点的个数,i=1,2,d
k
表示极点,k表示极点的个数,k=1,2,ω表示角频率,t表示采样周期,s表示拉普拉斯变换算子,r为大于零且小于1的常数;
[0047]
相应的,所述第一个零点的值为1,所述第二个零点的值为-1,第一个极点的值为re
jωts
,第二个极点的值为re-jωts
。
[0048]
具体的,当相位角为0(即ωt=0)时,通过零点公式可以计算得到c1=1;当相位角为π(即ωt=π)时,通过零点公式可以计算得到c2=-1;通过极点点公式可以计算得到d1=re
jωts
,由于第二个极点和第一个极点是共轭关系,由此可以得到d2=re-jωts
。
[0049]
新的传递函数为:
[0050][0051]
s120,将所述新的传递函数通过z反变换得到差分方程。
[0052]
在得到了新的传递函数之后,将新的传递函数通过z反变换可以得到差分方程。
[0053]
可选的,所述差分方程的表达式为:
[0054]
y(n 2)=ax(n 2)-ax(n) 2rcos(ωts)y(n 1)-r2y(n)
[0055]
其中,x表示自变量,y表示因变量,n表示大于或者等于1的整数。
[0056]
具体的,将进行z反变换就可以得到上面的差分方程表达式。
[0057]
s130,通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值。
[0058]
所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
[0059]
通过s120中的差分方程对继电保护过程中电力系统故障后的电流信号或者电压信号进行滤波,能够得到目标信号对应的基频分量提取值。
[0060]
进一步的,所述r的取值范围为:0.99<r<1。
[0061]
具体的,当r=0.9987时滤波效果较佳,幅值响应较大,能够较大的保留基频分量,且低频和高频幅值响应较小,能够对低频和高频分量进行良好的抑制。
[0062]
本实施例提供的技术方案,首先在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,然后将新的传递函数通过z反变换得到差分方程,最后通过差分方程对目标信号进行滤波,得到目标信号的基频分量提取值,能够有效滤除目标信号的直流分量和偶次谐波分量,满足继电保护的速动性要求,从而避免故障的误判断。
[0063]
实施例二
[0064]
图2为本发明实施例二提供的一种狭窄带通滤波方法的流程图。本发明实施例是在上述实施例的基础上进行优化。可选的,本实施例对得到目标信号的基频分量提取值的过程进行详细的解释说明。
[0065]
参见图2,本实施例的方法包括但不限于如下步骤:
[0066]
s210,在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数。
[0067]
s220,将所述新的传递函数通过z反变换得到差分方程。
[0068]
s230,对目标信号进行采样,得到预设个数的采样点以及与采样点个数相同的所述目标信号的采样值。
[0069]
其中,所述预设个数为偶数。
[0070]
假设采样周期t=20ms,对目标信号进行采样,可以得到预设个数(设为n)的采样点以及与采样点个数相同的目标信号的采样值,即x(1),x(2),
…
,x(n)。
[0071]
s240,将所述差分方程和半波傅里叶算法结合,对所述采样值进行处理,得到所述基频分量提取值。
[0072]
由于初始值会影响对采样值进行处理的速度,为了能使速度加快,此时可以将差分方程和半波傅里叶算法结合,对采样值进行处理,从而得到基频分量提取值。
[0073]
可选的,所述将所述差分方程和半波傅里叶算法结合,对所述采样值进行处理,得到所述基频分量提取值,可以具体包括:针对第一采样点进行半波傅里叶计算,得到第一初始值;针对第二采样点进行半波傅里叶计算,得到第二初始值,所述第一采样点和第二采样点位于所有采样点的中间,且所述第二采样点位于所述第一采样点之后;将第一初始值、第
二初始值以及所述目标信号的采样值代入所述差分方程进行计算,得到所述预设个数的点数据;对所述预设个数的点数据进行半波傅里叶计算,得到所述基频分量提取值。
[0074]
具体的,由于采样点个数为n,首先针对第一采样点(即第点)进行半波傅里叶计算,可以得到第一初始值,记为y1;接着针对第二采样点(即第 1点)进行半波傅里叶计算,可以得到第二初始值,记为y2;然后将第一初始值、第二初始值以及目标信号的采样值(x(1),x(2),
…
,x(n))代入差分方程的公式中进行计算,得到y3,y4,
…
,y
n
;最后对这n个数据进行半波傅里叶计算,就能够得到基频分量提取值。
[0075]
本发明实施例通过将差分方程和半波傅里叶算法结合,得到基频分量提取值,不仅收敛速度快,而且能够有效滤除目标信号的直流分量和偶次谐波分量,满足继电保护的速动性要求。
[0076]
进一步的,所述针对第一采样点进行半波傅里叶计算,得到第一初始值,可以具体包括:在确定采样时间大于采样周期的一半时,针对第一采样点进行半波傅里叶计算,得到第一初始值。
[0077]
具体的,由于半波傅里叶算法需要满足采样时间大于采样周期的一半,因此在确定采样时间大于采样周期的一半时,针对第一采样点进行半波傅里叶计算,得到第一初始值。
[0078]
本发明实施例通过在确定采样时间大于采样周期的一半时,计算第一初始值,使得得到的第一初始值更为准确,从而加快基频分量的提取过程。
[0079]
本实施例提供的技术方案,首先在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,接着将新的传递函数通过z反变换得到差分方程,然后对目标信号进行采样,得到预设个数的采样点以及与采样点个数相同的目标信号的采样值,最后将差分方程和半波傅里叶算法结合,对采样值进行处理,得到基频分量提取值,不仅能够加快基频分量的提取过程,而且能够有效滤除目标信号的直流分量和偶次谐波分量,满足继电保护的速动性要求,从而避免故障的误判断。
[0080]
实施例三
[0081]
图3为本发明实施例三提供的一种狭窄带通滤波装置的结构示意图,如图3所示,该装置可以包括:
[0082]
传递函数确定模块310,用于在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,其中,所述两个极点的值是共轭关系,所述传递函数是通过零极点配置法得到的;
[0083]
差分方程确定模块320,用于将所述新的传递函数通过z反变换得到差分方程;
[0084]
提取值确定模块330,用于通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值,所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
[0085]
本实施例提供的技术方案,首先在相位角为0时,通过零点公式确定出传递函数中
第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,然后将新的传递函数通过z反变换得到差分方程,最后通过差分方程对目标信号进行滤波,得到目标信号的基频分量提取值,能够有效滤除目标信号的直流分量和偶次谐波分量,满足继电保护的速动性要求,从而避免故障的误判断。
[0086]
进一步的,所述零点公式和极点公式为:
[0087]
c
i
=e
jωts
[0088]
d
k
=re
jωts
[0089]
其中,c
i
表示零点,i表示零点的个数,i=1,2,d
k
表示极点,k表示极点的个数,k=1,2,ω表示角频率,t表示采样周期,s表示拉普拉斯变换算子,r为大于零且小于1的常数;
[0090]
相应的,所述第一个零点的值为1,所述第二个零点的值为-1,第一个极点的值为re
jωts
,第二个极点的值为re-jωts
。
[0091]
进一步的,上述狭窄带通滤波装置,还可以包括:
[0092]
目标信号采样模块,用于对目标信号进行采样,得到预设个数的采样点以及与采样点个数相同的所述目标信号的采样值,所述预设个数为偶数;
[0093]
分量提取值确定模块,用于将所述差分方程和半波傅里叶算法结合,对所述采样值进行处理,得到所述基频分量提取值。
[0094]
进一步的,上述分量提取值确定模块,可以具体用于:
[0095]
针对第一采样点进行半波傅里叶计算,得到第一初始值;
[0096]
针对第二采样点进行半波傅里叶计算,得到第二初始值,所述第一采样点和第二采样点位于所有采样点的中间,且所述第二采样点位于所述第一采样点之后;
[0097]
将第一初始值、第二初始值以及所述目标信号的采样值代入所述差分方程进行计算,得到所述预设个数的点数据;
[0098]
对所述预设个数的点数据进行半波傅里叶计算,得到所述基频分量提取值。
[0099]
进一步的,上述狭窄带通滤波装置,还可以包括:
[0100]
第一初始值确定模块,用于在确定采样时间大于采样周期的一半时,针对第一采样点进行半波傅里叶计算,得到第一初始值。
[0101]
进一步的,所述差分方程的表达式为:
[0102]
y(n 2)=ax(n 2)-ax(n) 2rcos(ωts)y(n 1)-r2y(n)
[0103]
其中,x表示自变量,y表示因变量,n表示大于或者等于1的整数。
[0104]
进一步的,所述r的取值范围为:0.99<r<1。
[0105]
本实施例提供的狭窄带通滤波装置可适用于上述任意实施例提供的狭窄带通滤波方法,具备相应的功能和有益效果。
[0106]
实施例四
[0107]
图4为本发明实施例四提供的一种计算机设备的结构示意图,如图4所示,该计算机设备包括处理器410、存储装置420和通信装置430;计算机设备中处理器410的数量可以是一个或多个,图4中以一个处理器410为例;计算机设备中的处理器410、存储装置420和通信装置430可以通过总线或其他方式连接,图4中以通过总线连接为例。
[0108]
存储装置420作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行
程序以及模块,如本发明实施例中的狭窄带通滤波方法对应的模块(例如,用于狭窄带通滤波装置装置中的传递函数确定模块310、差分方程确定模块320和提取值确定模块330)。处理器410通过运行存储在存储装置420中的软件程序、指令以及模块,从而执行计算机设备的各种功能应用以及数据处理,即实现上述的狭窄带通滤波方法。
[0109]
存储装置420可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据终端的使用所创建的数据等。此外,存储装置420可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储装置420可进一步包括相对于处理器410远程设置的存储器,这些远程存储器可以通过网络连接至计算机设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
[0110]
通信装置430,用于实现服务器之间的网络连接或者移动数据连接。
[0111]
本实施例提供的一种计算机设备可用于执行上述任意实施例提供的狭窄带通滤波方法,具备相应的功能和有益效果。
[0112]
实施例五
[0113]
本发明实施例五还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现本发明任意实施例中的狭窄带通滤波方法,该方法具体包括:
[0114]
在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,其中,所述两个极点的值是共轭关系,所述传递函数是通过零极点配置法得到的;
[0115]
将所述新的传递函数通过z反变换得到差分方程;
[0116]
通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值,所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。
[0117]
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作,还可以执行本发明任意实施例所提供的狭窄带通滤波方法中的相关操作。
[0118]
通过以上关于实施方式的描述,所属领域的技术人员可以清楚地了解到,本发明可借助软件及必需的通用硬件来实现,当然也可以通过硬件实现,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如计算机的软盘、只读存储器(read-only memory,rom)、随机存取存储器(random access memory,ram)、闪存(flash)、硬盘或光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
[0119]
值得注意的是,上述狭窄带通滤波装置的实施例中,所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
[0120]
以上所述仅为本发明的优选实施例,并不用于限制本发明,对于本领域技术人员而言,本发明可以有各种改动和变化。凡在本发明的精神和原理之内所作的任何修改、等同
替换、改进等,均应包含在本发明的保护范围之内。
技术特征:
1.一种狭窄带通滤波方法,其特征在于,包括:在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,其中,所述两个极点的值是共轭关系,所述传递函数是通过零极点配置法得到的;将所述新的传递函数通过z反变换得到差分方程;通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值,所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。2.根据权利要求1所述的方法,其特征在于,所述零点公式和极点公式为:c
i
=e
jωts
d
k
=re
jωts
其中,c
i
表示零点,i表示零点的个数,i=1,2,d
k
表示极点,k表示极点的个数,k=1,2,ω表示角频率,t表示采样周期,s表示拉普拉斯变换算子,r为大于零且小于1的常数;相应的,所述第一个零点的值为1,所述第二个零点的值为-1,第一个极点的值为re
jωts
,第二个极点的值为re-jωts
。3.根据权利要求2所述的方法,其特征在于,所述通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值,包括:对目标信号进行采样,得到预设个数的采样点以及与采样点个数相同的所述目标信号的采样值,所述预设个数为偶数;将所述差分方程和半波傅里叶算法结合,对所述采样值进行处理,得到所述基频分量提取值。4.根据权利要求3所述的方法,其特征在于,所述将所述差分方程和半波傅里叶算法结合,对所述采样值进行处理,得到所述基频分量提取值,包括:针对第一采样点进行半波傅里叶计算,得到第一初始值;针对第二采样点进行半波傅里叶计算,得到第二初始值,所述第一采样点和第二采样点位于所有采样点的中间,且所述第二采样点位于所述第一采样点之后;将第一初始值、第二初始值以及所述目标信号的采样值代入所述差分方程进行计算,得到所述预设个数的点数据;对所述预设个数的点数据进行半波傅里叶计算,得到所述基频分量提取值。5.根据权利要求4所述的方法,其特征在于,所述针对第一采样点进行半波傅里叶计算,得到第一初始值,包括:在确定采样时间大于采样周期的一半时,针对第一采样点进行半波傅里叶计算,得到第一初始值。6.根据权利要求2所述的方法,其特征在于,所述差分方程的表达式为:y(n 2)=ax(n 2)-ax(n) 2r cos(ωts)y(n 1)-r2y(n)其中,x表示自变量,y表示因变量,n表示大于或者等于1的整数。7.根据权利要求6所述的方法,其特征在于,所述r的取值范围为:0.99<r<1。8.一种狭窄带通滤波装置,其特征在于,包括:传递函数确定模块,用于在相位角为0时,通过零点公式确定出传递函数中第一个零点
的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数,其中,所述两个极点的值是共轭关系,所述传递函数是通过零极点配置法得到的;差分方程确定模块,用于将所述新的传递函数通过z反变换得到差分方程;提取值确定模块,用于通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值,所述目标信号是继电保护过程中电力系统故障后的电流信号或者电压信号。9.一种计算机设备,其特征在于,所述计算机设备包括:一个或多个处理器;存储装置,用于存储一个或多个程序;当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-7中任一所述的狭窄带通滤波方法。10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-7中任一所述的狭窄带通滤波方法。
技术总结
本发明实施例公开了一种狭窄带通滤波方法、装置、设备及存储介质。其中,该方法包括:在相位角为0时,通过零点公式确定出传递函数中第一个零点的值,和在相位角为π时通过零点公式确定出传递函数中第二个零点的值,通过极点公式确定出传递函数中两个极点的值,得到新的传递函数;将所述新的传递函数通过Z反变换得到差分方程;通过所述差分方程对目标信号进行滤波,得到所述目标信号的基频分量提取值。本发明实施例提供的技术方案,能够有效滤除目标信号的直流分量和偶次谐波分量,满足继电保护的速动性要求,从而避免故障的误判断。从而避免故障的误判断。从而避免故障的误判断。
技术研发人员:叶飞
受保护的技术使用者:中船动力研究院有限公司
技术研发日:2020.11.17
技术公布日:2021/3/9
转载请注明原文地址:https://wp.8miu.com/read-54775.html