CN109067264B - Balance car system and control method thereof - Google Patents

Balance car system and control method thereof Download PDF

Info

Publication number
CN109067264B
CN109067264B CN201811083745.3A CN201811083745A CN109067264B CN 109067264 B CN109067264 B CN 109067264B CN 201811083745 A CN201811083745 A CN 201811083745A CN 109067264 B CN109067264 B CN 109067264B
Authority
CN
China
Prior art keywords
particle
value
balance car
follows
control
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811083745.3A
Other languages
Chinese (zh)
Other versions
CN109067264A (en
Inventor
林海
翟晋平
王正来
盛丹洁
张琳虎
曲正
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HUIJIAWANG (TIANJIN) TECHNOLOGY Co.,Ltd.
Original Assignee
Changan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changan University filed Critical Changan University
Priority to CN201811083745.3A priority Critical patent/CN109067264B/en
Publication of CN109067264A publication Critical patent/CN109067264A/en
Application granted granted Critical
Publication of CN109067264B publication Critical patent/CN109067264B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/04Arrangements for controlling or regulating the speed or torque of more than one motor
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B62LAND VEHICLES FOR TRAVELLING OTHERWISE THAN ON RAILS
    • B62KCYCLES; CYCLE FRAMES; CYCLE STEERING DEVICES; RIDER-OPERATED TERMINAL CONTROLS SPECIALLY ADAPTED FOR CYCLES; CYCLE AXLE SUSPENSIONS; CYCLE SIDE-CARS, FORECARS, OR THE LIKE
    • B62K11/00Motorcycles, engine-assisted cycles or motor scooters with one or two wheels
    • B62K11/007Automatic balancing machines with single main ground engaging wheel or coaxial wheels supporting a rider
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/006Controlling linear motors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/10Arrangements for controlling torque ripple, e.g. providing reduced torque ripple
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/34Modelling or simulation for control purposes

Abstract

The invention discloses a balance car system and a control method thereof.A state equation of extended Kalman filtering and a mathematical model of a permanent magnet linear motor are established according to a difference equation model and a drift error model of an attitude signal of the balance car system; driving the permanent magnet synchronous motor by using a second-order dynamic Terminal sliding film; and simultaneously optimizing parameters of a tracking controller, a synchronous controller, an anti-backlash controller and a switching function in the dual-motor system by using a particle swarm algorithm. The invention realizes that the signals acquired by the sensor of the balance car have larger noise and zero drift error; the robustness is poor; the brush motor has the advantages of complex structure, more faults, large maintenance workload, short service life, easy generation of electromagnetic interference by reversing sparks, and the problems of jitter, overlarge deviation, low safety and the like in the running process of the balance car.

Description

Balance car system and control method thereof
Technical Field
The invention belongs to the technical field of balance car control, and particularly relates to a balance car system and a control method thereof.
Background
The two-wheel balance car gradually enters the life of people due to the characteristics of low price, small volume, light weight, convenience in moving and carrying and the like. At present, the related technology of the balance car is continuously updated, the control of a balance system is more flexible and sensitive, and a human-computer interaction system tends to be delicate and humanized. However, in the driving mode, most of the balance vehicle brands in the current market are direct current brush motors, and electric energy enters an armature winding through a brush and a commutator to drive the motor. Due to the existence of the electric brush and the commutator, the brush motor has the characteristics of complex structure, more faults, large maintenance workload, short service life, easy generation of electromagnetic interference by commutation sparks and the like, and the production cost and the maintenance cost are increased to a certain extent. The traditional balance car has the defects of weak endurance, low precision of a control system, single function and the like.
Although balance cars on the market are various in product types, most of the product motors are insufficient in control design, and the traditional control method is easy to cause balance accidents due to the defects of the traditional control method. PID control, based on T _ S fuzzy control, self-adaptive control and the like have the following defects: the control variable armature voltage indirectly controls the balance of the vehicle body, so that the dynamic performance of the electric vehicle is difficult to ensure; the signal collected by the sensor has larger noise and zero drift error; the robustness is poor, and the method is not suitable for the application of a nonlinear system. The sliding mode control design is simple, the control precision is high, the sliding mode has complete robustness to perturbation and external disturbance of a system, and the method is an important method for solving the nonlinear problem. However, in general sliding mode control, a linear switching function is generally selected, and although system stability can be ensured, in theory, infinite time is required for a system state (or error) to converge to a balance origin. The Terminal sliding mode constructs a switching surface by a nonlinear function and can converge to a balance point within a limited time, which is inspired by an attractor concept in a neural network. But the problem of buffeting exists, and the application and the control effect of the sliding mode control in engineering are seriously influenced by the buffeting. Although effective in removing buffeting, methods such as saturation functions may degrade control performance.
In summary, the existing balance car driving control method does not comprehensively consider the tracking and monitoring of the controller on various performances of the permanent magnet synchronous motor, and the problems of jitter, overlarge deviation, low safety and the like often occur in the implementation process of the method. Therefore, the novel balance car control method has great practical engineering value.
Disclosure of Invention
The invention aims to solve the technical problem of providing a balance car system and a control method thereof aiming at the defects in the prior art, and the balance car system converges at a balance origin within a limited time on the basis of a permanent magnet synchronous motor based on a second-order dynamic Terminal sliding mode of Kalman filtering and a particle swarm algorithm, so that buffeting is eliminated, and the stability and durability of the balance car are improved.
The invention adopts the following technical scheme:
a balance car system comprises a control part, a sensor part, a man-machine interaction part, a power supply part, a main circuit part and an upper computer part; the control part adopts a microcontroller and is respectively connected with the sensor part, the man-machine interaction part, the power supply part and the main circuit part; the upper computer part is connected with the human-computer interaction part, and the human-computer interaction part is powered by the power supply part;
the sensor part comprises a gyroscope for acquiring the angular speed of the offset of the vehicle body, an accelerometer for acquiring the acceleration of the vehicle, a speed sensor for acquiring the actual speed, a corner sensor for acquiring the direction angle of the joystick, an infrared sensor for detecting the road condition, and a voltage and current sensor for acquiring related data of the inverter;
the man-machine interaction part comprises an L CD display screen used for displaying vehicle speed information, electric quantity information, battery temperature information and system running state feedback, a wireless Bluetooth used for transmitting data to a mobile phone APP and sending instructions to a controller, an upper computer communication module used for communicating with a PC, keys used for power switch control, mode selection and information inquiry, a L ED lamp and a buzzer used for prompting vehicle alarm information, and a power supply module used for supplying power to the mobile phone APP;
the power supply part comprises a charging module and a storage battery, wherein the charging module is used for charging batteries and supplying power to all parts of the system;
the main circuit part comprises two permanent magnet synchronous motors and a three-phase arm inverter capable of driving the two permanent magnet synchronous motors, and is used for realizing independent control and cooperative operation of the motors;
the upper computer part comprises a mobile phone for transmitting information and transmitting instructions through wireless Bluetooth and a PC (personal computer) end for communicating through a serial port.
A balance car system control method comprises the steps of establishing a state equation of extended Kalman filtering and a mathematical model of a permanent magnet linear motor according to a difference equation model and a drift error model of an attitude signal of a balance car system; driving the permanent magnet synchronous motor by using a second-order dynamic Terminal sliding film; and simultaneously optimizing parameters of a tracking controller, a synchronous controller, an anti-backlash controller and a switching function in the dual-motor system by using a particle swarm algorithm.
Specifically, the state equation of the extended kalman filter is as follows:
Figure BDA0001802586650000031
wherein the content of the first and second substances,
Figure BDA0001802586650000032
is the output value of the gyroscope,
Figure BDA0001802586650000033
the actual angular velocity measured by the gyroscope is represented by k, scale error and drift error, v (n-1) is white Gaussian noise with the mean value of zero, and T is a period.
Further, the difference equation model of the attitude signal of the balance car is as follows:
Figure BDA0001802586650000034
wherein, ω is0For initial angle values of the gyroscope, κ is the scale error, α is the drift error, and the first order autoregressive AR model for the gyroscope drift error is as follows:
Figure BDA0001802586650000035
wherein the content of the first and second substances,
Figure BDA0001802586650000036
is an autoregressive model parameter; v (n) is measurement white noise with a mean value of zero;
the variance of gaussian white noise is as follows:
Figure BDA0001802586650000037
Figure BDA0001802586650000041
wherein, Ci、CgInclinometer and gyroscope noise covariance, respectively;1is the standard deviation of the tilt angle;2is the gyroscope gaussian noise density standard deviation;3is the standard deviation of the noise in the gyro drift error AR model.
Specifically, the mechanical motion equation of the permanent magnet linear motor is as follows:
Figure BDA0001802586650000042
wherein s is mover displacement; v is the mover speed; m is the total mass of the rotor and the load carried by the rotor; b isvIs the viscous friction factor; fIs the total disturbance force.
Further, the total disturbance force FThe following were used:
F=Frip+Fload+Ffric
Figure BDA0001802586650000043
Figure BDA0001802586650000044
wherein, FloadIs the load resistance; fripEquivalent resistance to the end effect; fMThe end effect thrust fluctuation amplitude; tau is a polar distance;
Figure BDA0001802586650000045
is the initial phase electrical angle; ffricIs a friction force; f. ofcIs the coulomb friction coefficient; f. ofsIs the static friction coefficient; v is the mover speed; v. ofsIs the critical friction speed.
Specifically, driving the permanent magnet synchronous motor by using a second-order dynamic Terminal slip film specifically comprises the following steps:
let x be [ x ]1,x2]T=[s,v]TIs a change of state of the systemQuantity, input control quantity u ═ iqThe equation of state is as follows:
Figure BDA0001802586650000046
Figure BDA0001802586650000051
wherein k is1,k2,k3For unknown parameters, kfIs a coefficient of resistance, BvIs viscous friction factor, FM is total disturbance force, and M is total mass of the rotor and the load carried by the rotor;
let the system tracking error be e ═ x1 *-x1Wherein x is1 *Is x1Defining the second-order nonsingular fast terminal sliding mode variable of the system as follows:
Figure BDA0001802586650000052
wherein, 0<α<1,β∈R+P, q ∈ N are odd numbers, lambda is more than p/q, and p/q is more than 1 and less than 2;
when the system state is close to the balance point, the high-order term of the tracking error e (t) approaches to 0, and the convergence speed of the tracking error e (t) is approximate to that of a nonsingular terminal sliding mode; when the system is far away from a balance point, the high-order terms of the tracking error e (t) play a main role, and the convergence speed of the tracking error e (t) is higher than that of a nonsingular terminal sliding mode; the sliding mode control system meets the sliding mode variable sigma and the first derivative thereof
Figure BDA0001802586650000053
Converging to zero.
Further, the first derivative of the synovial variable
Figure BDA0001802586650000054
The following were used:
Figure BDA0001802586650000055
the synovial membrane control law is as follows:
Figure BDA0001802586650000056
wherein: r is a position command, ce is an error, f (t) is a function designed to achieve a global sliding mode surface, f (t) is s (0) exp (- λ t), λ >0, s (0) is an initial time s (t), sgn(s) is a step function, and B is a friction coefficient; k (t) is a switching gain, k (t) max (| e (t) |) + ρ, ρ > 0; when the system input ss <0, the synovium is present; if the system input ss >0, the switching gain K (t) should be increased; if the system input ss <0, the switching gain K (t) should be reduced.
Specifically, a performance index function F of a parameter set by using a particle swarm algorithm is as follows:
Figure BDA0001802586650000061
wherein, the smaller the first term of the performance index is, the smaller the steady-state error of the system is; the smaller the second term is, the smaller the energy consumption of the system is; t is time, e1(t) is a constraint term, u1,u2Is a comprehensive control law.
Further, the method comprises the following specific steps:
s401, initializing each parameter value of a particle swarm, calculating a fitness function of each particle, and calculating the performance index size of each group of parameters in a weighting mode;
updating the adaptation value with the new particle if the adaptation value of the new particle is smaller than that of the previous one; otherwise, the adaptation value remains unchanged.
Figure BDA0001802586650000062
Where pbest (t) is the best adaptation value at time t, f (-) is the objective function of the performance indicator, and X (t) is the position of each particle; pbest records the optimal solution searched by an individual, and gbest is used for recording the optimal solution searched by the whole group in one iteration;
the updated equations for velocity and particle position are as follows:
V[i]=w×v[i]+c1×rand()×(pbest[i]-present[i])+c2×rand()×(gbest-present[i])
wherein v [ i ] represents the velocity of the ith particle, w represents an inertia weight, c1 and c2 represent learning parameters, rand () represents a random number between 0 and 1, pbest [ i ] represents the optimal value searched by the ith particle, gbest represents the optimal value searched by the whole cluster, and present [ i ] represents the current position of the ith particle;
s402, when the minimum adaptive value in pbest is smaller than the global adaptive value, updating the global adaptive value by using the position of the corresponding minimum adaptive value; otherwise, the global adaptation value remains unchanged;
gbest(t+1)=argmin{f(pbest1(t)),f(pbest2(t)),...,f(pbestn(t))}
wherein gbest (t) is the global optimal adaptive value at the time t, and n is the total number of particles;
s403, updating the control parameter value according to the following formula
Figure BDA0001802586650000071
Xi'j(t+1)=Xi'j(t)+Vi'j(t+1)
Wherein Vi' j (t) is the speed of the ith population iteration particle of the jth dimension particle; xi ' j (t) is the position of the ith ' group iteration particle of the jth dimension particle, and i ' is taken as 20; omega is inertia weight, and is taken as 0.7; c. C1And c2To obtain the learning rate, take c1=c2=2;
S404, the current position of the particle exceeds the set maximum value and the set minimum value, and the particles beyond the range are re-assigned, namely
Figure BDA0001802586650000072
When the maximum velocity of the particle is exceeded, the velocity of the particle is reassigned to
Figure BDA0001802586650000073
Wherein, xmin (j) and xmax (j) are respectively the minimum position and the maximum position of the j dimension; vmin (j) and Vmax (j) are respectively the minimum speed and the maximum speed of the j dimension;
s405, when the iteration times are smaller than the maximum set times, turning to the step S403; otherwise, ending.
Compared with the prior art, the invention has at least the following beneficial effects:
the invention relates to a balance car system which comprises a control part, a sensor part, a man-machine interaction part, a power supply part, a main circuit part and an upper computer part, wherein a permanent magnet synchronous motor is adopted as a motor. The design of the five-bridge arm inverter improves the dynamic performance of the two-motor system, reduces torque fluctuation and can effectively reduce the energy loss rate. Under the condition of the same battery capacity, the endurance of a balance car system driven by the five-bridge-arm inverter is obviously stronger.
The invention also discloses a balance car system control method, which optimizes the parameters of the designed controller by utilizing the particle swarm algorithm, reduces the energy consumption of the system on the premise of ensuring the tracking performance, needs to minimize the value of the controlled variable to achieve the energy-saving effect, obtains the unknown parameters in the control law and the optimal value of β in the switching function by adopting the particle swarm algorithm to simultaneously meet the aims of tracking precision and energy saving, combines the second-order dynamic sliding mode with the Terminal sliding mode, effectively removes the sliding mode buffeting while quickly and highly accurately tracking the command signals, and simulates to show the effectiveness of the method.
Further, using kalman filtering, well-founded predictions are made of the next step of the system, which always indicates what is actually happening, even with various disturbances. The use of kalman filtering is highly desirable in continuously varying systems, has the advantage of small memory usage, (no need to retain historical data other than the previous state quantity) and is fast and well suited for real-time problem and embedded systems.
Furthermore, the permanent magnet synchronous linear motor has the advantages of no need of electromagnetism, high thrust density, high efficiency and the like. The simulation method overcomes the defects of the prior art, provides the running state of the linear motor, realizes the simulation of the permanent magnet linear synchronous motor, and provides great help for the research of the running state of the permanent magnet linear synchronous motor.
Further, when the second-order dynamic Terminal sliding mode traditional sliding mode is designed, in order to keep the system motion on the sliding mode surface, the system input is switched back and forth between different controllers, so that the system has buffeting. High order sliding modes and dynamic sliding modes can overcome this problem well. The high-order sliding mode selection controller not only enables the switching function to be equal to zero, but also enables the first-order derivative of the switching function to be equal to zero until the r-1-order derivative, and the control item only acts on the r-order derivative of the switching function, which is called r-order sliding mode control. The dynamic sliding mode embeds the discontinuous term into the derivative term of the controller.
Furthermore, the method for setting the parameters of the PID controller by using the particle swarm optimization algorithm gives specific design steps for each parameter and an adaptive function, selects a permanent magnet type direct current servo motor commonly used in control engineering design as a research object, and has certain practical value. Compared with the genetic algorithm, the simulation result has certain improvement. The optimal parameters obtained by the simulation result under different adaptive functions can be applied to actual servo motor speed control.
In conclusion, the invention realizes that the signals acquired by the sensor of the balance car have larger noise and zero drift error; the robustness is poor; the brush motor has the advantages of complex structure, more faults, large maintenance workload, short service life, easy generation of electromagnetic interference by reversing sparks, and the problems of jitter, overlarge deviation, low safety and the like in the running process of the balance car.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
FIG. 1 is a general block diagram of a novel balance car system;
FIG. 2 is a functional block diagram of the present invention;
FIG. 3 is a control flow diagram of the present invention;
FIG. 4 is a block diagram of a particle swarm algorithm of the present invention.
Detailed Description
Referring to fig. 1, the invention provides a balance car system, which adopts a second-order dynamic Terminal sliding mode and a particle swarm algorithm based on kalman filtering, and comprises a control part, a sensor part, a man-machine interaction part, a power supply part, a main circuit part and an upper computer part; the control part adopts a microcontroller and is respectively connected with the sensor part, the man-machine interaction part, the power supply part and the main circuit part; the upper computer part is connected with the human-computer interaction part, and the human-computer interaction part is powered by the power supply part.
The control part: STM32F103ZET6 control of selection ST company is as main control chip, and this chip has two senior timers and 512K's storage capacity, has 5 kinds of IO mouths, a plurality of serial ports communication interface, fine satisfying the design requirement.
A sensor portion: the gyroscope, the external sensor, the current sensor, the voltage sensor, the rotation angle sensor and the speed sensor are arranged. The gyroscope is used for acquiring angular velocity of vehicle body deviation, the accelerometer is used for acquiring acceleration of a vehicle, the speed sensor is used for acquiring actual speed, the corner sensor is used for acquiring a direction angle of the operating lever, the infrared sensor is used for detecting road conditions, and the voltage and current sensors are used for acquiring relevant data of the inverter.
The man-machine interaction part comprises an L CD display screen, a wireless Bluetooth, an upper computer communication module, a key, a L ED lamp and a buzzer.
The L CD display screen can display vehicle speed information, electric quantity information, battery temperature information and system running state feedback, the wireless Bluetooth is mainly used for transmitting data to the mobile phone APP and sending instructions to the controller, the upper computer communication module is used for communicating with the PC, the keys are used for power switch control, mode selection and information inquiry, and the L ED lamp and the buzzer are used for prompting vehicle alarm information.
A power supply section: the charging device comprises a charging module and a storage battery.
For battery charging and for powering the various parts of the system.
A main circuit part: there are power driver and two permanent magnet synchronous motors.
The power driver is a five-bridge arm inverter, can simultaneously drive two motors, and can realize independent control and cooperative operation of the motors.
The upper computer part: the system is divided into a mobile phone end and a PC end.
The mobile phone mainly transmits information and transmits instructions through wireless Bluetooth; the PC terminal is communicated through a serial port and is mainly used for programming and debugging a system.
The function module comprises a rotor interval calculation module, a speed calculation module, an angle subdivision module, a phase adjustment module, a speed control module, a current limiting control module, a waveform modulation module, a three-phase inverter and a motor starting module. As shown in fig. 2.
Referring to fig. 3, the method for controlling a balance car according to the present invention optimizes parameters of a designed controller by using a particle swarm algorithm, reduces energy consumption of a system on the premise of ensuring tracking performance, minimizes a value of a control quantity ui for achieving an energy saving effect, obtains optimal values of unknown parameters k1, k2, k3, k4 in a control law and β in a switching function for simultaneously satisfying tracking accuracy and energy saving, and specifically integrates tracking accuracy and energy consumption into a target function E, wherein the final objective of optimization is to obtain the minimum value of the target function E under the condition of satisfying constraint conditions.
The balance car works as follows:
the method comprises the following steps: starting, initializing the system;
step two: information interaction is carried out through a wireless communication Bluetooth module;
step three: acquiring a video shot by a vehicle bottom camera, and compressing, encoding and storing the video signal;
step four: keyboard key scanning;
step five: scanning obstacles; acquiring obstacle avoidance sensor signals and calculating a route through an obstacle avoidance algorithm;
step six: road condition scanning;
(1) the infrared sensor is used for detecting the distance from the vehicle bottom to the ground to acquire the ground flatness information, the road condition parameters are calculated through a road condition detection algorithm, and the road condition of the road which is running at present is acquired through further analyzing data, so that the road condition is judged.
(2) The method comprises the steps that the angle and the angular speed of a vehicle body are obtained through an accelerometer and a gyroscope sensor, the angle and the angular speed are input into a Kalman filter, the second-order dynamic Terminal sliding mode is used for calculating, and a group of output signals are output to a motor to control the adjustment of the angle of the vehicle body, so that the vehicle body can be stably kept upright;
(3) acquiring a signal of a speed sensor, acquiring an actual speed, and completing closed-loop control on the speed through a speed PID;
(4) sampling a voltage signal of the operating rod, comparing the voltage signal with a reference voltage to judge the direction and the angle of steering, generating a corresponding motor control signal to a control part, and controlling the left motor and the right motor to run at different speeds to complete steering;
(5) detecting the position of an obstacle by an obstacle avoidance sensor and acquiring the surrounding situation of the obstacle by an infrared sensor, firstly selecting all positions without the obstacle within a certain range, calculating and comparing routes corresponding to the positions according to the detected surrounding characteristics of the obstacle, and finding out an obstacle avoidance route with the shortest bypassing radius so as to avoid the obstacle;
(6) collecting motor voltage and winding current through a voltage and current sensor, calculating torque and flux linkage by using the motor winding current and voltage, generating PWM (pulse width modulation) of required 12 paths through a torque and flux linkage regulator and an existing and disclosed vector table, performing direct torque control, and outputting the PWM to a direct current brushless motor;
controlling a buzzer and L ED, displaying on a display screen, judging whether the operation is finished, if so, entering the step eight, if not, returning to the step two;
step eight: and (6) ending.
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of the embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The control method of the balance car comprises the following steps:
s1 Kalman signal fusion filtering
In the sensor system, the inclination angle and speed of the balance car body are measured by a gyroscope, the measurement result has errors which are mainly drift error and scale error, and the output characteristic is
Figure BDA0001802586650000121
In the formula:
Figure BDA0001802586650000122
is the output value of the gyroscope,
Figure BDA0001802586650000123
for the true angular velocity of the gyroscope measurements, κ is the scale error and α is the drift error.
If the equation (1) is subjected to differential operation, the model of the differential equation of the attitude signal of the balance car in one period T is
Figure BDA0001802586650000124
ω0For initial angle values of the gyroscope, a first-order Autoregressive (AR) model of the drift error of the gyroscope can be established according to an AR estimation method of the drift error of the gyroscope:
Figure BDA0001802586650000125
in the formula:
Figure BDA0001802586650000126
the values are 0.9185 for the autoregressive model parameters; v (n) is measurement white noise with a mean value of zero.
Establishing a state equation of the extended Kalman filtering according to a difference equation model and a drift error model of the attitude signal of the balance car:
Figure BDA0001802586650000131
in the formula: v (n-1) is white Gaussian noise with a mean value of zero.
Variance of Gaussian white noise of
Figure BDA0001802586650000132
Figure BDA0001802586650000133
In the formula: ci、CgInclinometer and gyroscope noise covariance, respectively;1is the standard deviation of the tilt angle;2is the gyroscope gaussian noise density standard deviation;3is the standard deviation of the noise in the model of the gyroscope drift error AR (1).
S2 mathematical model of permanent magnet linear motor
The electromagnetic thrust expression of the permanent magnet linear motor is
Figure BDA0001802586650000134
The mechanical equation of motion for PM L SM is:
Figure BDA0001802586650000135
wherein s is mover displacement; v is the mover speed; m is the total mass of the rotor and the load carried by the rotor; b isvIs the viscous friction factor; fIs the total disturbance force F=Frip+Fload+Ffric;FloadIs the load resistance; fripEquivalent resistance to the end effect;
Figure BDA0001802586650000141
wherein, FMThe end effect thrust fluctuation amplitude; tau is a polar distance;
Figure BDA0001802586650000142
is the initial phase electrical angle; ffricIs a friction force;
Figure BDA0001802586650000143
fcis the coulomb friction coefficient; f. ofsIs the static friction coefficient; v is the mover speed; v. ofsIs the critical friction speed.
S3, second-order dynamic Terminal sliding mode control
The permanent magnet synchronous motor is used as a control object, a magnetic circuit is not saturated, a space magnetic field is in sinusoidal distribution, hysteresis and eddy current loss influence are not counted, and a voltage equation is
Figure BDA0001802586650000144
The torque equation is
Tm=p(ψd*iqq*id) (10)
The mechanical equation of motion is
Figure BDA0001802586650000145
In the formula ud、uqIs the stator voltage dq axis component; i.e. id、iqIs stator current component, omega is rotor angular velocity Ld、LqEquivalent inductance of a dq axis of a stator winding; rsIs a stator resistor; psifFor each pair of magnetic pole flux linkage; t isLIs the load torque; j is moment of inertia; b is a friction coefficient; p is the rotor stage logarithm.
Let the speed given signal be ω, define the error state emThe rotation speed error system is obtained as
Figure BDA0001802586650000146
According to the properties of the high-order slip film, the system buffeting-free slip film control is realized through the slip film with two or more stages.
For convenience of description, let x ═ x1,x2]T=[s,v]TAs a state variable of the system, the input control quantity is u-iq
The equation of state for the system (3) is then in the form:
Figure BDA0001802586650000151
wherein the content of the first and second substances,
Figure BDA0001802586650000152
let the system tracking error be e ═ x1 *-x1Wherein x is1 *Is x1Given values of (a).
Defining a second-order nonsingular fast terminal sliding mode variable of the system as follows:
Figure BDA0001802586650000153
wherein, 0<α<1,β∈R+P, q ∈ N are odd numbers, λ > p/q, 1 < p/q <2 to satisfy the nonsingularity of the slide film surface.
According to the sliding mode surface, when the system state is close to a balance point, the high-order term of the tracking error e (t) approaches to 0 and can be ignored, and the convergence speed of the tracking error e (t) is similar to that of a nonsingular terminal sliding mode; when the system is far away from the balance point, the high-order terms of the tracking error e (t) play a main role, so the convergence speed of the system is faster than that of a nonsingular terminal sliding mode.
The sliding mode variable sigma and the first derivative thereof must be satisfied according to the design requirement of the sliding mode control system
Figure BDA0001802586650000154
Converging to zero.
First derivative of synovial variable
Figure BDA0001802586650000155
Is composed of
Figure BDA0001802586650000156
Assuming that the sliding mode variable σ designed in equation (14) is the output variable of the uncertain system, and the relative order of the second-order sliding mode control system is 1 as known from equation (15), u is selected to make the output variable σ and its derivative be appropriate
Figure BDA0001802586650000161
Convergence to zero within a finite time.
The equation of state for the sliding-mode variable of equation (15) can be written as
Figure BDA0001802586650000162
Wherein, r is the position instruction, through fuzzy synovial membrane algorithm, sets up error e to: e-r- θ.
The overall dynamic slide film surface is additionally arranged as follows:
s=e+ce-F(t) (17)
where f (t) is a function for achieving the global sliding mode surface design, s (0) exp (- λ t), λ >0, and s (0) is s (t) at the initial time.
Define L yapunov function as
Figure BDA0001802586650000163
Then:
Figure BDA0001802586650000164
the synovial membrane control law is therefore:
Figure BDA0001802586650000165
wherein: k (t) is a switching gain, and k (t) max (| e (t) |) + ρ, ρ > 0.
Bringing formula (19) into formula (18):
V=-K(t)|s|-E(t)s (20)
wherein V is less than or equal to- η | s | is less than or equal to 0.
From the above, the stability theorem of Liya Punuo Fu can be used
Figure BDA0001802586650000166
When ss<0, the slip film is present.
Fuzzy rules:
1. if ss >0, K (t) should be increased.
2. ss <0, K (t) should decrease.
Let ss be the system input and Δ K (t) be the system output, according to the fuzzy rule, the fuzzy is defined as follows:
ss={NB NM ZO PM PB}
ΔK={NB NM ZO PM PB} (21)
wherein NB, NM, ZO, PM, PB are respectively negative big, negative middle, zero, middle and positive big.
Estimating the upper limit of K (t) by adopting an integral method:
Figure BDA0001802586650000171
wherein: g is a proportionality coefficient, G > 0. Carrying in formula (19) to obtain:
Figure BDA0001802586650000172
wherein the content of the first and second substances,
Figure BDA0001802586650000173
r(t,x1,x2) Is an uncertain function and satisfies the following boundary conditions:
Figure BDA0001802586650000174
the state trajectory can be set to converge spirally around the origin to the origin within a finite time on the sigma-sigma plane using a supercoiling algorithm. The specific algorithm is as follows:
Figure BDA0001802586650000175
wherein M is1、M2Satisfies the following conditions:
Figure BDA0001802586650000181
for the system (13), a nonlinear sliding mode variable (14) is selected, a condition (24) is met, and under the action of a superspiral control law (25), if the condition (26) is met, the system converges in a limited time.
S4, simultaneously optimizing parameters of the tracking controller, the synchronous controller, the anti-backlash controller and the switching function in the dual-motor system by using a particle swarm optimization, as shown in figure 4.
And selecting a particle swarm algorithm-based setting parameter. In order to guarantee two conditions of minimum tracking performance and energy consumption, the following performance index functions are designed:
Figure BDA0001802586650000182
smaller first term of the performance index indicates smaller steady state error of the system; the smaller the second term, the smaller the power consumption of the system. The constraints are to ensure that the tracking accuracy is within 10% and that the total energy consumed by the system is minimal. And integrating the tracking precision and the energy consumption into an objective function E, wherein the final objective of the optimization is to obtain the minimum value of the objective function E under the condition of meeting the constraint condition.
S401, initializing each parameter value of the particle swarm, and calculating a fitness function of each particle. And calculating the size of the performance index (target function adaptive value) of each group of parameters in a weighted mode.
Updating the adaptation value with the new particle if the adaptation value of the new particle is smaller than that of the previous one; otherwise, the adaptation value remains unchanged.
Figure BDA0001802586650000183
Where pbest (t) is the best fit at time t, f (-) is usually an objective function of the performance metric, i.e., a weighting of multiple significant metrics, and x (t) is the position of each particle.
pbest is used for recording the optimal solution searched by the individual, and gbest is used for recording the optimal solution searched by the whole population in one iteration.
The updated equations for velocity and particle position are as follows:
V[i]=w×v[i]+c1×rand()×(pbest[i]-present[i])+c2×rand()×(gbest-present[i]) (29)
where v [ i ] represents the velocity of the ith particle, w represents the inertia weight, c1 and c2 represent the learning parameters, rand () represents a random number between 0 and 1, pbest [ i ] represents the optimum value searched for by the ith particle, gbest represents the optimum value searched for by the entire cluster, and present [ i ] represents the current position of the ith particle.
S402, when the minimum adaptive value in pbest is smaller than the global adaptive value, updating the global adaptive value by using the position of the corresponding minimum adaptive value; otherwise, the global adaptation value remains unchanged.
gbest(t+1)=argmin{f(pbest1(t)),f(pbest2(t)),...,f(pbestn(t))} (30)
Wherein, gbest (t) is the global optimum adaptive value at time t, and n is the total number of particles.
S403, updating the control parameter value according to the following formula
Figure BDA0001802586650000191
Xi'j(t+1)=Xi'j(t)+Vi'j(t+1) (31)
Wherein Vi' j (t) is the speed of the ith population iteration particle of the jth dimension particle; xi ' j (t) is the position of the ith ' group iteration particle of the jth dimension particle, and i ' is taken as 20; omega is inertia weight, and is taken as 0.7; c1 and c2 are learning rates, and c1 is equal to c2 is equal to 2.
S404, the current position of the particle exceeds the set maximum value and the set minimum value, and the particles beyond the range are re-assigned, namely
Figure BDA0001802586650000192
Similarly, for the current velocity of the particle, when the maximum velocity of the particle is exceeded, the velocity of the particle is reassigned to
Figure BDA0001802586650000201
Wherein, xmin (j) and xmax (j) are respectively the minimum position and the maximum position of the j dimension; vmin (j) and Vmax (j) are the minimum and maximum speeds, respectively, for dimension j.
S405, when the iteration times are smaller than the maximum set times, turning to the step S3; otherwise, the procedure is terminated.
In the control of the motor of the balance car, under the condition that the particle swarm algorithm finds the optimal parameters, not only the tracking error of the system does not generate tremble, but also the two motors only eliminate backlash when backlash occurs, and keep synchronous operation under other conditions, thereby realizing the comprehensive control target of tracking, synchronization and backlash elimination. And the two motors alternately pass through the tooth gaps, and when one motor passes through the tooth gaps, the other motor drags the load, so that the function of offsetting moment backlash elimination is achieved. Therefore, the integrated controller designed by the embodiment can eliminate the nonlinear influence of the backlash and ensure the synchronization of the two motors on the premise of ensuring the load tracking performance.
The above-mentioned contents are only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited thereby, and any modification made on the basis of the technical idea of the present invention falls within the protection scope of the claims of the present invention.

