119 results about "Residual sum of squares" patented technology

The invention discloses a high-order multi-stage auto-regressive distributed lag modeling method of thermal error compensation of a numerical control machine. The method is characterized by comprising the following steps of: providing a high-order multi-stage auto-regressive distributed lag modeling formula of a thermal error of a numerical control machine containing a coefficient to be solved, selecting numerical control machine thermal error lag phases and numerical control machine temperature lag phases as 1, 2, 3 and 4 respectively, and substituting experimental data to fit the coefficient to be solved in the formula according to a least square method so as to determine high-order multi-stage auto-regressive distributed lag models of thermal error compensation of the numerical control machine in different lag phases; substituting the experimental data into each model to obtain a residual sum of squares of each model; and substituting the residual sum of squares of each model into an akaike information criterion to determine an optimal lag phase to determine a high-order multi-stage auto-regressive distributed lag model of a thermal error of the numerical control machine. The invention discloses a modeling method of thermal error compensation of the numerical control machine, which has the advantages of convenience in application, easiness in modeling and high stability, and has higher accuracy than a traditional ADL (automatic data logger) model.

The invention discloses an output-only linear time-varying structure modal parameter identification method and belongs to the technical field of structural dynamics. Firstly, a cost function of a least squares support vector machine vector time-varying autoregressive model is deduced; secondly, a function space is built by means of a Wendland compactly supported radial basis function; a regular factor is determined through the non-parameter method based on Gamma testing, and a basis function width reduction coefficient is given on the basis of actual experiences; a time-varying autoregressive model order is determined according to the Bayesian information criterion and the Akaike information criterion; a function space order is determined according to the ratio of residual sum of squares to sequence sum of squares; finally, the matrix expression of the least squares support vector machine vector time-varying autoregressive model is solved according to the cost function, modal frequency of a system is solved according to a time freezing method, and linear time-varying structure modal parameter identification is finished. The method can improve calculation efficiency, improves system robustness, and is widely used in linear time-varying structure modal identification in structural dynamic engineering application.

The method for forming the optimal characteristic curve of a solar cell comprises the steps of: providing a first acceptable error; and determining a current-voltage polynomial regression equation, whose square root of the residual sum of square is less than the first acceptable error for a set of solar cell measured data. The order of the current-voltage polynomial regression equation is gradually increased until the square root of the residual sum of square of the current-voltage polynomial regression equation is less than the first acceptable error.

A method for tire parameter derivation, tire cornering characteristic calculation and tire design is used with a tire dynamic model constituted by using a plurality of tire dynamic element parameters including stiffness and friction coefficient and parameter defining a distribution of contact pressure of the tire. The parameters and tire cornering characteristic are derived by using the combined sum of squared residuals being obtained by weighted addition of a first sum of squared residuals of lateral force and a second sum of squared residuals of self-aligning torque. The tire dynamic model is a model for calculating a lateral force and for calculating a self-aligning torque separately as a lateral force-based torque component generated by the lateral force applied on a contact patch of the tire and a longitudinal force-based torque component generated by a longitudinal force applied on the contact patch of the tire.

A method of identifying endmember spectral values from multispectral image data, where each multispectral data value is equal to a sum of mixing proportions of each endmember spectrum. The method comprises the steps of processing the data to obtain a multidimensional simplex having a number of vertices equal to the number of endmembers. The position of each vertex represents a spectrum of one of the endmembers. Processing the data is conducted by providing starting estimates of each endmember spectrum for each image data value. The mixing proportions for each data value is estimated from estimates of the spectra of all the endmembers. The spectrum of each endmember is estimated from estimates of the mixing proportions of the spectra of all the endmembers for each image data value. The estimation steps are repeated until a relative change in the regularised residual sum of squares is sufficiently small. The regularised residual sum of squares includes a term which is a measure of the size of the simplex.

