66 results about "Functional expression" patented technology

The invention discloses a low-orbit constellation measurement and control resource scheduling scheme based on an improved genetic algorithm, and the scheme comprises the steps: firstly building a model through satellite simulation software STK, obtaining resource data, building a scheduling distribution model, converting measurement and control resources and measurement and control tasks into a mathematical form, and carrying out the measurement and control of the low-orbit constellation, thereby achieving the high-dimensional and multi-constraint problems of the measurement and control of the low-orbit constellation. A profit value is specified for each measurement and control task, an arc section is used as a measurement and control main body according to constraint conditions, a function expression for evaluating the advantages and disadvantages of a scheme is specified, finally, an improved genetic algorithm is applied, improved binary coding is adopted, the measurement and control arc section is used as a gene unit, a tournament selection mechanism is utilized, and data preprocessing and multiple constraint judgment correction are carried out. And solving an optimal scheme. The scheme can efficiently process complex scheduling, reduce invalid search and improve algorithm efficiency and scheme feasibility, meets actual operation requirements, and is suitable for scenes such as commercial satellite measurement and control.

The present application discloses a program compilation method and device, including: obtaining a first expression and a second expression of a program, the first expression and the second expression generating a dependency relationship based on a first variable, and the operation type of the first variable in the first expression is different from the operation type of the first variable in the second expression; updating the operation type of the first variable in the first expression and / or the operation type of the first variable in the second expression, the updated operation type of the first variable indicating that the first variable has no differential result or the result is zero after encountering a differential operation; processing the first expression and the second expression respectively according to the updated operation type of the first variable to obtain corresponding derivative function expressions; combining the derivative function expressions corresponding to the first expression and the second expression respectively based on the chain rule to obtain the compilation result of the program. By implementing the present application, the compilation result of the program can be simplified and the running performance of the compiled program can be improved.

The invention discloses a calculation method for a local axis of a spatial curved surface reticulated shell structural member, and the calculation method comprises the following steps: 1, obtaining a function of a spatial curved surface S, z = S (x, y); 2, for a component on the spatial curved surface reticulated shell, a starting point Pi and an end point Pj of the component are both located on a spatial curved surface S; 3, calculating the normal vector of the spatial curved surface S at the point Pi through a numerical method; and 4, averaging the normal vector of the spatial curved surface S at the starting point Pi position of the component and the normal vector of the spatial curved surface S at the end point Pj position of the component to obtain Vm. The spatial curved surface is an elementary spatial curved surface which can be expressed by an explicit function. A pure mathematical method is used for calculating the local axis of the spatial curved surface reticulated shell structural component, calculation is carried out through spatial analytic geometry, the local axis rotation angle, the local axial vector and the local axis plane of the component can be output, and the method is suitable for various local axis adjusting modes.

The invention provides a physical knowledge constrained intelligent aerodynamic prediction model modeling method. The physical knowledge constrained intelligent aerodynamic prediction model modeling method comprises the following five steps: (1) acquiring aerodynamic physical knowledge through data mining; (2) obtaining a function expression of aerodynamic force; (3) approximately replacing a function expression of aerodynamic force by adopting a depth operator neural network; (4) training a depth operator neural network; and (5) carrying out full-profile aerodynamic modeling based on the trained depth operator neural network.

The application discloses a Zanthoxylum bungeanum protein characteristic data-driven modeling method and system for functional expression, belongs to the technical field of food protein data modeling, obtains protein characteristic indexes, functional expression indexes and impurity interference indexes of Zanthoxylum bungeanum protein samples under different processing conditions, forms corresponding data and carries out standardization and abnormal data processing; based on the correlation degree of the protein characteristic indexes, the functional expression indexes and the impurity interference indexes, functional contribution factors and impurity domination factors are determined respectively, and comprehensive weights are generated; the protein characteristic data is weighted by using the comprehensive weights to obtain function-oriented protein characteristic data; then principal component extraction is carried out, and function-oriented principal components are screened according to the correlation between the principal components and the functional expression data and the impurity interference data; finally, a Zanthoxylum bungeanum protein functional expression prediction model is established; the application can weaken the influence of interference information such as pigments, turbidity, residual oil rate and non-protein particles on modeling, and improve the functional prediction stability.