Claims (2)

1. A control method of a balance car system is characterized in that the balance car system comprises a control part, a sensor part, a man-machine interaction part, a power supply part, a main circuit part and an upper computer part, wherein the control part adopts a microcontroller which is respectively connected with the sensor part, the man-machine interaction part, the power supply part and the main circuit part, the upper computer part is connected with the man-machine interaction part, the man-machine interaction part provides electric energy by the power supply part, the sensor part comprises a gyroscope for acquiring the angular velocity of the deviation of a car body, an accelerometer for acquiring the acceleration of the car, a speed sensor for acquiring the actual speed, a corner sensor for acquiring the direction angle of a control lever, an infrared sensor for detecting the road condition, a voltage sensor and a current sensor for acquiring relevant data of an inverter, the man-machine interaction part comprises an L CD display screen for displaying the speed information, the electric quantity information, the battery temperature information and the system running state feedback, a wireless Bluetooth module for transmitting the data to a mobile phone and sending instructions to the controller, the upper computer communication module for communicating with the PC, a power supply module for controlling the power switch, a mode selection, a key for prompting the vehicle alarm information, a L and a buzzer for prompting the vehicle alarm information, a battery, a power supply part, a battery, a charging system, a permanent magnet motor, a battery, a charging system, a permanent magnet motor and a permanent;
establishing a state equation of extended Kalman filtering and a mathematical model of the permanent magnet synchronous motor according to a difference equation model and a drift error model of an attitude signal of a balance car system; driving the permanent magnet synchronous motor by using a second-order dynamic Terminal sliding film; simultaneously optimizing parameters of a tracking controller, a synchronous controller, an anti-backlash controller and a switching function in a dual-motor system by using a particle swarm algorithm, wherein the state equation of the extended Kalman filtering is as follows:
Figure FDA0002520690540000011
wherein, k is a scale error and is a drift error, v (n-1) is white gaussian noise with a mean value of zero, T is a period, and a difference equation model of the attitude signal of the balance car is as follows:
Figure FDA0002520690540000021
wherein, ω is0The initial angle value of the gyroscope, κ is the scale error, and drift error, and the first-order Autoregressive (AR) model of the gyroscope drift error is as follows:
Figure FDA0002520690540000022
wherein the content of the first and second substances,
Figure FDA0002520690540000023
is an autoregressive model parameter; v (n) is white Gaussian noise with the mean value of zero;
the variance of gaussian white noise is as follows:
Figure FDA0002520690540000024
Figure FDA0002520690540000025
wherein, Ci、CgInclinometer and gyroscope noise covariance, respectively;1is the standard deviation of the tilt angle;2the mechanical motion equation of the permanent magnet synchronous motor is the standard deviation of the Gaussian noise density of the gyroscope as follows:
Figure FDA0002520690540000026
wherein s is mover displacement; v is the mover speed; m is the total mass of the rotor and the load carried by the rotor; b isvIs the viscous friction factor; fTotal disturbance force is total disturbance force FThe following were used:
F=Frip+Fload+Ffric
Figure FDA0002520690540000027
wherein, FloadIs the load resistance; fripEquivalent resistance to the end effect; fMThe end effect thrust fluctuation amplitude; tau is a polar distance;
Figure FDA0002520690540000028
is the initial phase electrical angle; ffricIs a friction force;
the driving of the permanent magnet synchronous motor by utilizing the second-order dynamic Terminal sliding film is specifically as follows:
let x be [ x ]1,x2]T=[s,v]TAs a state variable of the system, the input control quantity is u-iqThe equation of state is as follows:
Figure FDA0002520690540000031
Figure FDA0002520690540000032
wherein k is1,k2,k3For unknown parameters, kfIs a coefficient of resistance, BvIs viscous friction factor, FM is total disturbance force, and M is total mass of the rotor and the load carried by the rotor;
let the system tracking error be e ═ x1 *-x1Wherein x is1 *Is x1Defining the second-order nonsingular fast terminal sliding mode variable of the system as follows:
Figure FDA0002520690540000033
wherein, 0<α<1,β∈R+P, q ∈ N are odd numbers, lambda is more than p/q, and p/q is more than 1 and less than 2;
when the system state is close to the balance point, the high-order term of the tracking error e (t) approaches to 0, and the convergence speed of the tracking error e (t) is approximate to that of a nonsingular terminal sliding mode; when the system is far away from a balance point, the high-order terms of the tracking error e (t) play a main role, and the convergence speed of the tracking error e (t) is higher than that of a nonsingular terminal sliding mode; the sliding mode control system meets the sliding mode variable sigma and the first derivative thereof
Figure FDA0002520690540000034
Convergence to zero, first derivative of synovial variable
Figure FDA0002520690540000035
The following were used:
Figure FDA0002520690540000036
the synovial membrane control law is as follows:
Figure FDA0002520690540000037
wherein: r is a position command, ce is an error, f (t) is a function for achieving global sliding mode surface design, f (t) is s (0) exp (- λ t), λ >0, s (0) is an initial time of s (t), sgn(s) is a step function, and B is a friction coefficient; k (t) is a switching gain, k (t) max (| e (t) |) + ρ, ρ > 0; when the system input ss <0, the synovium is present; if the system input ss >0, the switching gain K (t) should be increased; if the system input ss <0, the switching gain K (t) should be reduced;
the performance index function F of the parameter setting by using the particle swarm optimization is as follows:
Figure FDA0002520690540000041
wherein, the smaller the first term of the performance index is, the smaller the steady-state error of the system is; the smaller the second term is, the smaller the energy consumption of the system is; t is time, e1(t) is a constraint term, u1,u2Is a comprehensive control law.
2. The balance car system control method according to claim 1, characterized by comprising the following specific steps:
s401, initializing each parameter value of a particle swarm, calculating a fitness function of each particle, and calculating the performance index size of each group of parameters in a weighting mode;
updating the adaptation value with the new particle if the adaptation value of the new particle is smaller than that of the previous one; otherwise, the adaptation value remains unchanged;
Figure FDA0002520690540000042
where pbest (t) is the best adaptation value at time t, f (-) is the objective function of the performance indicator, and X (t) is the position of each particle; pbest records the optimal solution searched by an individual, and gbest is used for recording the optimal solution searched by the whole group in one iteration;
the updated equations for velocity and particle position are as follows:
V[i]=w×v[i]+c1×rand()×(pbest[i]-present[i])+c2×rand()×(gbest-present[i])
wherein v [ i ] represents the velocity of the ith particle, w represents an inertia weight, c1 and c2 represent learning parameters, rand () represents a random number between 0 and 1, pbest [ i ] represents the optimal value searched by the ith particle, gbest represents the optimal value searched by the whole cluster, and present [ i ] represents the current position of the ith particle;
s402, when the minimum adaptive value in pbest is smaller than the global adaptive value, updating the global adaptive value by using the position of the corresponding minimum adaptive value; otherwise, the global adaptation value remains unchanged;
gbest(t+1)=argmin{f(pbest1(t)),f(pbest2(t)),...,f(pbestn(t))}
wherein gbest (t) is the global optimal adaptive value at the time t, and n is the total number of particles;
s403, updating the control parameter value according to the following formula
Figure FDA0002520690540000051
Xi'j(t+1)=Xi'j(t)+Vi'j(t+1)
Wherein Vi' j (t) is the speed of the ith population iteration particle of the jth dimension particle; xi ' j (t) is the position of the ith ' group iteration particle of the jth dimension particle, and i ' is taken as 20; omega is inertia weight, and is taken as 0.7; c. C1And c2To obtain the learning rate, take c1=c2=2;
S404, the current position of the particle exceeds the set maximum value and the set minimum value, and the particles beyond the range are re-assigned, namely
Figure FDA0002520690540000052
When the maximum velocity of the particle is exceeded, the velocity of the particle is reassigned to
Figure FDA0002520690540000053
Wherein, xmin (j) and xmax (j) are respectively the minimum position and the maximum position of the j dimension; vmin (j) and Vmax (j) are respectively the minimum speed and the maximum speed of the j dimension;
s405, when the iteration times are smaller than the maximum set times, turning to the step S403; otherwise, ending.
CN201811083745.3A 2018-09-17 2018-09-17 Balance car system and control method thereof Active CN109067264B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811083745.3A CN109067264B (en) 2018-09-17 2018-09-17 Balance car system and control method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811083745.3A CN109067264B (en) 2018-09-17 2018-09-17 Balance car system and control method thereof