The invention relates to an equivocation determining method, in particular to a dynamic equivocation determining method. The invention adopts an LAMBDA method to search a plurality of groups of equivocations, and then utilizes the square sum of residuals which is alternatively combined to detect, thereby the reliability of the equivocation determination can be increased. The invention utilizes the known wide-lane equivocation and L1 equivocation to express the non-ionosphere equivocation combination for carrying out the equivocation search so that the influence on a first-order phase of an ionosphere is eliminated, and the integer performance of the searched equivocation can be caused, thereby the distance of the equivocation search is enlarged, various conditions for rapidly and dynamically determining the long-distance equivocation are fully utilized, and the equivocation search is accelerated.

The invention discloses an optimum eight-place standardization method of a flexible gyroscope, the flexible gyroscope is arranged on a three-shaft place rate rotary table, wherein, the flexible gyroscope is connected with data collecting equipment which is communicated with a computer. The invention discloses a specific position in the optimum eight-place standardization, and the output of the flexible gyroscope is effectively improved in the way that the invention compensates a measurement value of a compensation model G0 of obtained optimum eight-place drifting coefficients and vibrating reed gyro quiet errors. In a testing process of the vibrating reed gyro in an inertial navigation testing center, a traditional eight-place method and the optimum eight-place method are respectively adopted to obtain the drifting coefficients which respectively and objectively output compensated evaluation results of another three-place top in space. By residual sum of square of a top measurement value, compared with the traditional eight-place method, the results of compensated drifting coefficients solved by the optimum eight-place experiment design method of the flexible gyroscope are increased by 4 to 8 times.

The invention discloses a method for forecasting the acute toxicity of organic compounds by building a quantitative structure-activity relationship model with a quantum chemistry method. The method fully geometrically optimizes compound structures by using a Gaussian procedure so as to obtain quantum chemistry parameters including molecular volume, relative molecular mass, highest occupied molecular orbital energy, lowest unoccupied molecular orbital energy, energy gaps of frontier molecular orbital, dipole moment, solvation energy, electron energy and the like; using the quantum chemistry parameters and a hydrophobicity parameter as structural descriptors; in combination with toxicity data, quantitative relationship equations between various structural descriptors and toxicity are established according to a written procedure based on partial least square stepwise linear regression to obtain the multiple correlation coefficient, F-test value and sum of squared residuals, and then the model is verified so as to guarantee the external predictive ability. Therefore, the method can quickly and effectively forecast the toxicity of organic compounds to be studied, and provide necessary basic data for risk assessment and supervision of chemicals.

The invention discloses an optimization solution method for locating parameters of a micro-earthquake source, which takes the minimum residual square sum of the difference between the arrival time picked up by two vibration pickups and the calculated arrival time as an objective, calculates the source position and the equivalent wave velocity by using particle swarm algorithm, and directly calculates the starting time according to the source position and the equivalent wave velocity. The invention provides constraint rules of parameters such as inertia weight, acceleration constant, maximum flying velocity of particles and the like in the PSO solution algorithm in combination with engineering practice, avoids algorithm failure due to improper selection of these parameters, improves the robustness of the solution algorithm, achieves the purpose of reducing iteration times, improving positioning accuracy and algorithm speed.

In a crystal structure analysis using an X-ray diffraction method or the like, using a measured value y0i of a sample and an expected crystal structure parameter, a vector F having the logarithmic conversion value of the measured value as a matrix element is determined by fi=kxlog(yi+ delta -bi) from a value yi of the measured value y0i after count missing correction of a detector, a background strength bi and a positive value delta of less than 1, a vector Fc having the logarithmic conversion value of the calculated value as a matrix element is obtained by Fci=kxlog(yci) from the value vector obtained by calculation from the crystal structure parameter, a weight matrix W is obtained from device function matrix, systematic error, and accidental error, and the calculated value vector Fc is determined so that a residual square sum s (=(F-Fc)tW(F-Fc)) obtained by multiplying a difference between F and Fc by the weight matrix is minimized and converged to obtain a crystal structure parameter.

An optimum measurement subset with a specified number of elements is generated from a set of input global navigation satellite system (GNSS) measurements. A design matrix and a weight matrix are generated. Values of a set of coefficients corresponding to the set of input GNSS measurements are calculated. The value of a specific coefficient is calculated as the ratio of the change in value of the at least one target parameter resulting from the change in value of the specific input GNSS measurement to the change in value of the sum of squared residuals resulting from the change in value of the specific input GNSS measurement. The optimum measurement subset is selected based at least in part on the values of the set of coefficients.