The present invention discloses a topology optimization method and system for precise microstructure filling. The method includes two aspects: a synthetic microstructure interpolation method based on an initial microstructure and an update method for the unit filling microstructure during the optimization process. The synthetic microstructure based on the initial microstructure interpolates the implicit function expressions of different microstructures to obtain the implicit function expression of the synthetic microstructure, and a corresponding synthetic microstructure geometric model can be established. The update method for the unit filling microstructure during the optimization process combines the distribution of different initial microstructure filling domains within the unit with the synthetic microstructure interpolation method to obtain the microstructure of each unit, and updates are performed in each iteration to clarify the filling microstructure of each unit after each iteration. The geometric configuration of the optimization result of the present invention is clear, without intermediate density units and unclear unit filling microstructures, and has good manufacturability.

The present application is a kind of geotechnical particle contact creep deformation discrete element numerical simulation method, the method is used to simulate long time contact creep between coarse particles, respectively determine four parameters in burgers model of normal coarse particle contact creep and tangential coarse particle contact creep;Obtain the corresponding four parameters under all different normal loads, the corresponding four parameters under all different normal loads and different tangential loads, respectively function fitting with load as independent variable and parameter as dependent variable, respectively obtain the function expression of eight parameters and load, and embed burgers model in PFC3D software;When the input normal load and tangential load are determined, the parameters in burgers model can be dynamically adjusted, and then the particle contact creep displacement under any load condition and time can be obtained. Break through the limitation of laboratory test, realize the quantitative prediction of rheological deformation of coarse aggregate in actual engineering.

The invention provides a multi-port system self-adaptive macro modeling method based on a Roca matrix and vector fitting, and the method comprises the following steps: obtaining S parameters of a multi-port system under each frequency point, and selecting an initial frequency point and a corresponding S parameter matrix; after the frequency is converted to an S domain, initial vector fitting is carried out based on a Loa matrix, and a system initial rational function expression is generated; adopting a dichotomy for interpolation in the frequency band with the maximum change of the rational function, inserting a new frequency point in the midpoint of the frequency band, and carrying out the next vector fitting based on the Roca matrix; and quantizing the sum of the port errors of each frequency band by using Gaussian integral, sorting according to the size, carrying out dichotomy interpolation again in the frequency band with the larger error, and iteratively executing the process until the fitting result reaches the convergence condition, thereby finally obtaining the macro modeling of the multi-port system. On the basis of the interpolation method and the Rocner vector fitting algorithm, the calculation memory requirement is remarkably reduced and the solving efficiency is improved while the solving precision is ensured.

This invention provides an FPGA-accelerated database query method, apparatus, device, and medium, comprising: inputting a query SQL into a row parsing module in the FPGA hardware core to obtain a first data stream; inputting the first data stream into a filtering module to obtain a second data stream; inputting the second data stream into a calculation module to obtain a third data stream; and inputting the third data stream into a projection output module to obtain a target data table. The row parsing module is used to split the data table into column data based on the configuration corresponding to the query SQL and put it into the stream; the filtering module is used to execute the filtering condition clause corresponding to the query SQL on the input data stream; the calculation module is used to execute the function expression calculation corresponding to the query SQL on the input data stream; and the projection output module is used to encapsulate the input data stream into row data. This invention aims to improve the versatility of FPGA-accelerated database SQL queries.

The invention discloses a spatial geometric model display method and system based on advanced mathematics, which can realize full-automatic and high-fidelity conversion from an original function expression to a three-dimensional rotating curved surface by constructing an advanced mathematics symbol analysis method. A dynamic sampling strategy based on a function change rate and a quaternion space rotation mechanism are adopted, and geometric accuracy and calculation stability of complex functions, especially piecewise functions and implicit functions under any rotation axis are guaranteed; real-time rendering and mathematical formula synchronous deduction are combined, an intuitive mapping relation from algebraic expression to spatial form is established, and the understanding depth of students on the rotating curved surface generation principle is remarkably enhanced; interactive operation and a multi-level teaching auxiliary function are introduced, the sense of participation and teaching applicability of the demonstration process are improved, the method is suitable for university mathematics public courses and engineering professional basic teaching scenes, and the structural defects of an existing visualization means in the aspects of semantic comprehension, adaptation flexibility and teaching fusion degree are effectively overcome.

