CN115165270B - Method and device for processing multi-channel flutter test signals - Google Patents

Method and device for processing multi-channel flutter test signals Download PDF

Info

Publication number
CN115165270B
CN115165270B CN202210674150.5A CN202210674150A CN115165270B CN 115165270 B CN115165270 B CN 115165270B CN 202210674150 A CN202210674150 A CN 202210674150A CN 115165270 B CN115165270 B CN 115165270B
Authority
CN
China
Prior art keywords
density function
power spectrum
singular value
spectrum density
function matrix
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
CN202210674150.5A
Other languages
Chinese (zh)
Other versions
CN115165270A (en
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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202210674150.5A priority Critical patent/CN115165270B/en
Publication of CN115165270A publication Critical patent/CN115165270A/en
Application granted granted Critical
Publication of CN115165270B publication Critical patent/CN115165270B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00Vibration-testing of structures; Shock-testing of structures
    • G01M7/02Vibration-testing by means of a shake table
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60Testing or inspecting aircraft components or systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明实施例公开了一种信号处理的方法和装置。该信号处理的方法包括:依据获取到的多通道紊流响应信号计算功率谱密度函数矩阵;依据功率谱密度函数矩阵获取最大奇异值曲线;依据最大奇异值曲线获取单自由度系统的自功率谱密度函数;依据自功率谱密度函数得到单模态系统模态参数。通过本发明,解决相关技术中由于现有技术在颤振试飞试验中紊流响应信号分析问题,达到了提高模态参数参数估计的精度,同时避免由于随机子空间方法带来的计算实时性不足的技术效果。

The embodiment of the present invention discloses a method and device for signal processing. The signal processing method includes: calculating a power spectrum density function matrix based on the acquired multi-channel turbulence response signal; obtaining a maximum singular value curve based on the power spectrum density function matrix; obtaining the self-power spectrum density function of the single-degree-of-freedom system based on the maximum singular value curve; and obtaining the modal parameters of the single-mode system based on the self-power spectrum density function. Through the present invention, the problem of turbulence response signal analysis in flutter flight tests in the related technology due to the existing technology is solved, and the accuracy of modal parameter estimation is improved, while avoiding the technical effect of insufficient calculation real-time performance caused by the random subspace method.

Description

Method and device for processing multichannel flutter test signals
Technical Field
The present invention relates to the field of aeronautical technology, and in particular, to a method and apparatus for signal processing.
Background
The flutter test is used as the verification of flutter design, is an unavoidable important link in the development process of the aircraft, and is an important task for accurately and effectively carrying out characteristic analysis on flutter response signals and further processing test data. Generally, the artificial excitation source (such as control surface sweep frequency, eccentric wheel excitation, rudder pedaling, small rocket excitation and the like) is implemented in the test, so that the quality of the flutter response signal can be effectively improved, the reliability of subsequent data processing conclusion is improved, however, the active excitation mode is often severely limited by the condition of the flutter test itself, such as the wind tunnel test is limited by the size and structure of a model, the active excitation is inconvenient to implement, and the aircraft needs to be modified by installing excitation equipment in test flight, so that the risk of test flight subjects such as dive, low-altitude large gauge speed and the like is further increased. Meanwhile, due to the general narrow-band characteristic of the active excitation signal, high-frequency mode information in the structure is not easy to sufficiently excite, and the analysis effect of a higher-frequency mode is further affected.
The turbulent flow excitation is an unavoidable excitation form in the flutter test flight/wind tunnel test, meanwhile, no additional modification is needed to be carried out on an aircraft/aeroelastic model, and the turbulent flow excitation belongs to broadband noise, so that the turbulent flow excitation becomes a necessary and effective excitation mode in high-risk test flight subjects, wind tunnel tests and high-frequency modal analysis. Because the atmospheric turbulence energy is dispersed and not measurable, although the turbulence response signal contains abundant and precious flutter test information, the randomness and low quality of the data often cause difficulty in modal analysis, thereby influencing the accurate acquisition of test conclusion.
Disclosure of Invention
The embodiment of the invention provides a signal processing method and device, which at least solve the problem of turbulent flow response signal analysis in a flutter test flight test due to the prior art in the related art.
According to one aspect of the embodiment of the invention, a signal processing method is provided, which comprises the steps of calculating a power spectrum density function matrix according to an acquired multi-channel turbulence response signal, acquiring a maximum singular value curve according to the power spectrum density function matrix, acquiring a self-power spectrum density function of a single-degree-of-freedom system according to the maximum singular value curve, and acquiring modal parameters of the single-mode system according to the self-power spectrum density function.
Optionally, obtaining the power spectrum density function matrix according to the obtained multi-channel turbulence response signal comprises performing power spectrum analysis on the multi-channel turbulence response signal to obtain a power spectrum density function matrix corresponding to each signal in the multi-channel turbulence response signal, wherein the power spectrum density function matrix is an autocorrelation power spectrum density function matrix or a cross correlation power spectrum density function matrix corresponding to each signal.
Further, optionally, obtaining the maximum singular value curve according to the power spectral density function matrix includes performing singular value decomposition on the power spectral density function matrix at each frequency point based on frequency domain orthogonality to obtain the maximum singular value curve of each frequency band.
Optionally, the method further comprises the step of carrying out singular value decomposition on the power spectrum density function matrix at each frequency point based on frequency domain orthogonality to obtain a left singular value vector corresponding to each frequency point, wherein the left singular value vector represents the mode shape of the system.
Further, optionally, obtaining the self-power spectral density function of the single degree of freedom system according to the maximum singular value curve comprises analyzing the left singular value vector based on a modal shape coherence criterion to obtain the self-power spectral density function of the single degree of freedom system of the corresponding frequency.
Optionally, obtaining the modal parameters of the single-mode system according to the self-power spectral density function comprises performing frequency domain fitting on the self-power spectral density function through an orthogonal polynomial, and performing modal parameter estimation on the single-degree-of-freedom system through polynomial fitting to obtain the modal parameters of the single-mode system, wherein the modal parameters of the single-mode system comprise frequency and damping of the single-degree-of-freedom system.
According to another aspect of the embodiment of the invention, a signal processing device is provided, which comprises a signal processing module, a first acquisition module, a second acquisition module and a third acquisition module, wherein the signal processing module is used for calculating a power spectrum density function matrix according to an acquired multi-channel turbulence response signal, the first acquisition module is used for acquiring a maximum singular value curve according to the power spectrum density function matrix, the second acquisition module is used for acquiring a self-power spectrum density function of a single-degree-of-freedom system according to the maximum singular value curve, and the third acquisition module is used for acquiring modal parameters of the single-mode system according to the self-power spectrum density function.
Optionally, the signal processing module comprises a signal processing unit, a power spectrum analysis unit and a power spectrum analysis unit, wherein the signal processing unit is used for carrying out power spectrum analysis on the multi-channel turbulence response signals to obtain a power spectrum density function matrix corresponding to each signal in the multi-channel turbulence response signals, and the power spectrum density function matrix is an autocorrelation power spectrum density function matrix or a cross correlation power spectrum density function matrix corresponding to each signal.
Further, the first acquisition module comprises a first acquisition unit, a second acquisition unit and a third acquisition unit, wherein the first acquisition unit is used for carrying out singular value decomposition on the power spectrum density function matrix at each frequency point based on frequency domain orthogonality to obtain a maximum singular value curve of each frequency band.
Optionally, the first acquisition module further comprises a decomposition unit, which is used for carrying out singular value decomposition on the power spectrum density function matrix at each frequency point based on frequency domain orthogonality to obtain a left singular value vector corresponding to each frequency point, wherein the left singular value vector represents the mode shape of the system.
According to the embodiment of the invention, a power spectrum density function matrix is calculated according to the acquired multi-channel turbulence response signals, a maximum singular value curve is acquired according to the power spectrum density function matrix, a self-power spectrum density function of a single-degree-of-freedom system is acquired according to the maximum singular value curve, and a modal parameter of the single-mode system is acquired according to the self-power spectrum density function. That is, the embodiment of the invention can solve the problem of turbulent flow response signal analysis in the flutter test flight test in the related technology, achieves the technical effects of improving the accuracy of modal parameter estimation and avoiding the defect of calculation instantaneity caused by a random subspace method.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a limitation on the application. In the drawings:
fig. 1 is a flow chart of a signal processing method according to an embodiment of the present invention;
FIG. 2 is a flowchart of another signal processing method according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a signal processing apparatus according to an embodiment of the present invention.
Detailed Description
In order that those skilled in the art will better understand the present invention, a technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in which it is apparent that the described embodiments are only some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
It should be noted that the terms "first," "second," and the like in the description and the claims of the present invention and in the drawings are used for distinguishing between different objects and not for limiting a particular order.
In the related art, for the frequency domain decomposition method and the processing of the turbulent flow response signal of the flutter test, the following signal processing methods mainly exist:
According to a traditional frequency domain decomposition method, a self-power spectrum density function of a single-degree-of-freedom system is calculated by carrying out frequency domain decomposition based on a modal coherence criterion, a single-degree-of-freedom time domain impulse response is calculated by interpolation and zero padding and then carrying out inverse Fourier transformation, and modal parameter estimation is carried out;
Aiming at the turbulent flow response signal analysis of the flutter test, generally, in the wind tunnel test, the turbulent flow response signal is processed through a random decrement technology because the test speed is closer to the flutter boundary, so as to obtain free attenuation response, and then time domain modal parameter estimation is carried out;
aiming at the situation that test data in flutter test flight is far lower than flutter speed and the interference of atmospheric turbulence is large, the modal parameter estimation of the turbulence response signal is carried out by a random subspace method.
However, the above signal processing schemes have better effects under certain technical conditions, but have certain technical defects as well:
in the traditional frequency domain decomposition method, the effect of calculating the time domain impulse response signal by the inverse Fourier transform is poor aiming at the inverse Fourier transform of the self-power spectrum of the single-degree-of-freedom system due to the problem of spectral line density, and the influence on the fitting result is larger;
aiming at the problems of random decrement technology commonly used in flutter wind tunnel test, such as atmospheric turbulence, principle flutter boundary and the like in flutter test flight, the error of the free attenuation response signal calculated by the technology is larger, and the effect of estimating the modal parameters of the flutter test flight is poorer;
The random subspace method has better performance on the modal parameter estimation of the flutter test flight, but in order to obtain a better estimation result, the subspace method obtains a large number of modal parameters through the modal parameter estimation of different orders, and the steady-state pattern method is based on the steady-state pattern method for identifying the steady mode, and the signal processing method has higher time requirements and is a larger defect in practical engineering application.
Therefore, the signal processing method provided by the embodiment of the application aims at the turbulent flow response signal of the flutter test, improves the traditional frequency domain decomposition method, directly carries out the modal parameter estimation of the orthogonal polynomial in the frequency domain, directly carries out the modal parameter estimation in the frequency domain, can better avoid estimation errors caused by insufficient spectral line density, improves the result of the modal parameter estimation, and only needs to obtain the self-power spectral density function of the system according to the technical route of the modal parameter estimation of the turbulent flow response signal of the flutter test, thereby obtaining the modal parameter result of the system by fitting in the frequency domain based on the scheme.
Specifically, according to an aspect of the embodiment of the present application, a signal processing method is provided, and fig. 1 is a schematic flow chart of the signal processing method provided by the embodiment of the present application. As shown in fig. 1, the signal processing method provided by the embodiment of the application includes:
step S102, calculating a power spectral density function matrix according to the acquired multi-channel turbulence response signals;
Optionally, the step S102 of calculating the power spectrum density function matrix according to the acquired multi-channel turbulence response signals includes performing power spectrum analysis on the multi-channel turbulence response signals to obtain a power spectrum density function matrix corresponding to each signal in the multi-channel turbulence response signals, where the power spectrum density function matrix is an autocorrelation power spectrum density function matrix or a cross correlation power spectrum density function matrix corresponding to each signal.
Specifically, fig. 2 is a schematic flow chart of another signal processing method provided by the embodiment of the present application, as shown in fig. 2, in the embodiment of the present application, power spectrum analysis is performed on a multi-channel turbulence response signal to obtain a self (mutual) correlation power spectrum density function matrix of the corresponding signal, and the multi-channel turbulence response signal obtains the power spectrum density function matrix through a periodic chart method.
Step S104, obtaining a maximum singular value curve according to a power spectrum density function matrix;
Optionally, the step S104 of obtaining the maximum singular value curve according to the power spectrum density function matrix includes performing singular value decomposition on the power spectrum density function matrix at each frequency point based on frequency domain orthogonality to obtain the maximum singular value curve of each frequency band.
Optionally, the signal processing method provided by the embodiment of the application further comprises the step of carrying out singular value decomposition on the power spectrum density function matrix at each frequency point based on frequency domain orthogonality to obtain a left singular value vector corresponding to each frequency point, wherein the left singular value vector represents the mode shape of the system.
Specifically, as shown in fig. 2, based on frequency domain orthogonality, singular value decomposition is performed on the power spectrum density function matrix at each frequency point to obtain a maximum singular value curve of each frequency band and a left singular value vector corresponding to each frequency point, and based on the power spectrum density function matrix obtained in step S102, the maximum singular value curve is obtained through singular value decomposition (Singular Value Decomposition).
Step S106, obtaining a self-power spectral density function of the single degree of freedom system according to the maximum singular value curve;
Optionally, the step S106 of obtaining the self-power spectral density function of the single degree of freedom system according to the maximum singular value curve comprises analyzing the left singular value vector based on the mode shape coherence criterion to obtain the self-power spectral density function of the single degree of freedom system with corresponding frequency.
Specifically, as shown in fig. 2, the left singular value vector in step S104 represents the mode shape of the system, the left singular value vector is analyzed based on the mode shape coherence criterion to obtain a self-power spectral density function of the single degree of freedom system with corresponding frequency, and the maximum singular value curve obtained in step S104 is based on the mode frequency to be analyzed to obtain the self-power spectral density function of the single degree of freedom system according to the mode amplitude coherence condition.
Step S108, obtaining the modal parameters of the single-mode system according to the self-power spectral density function.
Optionally, obtaining the modal parameters of the single-mode system according to the self-power spectral density function in step S108 includes performing frequency domain fitting on the self-power spectral density function through an orthogonal polynomial, and performing modal parameter estimation on the single-degree-of-freedom system through polynomial fitting to obtain the modal parameters of the single-mode system, wherein the modal parameters of the single-mode system include frequency and damping of the single-degree-of-freedom system.
Specifically, as shown in fig. 2, the frequency domain fitting is performed on the self-power spectral density function of the single-degree-of-freedom system obtained in step S106 through an orthogonal polynomial, and the modal parameter estimation is directly performed on the single-degree-of-freedom system through polynomial fitting, so that the frequency and the damping of the single-degree-of-freedom system are finally obtained.
In summary, the signal processing method provided by the embodiment of the application aims at the turbulent flow response signal analysis in the flutter test flight, can improve the accuracy of modal parameter estimation, and simultaneously avoids the problem of insufficient calculation instantaneity caused by a random subspace method.
According to the embodiment of the invention, a power spectrum density function matrix is calculated according to the acquired multi-channel turbulence response signals, a maximum singular value curve is acquired according to the power spectrum density function matrix, a self-power spectrum density function of a single-degree-of-freedom system is acquired according to the maximum singular value curve, and a modal parameter of the single-mode system is acquired according to the self-power spectrum density function. That is, the embodiment of the invention can solve the problem of turbulent flow response signal analysis in the flutter test flight test in the related technology, achieves the technical effects of improving the accuracy of modal parameter estimation and avoiding the defect of calculation instantaneity caused by a random subspace method.
According to another aspect of the embodiment of the present invention, a signal processing apparatus is provided, and fig. 3 is a schematic diagram of the signal processing apparatus provided in the embodiment of the present invention. As shown in fig. 3, the system comprises a signal processing module 32 for calculating a power spectrum density function matrix according to the acquired multi-channel turbulence response signal, a first acquisition module 34 for acquiring a maximum singular value curve according to the power spectrum density function matrix, a second acquisition module 36 for acquiring a self-power spectrum density function of the single-degree-of-freedom system according to the maximum singular value curve, and a third acquisition module 38 for obtaining a modal parameter of the single-mode system according to the self-power spectrum density function.
Optionally, the signal processing module 32 includes a signal processing unit, configured to perform power spectrum analysis on the multi-channel turbulence response signal, to obtain a power spectrum density function matrix corresponding to each signal in the multi-channel turbulence response signal, where the power spectrum density function matrix is an autocorrelation power spectrum density function matrix or a cross correlation power spectrum density function matrix corresponding to each signal.
Further, the first obtaining module 34 may optionally include a first obtaining unit, configured to perform singular value decomposition on the power spectrum density function matrix at each frequency point based on the orthogonality of the frequency domains, so as to obtain a maximum singular value curve of each frequency band.
Optionally, the first obtaining module 34 further includes a decomposition unit, configured to perform singular value decomposition on the power spectrum density function matrix at each frequency point based on the frequency domain orthogonality, to obtain a left singular value vector corresponding to each frequency point, where the left singular value vector represents a mode shape of the system.
Further, optionally, obtaining the self-power spectral density function of the single degree of freedom system according to the maximum singular value curve comprises analyzing the left singular value vector based on a modal shape coherence criterion to obtain the self-power spectral density function of the single degree of freedom system of the corresponding frequency.
Optionally, obtaining the modal parameters of the single-mode system according to the self-power spectral density function comprises performing frequency domain fitting on the self-power spectral density function through an orthogonal polynomial, and performing modal parameter estimation on the single-degree-of-freedom system through polynomial fitting to obtain the modal parameters of the single-mode system, wherein the modal parameters of the single-mode system comprise frequency and damping of the single-degree-of-freedom system.
The foregoing description is only of the preferred embodiments of the present invention, and is not intended to limit the scope of the present invention.

Claims (5)

1.一种信号处理的方法,其特征在于,包括:1. A signal processing method, comprising: 依据获取到的多通道紊流响应信号计算功率谱密度函数矩阵;包括:Calculate the power spectrum density function matrix based on the acquired multi-channel turbulence response signal; including: 对所述多通道紊流响应信号进行功率谱分析,得到所述多通道紊流响应信号中各信号对应的所述功率谱密度函数矩阵,其中,所述功率谱密度函数矩阵为所述各信号对应的自相关功率谱密度函数矩阵或互相关功率谱密度函数矩阵;Performing power spectrum analysis on the multi-channel turbulence response signal to obtain the power spectrum density function matrix corresponding to each signal in the multi-channel turbulence response signal, wherein the power spectrum density function matrix is an autocorrelation power spectrum density function matrix or a cross-correlation power spectrum density function matrix corresponding to each signal; 依据所述功率谱密度函数矩阵获取最大奇异值曲线;包括:Obtaining a maximum singular value curve according to the power spectrum density function matrix; comprising: 基于频域正交性,在各频点对所述功率谱密度函数矩阵进行奇异值分解,得到各频段的所述最大奇异值曲线;Based on frequency domain orthogonality, singular value decomposition is performed on the power spectrum density function matrix at each frequency point to obtain the maximum singular value curve of each frequency band; 基于频域正交性,在所述各频点对所述功率谱密度函数矩阵进行奇异值分解,得到所述各频点对应的左奇异值向量;其中,所述左奇异值向量表示系统的模态振型;Based on frequency domain orthogonality, singular value decomposition is performed on the power spectrum density function matrix at each frequency point to obtain a left singular value vector corresponding to each frequency point; wherein the left singular value vector represents a mode vibration shape of the system; 依据所述最大奇异值曲线获取单自由度系统的自功率谱密度函数;包括:Obtaining the auto-power spectral density function of the single degree of freedom system according to the maximum singular value curve; comprising: 基于模态振型相干准则对所述左奇异值向量进行分析,得到对应频率的所述单自由度系统的自功率谱密度函数;Analyzing the left singular value vector based on the modal vibration shape coherence criterion to obtain the autopower spectral density function of the single degree of freedom system at the corresponding frequency; 依据所述自功率谱密度函数得到单模态系统模态参数;包括:Obtaining the modal parameters of the single-mode system according to the auto-power spectral density function; including: 对所述自功率谱密度函数,通过正交多项式进行频域拟合,并通过多项式拟合的单自由度系统进行模态参数估计,得到所述单模态系统模态参数,其中,所述单模态系统模态参数包括:单自由度系统的频率和阻尼。The autopower spectral density function is fitted in the frequency domain by orthogonal polynomials, and modal parameters of the single-degree-of-freedom system fitted by the polynomials are estimated to obtain the modal parameters of the single-mode system, wherein the modal parameters of the single-mode system include: frequency and damping of the single-degree-of-freedom system. 2.一种实现权利要求1所述信号处理方法的装置,其特征在于,包括:2. A device for implementing the signal processing method according to claim 1, characterized in that it comprises: 信号处理模块,用于依据获取到的多通道紊流响应信号计算功率谱密度函数矩阵;A signal processing module, used for calculating a power spectrum density function matrix based on the acquired multi-channel turbulence response signal; 第一获取模块,用于依据所述功率谱密度函数矩阵获取最大奇异值曲线;A first acquisition module, used for acquiring a maximum singular value curve according to the power spectrum density function matrix; 第二获取模块,用于依据所述最大奇异值曲线获取单自由度系统的自功率谱密度函数;A second acquisition module is used to acquire the autopower spectral density function of the single degree of freedom system according to the maximum singular value curve; 第三获取模块,用于依据所述自功率谱密度函数得到单模态系统模态参数。The third acquisition module is used to obtain the modal parameters of the single-mode system according to the autopower spectral density function. 3.根据权利要求2所述的装置,其特征在于,所述信号处理模块包括:3. The device according to claim 2, characterized in that the signal processing module comprises: 信号处理单元,用于对所述多通道紊流响应信号进行功率谱分析,得到所述多通道紊流响应信号中各信号对应的所述功率谱密度函数矩阵,其中,所述功率谱密度函数矩阵为所述各信号对应的自相关功率谱密度函数矩阵或互相关功率谱密度函数矩阵。A signal processing unit is used to perform power spectrum analysis on the multi-channel turbulence response signal to obtain the power spectrum density function matrix corresponding to each signal in the multi-channel turbulence response signal, wherein the power spectrum density function matrix is an autocorrelation power spectrum density function matrix or a cross-correlation power spectrum density function matrix corresponding to each signal. 4.根据权利要求3所述的装置,其特征在于,所述第一获取模块包括:4. The device according to claim 3, wherein the first acquisition module comprises: 第一获取单元,用于基于频域正交性,在各频点对所述功率谱密度函数矩阵进行奇异值分解,得到各频段的所述最大奇异值曲线。The first acquisition unit is used to perform singular value decomposition on the power spectrum density function matrix at each frequency point based on frequency domain orthogonality to obtain the maximum singular value curve of each frequency band. 5.根据权利要求4所述的装置,其特征在于,所述第一获取模块还包括:5. The device according to claim 4, characterized in that the first acquisition module further comprises: 分解单元,用于基于频域正交性,在所述各频点对所述功率谱密度函数矩阵进行奇异值分解,得到所述各频点对应的左奇异值向量;其中,所述左奇异值向量表示系统的模态振型。A decomposition unit is used to perform singular value decomposition on the power spectrum density function matrix at each frequency point based on frequency domain orthogonality to obtain a left singular value vector corresponding to each frequency point; wherein the left singular value vector represents a modal vibration shape of the system.
CN202210674150.5A 2022-06-15 2022-06-15 Method and device for processing multi-channel flutter test signals Active CN115165270B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210674150.5A CN115165270B (en) 2022-06-15 2022-06-15 Method and device for processing multi-channel flutter test signals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210674150.5A CN115165270B (en) 2022-06-15 2022-06-15 Method and device for processing multi-channel flutter test signals

Publications (2)

Publication Number Publication Date
CN115165270A CN115165270A (en) 2022-10-11
CN115165270B true CN115165270B (en) 2025-03-28

Family

ID=83484480

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210674150.5A Active CN115165270B (en) 2022-06-15 2022-06-15 Method and device for processing multi-channel flutter test signals

Country Status (1)

Country Link
CN (1) CN115165270B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4811404A (en) * 1987-10-01 1989-03-07 Motorola, Inc. Noise suppression system
US6779404B1 (en) * 1999-11-03 2004-08-24 Rune Brincker Method for vibration analysis

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8010069B2 (en) * 2008-04-09 2011-08-30 Mstar Semiconductor, Inc. Method and apparatus for processing radio signals to identify an active system in a coexisting radio network
KR101786856B1 (en) * 2013-11-12 2017-10-18 울산과학기술원 Method for real time frequency domain decomposition
US9983776B1 (en) * 2016-06-30 2018-05-29 Bentley Systems, Incorporated Software system for dynamic feature extraction for structural health monitoring
CN107133195B (en) * 2017-04-14 2019-08-09 大连理工大学 A Model Ordering Method for Modal Identification of Engineering Structures
CN107238480B (en) * 2017-06-20 2019-02-12 西北工业大学 Damping calibration method in milling process based on working modal analysis
KR102033927B1 (en) * 2017-12-29 2019-11-08 연세대학교 산학협력단 Seismic Analysis Method considering refrigerant flow according to operating conditions of Nuclear power plant
US11170070B2 (en) * 2018-03-06 2021-11-09 Dalian University Of Technology Sparse component analysis method for structural modal identification when the number of sensors is incomplete
GB2577276A (en) * 2018-09-19 2020-03-25 Guided Ultrasonics Ltd Signal processing
CN110470450A (en) * 2019-08-27 2019-11-19 中国空气动力研究与发展中心高速空气动力研究所 Wind tunnel test flutter stability parameter prediction method and device
CN110470451A (en) * 2019-08-27 2019-11-19 中国空气动力研究与发展中心高速空气动力研究所 Wind tunnel test data processing method and device
CN114354170B (en) * 2022-01-07 2022-10-25 大连理工大学 A Structural Damping Ratio Identification Method Based on Unknown Impulse Response

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4811404A (en) * 1987-10-01 1989-03-07 Motorola, Inc. Noise suppression system
US6779404B1 (en) * 1999-11-03 2004-08-24 Rune Brincker Method for vibration analysis

Also Published As

Publication number Publication date
CN115165270A (en) 2022-10-11

Similar Documents

Publication Publication Date Title
CN106052849B (en) A method for identifying non-stationary abnormal noise sources in automobiles
CN107860548B (en) Approximate online flutter boundary prediction method
CN109459745A (en) A method of moving acoustic sources speed is estimated using radiated noise
US11832067B2 (en) Open-loop multichannel audio impulse response measurement and capture path evaluation
CN113642417B (en) Denoising method for partial discharge signal of insulated overhead conductor based on improved wavelet algorithm
CN110596247B (en) Ultrasonic structure health monitoring method in temperature change environment
CN104462862A (en) Mechanical structure dynamic load identification method based on cubic b-spline scaling function
CN108398260A (en) The fast evaluation method of gear-box instantaneous angular velocity based on mixing probabilistic method
CN113901379A (en) A dynamic online fast processing method for real-time data at the edge
CN109001672B (en) Time difference and frequency difference parameter estimation method and device
CN110263762B (en) An output-based energy transfer path analysis method for semi-submersible offshore platforms
US20220412793A1 (en) Method, device and computer program for monitoring a rotating machine of an aircraft
CN115165270B (en) Method and device for processing multi-channel flutter test signals
CN115902773A (en) A sound source localization method and system applicable to different microphone arrays
CN106407955B (en) A kind of low-frequency oscillation of electric power system mode on-line identification method
Janssens et al. Time-domain source contribution analysis method for in-room pass-by noise
CN112975574B (en) An online detection system for the surface quality of aluminum alloy thin-walled parts during milling
CN112688324B (en) Mode identification method for low frequency oscillation of power system based on FastICA and TLS-ESPRIT
CN114813129A (en) Acoustic signal fault diagnosis method of rolling bearing based on WPE and EMD
CN115014695B (en) Method and device for obtaining modal parameters
CN118836964A (en) Noise testing method, sound pressure signal testing system, analysis device and storage medium
Cheng et al. AR model-based crosstalk cancellation method for operational transfer path analysis
Wang et al. Time-frequency post-processing method based on generalized S-transform and its application to aero engine rotor system fault diagnosis
CN111353415A (en) Detection method of harmonic component in impulse response
CN104156509A (en) Method for synthesizing noise

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