Publications (2)

Publication Number Publication Date
CN109067264A CN109067264A (en) 2018-12-21
CN109067264B true CN109067264B (en) 2020-07-28

Family

ID=64762912

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811083745.3A Active CN109067264B (en) 2018-09-17 2018-09-17 Balance car system and control method thereof

Country Status (1)

Country Link
CN (1) CN109067264B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3611587A1 (en) * 2018-08-16 2020-02-19 Siemens Aktiengesellschaft System for controlling and monitoring of adaptive cyber physical systems
CN109861618B (en) * 2019-01-11 2020-11-20 江苏大学 Construction method of anti-interference composite controller of BSG alternating current motor for hybrid electric vehicle
CN110109353B (en) * 2019-04-17 2022-01-11 杭州电子科技大学 Fuzzy self-adaptive sliding-mode control system of counteractive wheel balance bicycle robot
CN110275439B (en) * 2019-06-28 2020-05-26 四川大学 Control method of self-balancing trolley and design method and device of controller
CN110504880B (en) * 2019-07-24 2021-01-26 东南大学盐城新能源汽车研究院 Feedforward compensation control method for interference observation of flux switching permanent magnet linear motor
CN110350841B (en) * 2019-07-31 2021-01-01 沈阳工业大学 Permanent magnet linear synchronous motor control device and method for numerical control machine tool
CN110794858A (en) * 2019-11-01 2020-02-14 浙江阿尔郎科技有限公司 Balance car main control system, control system and balance car
CN110794877B (en) * 2019-11-22 2020-10-13 北京理工大学 Vehicle-mounted camera holder servo system and control method
CN111224593B (en) * 2020-02-20 2021-12-31 浙江理工大学 Fuzzy self-adaptive sliding mode control method and system based on differential evolution algorithm optimization
CN113359434A (en) * 2021-04-15 2021-09-07 山东师范大学 Finite time tracking control method and system for electric balance car

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5325026A (en) * 1992-06-29 1994-06-28 General Electric Company Microprocessor-based commutator for electronically commutated motors

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316055B (en) * 2014-09-19 2018-08-17 南京航空航天大学 A kind of double-wheel self-balancing robot attitude algorithm method based on improved expanded Kalman filtration algorithm
CN105116729B (en) * 2015-08-17 2017-11-07 杭州电子科技大学 A kind of double-wheel self-balancing robot adaptive sliding mode variable structure control method
CN105425802B (en) * 2015-12-10 2018-04-06 长安大学 A kind of two-wheeled intelligent balance car and its control method
CN105785795B (en) * 2016-05-05 2018-11-23 北京交通大学 A kind of train running speed curve energy conservation optimizing method based on particle swarm algorithm

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5325026A (en) * 1992-06-29 1994-06-28 General Electric Company Microprocessor-based commutator for electronically commutated motors