The invention provides a database access method and system supporting unified authentication and business model mapping, and the method comprises the steps: uniformly converting a query request into a standard abstract syntax tree, and intercepting an identity certificate in the query request to construct a complete authority context object; mapping the permission context object into a temporary role name with the minimum required permission, and sending the temporary role name to a key management system to generate a temporary database connection voucher; the permission context object is received, and a specific RLS predicate character string and a CLS function expression are generated based on a preset RBAC strategy rule, an RLS constraint predicate template and a CLS constraint desensitization function template; based on the business model, mapping the standard abstract syntax tree to generate an initial logic query plan; and according to the RLS predicate character string and the CLS function expression, forcibly solidifying the security policy into the initial logic query plan to generate a target logic query plan. The technical problem that authentication and authorization mechanisms are separated in the prior art is solved.

The invention provides a rice heat-resistant candidate gene screening method based on population inheritance and transcriptome integration, which comprises the following steps: acquiring a rice population, performing population differentiation analysis on the rice population, and screening out highly differentiated sites. And carrying out heterozygosity screening on the rice population to obtain a target SNP site. And taking a target SNP site which is obviously differentiated among the subgroups and is homozygous in the subgroups as a candidate SNP site. False positive signals caused by genetic background differences can be effectively eliminated through combined screening of genetic differentiation between subgroups and homozygosity in groups. And performing function annotation association on the candidate SNP sites, and performing transcriptome differential expression analysis on the rice population. And determining a final candidate gene according to the potential regulatory gene and the differential expression gene. Through collaborative analysis of population genetics and transcriptomics, genetic differentiation and functional expression are verified at the same time, the accuracy and reliability of candidate genes can be remarkably improved, and limitation of single-dimension screening is avoided.

This invention proposes an online correction method and system for sine and cosine signals of a rotary transformer, relating to the technical field of servo control. The method involves fitting the demodulated original sine and cosine signals using second-order polynomial functions to obtain second-order sine and cosine polynomial expressions, respectively. The term matrix and the parameter matrix to be optimized are then extracted to obtain the sine and cosine matrices. Based on these matrices, a sine and cosine correction error function is constructed. With the goal of minimizing the sine and cosine correction error function, the method uses the parameter matrix to perform online iterative solutions, thereby completing the online correction of the rotary transformer's sine and cosine signals. This method has low computational complexity and short processing time. Compared to directly decoding uncorrected sine and cosine signals, it improves the linearity of the subsequent decoding angle and reduces decoding speed fluctuations.

The invention provides an assembling mode of RNA liquid drops deformed through photostimulation and application of the RNA liquid drops, and belongs to the technical field of preparation of nucleic acid for analysis. The method comprises the following steps: mixing an RNA single strand with a buffer solution through structural design to obtain a mixed system with a specific concentration; after the mixed system is subjected to heat shock to complete base pairing and assembling, sealing and packaging are performed, the RNA liquid drops deformed through photostimulation are obtained, and the diameter of the liquid drops is 7.0-13.0 microns. The RNA liquid drop obtained in the embodiment can deform under light stimulation, behavior regulation and control can be achieved through group modification, and feasibility is provided for construction and functional expression of artificial cells.

The application discloses a water adding flow calculation method for a moisture regaining machine. The moisture regaining machine is provided with a motor, a rotating shaft and helical blades. The motor is used to drive the rotating shaft and the helical blades to rotate in the moisture regaining machine, so as to transport materials from the inlet of the moisture regaining machine to the outlet of the moisture regaining machine. The calculation method comprises the following steps: obtaining an outlet moisture content target value, an inlet moisture content measured value and an inlet flow measured value of the materials; calculating a water adding flow basic value according to the outlet moisture content target value, the inlet moisture content measured value, the inlet flow measured value and a first functional expression; obtaining a frequency of the motor and a steam pressure value in the moisture regaining machine; calculating a water adding flow adjustment value according to the frequency of the motor, the steam pressure value and a second functional expression; and calculating a total water adding flow according to the water adding flow basic value, the water adding flow adjustment value and a third functional expression. The water adding flow of the moisture regaining machine can be calculated according to the relevant parameters of the materials at the inlet of the moisture regaining machine. The application further discloses a computer readable storage medium.

