50 results about "Transfer function matrix" patented technology

The invention relates to the technical field of multi-axis servo systems, and discloses a synchronous control method, device and equipment for a multi-axis servo system, and the method comprises the steps: measuring a mechanical coupling relation parameter, a motor dynamic characteristic parameter and a sampling delay time parameter of the multi-axis servo system, and substituting the parameters into a state equation modeling formula, calculating to obtain a discrete time system state matrix; performing Z transformation processing on the discrete time system state matrix and calculating frequency response characteristics to obtain a transfer function matrix and a target coupling axis pair set; executing hierarchical optimization to obtain a preliminary control instruction; a trajectory tracking error and an inter-axis synchronization error are calculated and substituted into the optimization objective function for sampling delay compensation processing, and a compensation control instruction is obtained; the compensation control instruction is processed in real time, synchronous control signals are output to the servo shafts to achieve multi-shaft coordinated movement, and the synchronous control precision in the multi-shaft servo system is improved by accurately recognizing the key coupling shaft pairs and the sensitive frequency points.

The invention relates to the technical field of production debugging control, and discloses a production debugging control method and system for a plastic container. The method comprises the following steps: arranging a plurality of different sensors in plastic container blow molding equipment, and simultaneously collecting blow molding process parameter data; performing wavelet threshold denoising and anomaly identification on the blow molding process parameter data to obtain a process parameter anomaly identification result; establishing a transfer function matrix between the process parameters and quality indexes according to the process parameter anomaly identification result, and obtaining a target process parameter set through multi-target optimization; and inputting the target process parameter set into a double-integral enhanced recurrent neural network for parameter regulation and control calculation to obtain a process parameter adjustment amount. According to the method, early detection and accurate positioning of process abnormity in the blow molding process are realized, the quality fluctuation risk and the defective product rate are greatly reduced, and the technical problem that different response characteristic parameters are difficult to coordinate is solved.

The invention relates to an impedance modeling method and device of a grid-connected converter system, a grid-connected system, computer equipment and a computer readable storage medium. The grid-connected converter system comprises a converter and a filtering module which are connected, the converter is used for connecting a direct-current power supply end, and the filtering module is used for connecting a power grid; the method comprises the steps of injecting a disturbance signal corresponding to a target operation parameter into a grid-connected point in a process that a grid-connected converter system operates at the target operation parameter; the grid-connected point is a connection point of the filtering module and a power grid; obtaining measurement electrical parameters of the grid-connected converter system; determining a transfer function matrix of the converter according to the measured electrical parameters; and determining an impedance model of the grid-connected converter system based on the transfer function matrix of the converter. In the process of constructing the impedance model by adopting the method, the transfer function matrix of the converter is obtained by collecting the real electrical parameters of the grid-connected converter system during operation, and the impedance model of the grid-connected converter system obtained based on the transfer function matrix is more accurate and reliable.

The invention provides an engineering structure random load power spectral density determination method and related equipment, and the method comprises the steps: obtaining a first power spectral density corresponding to an acceleration associated with an engineering structure, and carrying out the spectral decomposition of the first power spectral density, and obtaining a virtual acceleration response value; acquiring a transfer function matrix corresponding to each frequency point in the first power spectral density based on engineering parameters of the engineering structure; performing bilateral balance transformation on each transfer function matrix to obtain an improved transfer function matrix corresponding to each transfer function matrix; determining a virtual load based on the virtual acceleration response value and each improved transfer function matrix, and determining a second power spectral density corresponding to the random load based on the virtual load, the technical problem in the prior art that the identification error of the power spectrum density corresponding to the random load of the engineering structure under part of frequency points (especially the frequency points near the inherent frequency of the engineering structure) is large is solved, and the accuracy of determining the random load borne by the engineering structure is improved.

The invention discloses a single-loop droop control method and system for a wind generating set, and the method comprises the following steps: constructing a doubly-fed fan-power converter-power grid three-module coupling model, and completing the dynamic adaptation of the characteristics of a doubly-fed fan of a single-loop droop control architecture without inner-loop current control; based on the input signal-model collaborative integration body and the adaptive architecture, MIMO-GBC is adopted to analyze a phase margin and a resonance peak value of an open-loop transfer function matrix eigenvalue of a dynamic small signal model, and double risks are identified in a full slip rate interval; aiming at the resonance risk, a current induction type virtual impedance strategy is adopted; aiming at insufficient phase margin, adopting a transient phase calibration strategy; and according to the real-time slip ratio of the DFIG and the stator power change, inter-partition linkage optimization is carried out on the virtual resistance and the phase calibration angle so as to minimize a dual-power step response fluctuation square mean value, a parameter setting period and a rotor side converter control period are synchronously matched, and all-working-condition stable control is realized.