The invention provides a method for hotspot discovery and evolution situation analysis in a network forum, using a crawler program to extract the HTML text of all posts in a forum; extracting pageviews and replies from each post text to form a binary group, assuming the first The binary group of the i post is (xi, yi); use the formula to calculate the score of the i-th post; obtain a popular post, and count the number of replies that occur every day from the topic posting to the capture time according to the reply record to get the binary group . The fitting effect and evaluation index are listed below: SSE=1.548e+07. SSE is the sum of squares of the error term, which reflects the discrete state of each observation value of each sample, and is also called the sum of squares within the group or the sum of squares of the residual error. R-square = 0.8339. R-square is the fitting coefficient, the larger the value, the better the fitting degree. RMSE = 525.7. RMSE is the root mean square error, which can be used as a numerical indicator to measure the measurement accuracy. After analyzing these indicators, it can be seen that the fitting effect is relatively ideal. Find the largest extreme point as x0=14; take x1=15>x0, y'(x0)<0, so the popularity of this post is declining.

The invention discloses a method for forecasting data. The method comprises the following steps of: acquiring a trend value and a residual value of each observation point; determining the principal direction of data forecasting; calculating the distances between one observation point and other observation points one by one, and included angles between the connection lines between one observation point and other observation points and the principal direction; classifying the observation points; calculating the average residual sum of squares of each type of the observation points respectively according to the residual values respectively; and calculating the distance between a target point and each observation point and an included angle between the connection line between the target point and each observation point and the principal direction, selecting the average residual sums of squares corresponding to the distances and the included angles, and forecasting the observed value of the target point. By applying the technical scheme, the distance correlation and the direction correlation are considered comprehensively in the data forecasting process, the accuracy of the average residual sum of squares is improved relative to Kriging interpolation in which only the distance correlation is considered, so that accuracy of the observed value forecasting of the target point is improved.

Provided is a cable running state assessment method and system. The method is characterized by obtaining characteristic quantity of a cable, wherein the characteristic quantity comprises insulation resistance value, insulation resistance absorption ratio, average value of dielectric loss angle tangent values, linear fitting slope, residual sum of squares, partial discharge quantity, partial discharge quantity increment, partial discharge frequency, concentrated partial discharge point positions and harmful discharge information; judging grade of each characteristic quantity according to a preset standard; and obtaining a state assessment result according to the weight corresponding to the grade of each characteristic quantity. According to the cable running state assessment method and system, the obtained characteristic quantity comprises the insulation resistance value, the insulation resistance absorption ratio, the average value of dielectric loss angle tangent values, the linear fitting slope, the residual sum of squares, the partial discharge quantity, the partial discharge quantity increment, the partial discharge frequency, the concentrated partial discharge point positions and the harmful discharge information, so that the running state of the cable cam be evaluated more completely and accurately.

The invention provides a method and device for predicting the output of an oil and gas well, and the method comprises the steps: obtaining the actual effective output data of the oil and gas well in a plurality of historical production time periods; respectively predicting the output of the oil and gas well in at least one historical production time period and one future production time period through employing a hyperbolic decline method, an exponential decrease method and a curve fitting method and the actual effective output data of the oil and gas well; calculating the residual sums of squares of prediction results of the above three methods according to various types of obtained output prediction data of the oil and gas well and the actual effective output data of the oil and gas well; selecting the minimum residual sum of squares, and enabling the data, obtained through the method which obtains the minimum residual sum of squares, of the oil and gas well in the future production time period to serve as the final output prediction result of the oil and gas well. The method selects the oil and gas output prediction result with the highest accuracy from the results of three prediction methods through residual error evaluation.