The invention discloses a discrete heterogeneous multi-agent system grouping consistency control method based on a cooperation-competition relation. Belongs to the technical field of artificial intelligence control. The method comprises the following steps: 1, establishing a heterogeneous multi-agent system model composed of m second-order agents and n first-order agents; 2, establishing a communication topology of the heterogeneous multi-agent system; 3, establishing a cooperation-competition relation functional expression between the intelligent agents; 4, establishing a heterogeneous multi-agent system consistency controller based on a competition-cooperation relationship; and 5, applying the established controller to the established heterogeneous multi-agent system model to realize the consistency control of the heterogeneous multi-agent system based on the competition-cooperation relationship. According to the method, communication between heterogeneous multi-agent systems is carried out based on the competition-cooperation relationship between the agents, and under the condition that packet loss and the interaction strength between the agents are considered, sufficient conditions for ensuring the agent grouping consistency are obtained.

The present disclosure relates to a distributed dynamic performance improvement control method and device. The method comprises: establishing and generating an inductor current prediction model for a parallel Buck converter system model based on droop control; converting the inductor current prediction model into a linear function expression based on external disturbance and model uncertainty analysis; establishing a residual-based performance improvement controller, and solving it based on a gradient descent algorithm to improve distributed dynamic performance control. The present disclosure can address external disturbances and model uncertainty issues in the system, while also being plug-and-play and applicable to a wide range of applications.

The invention relates to a valve side dry-type bushing insulation damp state evaluation method, and belongs to the technical field of power equipment performance evaluation. The method comprises the following steps: S1, carrying out high-voltage impedance spectroscopy test on dry-type bushing samples with different moisture contents to obtain a function expression of frequency and moisture so as to obtain a frequency domain impedance spectrum of the valve-side dry-type bushing; s2, converting the frequency domain impedance spectrum into a multi-time-scale relaxation dynamic spectrum of a time domain by using an improved inverse Laplace transformation algorithm and combining a regularization technology and a machine learning algorithm; s3, extracting multi-time scale relaxation dynamics behavior spectrum characteristics strongly related to the moisture content; and S4, establishing a nonlinear mapping relationship with the moisture content characteristic parameters, and constructing a moisture content evaluation model. According to the method, the applicability and robustness of evaluation are remarkably improved, and technical support is provided for guaranteeing safe and stable operation of the valve side dry type casing pipe.

The invention discloses a dynamic load prediction method and device for a double-layer vibration isolation system under the variable stiffness condition. The method comprises the steps of firstly establishing a multi-degree-of-freedom dynamic model of the double-layer vibration isolation system, then describing structural displacement based on a frequency domain transfer function, further deriving a function expression of a to-be-identified dynamic load about dynamic characteristics of the double-layer vibration isolation system, and finally obtaining a relational expression between the to-be-identified dynamic load before and after stiffness change. When the rigidity of a lower-layer structure in a double-layer vibration isolation system of mechanical equipment changes, the dynamic load borne by the lower-layer structure needs to be accurately predicted so as to guarantee safe and stable operation of the equipment. The method is simple and easy to implement, the dynamic load can be recognized only through the dynamic characteristics of the structure, inverse operation in a traditional dynamic load recognition method is avoided, and the recognition precision is improved. Besides, under the variable stiffness condition, the method can quickly predict the load after the stiffness is changed on the basis of the load identified before the stiffness is changed, and compared with a traditional dynamic load identification method, the identification efficiency is improved.

The invention discloses a distributed N + 1 salary settlement system based on a block chain, and relates to the technical field of computer software. The system comprises a configuration center module for constructing salary item configuration and rule binding data, a salary item management module for generating a structured salary item model, a rule engine module for analyzing rule expression and field dependence, a data adapter module for completing field normalization and standardization processing, a dynamic formula analyzer for generating a rule function set, and a rule function module for generating rule function data. The multi-source data access module constructs a callable field dictionary, the 1 + N salary model module generates aggregation configuration data, the settlement core module executes a distributed calculation task, and the salary report module completes multi-role report generation and chain registration. By constructing an extensible salary structure, field normalization and a dynamic rule calculation mechanism, multi-source data adaptation, rule functional expression and on-chain right confirmation processing of a distributed settlement result are achieved, and the configuration flexibility, calculation concurrency and data credibility of a salary system are improved.

The application discloses a processing method, device and equipment of a flicker pulse and a storage medium. The processing method comprises the following steps: presetting a plurality of reference threshold values for a working device used for multi-threshold sampling, and acquiring a threshold correction function model; presetting a plurality of reference working temperatures, and performing multi-threshold sampling on an auxiliary correction pulse at each reference working temperature to determine a plurality of reference pulse widths corresponding to the plurality of reference threshold values at the reference working temperatures; performing contrast sampling on the auxiliary correction pulse to determine a conversion relationship between the threshold value and the pulse width of the auxiliary correction pulse; determining a plurality of actual working threshold values related to the plurality of reference pulse widths at each reference working temperature based on the conversion relationship; determining a function expression of the threshold correction function model; and setting a plurality of actual preset threshold values for the working device based on the function expression. The application is especially suitable for threshold correction of a threshold comparison working device in a wide temperature range.