The application discloses a method for analyzing the interaction of multiple grid-connected inverters based on matrix norm theory, and aims at the interaction between grid-connected inverters, takes a power distribution system containing two grid-connected inverters as an analysis object, establishes an equivalent grid-connected model containing two inverters based on the Norton theorem, derives a transfer function matrix of a control system, and quantitatively analyzes the relationship between the interaction between inverters and the system frequency, equivalent grid-connected inductance and controller combination mode based on the spectral norm and condition number of the transfer function matrix, so as to provide a reference for the coordinated control of the inverters. A simulation model is built in Matlab / Simulink for analysis. The application provides a simple and effective method for analyzing the interaction of inverters, and provides a theoretical method for the coordinated control of the inverters.

The invention discloses an RTP equipment temperature control method and device based on feedforward decoupling. A heating model is obtained based on input and output response data fitting of a heating area of RTP equipment; selecting an element with the maximum lag time in diagonal elements of the heating model as a reference basis, constructing a diagonal target transfer function matrix containing time delay characteristics, and ensuring that the lag time of each diagonal element is not less than the actual physical lag time of a corresponding channel; a feedforward decoupling controller model is solved based on the heating model, a direct control loop matrix and a feedback control loop matrix are obtained through calculation, the direct control loop matrix and the feedback control loop matrix are deployed in a temperature control system, and closed-loop temperature control is conducted on the multiple heating areas of the RTP equipment in combination with a PID controller. According to the method, the advantages of low complexity and quick response of the feedforward decoupling algorithm are combined, the physical realizability problem of the decoupler is solved by introducing the target transfer function matrix, and high-precision multi-region temperature synchronous control is realized.

The application provides a boiler-turbine coupling process modeling method based on an improved particle swarm optimization algorithm, and relates to the field of thermal power generation technology. The method first establishes a three-input three-output boiler-turbine coupling process transfer function matrix model; collects historical operation data of a thermal power generating unit and performs zero-mean preprocessing to construct a training data set and a test data set for verification; then uses the improved particle swarm optimization algorithm to identify the parameters to be identified in the boiler-turbine coupling process transfer function model; substitutes the identified parameters into the boiler-turbine coupling process transfer function matrix model to form a complete boiler-turbine coupling process transfer function matrix model, and uses the test set to verify the model precision. The boiler-turbine coupling process transfer function matrix model established based on actual operation data provides a reliable theoretical basis and technical support for the subsequent advanced control strategy design and optimized operation of ultra-supercritical units.

The invention provides a symmetrical current three-wire common inductance direct current power flow control transient analysis method and system. The method comprises the following steps: S1, establishing a modal space model of a symmetrical current control three-wire common inductance type direct current power flow controller; s2, a dynamic segmentation method is adopted for the modal space model, and a transfer function matrix of the autonomous system is obtained; and S3, performing transient oscillation characteristic analysis on the autonomous system by using a singular value analysis method, and outputting an analysis result. And S4, analyzing the influence of different types of parameters on the transient oscillation contribution degree. According to the direct current power flow controller transient state analysis method based on singular value decomposition, transient state oscillation characteristics closer to multi-end direct current transmission are obtained in the three-wire direct current power flow controller topology, and reference is provided for more comprehensive and more-dimensional direct current power flow transient state oscillation suppression.

The invention relates to the technical field of power equipment state monitoring, and discloses a main transformer vibration characteristic intelligent analysis method and system, and the method comprises the steps: collecting vibration and current signals of a multi-load level in a health state, carrying out the decoupling of a vibration source, calculating a transfer function, and constructing a time-varying transfer function tensor representing the health state; during online monitoring, real-time signals are obtained, and a real-time transfer function matrix is calculated; interpolating the tensor according to a real-time load to obtain a health reference transfer function matrix matched with the current working condition; calculating a difference value between the real-time matrix and the health reference matrix, and obtaining a differential response field capable of stripping the influence of the working condition; and finally, on the basis of the differential response field, establishing and solving a physical inverse problem to realize reverse inversion and positioning of the abnormal state in the main transformer. According to the method, the interference of working condition change on vibration analysis is effectively solved, the decoupling of a fault source and the accurate positioning of internal abnormity are realized, and the accuracy and reliability of diagnosis are improved.