The invention relates to an adaptive data driving fault diagnosis method and a device in a complex refining process. The method comprises steps that, a relation graph among relevant process parameters causing a fault is established, and the process parameters causing the fault are determined; time sequence of all the process parameters and time sequence of the fault are acquired through the process parameters; whether the time sequence of the process parameters is covariance stationary is determined; if not, first order difference processing on the time sequence of the process parameters is carried out; a regression equation is constructed by utilizing the time sequence of the process parameters and the time sequence of the fault after covariance stationary, and a first residual error quadratic sum and a second residual error quadratic sum are acquired according to the regression equation; F examination is carried out; a largest path of a cause and effect relation magnitude is searched according to the relation graph; the time sequence of the process parameter which corresponds to a start point of the path is the root reason of the fault, and the terminal point of the path corresponds to the time sequence of the fault; and the largest path of the cause and effect relation magnitude is the propagation path of the fault in a system.

The invention discloses a Fresnel matrix THz wave propagation model based high-precision thickness detection method for coatings. Based on a electromagnetic field theory model, a reflecting propagation one-dimensional electric field model of normal incidence THz wave in a multi-layer medium structure is established according to a reflecting type THz time domain spectrum system used in coating thickness detection. On this basis, a Fresnel coefficient matrix and a phase matrix are introduced and a Fresnel matrix THz wave propagation model is proposed and validated. Quantitative detection of coating thickness is realized by utilizing the Fresnel matrix THz wave propagation model based LS (Least Squares) thickness optimization method and the coating thickness is solved through iterative optimization so as to minimize the minimal residual sum of squares of a simulation result and a measurement result and to realize the high precision optimization of thickness parameters. A distance parameter DTS between a THz emitter and a to-be-measured piece is introduced and visible treatment is performed on distance errors between the THz emitter and the to-be-measured piece, so that the accuracy ofcoating thickness detection is improved further.

The invention discloses a multiple linear regression model-based belt weigher main error factor analysis method. The method includes the following steps that: a belt weigher experiment platform is set up, main error factors such as tensions, temperatures and equivalent flow rates are changed according to actual situations, the tensions, temperatures, equivalent flow rates, calibration values and hopper weighers of various test points are recorded, and the relative errors of the calibration values and hopper weighers are calculated; the correlation coefficients of the main error factors of a belt weigher are calculated, and the correlation of the main error factors is judged; the relations of the factors are determined according to the correlation, and a multiple linear model is set and solved; residual sum of squares, determination coefficients and MS variance of residuals are adopted as check indexes, and the regression effect of the established model is evaluated; and the test values of the test points are predicted through using a fitting model, and the test values are compared with practical values, so that prediction errors can be calculated, and the accuracy of a fitting result is determined. With the method of the invention adopted, a theoretical basis is provided for quantitative analysis on the degree of influence of the main error factors on the accuracy of the belt weigher.

The invention provides an estimation method of external network Thevenin equivalent parameters based on single point measurement information, relating to the power system external network Thevenin equivalent parameters estimation method technology field. According to the invention, a computer and a data acquisition apparatus of a line terminal are utilized, through a program, firstly a parameter of a monitored line is inputted, then according to SCADA measurement data of a plurality of periods of a terminal of the monitored line, a least square equivalence model with minimum of square sum of residuals as an object is established, the model is solved, and equivalent built-in potential and an equivalent impedance parameter of a network outside of the monitored line are obtained. According to the invention, a physical meaning of the external network Thevenin equivalent parameters is fully considered, rationality of the external network equivalent parameters can be effectively ensured, influence of a measurement error on an equivalent parameters calculation result can be effectively avoided, accuracy is high, the method is simple, estimation precision is high, a calculation speed is fast, and the method can be widely applied to the estimation of the external network Thevenin equivalent parameters of any power line whose terminal is installed with the data acquisition apparatus.

The invention relates to a method for processing the reflector mirror shape of a high-precision space camera. According to the method, the node coordinates of a deformed mirror are rotated, adjusted in scale and translated to realize the coordinate transformation. After the coordinate transformation process, the residual sum of squares between the node of the deformed mirror and the curved surface of an ideal mirror is minimized. After that, the approximate values of a scaling factor, a coordinate translation amount and a coordinate rotation amount are solved in the least-squares sense. The scaling factor is adopted as an estimated parameter and a fitting equation approximate to an ideal equation is fitted based on known data points. Meanwhile, based on the scaling factor obtained through calculation, a curved surface equation fitted to deformed nodes can be determined. The above processing method lowers the data of the residual error and improves the calculation accuracy. Meanwhile, the credibility and the reliability of fitting results are improved.