The application provides a hyperspectral image classification method based on function expression dimension reduction; first, the pixel spectral measure of the hyperspectral image is sequentially accumulated to construct an accumulated spectral curve with a monotone increasing characteristic; the accumulated spectral curve is fitted by a polynomial to obtain an optimal function expression; finally, the polynomial coefficients are taken as features, and SVM is used to realize hyperspectral image classification; in order to verify the effectiveness of the application, three hyperspectral image data are respectively experimented, and the application is compared with three related comparative algorithms for comparative analysis, and the results show that the overall accuracy of the application on the three hyperspectral images is 97.83%, 96.22% and 98.15% respectively, the overall accuracy is high, it is proved that the polynomial function expression coefficients are taken as features, the purpose of hyperspectral image feature dimension reduction can be achieved, the class separability is high, and high classification accuracy can be achieved.

The present invention provides a method for obtaining rock mass structural information based on rock mass surface and borehole scanning, comprising: collecting an original digital image of the exposed surface of the rock mass excavation; collecting an original digital image of the inner wall of the borehole and recording the borehole position information; processing the original digital image to obtain a binary image containing only crack patterns; obtaining the coordinates of all crack patterns in the binary image, and obtaining a function expression of the crack patterns in the binary image of the inner wall of the borehole by fitting; analyzing and calculating to obtain a crack information set in the rock mass at the engineering site; inputting the crack information set into an information induction model for classification and organization to obtain a structural surface information set; constructing a three-dimensional model of the structural surface; and building a rock mass structural model to obtain rock mass structural information. The present invention makes the process of obtaining structural surface parameters more reliable and complete, and provides a more efficient, convenient and reliable method for obtaining rock mass structural information in actual engineering projects, which has a good guiding significance for engineering sites.

The invention discloses a design method for a non-linear variable air damper regulating device of a vehicle-mounted air conditioner, which comprises the following steps: First, establish a cross coordinate system, and take the rotation center of the track disc as the zero point of the track for controlling the swing of the control rocker; Then, establish the relational expression of the track in the cross coordinate system; Finally, use the "LAW" command in the CATIA software to convert the relational expression of the track into a relational expression about the abscissa and ordinate of any point on the track with α n , that is, the included angle between the track-driven rocker and the abscissa of the coordinate system as the independent variable; Then use the "Space Blend" command in the surface module of the CATIA software to obtain the function curves of the abscissa and ordinate of any point on the track with α n as the independent variable, that is, the track according to α n function expression. The invention adopts a computer-aided design method to design the track disc for driving in a vehicle-mounted air conditioner with non-linear heating, and can design the non-linear regulating device quickly and pertinently.

The present invention discloses a method and system for detecting Ponzi scheme source code in smart contracts based on a large model. The method first compiles the contract to obtain an abstract syntax tree (AST) of the smart contract code. The AST is then parsed to generate an intermediate representation (IR) of the contract code and a contract object containing the contract's structured information. The IR of the smart contract code is used to model the smart contract as a hypernode graph. A mapping is established between each node and the contract source code. The funds sent to the contract and the sender's address are used as taint sources to perform static taint analysis on the contract and track the flow of incoming funds within the contract. A taint variable propagation graph is generated. The vertices in the IR graph are then divided into three levels: contracts, functions, and expressions. Based on the smart contract IR graph and its subgraphs, the code corresponding to the function nodes involved in tainted variable propagation is sliced. Finally, the IR is used to detect Ponzi schemes. The present invention achieves rapid and accurate Ponzi scheme contract detection.