This invention proposes a secondary channel modeling method for a vibration active control system based on a top-mounted rotor hub. The method includes: activating the actuator of the vibration active control system when the rotor hub is not rotating, outputting actuating forces in both directions. These forces serve as load inputs to the test bench. The actuating forces are transmitted via a path that sequentially passes through the rotor hub actuator mounting location, the collector ring, the shaft, the spring, and the base, ultimately reaching the vibration sensor mounting location. The vibration sensor collects the vibration response at this location, identifying the secondary channel transfer function matrix in the fixed state. The secondary channel transfer function matrix is ​​then substituted into the output force transmission relationship formulas for the fixed and rotating states to obtain the initial value of the secondary channel model in the rotating state. Using this initial value, an adaptive algorithm is used as the physical model for online identification of the secondary channel, enabling online modeling of the system's secondary channels under rotor hub rotation conditions. This method allows for online iterative construction of time-varying secondary channel models, resulting in high modeling efficiency.

The application provides a decentralized compensation active disturbance rejection control method for an isobaric compressed air energy storage system, and is characterized by comprising the following steps: constructing an input-output relationship model based on a system transfer function matrix for a high-order multivariable isobaric compressed air energy storage system; determining a control structure based on different working stages; designing a compensation ADRC controller for a diagonal transfer function corresponding to each loop of the isobaric compressed air energy storage system; parameter tuning of the decentralized compensation ADRC controller is carried out based on a bandwidth parameterization method, gain parameters of an extended state observer and a state feedback control law in each loop are determined, and a designed decentralized compensation ADRC controller is obtained; and the designed decentralized compensation ADRC controller is applied to control of a compression stage and an expansion stage of the isobaric compressed air energy storage system. The application improves a traditional decentralized ADRC structure by introducing a high-order compensation link, and realizes effective control of a low-order controller on a high-order multivariable system.

An estimation device (101) performs a prescribed process for each of N×M×S elements in a first complex transfer function matrix that includes, as elements of a three-dimensional array of N×M×S, complex transfer functions respectively obtained for S subcarriers and N×M combinations, and thereby calculates a second complex transfer function matrix. In the prescribed process, an amplitude average of a plurality of first elements including an element to be processed is calculated, and the element to be processed is divided by the amplitude average. The plurality of first elements are included in S×M elements obtained for one reception antenna element which corresponds to the element to be processed.

The invention relates to the technical field of electrical engineering, and discloses a transient oscillation characteristic analysis method, system and device and a medium, and the method comprises the steps: obtaining the inertia characteristic of a receiving-end converter station, and building a small-signal control model of a DC power flow controller based on the inertia characteristic of the receiving-end converter station; performing first conversion operation on the small signal control model; solving the transfer function matrix, and obtaining the transient characteristics of the DC power flow controller according to the solving result; the method effectively solves the problem that transient oscillation analysis of the DC power flow controller is inaccurate due to the fact that inertia characteristics of a receiving end converter station are ignored in traditional analysis. Through accurate modeling and analysis, the influence of the parameters of the direct current power flow controller on transient oscillation under different system working conditions and parameters can be comprehensively evaluated, and powerful support is provided for stable operation of a power system.

The invention provides a networked control system tracking performance optimization method and system, and relates to communication signal tracking control, and the method comprises the steps: building a two-degree-of-freedom networked control system model; according to the two-degree-of-freedom networked control system model, obtaining an output signal and an error signal of the two-degree-of-freedom networked control system; according to the output signal, co-prime decomposition based on a rational transfer function matrix, a Bezout equation and a Youla parameterization form of a two-degree-of-freedom controller are calculated; the error signal is simplified; calculating an optimal performance expression of the two-degree-of-freedom networked control system model according to the simplified error signal; the optimal performance expression is simplified; calculating an optimal tracking performance expression of the two-degree-of-freedom networked control system model according to the simplified optimal performance expression; and optimizing the tracking performance of the networked control system according to the optimal tracking performance expression of the two-degree-of-freedom networked control system model.