An optimum measurement subset with a specified number of elements is generated from a set of input global navigation satellite system (GNSS) measurements. A design matrix and a weight matrix are generated. Values of a set of coefficients corresponding to the set of input GNSS measurements are calculated. The value of a specific coefficient is calculated as the ratio of the change in value of the at least one target parameter resulting from the change in value of the specific input GNSS measurement to the change in value of the sum of squared residuals resulting from the change in value of the specific input GNSS measurement. The optimum measurement subset is selected based at least in part on the values of the set of coefficients.

The method for forming the optimal characteristic curve of a solar cell comprises the steps of: providing a first acceptable error; and determining a current-voltage polynomial regression equation, whose square root of the residual sum of square is less than the first acceptable error for a set of solar cell measured data. The order of the current-voltage polynomial regression equation is gradually increased until the square root of the residual sum of square of the current-voltage polynomial regression equation is less than the first acceptable error.

The invention provides a method for confirming yield through spectrum unfolding analysis of a stratum element logging instrument. The method comprises the following steps of confirming a segment area number n and left and right boundaries kn and kn+1 of each segment area; calculating a standard trap spectrum normalization response matrix A; calculating a total numerical value ci of a spectrum to be unfolded in each segment area; calculating an opposite angle weighting matrix W; calculating a weighted least square method spectrum unfolding to obtain the yield yj; utilizing a mean value delta 2 of residual sum of squares to judge an imitative effect. The invention adopts a segment clamping manner spectrum unfolding method. The segment clamping manner is characterized by selectively selecting a channel area, in which an characteristic energy spectrum of an element is located, to unfold the spectrum. Through MC numerical simulation, the method has no large effect on spectrum unfolding.

The invention discloses an evaluation method and system of water quality based on fish electrocardio indexes. The method comprises the following steps: receiving fish body electrocardio signals beforeand after fish contact with a solution which contains specific organic matter, and respectively extracting multiple kinds of electrocardio indexes; analyzing the relativity between environmental stress and each extracted electrocardio index according to a Pearson correlational analysis method, and screening the electrocardio indexes for evaluating the specific organic matter in the water quality;analyzing and evaluating environmental stress and linear regression models between various electrocardio indexes according to Akaike information criteria and the residual sum of squares so as to obtain the electrocardio indexes used for finally evaluating the organic matter in the water quality. Receive the fish body electrocardio signals before and after the fish contact with to-be-determined water, and respectively extracting an electrocardio index which is corresponding to a certain organic matter in determined evaluation water quality; contrast the extracted electrocardio indexes of the fish body electrocardio signals before and after the fish contact with the to-be-determined water so as to obtain electrocardio index changes for evaluating the corresponding organic matter in the to-be-determined water quality.

The invention discloses a statistical method of the cumulative probability distribution of multi-section lightning current amplitude. Flashes of different frequencies of the same line are compared to fit a statistical result, a fit curve and a fitting residual error quadratic sum of the cumulative probability of the lightning current amplitude, and 2000 frequencies of lightning flashes are selected as threshold values to segment a line corridor. After the line corridor finishes being divided, curve fitting is carried out by sections to obtain the cumulative probability distribution formula of the lightning current amplitude among the intervals of a power transmission line corridor. Therefore, cumulative probability distribution parameter, which is obtained by statistics, of the lightning current amplitude of the power transmission line corridor can more accurately reflect the lightning activity routines of the line, and guide is provided for the anti-lightning calculation and the engineering design of the line.

