CN202339381U - 一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量*** - Google Patents
一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量*** Download PDFInfo
- Publication number
- CN202339381U CN202339381U CN201120405403.6U CN201120405403U CN202339381U CN 202339381 U CN202339381 U CN 202339381U CN 201120405403 U CN201120405403 U CN 201120405403U CN 202339381 U CN202339381 U CN 202339381U
- Authority
- CN
- China
- Prior art keywords
- nuttall
- unit
- self
- harmonic
- electric energy
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Complex Calculations (AREA)
Abstract
一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量***,该***的加纳托尔自卷积窗单元连接信号采样单元与FFT运算单元,第h次谐波的频率范围计算单元连接FFT运算单元与幅度最大和次大谱线的局部搜索单元,频谱插值计算单元连接幅度最大和次大谱线的局部搜索单元与频率偏移量计算单元,谐波参数计算单元连接频率偏移量计算单元与谐波电能输出单元。采用该***能够实现复杂电力信号中的整数次谐波电能计量。
Description
技术领域
本实用新型涉及一种谐波电能计量***,具体是一种纳托尔自卷积窗加权傅里叶变换的电力谐波电能计量***。属于信号处理领域。
背景技术
高精度信号分析与处理为谐波潮流计算、设备检测、电力***谐波补偿与抑制、振动信号分析和随机故障处理等提供科学依据。采用傅立叶变换理论进行信号处理时,由于信号频率存在波动或干扰,严格的同步采样无法实现,频谱泄漏和栅栏效应会引入较大误差。
经文献检索发现,采用窗函数对信号加权可以有效减少频谱泄漏,窗函数频谱的主瓣与频率分辨力直接有关,主瓣宽、频率分辨力低;旁瓣与泄漏直接有关,旁瓣大,泄漏多;旁瓣衰减斜率反映旁瓣衰减的速度,旁瓣衰减越快,对泄漏抑制越强。国内外学者提出了一系列窗函数,如Hanning窗、Blackman-Harris窗、Rife-Vincent(I)窗、纳托尔窗和矩形卷积窗等,并将它们运用到信号加权处理算法中,提高了信号处理精度。
基于经典窗函数(Hanning窗、Blackman-Harris窗、Rife-Vincent(I)窗、纳托尔窗和矩形卷积窗等)的信号加权分析算法利用窗函数的频谱旁瓣下降特性可以在一定程度上可以减少频谱泄漏。但是,由于经典窗的旁瓣特性仍不够理想,对频谱泄漏的抑制作用仍然有限,谐波间的相互干扰不可忽略,信号分析处理精度受到限制。
发明内容
本实用新型目的在于克服现有加窗傅里叶变换谐波电能计量***的缺陷,提供一种简便的纳托尔自卷积窗加权傅里叶变换谐波电能计量***,解决信号加权处理过程中频谱泄漏过大,谐波间相互干扰的问题,从而提高信号谐波分析的准确度和实用性,为进一步进行的信号参数识别与谐波电能计量提供可靠依据。
本实用新型通过以下技术方案实现上述目的:一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量***,由信号采样单元、加纳托尔自卷积窗单元、FFT运算单元、第h次谐波的频率范围计算单元、幅度最大和次大谱线的局部搜索单元、频谱插值计算单元、频率偏移量计算单元、谐波参数计算单元以及谐波电能输出单元组成,加纳托尔自卷积窗单元连接信号采样单元与FFT运算单元,第h次谐波的频率范围计算单元连接FFT运算单元与幅度最大和次大谱线的局部搜索单元,频谱插值计算单元连接幅度最大和次大谱线的局部搜索单元与频率偏移量计算单元,谐波参数计算单元连接频率偏移量计算单元与谐波电能输出单元。
所述加纳托尔自卷积窗单元是先构建长度为M的离散纳托尔窗序列,然后将p个离散纳托尔窗序列作p-1次卷积运算,得到卷积序列,在卷积序列的首或尾补p-1个零,即得到长度为N=pM的p阶离散纳托尔自卷积窗。
所述频谱插值计算单元是运用最小二乘法对谐波频率附近的幅度最大和次大的谱线进行拟合,获得该谐波频率的准确频谱。
该***实现步骤如下:首先构建长度为离散纳托尔窗序列,然后将离散纳托尔窗卷积运算,得到卷积序列,在卷积序列的首或尾补零,得到离散纳托尔 自卷积窗,采用纳托尔自卷积窗序列对离散电力信号进行加权运算,再采用傅里叶变换获得电力信号的离散频谱,在离散频谱对应的谐波频率附近寻找幅度最大和次大的谱线,运用最小二乘法对峰值谱线进行拟合,获得电力谐波的幅值、频率和初相角,从而实现复杂电力信号中的整数次谐波电能计量。
本实用新型的原理是:
纳托尔窗是余弦组合窗,其时域表达式为
式中,M为窗函数的项数;n=0,1,…,N-1;bg满足约束条件
表1给出了典型的纳托尔窗函数的系数
表1纳托尔窗函数的系数
时域离散的纳托尔自卷积窗为
其中,下标p表示参与卷积运算的纳托尔自卷积窗的个数,下标N表示参与卷 积的纳托尔窗,为便于快速傅里叶变换的实现,对卷积结果进行补零操作,则离散纳托尔自卷积窗的长度为N=pM。
设仅含单一频率成分的信号x(t),通过采样频率为fs的模数转换后,得到离散序列:
式中,A0、f0、φ0分别为信号的幅值、频率和初相角。
对离散化后的信号加长度为N的p阶纳托尔自卷积窗wN-p(n)(n=0,1,…,N-1),得到截短后的序列为
x(n)=x(m)wN-p(n),n=0,1,…,N-1
对序列x(n)进行离散傅里叶变换后,得到其离散频谱为
式中,k0=f0N/fs。
设离散频谱中,频率f0附近的局部幅值最大和次最大的谱线分别为第k1和k2根,满足k1≤k0≤k2=k1+1。在频率f0附近采用局部峰值搜索策略找到这两条谱线,即可确定k1和k2。设这两条谱线幅值分别是y1和y2,即
考虑到0≤k0-k1≤1,定义系数α为
α=k0-k1-0.5,-0.5≤α≤0.5
则y1和y2可改写为关于α的函数:
定义系数β为
将式上式中的y1和y2进行置换,则β可写为关于α的函数:
在α∈[-0.5,0.5]范围内,以等间隔取样若干个点,利用每个αj值计算出相应的βi,确定待多项式拟合次数K,利用式上求解系数ak(k=0,1,…,K-1),从而建立频谱插值多项式SK(x),可记为
α=g-1(β)≈a1β+a2β2+a3β3+…+aKβK
式中,a1,a2,...,aK为拟合多项式的系数。
频率f0的计算式为
幅值为
或
初相角为
对于各次谐波分量,均可按上述方法求解参数。
附图说明
图1为本实用新型所述的纳托尔自卷积窗加权傅里叶变换的谐波电能计量***的方框图。
图2是本实用新型的离散频谱插值处理示意图。
具体实施方式
以下通过实施例对本实用新型的技术方案进一步说明。
如图1所示,一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量***,由信号采样单元、加纳托尔自卷积窗单元、FFT运算单元、第h次谐波的频率范围计算单元、幅度最大和次大谱线的局部搜索单元、频谱插值计算单元、频率偏移量计算单元、谐波参数计算单元以及谐波电能输出单元组成,加纳托尔自卷积窗单元连接信号采样单元与FFT运算单元,第h次谐波的频率范围计算单元连接FFT运算单元与幅度最大和次大谱线的局部搜索单元,频谱插值计算单元连接幅度最大和次大谱线的局部搜索单元与频率偏移量计算单元,谐波参数计算单元连接频率偏移量计算单元与谐波电能输出单元。
本实用新型所述的纳托尔自卷积窗加权傅里叶变换的谐波电能计量***的实现方法按如下步骤进行:
首先建立M=64、G=3的离散4项最小旁瓣纳托尔窗,即选择b0=0.3635819、b1=0.4891775、b2=0.1365995、b3=0.0106411,通过下式计算4项最小旁瓣纳托尔窗的序列:
其次,对w(n)进行2~4阶自卷积运算,即p分别取值为2、3和4,再进行对应 的补零操作,使得2~4阶4项最小旁瓣纳托尔自卷积窗的长度分别为128、192和256。
设有一多频率成分的电压信号为
式中,A0=0.2V;A1=6V;f1=20.2Hz;φ1=0.1rad;A3=1;f3=60.6Hz;φ3=0rad。
采用所建立的4阶4项最小旁瓣纳托尔自卷积窗对离散序列x(n)进行加权,再进行离散傅里叶变换,得到电压信号加权后的离散频谱。对离散频谱进行归一化处理,频率f1附近的频谱幅度分布如图2所示。信号频率f1不与离散谱线重合,且位于离散频谱幅值最大谱线k1=4和次大谱线k2=5之间。确定k1和k2分别为最大和次大峰值谱线后,得到对应的幅度值y1和y2,根据前述发明内容中的最小二乘峰值频谱拟合过程,计算得到多项式最高次数设置为K=5时,加4阶4项最小旁瓣纳托尔自卷积窗后,离散频谱插值多项式为
α4th=0.4026β5+0.7644β3+11.9048β
将 带入上式,即计算出f1为20.1999999996Hz,A0为0.20000000007V,A1为5.9999999999V,φ1为0.10000008rad。依次类推,可计算出A3为0.99999999997V,f3为60.60000000001Hz,φ3为-0.00000002rad。
如图2所示,离散频谱插值处理原理为:
设在离散频谱中,频率f0附近的局部幅值最大和次最大的谱线分别为第k1和k2根,设这两条谱线幅值分别是y1和y2。考虑到0≤k0-k1≤1,定义系数α为
α=k0-k1-0.5,-0.5≤α≤0.5
则y1和y2可改写为关于α的函数,定义系数β为
将上式中的y1和y2进行置换,则β可写为关于α的函数:
在α∈[-0.5,0.5]范围内,以等间隔取样若干个点,利用每个αi值计算出相应的βi,确定待多项式拟合次数K,利用式上求解系数ak(k=0,1,…,K-1),从而建立频谱插值多项式SK(x),可记为
α=g-1(β)≈a1β+a2β2+a3β3+…+aKβK
式中,a1,a2,...,aK为拟合多项式的系数。从而实现离散频谱插值处理。
Claims (3)
1.一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量***,其特征在于,该***由信号采样单元、加纳托尔自卷积窗单元、FFT运算单元、第h次谐波的频率范围计算单元、幅度最大和次大谱线的局部搜索单元、频谱插值计算单元、频率偏移量计算单元、谐波参数计算单元以及谐波电能输出单元组成,加纳托尔自卷积窗单元连接信号采样单元与FFT运算单元,第h次谐波的频率范围计算单元连接FFT运算单元与幅度最大和次大谱线的局部搜索单元,频谱插值计算单元连接幅度最大和次大谱线的局部搜索单元与频率偏移量计算单元,谐波参数计算单元连接频率偏移量计算单元与谐波电能输出单元。
2.根据权利要求1所述的纳托尔自卷积窗加权傅里叶变换的谐波电能计量***,其特征在于,所述加纳托尔自卷积窗单元是先构建长度为M的离散纳托尔窗序列,然后将p个离散纳托尔窗序列作p-1次卷积运算,得到卷积序列,在卷积序列的首或尾补p-1个零,即得到长度为N=pM的p阶离散纳托尔自卷积窗。
3.根据权利要求1所述的纳托尔自卷积窗加权傅里叶变换的谐波电能计量***,其特征在于,所述频谱插值计算单元是运用最小二乘法对谐波频率附近的幅度最大和次大的谱线进行拟合,获得该谐波频率的准确频谱。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103237195A CN102435845A (zh) | 2011-10-19 | 2011-10-19 | 一种Nuttall自卷积窗加权傅里叶变换的谐波电能计量方法 |
CN201120405403.6U CN202339381U (zh) | 2011-10-19 | 2011-10-19 | 一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量*** |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011103237195A CN102435845A (zh) | 2011-10-19 | 2011-10-19 | 一种Nuttall自卷积窗加权傅里叶变换的谐波电能计量方法 |
CN201120405403.6U CN202339381U (zh) | 2011-10-19 | 2011-10-19 | 一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN202339381U true CN202339381U (zh) | 2012-07-18 |
Family
ID=51302856
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201120405403.6U Expired - Fee Related CN202339381U (zh) | 2011-10-19 | 2011-10-19 | 一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量*** |
CN2011103237195A Pending CN102435845A (zh) | 2011-10-19 | 2011-10-19 | 一种Nuttall自卷积窗加权傅里叶变换的谐波电能计量方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011103237195A Pending CN102435845A (zh) | 2011-10-19 | 2011-10-19 | 一种Nuttall自卷积窗加权傅里叶变换的谐波电能计量方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN202339381U (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103207319A (zh) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | 数字化变电站电力信号非同步采样条件下的谐波测量方法 |
CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
CN104062528A (zh) * | 2014-07-04 | 2014-09-24 | 武汉大学 | 基于汉宁乘积窗的信号谐波分析方法及*** |
CN105137185A (zh) * | 2015-07-23 | 2015-12-09 | 河海大学 | 一种基于离散傅里叶变换的频域插值电力谐波分析方法 |
CN105137181A (zh) * | 2015-09-15 | 2015-12-09 | 福州大学 | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 |
CN106018956A (zh) * | 2016-08-10 | 2016-10-12 | 北京妙微科技有限公司 | 一种加窗谱线插值的电力***频率计算方法 |
CN106771586A (zh) * | 2016-11-24 | 2017-05-31 | 云南电网有限责任公司电力科学研究院 | 一种直流控制保护板卡的回路信号分析方法及装置 |
CN109782063A (zh) * | 2018-10-23 | 2019-05-21 | 国网安徽省电力有限公司芜湖供电公司 | 一种基于纳托尔自卷积窗三谱线插值fft的动态间谐波分析方法 |
CN110763930A (zh) * | 2019-09-28 | 2020-02-07 | 沈阳工程学院 | 基于布莱克曼双峰插值谐波分析法的避雷器阻性电流在线监测*** |
CN117076833A (zh) * | 2023-10-16 | 2023-11-17 | 深圳大学 | 一种基于fpga的单比特降维快速傅里叶变换方法 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102998528B (zh) * | 2012-12-12 | 2014-12-10 | 湖南大学 | 动态正弦畸变信号中频率分量的提取方法及其实现装置 |
CN103105529A (zh) * | 2013-01-22 | 2013-05-15 | 广西电网公司电力科学研究院 | 一种基于参数分析的谐波电能计量***及其控制方法 |
CN103105532A (zh) * | 2013-01-22 | 2013-05-15 | 广西电网公司电力科学研究院 | 一种改进傅里叶变换的谐波电能计量***及其控制方法 |
CN104111373B (zh) * | 2014-07-23 | 2016-11-09 | 国家电网公司 | 一种用于智能变电站数字化电能的计量方法 |
CN105372492B (zh) * | 2014-08-31 | 2018-06-22 | 常州昊云工控科技有限公司 | 基于三条dft复数谱线的信号频率测量方法 |
CN104749433A (zh) * | 2015-01-04 | 2015-07-01 | 云南电网公司电力科学研究院 | 一种准确计算分析谐波及电能质量的方法 |
CN109633262A (zh) * | 2019-01-29 | 2019-04-16 | 国网湖南省电力有限公司 | 基于组合窗多谱线fft的三相谐波电能计量方法、装置 |
CN109900959B (zh) * | 2019-04-17 | 2020-10-02 | 贵州电网有限责任公司 | 一种动态正弦畸变信号中谐波成分的提取方法 |
CN110133564B (zh) * | 2019-06-02 | 2021-04-09 | 国网江西省电力有限公司电力科学研究院 | 一种电流传感器相位特性的同步测试方法及装置 |
CN110221125B (zh) * | 2019-07-15 | 2020-06-12 | 福州大学 | 一种采用直流分量恢复的直流***谐波分析方法 |
CN114258495B (zh) * | 2020-07-23 | 2023-08-15 | 刘保国 | 一种有限复数信号测量***与高精度分解方法 |
CN111929532A (zh) * | 2020-08-06 | 2020-11-13 | 国网北京市电力公司 | 一种10kV配电电缆中间接头定位方法 |
CN111999635A (zh) * | 2020-08-23 | 2020-11-27 | 国网山东省电力公司莱芜供电公司 | 一种基于4项5阶Nuttall窗的板卡故障信号分析方法及终端 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100454026C (zh) * | 2005-07-28 | 2009-01-21 | 泰州苏源集团科电有限公司 | 电网电能质量中的周波采样方法 |
CN100516897C (zh) * | 2007-08-03 | 2009-07-22 | 杭州万工科技有限公司 | 一种谐波电能的计量方法 |
JP2009074941A (ja) * | 2007-09-20 | 2009-04-09 | Fujitsu Ltd | 周期信号の振幅測定方法および装置並びに磁気ヘッドの試験方法および装置 |
CN101261618A (zh) * | 2008-04-14 | 2008-09-10 | 湖南大学 | 三角自卷积窗加权算法 |
CN101852826B (zh) * | 2009-03-30 | 2012-09-05 | 西门子公司 | 一种电力***的谐波分析方法及其装置 |
-
2011
- 2011-10-19 CN CN201120405403.6U patent/CN202339381U/zh not_active Expired - Fee Related
- 2011-10-19 CN CN2011103237195A patent/CN102435845A/zh active Pending
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103207319A (zh) * | 2013-03-12 | 2013-07-17 | 广东电网公司电力科学研究院 | 数字化变电站电力信号非同步采样条件下的谐波测量方法 |
CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
CN104062528A (zh) * | 2014-07-04 | 2014-09-24 | 武汉大学 | 基于汉宁乘积窗的信号谐波分析方法及*** |
CN105137185A (zh) * | 2015-07-23 | 2015-12-09 | 河海大学 | 一种基于离散傅里叶变换的频域插值电力谐波分析方法 |
CN105137185B (zh) * | 2015-07-23 | 2017-09-22 | 河海大学 | 一种基于离散傅里叶变换的频域插值电力谐波分析方法 |
CN105137181A (zh) * | 2015-09-15 | 2015-12-09 | 福州大学 | 基于Nuttall-Kaiser组合窗双谱线插值的谐波分析方法 |
CN106018956B (zh) * | 2016-08-10 | 2018-10-16 | 北京妙微科技有限公司 | 一种加窗谱线插值的电力***频率计算方法 |
CN106018956A (zh) * | 2016-08-10 | 2016-10-12 | 北京妙微科技有限公司 | 一种加窗谱线插值的电力***频率计算方法 |
CN106771586A (zh) * | 2016-11-24 | 2017-05-31 | 云南电网有限责任公司电力科学研究院 | 一种直流控制保护板卡的回路信号分析方法及装置 |
CN106771586B (zh) * | 2016-11-24 | 2019-07-26 | 云南电网有限责任公司电力科学研究院 | 一种直流控制保护板卡的回路信号分析方法及装置 |
CN109782063A (zh) * | 2018-10-23 | 2019-05-21 | 国网安徽省电力有限公司芜湖供电公司 | 一种基于纳托尔自卷积窗三谱线插值fft的动态间谐波分析方法 |
CN110763930A (zh) * | 2019-09-28 | 2020-02-07 | 沈阳工程学院 | 基于布莱克曼双峰插值谐波分析法的避雷器阻性电流在线监测*** |
CN117076833A (zh) * | 2023-10-16 | 2023-11-17 | 深圳大学 | 一种基于fpga的单比特降维快速傅里叶变换方法 |
CN117076833B (zh) * | 2023-10-16 | 2023-12-12 | 深圳大学 | 一种基于fpga的单比特降维快速傅里叶变换方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102435845A (zh) | 2012-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN202339381U (zh) | 一种纳托尔自卷积窗加权傅里叶变换的谐波电能计量*** | |
CN103454497B (zh) | 基于改进加窗离散傅立叶变换的相位差测量方法 | |
CN101603985B (zh) | 高准确度正弦信号测量方法 | |
CN106483374B (zh) | 一种基于Nuttall双窗全相位FFT的谐波间谐波检测方法 | |
CN102435844B (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN101701984B (zh) | 基于三项系数Nuttall窗插值FFT的基波与谐波检测方法 | |
CN104897960B (zh) | 基于加窗四谱线插值fft的谐波快速分析方法及*** | |
CN103576002B (zh) | 一种容性绝缘设备介质损耗角的计算方法 | |
CN101261292A (zh) | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
CN104062528A (zh) | 基于汉宁乘积窗的信号谐波分析方法及*** | |
CN102890190B (zh) | 基于多小波的随机非线性负荷有功电能计量方法 | |
CN101441233A (zh) | 基于Kaiser窗双谱线插值FFT的基波与谐波检测方法 | |
CN103207319A (zh) | 数字化变电站电力信号非同步采样条件下的谐波测量方法 | |
CN103353550A (zh) | 一种测量电力***信号频率及谐波参数的方法 | |
CN104062500A (zh) | 基于汉明乘积窗的信号谐波分析方法及*** | |
CN104391178A (zh) | 一种基于Nuttall窗的时移相位差稳态谐波信号校正方法 | |
CN102662106A (zh) | 谐波电网电能计量方法 | |
CN103675447B (zh) | 一种电气化铁路的高精度实时谐波分析方法 | |
CN104914408A (zh) | 基于中国余数定理的频率、doa联合测量方法以及装置 | |
CN110542828A (zh) | 一种配网故障定位方法及*** | |
CN104714075A (zh) | 一种电网电压闪变包络参数提取方法 | |
CN108776263B (zh) | 基于高阶汉宁自卷积窗及改进插值算法的谐波检测方法 | |
CN102135552A (zh) | 电网有功功率和无功功率的数字化实时检测方法 | |
CN103795411A (zh) | 基于五项最大旁瓣衰减窗三谱线插值测试sfdr的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120718 Termination date: 20151019 |
|
EXPY | Termination of patent right or utility model |