The invention provides a transient oscillation characteristic analysis method and system for a single-loop PI control direct current power flow controller. The method comprises the following steps: constructing a transfer function matrix of a two-line common inductance type direct current power flow controller; dynamically dividing the direct current power flow controller into a full-control system and a semi-autonomous system based on the transfer function matrix; solving the maximum singular value of the semi-autonomous system by using the singular value decomposition technology, and researching the transient oscillation characteristic of the direct current power flow controller; the singular value decomposition result of the direct-current power flow controller semi-autonomous system is analyzed through frequency sweeping, the line resistance and the control reference value are comprehensively considered, and the influence rule of the key parameters of the direct-current power flow controller on the transient oscillation of the system is analyzed. According to the method, the active response characteristic of the direct current power flow controller to direct current micro-grid disturbance is fully utilized, reference is provided for transient oscillation suppression of the direct current micro-grid based on system dynamic, and direct current power flow control adapting to a direct current power grid complex scene can be achieved.

A method of recovering quantum error induced by non-markovian noise is disclosed. The method of recovering quantum error induced by non-markovian noise according to an exemplary embodiment of the present invention, includes obtaining a transfer function matrix including non-Markovian noise, obtaining an ideal operator for a requested gate operator, generating a recovery matrix by multiplying the ideal operator by an inverse matrix of a transfer function matrix including non-markovian noise, and multiplying the requested gate operator by a recovery matrix.

The present invention relates to a method for inverting the particle size spectrum of ultrafine particles under non-negativity constraints, comprising: obtaining a transfer function matrix and a corresponding response matrix; establishing a relationship between the transfer function matrix and the response matrix; using a logarithmic regularization term matrix to provide a fixed positive size distribution; fitting the logarithmic regularization term matrix, substituting an initial number concentration matrix and an initial coefficient matrix into a regularization solution formula, iteratively obtaining N regularization solutions, residual terms, and regularization terms; obtaining an optimal regularization parameter λ, and then obtaining an optimal coefficient matrix, and plotting a particle size spectrum. During the particle size spectrum inversion process, an L-curve and regularization method are used to address the morbidity of the transfer function matrix during the inversion process, reduce the particle size spectrum inversion error caused by response errors, address the non-negativity problem existing in the inversion process, and restore a sharp size distribution.

The invention discloses a multi-channel iteration method based on a multi-degree-of-freedom electro-hydraulic servo vibration test stand servo system decoupling method and a storage medium, and belongs to the technical field of electro-hydraulic servo control. In order to solve the problems of multi-degree-of-freedom coupling interference and insufficient waveform reproduction precision of a multi-channel electro-hydraulic servo vibration test stand, an H1 non-parameter identification method is adopted to carry out system identification to obtain a multi-channel transfer function matrix, and then a decoupling controller is designed based on the transfer function matrix to eliminate coupling interference among multiple channels. And finally, through correcting iteration dynamic optimization control input, the purposes of improving the control precision of the multi-degree-of-freedom electro-hydraulic servo vibration table and accurately simulating the structure reaction under the earthquake action are achieved. The test data is accurate and reliable. The invention is suitable for the civil engineering field, the traffic field, the mechanical field and the like.

The invention provides a dynamic reconstruction method for enabling an MIMO impedance matrix to be equivalent to an SISO complex impedance vector, and the method comprises the following steps: building a 2 * 2 linear impedance matrix of a converter subsystem and a power grid subsystem of a key dynamic state based on a control architecture and power grid parameters; a closed-loop transfer function matrix of the converter-power grid interconnection system is deduced, and an accurate characteristic value of the system is obtained and serves as a benchmark reference for stability analysis; a dynamic reconstruction method is applied, the 2 * 2 impedance model is converted into a complex vector expression, and finally an SISO impedance transfer function with the complete dynamic state reserved is obtained; and establishing an SISO open-loop transfer function of the interconnection system, carrying out cross validation on the SISO open-loop transfer function and a characteristic value analysis result, and realizing dynamic characteristic analysis through characteristic values. According to the method, the MIMO impedance matrix of the converter is accurately equivalent to the complex vector of SISO, and the dynamic characteristics of the system are completely reserved while dimension reduction mapping is realized, so that a more accurate dynamic analysis result is provided.