The invention provides a method for determining the yield of a stratum element logging instrument by virtue of spectrum unfolding. The method comprises the following steps: calculating a standard capture spectrum normalization response matrix A, calculating actual measurement stratum capture spectrum count ci, calculating an opposite angle weighting matrix W, calculating the yield yj obtained by virtue of spectrum unfolding based on a weighted least squares method, and judging the fitting effect by use of residual sum of squares and the average number delta 2. According to the method, a one-by-one spectrum unfolding method is adopted, complete 256 energy windows are adopted for spectrum unfolding. According to MC numerical simulation, the method has slight influence on spectrum unfolding; a one-by-one method is adopted for performing spectrum unfolding on an unknown stratum, resolution result errors are large but can basically reflect the types and relative contents of elements in the stratum.

A method of identifying endmember spectral values from multispectral image data, where each multispectral data value is equal to a sum of mixing proportions of each endmember spectrum. The method comprises the steps of processing the data to obtain a multidimensional simplex having a number of vertices equal to the number of endmembers. The position of each vertex represents a spectrum of one of the endmembers. Processing the data is conducted by providing starting estimates of each endmember spectrum for each image data value. The mixing proportions for each data value is estimated from estimates of the spectra of all the endmembers. The spectrum of each endmember is estimated from estimates of the mixing proportions of the spectra of all the endmembers for each image data value. The estimation steps are repeated until a relative change in the regularised residual sum of squares is sufficiently small. The regularised residual sum of squares includes a term which is a measure of the size of the simplex.

The invention discloses a Spectralon diffuse reflection plate correcting method. The method comprises the steps: 1) illuminating the Spectralon diffuse reflection plate from different zenith angles by a light source, and using the spectrograph to collect the corresponding reflection spectroscopic data from the different zenith angles and the azimuth angles; 2) intercepting the spectroscopic data within the valid wave band range from the reflection spectroscopic data, calculating the average value, and eliminating the spectroscopic data with relatively large difference; 3) calculating the residual sum of squares Rss of the residual spectroscopic data, determining the value of the most stable probe receiving zenith angle (thetar, degree), and establishing the relational expression of the random probe receiving zenith angle thetar and (thetar, degree); and 4) substituting the relational expression to the BRDF expression of the Spectralon diffuse reflection plate, and finishing the Spectralon diffuse reflection plate correction. The method is capable of using the known data of the Spectralon diffuse reflection plate to correct, calculating the universal expression of the probe receiving zenith angle orientation, and shortening the time for the follow-up test process.

An LS (least square)-based multi- model adaptive state estimation method includes the following steps: S1, same moving object observation sequences of sensors are acquired to form moving object X, Y and Z direction observation sequences in the space rectangular coordinate system at different times; S2, the degrees and the number of time-variable polynomial equations are determined, and LS (least square)-based polynomial equation coefficient estimation of the moving object X, Y and Z direction observation sequences is performed respectively; S3, according to the polynomial equations in the step S2, regression sum of squares and residual sum of squares are respectively calculated for F-test to determine a best fit polynomial; S4, a corrected best fit polynomial is determined; S5, a first time partial derivative and a second time partial derivative are sought respectively to obtain moving object corresponding velocity and acceleration polynomials; and S6, object position state estimator coordinate conversion is performed to obtain corresponding position coordinates, and then by combining with velocity estimated values, object state estimation can be completed.

The present invention discloses an optimal weight combination prediction model of tourism data. The method comprises: firstly, analyzing tourism data influential factors, and classifying the tourism data influential factors according to tourists, destinations and external environments; analyzing the tourism data influential factors that influence tourism revenues and international tourist arrivals by using a grey correlation; sequencing the correlation according to weight values, so as to obtain importance level distribution of the influential factors; then constructing an optimal weight combination prediction model based on a least squares criterion, wherein the model uses the least squares criterion, and an optimal weight value of each model is obtained under a condition of a minimum residual square-sum; and combining a GM(1,1) model and a support vector regression model by using an optimal weight; and finally, obtaining a final prediction value by using the optimal weight combination prediction model. The present invention has the advantageous effects of giving the optimal weight combination prediction model of tourism data and obtaining efficient prediction, so that scenic spot control and tourism-related industry development are supported by efficient decision data, thereby promoting economic benefit improvement in scenic spots and optimization of tourism environment.