Also Published As

Publication number Publication date
CN109067264A (en) 2018-12-21

Similar Documents

Publication Publication Date Title
CN109067264B (en) Balance car system and control method thereof
CN107368081B (en) A kind of double-wheel self-balancing robot adaptive sliding mode variable structure control system
CN110022109B (en) Method and system for controlling torque ripple of SRM (sequence-derived minimum-mean-square) of torque-current neural network model
CN109143868B (en) Nonlinear anti-interference control method and device for electronic throttle system
CN105785999A (en) Unmanned surface vehicle course motion control method
CN109343351B (en) Switched reluctance motor torque control system with improved PID control
CN105928514A (en) AGV composite guiding system based on image and inertia technology
CN108177693A (en) Wheel hub drives the electronic differential control system of electric vehicle
Prasetyo et al. Controls of BLDC motors in electric vehicle Testing Simulator
CN102198032A (en) A controller and control method for a motorised vehicle
CN103513574B (en) The building method of axial mixed magnetic bearing fish-swarm algorithm Fractional Order PID Controller
CN106877774B (en) Input supersonic motor servo adaptive control system and method under saturation conditions
CN113852305B (en) DC motor terminal sliding mode control method, system, equipment and medium
CN111273542B (en) Cubic robot, control system and method, and design method
CN108710380A (en) A kind of control system of Double-wheel self-balancing trolley
CN109884881A (en) A kind of design for surely taking aim at servo controller based on nonlinear PID controller technology
CN106950823A (en) A kind of Neural network PID System with Sliding Mode Controller for intelligent wheel chair
Zaheer et al. Performance analysis of path planning techniques for autonomous mobile robots
CN108388177B (en) Half-size computer mouse motion control system
CN113343592B (en) Intelligent control method for DQN (direct current) of permanent magnet synchronous motor of new energy aircraft
Deshmukh et al. An improvement to the conventional PD scheme for the speed control of a two-wheeled differential steering mobile robot
CN106842957B (en) Ultrasonic motor servo control method based on output feedback controller
Antar et al. Sensorless speed control of separately excited DC motor using neuro-fuzzy controller
Atay et al. Real-time power optimization of a bi-modal rolling-flying vehicle via extremum seeking control
Zhang et al. Nonlinear Model Predictive Control of an Omnidirectional Mobile Robot with Self-tuned Prediction Horizon

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210114

Address after: 1-1-505-1, Xindu building, southwest of the intersection of Weijin road and wandezhuang street, Nankai District, Tianjin

Patentee after: HUIJIAWANG (TIANJIN) TECHNOLOGY Co.,Ltd.

Address before: 710064 middle section of South Second Ring Road, Beilin District, Xi'an City, Shaanxi Province

Patentee before: CHANG'AN University

TR01 Transfer of patent right
OL01 Intention to license declared
OL01 Intention to license declared