The invention discloses a plane symmetry rocket multivariable PD control parameter adjustment method based on frequency response, and belongs to the field of automatic control. Firstly, a transfer function matrix model is obtained according to a state-space equation of a plane symmetry rocket multivariable system; then, a reference model structure is determined, and an unexpected coupling effect is suppressed by enabling a closed-loop system to track a diagonal reference model matrix; then, obtaining a frequency domain model of a compensation controller required by the tracking reference model matrix by utilizing a sensitivity function matching condition of a closed-loop system and an expected closed-loop system; finally, a control gain calculation formula in a compensation controller matrix under PD control is comprehensively obtained based on the frequency response data; and a selection criterion of key variables in the frequency domain parameter adjustment method is designed. According to the method, mutual coupling among multiple channels of the plane-symmetric rocket can be processed, so that the response performance of a rocket multivariable system is improved, and the flight reliability is improved.

The invention discloses a self-adaptive control method and system for a lithium battery pole piece roller press. The method comprises the steps that controlled variable real-time data and online quality index real-time data in the rolling production process are obtained; wherein the controlled variable comprises at least one of main cylinder pressure, bending cylinder pressure, roller temperature and rolling speed, and the online quality index comprises pole piece thickness or surface density; based on the historical data of the controlled variables and the online quality indexes, a transfer function matrix model is established, and the transfer function matrix model is used for describing the dynamic response relation of the controlled variables to the online quality indexes; based on the transfer function matrix model and the online quality index real-time data, calculating the real-time adjusting quantity of the controlled variable through a model prediction control algorithm; a control instruction is generated according to the real-time adjusting quantity, and the control instruction is used for controlling the lithium battery pole piece roller press to adjust the controlled variable according to the real-time adjusting quantity. According to the invention, the consistency of mass production of the lithium battery pole pieces can be ensured, and the product quality stability is obviously improved.

The invention relates to the technical field of power electronics, and particularly provides a harmonic prediction and suppression method for a droop control converter based on a multi-frequency estimator, and the method comprises the steps: obtaining a main circuit parameter, a control parameter and a steady-state working point parameter of the droop control converter; determining a main circuit matrix of a harmonic state space according to the main circuit parameters based on the steady-state working point parameters, and determining a harmonic state space control loop matrix according to the control parameters; splitting by columns according to a target parameter matrix of the main circuit matrix, and separating components related to the switching disturbance to obtain a plurality of parameter sub-matrixes corresponding to the switching disturbance; the control loop matrix is substituted into the parameter sub-matrix, a harmonic state space equation is obtained, a harmonic transfer function matrix is calculated, and the harmonic transfer function matrix is used for conducting harmonic prediction according to the main circuit parameters, the control parameters and the steady-state working point parameters. According to the invention, the harmonic prediction accuracy is effectively improved.

The invention relates to a sound field partition control method, which comprises the following steps: determining bright areas and dark areas in a sound field of a vehicle, and obtaining a transfer function matrix corresponding to at least one bright area and a transfer function matrix corresponding to at least one dark area; based on the transfer function matrix corresponding to the bright area and the transfer function matrix corresponding to the dark area, determining an error sample corresponding to the bright area and an error sample corresponding to the dark area; based on the error sample corresponding to the bright area and the error sample corresponding to the dark area, correcting the correlation matrix corresponding to the bright area and the correlation matrix corresponding to the dark area; and establishing a sound field partition control model based on the obtained target bright region correlation matrix and the target dark region correlation matrix, and adjusting a driving signal of a loudspeaker array based on the sound field partition control model so as to perform partition control on the sound field of the vehicle. According to the invention, the robustness of the sound field partition control model can be improved, and the partition control of the sound field can be realized under different temperatures and humidity.

The invention provides a machine-furnace coupling process modeling method based on an improved particle swarm optimization algorithm, and relates to the technical field of thermal power generation. The method comprises the following steps: firstly, establishing a three-input three-output boiler-steam turbine coupling process transfer function matrix model; historical operation data of the thermal generator set are collected, zero-mean preprocessing is carried out, and a training data set and a test data set for verification are constructed; identifying parameters to be identified in the boiler-steam turbine coupling process transfer function model by adopting an improved particle swarm algorithm; and substituting the identified parameters into the boiler-steam turbine coupling process transfer function matrix model to form a complete turbine-steam turbine coupling process transfer function matrix model, and verifying the model precision by using a test set. According to the method, a boiler-steam turbine coupling process transfer function matrix model is established based on actual operation data, and reliable theoretical basis and technical support are provided for advanced control strategy design and optimized operation of a subsequent ultra-supercritical unit.