The invention discloses a function expression method for milling an aviation aluminum alloy stress wave in a real number field, relates to the technical field of aviation aluminum alloy milling, and aims to solve the problems that an existing three-dimensional stress wave model is complex in solution, time-consuming and difficult to guarantee in accuracy. The method comprises the following steps: firstly, judging the contact condition of an end mill and a workpiece, and calculating the instantaneous cutting layer thickness of a milling infinitesimal by utilizing a cutter tooth instantaneous contact relationship; on the basis of a one-dimensional stress wave theory, the size of a milling stress wave excitation source is calculated by combining the obtained thickness of the cutting layer; and finally, constructing a longitudinal wave fluctuation equation by using a Lame-Navier equation, and solving a partial differential equation by using a Laplace transformation method to obtain a function expression of stress waves in a real number field. According to the method, the stress wave solving process is simplified, the calculation precision and efficiency are improved, and a basic calculation method is provided for high-stability machining of the milling aviation aluminum alloy and stress wave related research.

The present invention discloses a rapid search method and system for three-component second-order threshold implementation solutions, comprising: determining a component coordinate combination based on incompleteness in second-order safety threshold implementation characteristics; determining initial component function expressions of a vector Boolean function based on the component coordinate combination and the correctness of the safety threshold implementation; determining candidate linear additions for each component function based on the component function expressions and the combination probability identical distribution and uniformity in the safety threshold implementation characteristics; placing the component functions corresponding to the coordinate functions that satisfy the probability identical distribution and uniformity and to which the candidate linear additions are added in the memory unit space of the coordinate function; extracting a component function from the memory unit space corresponding to each coordinate function, combining all the extracted component functions into a whole, and taking the whole that satisfies the probability identical distribution, uniformity, and correctness as the final threshold implementation solution. Rapid search for threshold implementation solutions is achieved.

The invention provides a quick modeling method for a simple geometry of a ceramic matrix composite, which comprises the following steps of: acquiring real braided structure data of an actual simple geometry in a CT (Computed Tomography) scanning mode, and extracting characteristic parameters of the real braided structure data; inputting the characteristic parameters and preset modeling auxiliary parameters to a Matlab end, performing function expression calculation of the number of cycles and a sweeping curve to generate discrete point coordinates of a yarn sweeping line, and storing the discrete point coordinates to a control file; and on the basis of the discrete point coordinates, constructing a yarn entity array of the flat plate geometry in a UG GRIP program, outputting a flat plate weaving entity, carrying out modeling on the revolving body geometry, obtaining a yarn entity array of a revolving body weaving structure, and outputting a revolving body simple geometry weaving entity. According to the method, the matlab and UG secondary development program is adopted, rapid modeling of the simple geometry is achieved, the modeling time is saved, program parameterization is achieved, the size of parameters can be changed at will, and the practicability and applicability of modeling are greatly improved.

The invention discloses a calculation method of an integral transformation fracture phase field model of an exponential cohesion law, and belongs to the technical field of finite element analysis of phase field models. The method comprises the following steps: deducing a cohesion eigenfunction of an index cohesion law by using a subprogram, deducing an analytical expression of an index energy density eigenfunction, and solving the expression of the index energy density eigenfunction and first two-order derivatives thereof by adopting multiple integral transformations, the Maclaurin series is adopted to determine an index energy density eigenfunction expression and an expansion expression of the first two-order derivatives at the zero point of the index energy density eigenfunction expression; and inputting the index energy density eigenfunction expression, the first two-order derivatives of the index energy density eigenfunction expression, and the expansion expression of the index energy density eigenfunction expression and the first two-order derivatives of the index energy density eigenfunction expression at the zero point into a solver for performing an integral transformation fracture phase field model to perform phase field model simulation calculation. According to the method, the simulation calculation of the exponential integral transformation fracture phase field model can be realized under a universal finite element framework, and the accuracy and the calculation efficiency are greatly improved.

The invention discloses a non-circular gear pair pitch curve design method without a function expression, and the method comprises the following steps: S1, drawing a driving wheel rotation angle-radial curve; s2, extracting point set coordinates; s3, setting a center distance; s4, calculating the rotating angle of the driven wheel; s5, driven wheel pitch curve closure judgment; s6, a driving wheel pitch curve is constructed through spiral lines; s7, a planet carrier is constructed; s8, creating a pitch curve contact point tracking circle; s9, tracking a pitch curve contact point; and S10, connecting the pitch curve contact points in series to form a driven wheel pitch curve. According to the method, any spline curve is drawn to serve as a driving wheel rotation angle-radial curve, data calculation is conducted after point set coordinates are extracted, the driving wheel pitch curve and the driven wheel pitch curve are obtained through spiral line fitting and motion simulation means, complex pitch curve function expression derivation is avoided, design efficiency is high, and the method is easy to implement. And the whole design process is more free, and the method has a good application prospect.