The invention discloses a bushing dynamic stiffness test method, device, test equipment, system and medium. The method comprises the following steps: exciting two subsystems, and determining an initial transfer function matrix; reducing the initial transfer function matrix, and determining an original point transfer function matrix corresponding to the two subsystems and an original point transfer function matrix corresponding to the coupling system; and determining the dynamic stiffness of the bushing based on the original point transfer function matrixes corresponding to the two subsystems and the original point transfer function matrix corresponding to the coupling system. According to the method, the dynamic stiffness of the bushing is tested through the coupling system connected with the bushing, so that the problem that the dynamic stiffness of the bushing is inaccurately tested in a high-frequency range due to resonance and other problems of existing test equipment can be effectively avoided, and the relatively accurate and stable dynamic stiffness test can be performed on the bushing in a relatively large frequency range.

The invention relates to a power system stability evaluation method, device and equipment, a medium and a program product. The method comprises the following steps: determining an output differential passivity index of to-be-evaluated equipment according to the equipment type and the equipment parameter by acquiring the equipment type and the equipment parameter of the to-be-evaluated equipment in the power system, and determining the output differential passivity index of the to-be-evaluated equipment according to a power flow model of a transmission network in the power system. And determining a transfer function matrix of a transmission network for representing a relationship between the to-be-evaluated equipment input and the to-be-evaluated equipment output, and determining a stability evaluation result of the power system based on the output differential passivity index and the transfer function matrix, the stability of the power system is evaluated in a distributed manner based on the output differential passivity index of the to-be-evaluated device, that is, the stability of the power system can be evaluated only by providing the passivity index of each distributed device, so that the stability of the novel power system can be better evaluated.

The invention discloses an ultra-supercritical unit flexible operation scheme based on improved error active disturbance rejection control and an improved northern eagle optimizer. Firstly, a coordinated control system of an ultra-supercritical unit is used as a research object, and modeling and control difficulties of the coordinated control system are analyzed; then, a multi-working-condition transfer function matrix model of the unit is identified based on an improved northern eagle optimizer; and then, an improved error active disturbance rejection control strategy is designed. And finally, verifying the feasibility of the flexible operation scheme based on a simulation platform, and carrying out quantitative statistic analysis on the effectiveness of the control strategy by adopting a performance index. According to the invention, the reduced-order extended state observer is applied to the error active disturbance rejection controller, the complexity of the controller structure is reduced, the rapid load control performance of the system is effectively improved, and the flexible operation capability of the unit under the wide load working condition is improved.

The application relates to the technical field of power system modeling and analysis, in particular to a method and device for establishing an MMC wide-frequency three-port frequency coupling impedance model, wherein the method comprises the following steps: a three-dimensional transfer function matrix of an MMC internal control link is established to determine a reference voltage generated by the MMC internal control link; based on the reference voltage and disturbance components of AC and DC sides of the MMC, dynamic quadratic iteration solving of MMC internal bridge arms is carried out to obtain actual output voltages of upper and lower bridge arms of the MMC; the actual output voltages of the upper and lower bridge arms are used to calculate a differential mode voltage and a common mode voltage of the MMC; and a wide-frequency three-port frequency coupling impedance model of the MMC for outputting stability of wide-frequency oscillation of a flexible AC / DC transmission system is established by using a differential mode voltage loop and a common mode voltage loop. Therefore, the problems that the impedance model established in the related art has a narrow applicable frequency range and low precision and cannot meet the requirements of wide-frequency oscillation stability analysis are solved.

The application relates to the technical field of vehicles, in particular to a noise identification method and device of a vehicle, a vehicle and a storage medium, wherein the method comprises the following steps: acquiring a vibration acceleration signal and a vehicle interior noise signal when the vehicle drives on a preset impact road and performing time domain preprocessing; constructing a cross spectrum matrix based on the vibration acceleration signal after the time domain preprocessing and performing target principal component analysis to obtain a principal component signal of the cross spectrum matrix; then obtaining a transfer function matrix from each wheel center excitation point to an interior response point by using a hammering method test; calculating and identifying road impact loads at each wheel center of the vehicle based on the principal component signal and the transfer function matrix by using an inverse matrix method; and calculating contribution amount quantization results of each noise transmission path to the vehicle interior noise signal according to the road impact loads and the transfer function matrix. Therefore, the problems of noise identification result distortion and instability caused by the ill-conditioned transfer function matrix and the coupling of multiple reference point signals in the related art are solved.