CN106604229A - Indoor positioning method based on manifold learning and improved support vector machine - Google Patents

Indoor positioning method based on manifold learning and improved support vector machine Download PDF

Info

Publication number
CN106604229A
CN106604229A CN201611225272.7A CN201611225272A CN106604229A CN 106604229 A CN106604229 A CN 106604229A CN 201611225272 A CN201611225272 A CN 201611225272A CN 106604229 A CN106604229 A CN 106604229A
Authority
CN
China
Prior art keywords
support vector
data
vector machine
value
algorithm
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.)
Granted
Application number
CN201611225272.7A
Other languages
Chinese (zh)
Other versions
CN106604229B (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.)
Southeast University
Original Assignee
Southeast 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 Southeast University filed Critical Southeast University
Priority to CN201611225272.7A priority Critical patent/CN106604229B/en
Publication of CN106604229A publication Critical patent/CN106604229A/en
Application granted granted Critical
Publication of CN106604229B publication Critical patent/CN106604229B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • H04W4/021Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W64/00Locating users or terminals or network equipment for network management purposes, e.g. mobility management

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

The invention discloses an indoor positioning method based on manifold learning and an improved support vector machine. The method comprises a step of determining a positioning area, dividing the positioning area according to an indoor structural characteristic and a layout characteristic, and obtaining a classification result, a step of obtaining offline training data, and collecting hotspot RSS signal values which can be received by the reference points in different classification area as a training data set, a step of using an isometric mapping algorithm to carry out training data characteristic extraction, a step of using the training data to carry out support vector machine classified training, using a taboo search algorithm to carry out support vector machine classification hyper parameter searching, and establishing the support vector regression model of each category at the same time, a step of carrying out online positioning, collecting the RSS signal value of each hotspot of a target, using a support vector machine classification model to carry out classification, and obtaining the rough positioning area of the target, and a step of carrying out the accurate positioning of the target by using the support vector regression model according to the classification result. According to the method, the time-varying characteristic of the wireless signal intensity is effectively suppressed, and the precision is obviously improved.

Description

It is a kind of based on manifold learning with improve support vector machine indoor orientation method
Technical field
Indoor positioning field of the present invention, more particularly to it is a kind of based on manifold learning with improve support vector machine indoor positioning Method.
Background technology
It is accurately positioned in problem indoors, Chinese scholars have carried out correlational study, is constantly proposed numerous effective and feasible Location algorithm.Wherein as WLAN (WLAN) can accomplish all standing, low cost and other advantages, the scene based on WLAN point Analysis indoor positioning algorithms become the main direction of studying of domestic and international colleges and universities and research institution.This algorithm is referred to as fingerprint location, Its essence is to set up the mapping relations one by one of RSS (signal intensity) and particular location, so as to be positioned.The difficult point of this method It is the foundation of mapping relations data base (fingerprint database).And the RSS values of each signal access point (AP) have time-varying characteristics, Cause positioning precision to reduce, therefore noise reduction and feature extraction are carried out to RSS signal demands, so as to reduce the introducing of signal time-varying characteristics Error.Meanwhile, the rise with machine Learning Theory with develop rapidly, researchers are by neutral net, support vector machine etc. Machine learning algorithm is fused in wireless indoor fingerprinting localization algorithm, improves the precision of wireless indoor positioning.Numerous studies are tied Fruit shows that support vector machine (SVM) are strong with generalization ability, compared to neural network algorithm, can effectively overcome Expired Drugs. Thus SVM has the advantages that other machines learning algorithm is incomparable in terms of positioning indoors.During SVM sets up, which surpasses The selection of parameter directly affects result of calculation, and traditional grid-search algorithms precision is low, and search time is long, affects indoor positioning Real-time, thus it is also particularly important to the choosing method of SVM hyper parameters.
The content of the invention
Goal of the invention:The invention aims to improve indoor position accuracy, it is proposed that it is a kind of based on manifold learning with Improve the algorithm of support vector machine.Invention carries out being accurately positioned for indoor objects using the algorithm, using the method for manifold learning Data Dimensionality Reduction is carried out except making an uproar, indoor wireless positioning precision is improved.
Technical scheme:It is a kind of based on manifold learning and the indoor orientation method for improving support vector machine, it is including following Step:
Step 1:Determine positioning region, root doors structure feature is classified to positioning region, obtains specification area;
Step 2:Offline training data is obtained, each focus that the reference point in different classifications region can be received is gathered RSS signal values are used as training dataset;
Step 3:Offline training data feature extraction is carried out using Isometric Maps algorithm to the training dataset in step 2, Carry out dimensionality reduction and denoising;
Step 4:Vector machine classification based training is supported using the Offline training data in step 2, is calculated using TABU search Method is supported the search of vector machine hyper parameter, obtains the support vector cassification model for training, while setting up of all categories propping up Hold vector regression model;
Step 5:Tuning on-line is carried out, the RSS signal values of each focus at target are gathered, using propping up for training in step 4 Hold vector machine disaggregated model to be classified, obtain target positioning region substantially;Wherein, carry out when gathering echo signal value multi-party To collection;
Step 6:According to classification results, the accurately fixed of target is carried out using the support vector regression model set up in step 4 Position.
Further, carry out to positioning region in the step 1 classification is whether to carry out mutation to carry out according to signal Classification.
Further, the step 2 is specially:According to the size cases of each specification area, arrange with reference to points, every Multi-direction gathered data in individual reference point, it is ensured that the data multidirectional of reference point;Signal strength threshold is set simultaneously, works as signal When intensity is less than threshold value, the focus is 0 for the positioning point of reference of the reference point, and its RSS value is set to 0, obtains training number According to collection.
Further, data are trained using Isometric Maps algorithm to the training dataset in step 2 in the step 3 Feature extraction, carries out dimensionality reduction and denoising;Specially:
To the RSS signal data collection X for givingi, i=1,2 ..., k, k are sample total number:
4) the field figure of each point in data sample X is set up by K nearest neighbor algorithms,
5) geodesic curve distance is calculated, sets up X*X distance relation matrixes;If sample point and each other neighborhood sample point, the two Between geodesic curve distance be Euclidean distance, otherwise, using heap optimization dijkstra's algorithm calculate shortest path, use shortest path To represent geodesic distance;
6) Nonlinear Dimension Reduction is carried out by multidimensional scaling algorithm, obtains the low-dimensional character representation of high dimensional data;Adjust the distance square Battle array centralization simultaneously carries out singular value decomposition, calculates d eigenvalue of maximum (λ of distance matrix1,...,λd) and corresponding feature Vectorial U, constructs diagonal matrix Λ=diag (λ1,...,λd), obtain the low-dimensional character representation of sample data X.
Further, the step 4 is specially:
In the case of linear separability, an optimum Optimal Separating Hyperplane is found, be spaced geometry and reach maximum, hyperplane For:
ωTX+b=0 formulas (1)
Wherein, w represents weight vector, and x is input vector, and b represents bias vector;
If set of data samples and its being expressed as (xi,yi), i=1 ..., k, x ∈ Rd, y ∈ { -1 ,+1 }, wherein Inputs of the x for data, y are exported for result;In order to find this hyperplane, the quadratic programming problem of following formula is solved:
s.t.,yiTxi+ b) >=1, i=1 ..., n
Above-mentioned quadratic programming problem is solved by Lagrange duality method, abbreviation obtains following formula:
Wherein α is Lagrange multiplier, and the α in formula (3) is fixed, partial derivative is sought to w and b respectively, their local derviation is made Number is zero, and result is rewinded formula (3), and abbreviation is obtained:
s.t.,αi>=0, i=1 ..., n formulas (5)
α is obtained according to formula (5)iValue, just can rewind the value that formula (4) obtains w, the value of b is solved by supporting vector, and Suzanne Lenglen day multiplier αiThe corresponding sample that is not zero is supporting vector;
Under nonlinear situation, a mapping phi () is introduced, nonlinear data is mapped to into a higher dimensional space so as to Become linear separability, the Nonlinear Classification of indoor positioning is solved the problems, such as by introducing gaussian kernel function, introduce slack variable and come This problem is solved, the abbreviation form of non-linear Lagrange function is obtained:
s.t.,0≤αi≤ C, i=1 ..., n formulas (6)
In formula:C be introduce slack variable after parameter, referred to as penalty factor, and support vector machine need optimization amount, Gaussian kernel functionNonlinear Classification function is accordingly:
Wherein, it is as follows come the Nonlinear Classification flow process for solving the problems, such as indoor positioning by introducing gaussian kernel function:
(5) Tabu search algorithm initialization;According to parameter searching step sizes, the initialization of taboo list length is carried out, generated super The search grid of parameter (C, σ), designs the producing method and maximum iteration time of neighborhood;Neighborhood generating mode is:According to number Radius, the square neighborhood centered on current solution are set according to size;
(6) fitness is calculated, checks estimated performance;According to each solution design support vector machine, tested using k- folding intersections The fitness of card method inspection solution, k value determine according to training data size, limit minimum and must not be less than 4, if k- folding cross validations Prediction probability inverse as fitness value;
(7) update the taboo element of taboo list and record current optimal value, it is first as taboo using changing value is solved Element;
(8) cyclic search, checks whether and reaches termination condition or reach maximum iteration time, reach, stop search, Output optimal solution.
Beneficial effect:The characteristics of present invention is using according to building structure and indoor arrangement carries out region division, does not rely on Hotspot's distribution situation, strong adaptability.The dimensionality reduction of data is carried out except making an uproar using the method for manifold learning, RSS signals time-varying is reduced special Property the impact that brings, improve indoor position accuracy.Combined with returning using support vector cassification, entered using tabu search algorithm Row hyper parameter is searched for, and improves positioning precision.The present invention effectively inhibits the time-varying characteristics of wireless signal strength, in off-line training number In the case of less, compared to traditional location algorithm, precision is significantly improved.
Description of the drawings
Fig. 1 is the inventive method schematic flow sheet;
Fig. 2 is equidistant mapping algorithm flow chart in the present invention;
Fig. 3 is the residual variance schematic diagram in the present invention;
Fig. 4 is the simulated environment schematic diagram of the present invention;
Fig. 5 is the simulation result comparison diagram of the present invention.
Specific embodiment
The present invention is described in further detail with implementation example below in conjunction with the accompanying drawings:
As shown in figure 1, the present invention be it is a kind of based on manifold learning with improve support vector machine indoor orientation method, including Following step:
Step 1:Determine positioning region, positioning region is classified.Any one piece of region in building has oneself Construction characteristic, any one room have its spatial layout feature, these features be all to maintain in most of time it is constant, or even It is changeless.And the RSS signals of focus (AP) also can be distributed according to such structure, RSS signal intensitys are general not Can be mutated, but if running into obstacle or certain reason can cause mutation, so whether we be able to can enter according to signal Being classified, the signal intensity feature of such classification can be especially apparent for row mutation, to classification based training with predict it is also further accurate. And judge that the mutation of signal is determined according to the construction featuress and layout of building oneself itself here.For example we can will be every Individual room is used as a region class, or is have a bookshelf in room, or it is other decorate to being separated by, such room Between can be divided into two classes according to separator.This is a kind of sorting technique for not relying on AP distributions, compared to according to AP distributions To be classified, the suitability is higher.
Step 2:Offline training data is obtained, each focus that the reference point in different classifications region can be received is gathered (AP) RSS signal values are used as training dataset.According to the size cases of each specification area, appropriate reference points are set, In time allowed band, reference point is as more as possible.The multi-direction gathered data in each reference point, it is ensured that the data of reference point are more Directivity.Signal strength threshold is set simultaneously, when signal intensity is less than threshold value, it is believed that the AP is determined for the reference point Position point of reference is 0, and its RSS value is set to 0, simplifies and calculates.
Step 3:Data characteristicses extraction is trained using Isometric Maps algorithm, dimensionality reduction and denoising effect is reached.It is initial to instruct Practice data and be usually associated with numerous noises, and RSS signals also have the change such as time-varying characteristics, temperature, humidity affect RSS letters Number.Therefore, in order to adapt to this change, need to find the dependency between signal.These changes overall can change RSS signals, But the dependency between each reference point still can be remained.Meanwhile, present wireless network application is ripe, and one piece of region has More than ten or even tens AP access points, this can cause the training data of higher-dimension, cause the difficult of calculating.And manifold learning This two point are solved the problems, such as simultaneously can.
Isomap (Isometric Maps algorithm) is that Josh Tenenbaum exist《Science》On the manifold learning delivered it is initiative Algorithm, it has started the new battlefield of a data processing, and its advantage is to maintain the geometry feature between all data points pair, Noise resisting ability is strong, and learning performance is stable, it is ensured that the robustness and Global Optimality of result, has to processing RSS signal noises Good rejection ability.The theoretical frame of Isomap is MDS (Multidimensional Scaling), but is placed on manifold In theoretical frame.MDS is a kind of Method of Data with Adding Windows, and its principle is just so that distance keeps not the point after dimensionality reduction between any two Become.But MDS is designed just for theorem in Euclid space, for the calculating of distance is also what is completed using Euclidean distance.And Isomap Euclidean distance has simply been changed into algorithm geodesic curve (geodesic) distance in manifold.In Isomap, geodesic curve distance is just It is that beeline algorithm can use the classic algorithm in computer graph theory, such as come approximate with the beeline on figure between 2 points Dijkstra's algorithm etc..If to the RSS signal data collection X, i=1,2 ... that give, k, k are sample total number, then the algorithm stream Journey is as shown in Figure 2:
1) field of each point in data sample X is set up by KNN (kNN, k-NearestNeighbor) nearest neighbor algorithm Figure, the value of Neighbourhood parameter k is too small here, and precision can be caused to reduce, while influence of noise is larger, k value is excessive, can bring not phase into The point of pass, reduces efficiency of algorithm, and the principle that we choose k value here is
2) geodesic curve distance is calculated, sets up X*X distance relation matrixes.If sample point and each other neighborhood sample point, they it Between geodesic curve distance be the Euclidean distance between them, otherwise, using heap optimization dijkstra's algorithm calculate shortest path, Geodesic distance is represented with shortest path.
3) Nonlinear Dimension Reduction is carried out by MDS (multidimensional scaling algorithm), obtains the low-dimensional character representation of high dimensional data.To away from From matrix centralization and singular value decomposition is carried out, calculate d eigenvalue of maximum (λ of distance matrix1,...,λd) and it is corresponding Characteristic vector U, constructs diagonal matrix Λ=diag (λ1,...,λd), so just can calculate the low-dimensional feature of sample data X Represent.
By the algorithm, we can obtain the dimensionality reduction result of any low-dimensional, and Fig. 3 is the data redundancy after 1 to 5 dimension is reduced to Variogram.
Step 4:Vector machine (SVM) classification based training is supported using training data, SVM is carried out using tabu search algorithm Hyper parameter is searched for, and gives up traditional grid search, improves search speed.Set up support vector regression model of all categories simultaneously.
Support vector machine (Support Vector Machine, SVM) are a kind of engineerings based on Statistical Learning Theory Learning method, improves the generalization ability of learning machine by seeking structuring least risk, realizes empiric risk and fiducial range Minimize, so as to reach in the case where statistical sample amount is less, can also obtain the purpose of good statistical law.The essence of SVM It is two graders, in the case of linear separability, SVM attempts to look for an optimum Optimal Separating Hyperplane, is spaced geometry Maximum is reached, hyperplane is as shown in Equation 1:
ωTX+b=0 (1)
Wherein, w represents weight vector, and b represents bias vector.
If set of data samples and its being expressed as (xi,yi), i=1 ..., k, x ∈ Rd, y ∈ { -1 ,+1 }, wherein Inputs of the x for data, y are exported for result.In order to find this hyperplane, can pass through to solve the quadratic programming problem of following formula:
s.t.,yiTxi+ b) >=1, i=1 ..., n
Above-mentioned quadratic programming problem can be solved by Lagrange duality method, abbreviation obtains following formula:
Wherein α is Lagrange multiplier, and the α in (3) formula is fixed, partial derivative is sought to w and b respectively, their local derviation is made Number is zero, and result is rewinded (3) formula, and abbreviation is obtained:
s.t.,αi>=0, i=1 ..., n (5)
By (4), (5) are as can be seen that no w, two variables of b, it is only necessary to according to (5) in Lagrangian formula Obtain αiValue, just can rewind the value that (4) obtain w, the value of b can be solved by supporting vector, and Suzanne Lenglen day multiplier αiNo Be sample corresponding to zero be exactly supporting vector required for us.
Under nonlinear situation, it would be desirable to introduce a mapping phi (), nonlinear data is mapped to into a higher-dimension Space so as to become linear separability, but the nonlinear mapping for introducing can cause the dimension disaster for calculating, and kernel function exactly may be used To avoid this problem.Therefore, herein by introducing gaussian kernel function solving the problems, such as the Nonlinear Classification of indoor positioning.Due to There is noise jamming in RSS signals, indivedual supporting vectors can be caused to deviate just datas, it is therefore desirable to by introduce slack variable come Solve this problem.(6) abbreviation form of the formula for non-linear Lagrange function.
s.t.,0≤αi≤ C, i=1 ..., n (6)
In formula:C be introduce slack variable after parameter, referred to as penalty factor, and support vector machine need optimization amount, Gaussian kernel functionNonlinear Classification function is accordingly:
The quality of the prediction effect of SVM has vital effect to positioning precision, finds in research, different types of Kernel function is less to SVM performance impacts, and the parameter of kernel function and penalty factor be affect SVM prediction effects key because Element.TABU search (TS) algorithm is a kind of heuristic global Stepwise optimization algorithm, is one kind simulation to human mind's process.Prohibit Avoid fully demonstrate in searching algorithm collection neutralization diffusion two it is tactful, embody a concentrated reflection of in the local search ability of TS, can from one Row solution is set out, and seeks preferably solution, terminated with reaching local optimum in the neighborhood of this feasible solution.Meanwhile, in order to jump out Locally optimal solution, algorithm devise Diffusion Strategy.Function of the Diffusion Strategy by taboo list, in taboo list, record has arrived at Some information, algorithm is by the taboo to these table midpoints, and reaches some points do not searched for, so as to realize bigger region Search.
Gaussian kernel function is selected herein, and hyper parameter (C, σ) is the parameter that we need to optimize.Algorithm flow is as follows:
1) Tabu search algorithm initialization.According to parameter searching step sizes, the initialization of taboo list length is carried out, generate super ginseng The search grid of number (C, σ), designs the producing method and maximum iteration time of neighborhood.The neighborhood generating mode for adopting herein For:Radius, the square neighborhood centered on current solution are set according to size of data.
2) fitness is calculated, checks estimated performance.According to each solution design SVM, using k- folding cross validation method inspections The fitness of solution is looked into, k value is determined according to training data size, but in order to ensure precision, limit minimum and must not be less than 4.This literary grace The Prediction sum squares of cross validation are rolled over as fitness value with k-.
3) update the taboo element of taboo list and record current optimal value, adopt herein and changing value is solved as taboo Element.
4) cyclic search, checks whether and reaches termination condition or reach maximum iteration time, reach, stop search, defeated Go out optimal solution.
Step 5:Tuning on-line is carried out, the RSS signal values of each AP at target is gathered, using the svm classifier model for training Classified, obtained target positioning region substantially.Equally, need to carry out multi-direction collection when gathering echo signal value.
Step 6:According to classification results, being accurately positioned for target is carried out using support vector regression model.
A kind of manifold learning and improvement support vector machine method of being based on proposed by the present invention is by MATLAB simulation softwares Carry out emulation experiment, it is contemplated that many problems of indoor AP receiving points, we propose 10, and simulated environment is as shown in Figure 4.Simulation hardware Environment is Intel (R) Core (TM) i5-2410M CPU 2.30GHz, 10 operating system of 4G RAM, Windows.Such as Fig. 5 It is shown, it is IMTS-SVM algorithms and traditional SVM algorithm and neural network algorithm accuracy comparison figure, it can be seen that neural Network algorithm due to less in data volume, positioning it is interval it is big in the case of, study generalization ability is poor, and positioning precision is relatively low, is determining Precision in the range of error 2m of position only has 48.4%., because of the good advantage of its generalization, positioning precision is higher than ANN for SVM algorithm, fixed Precision in the range of error 2m of position is 53.2%, and the precision within 4m can reach 78.0%.But also due to positioning Interval larger, the impact of RSS signals easily fluctuation, positioning precision are extremely difficult to satisfied degree.And set forth herein IMTS-SVM Precision of the algorithm in the range of position error 2m reaches 63.6%, and the precision in the range of 4m has been up to 90.8%, compared to In traditional SVM algorithm 2m, precision improves 10.4%, and the precision in 4m improves 12.8%, and positioning precision obtains higher lifting.
The present invention effectively inhibits the time-varying characteristics of wireless signal strength (RSS), in the less situation of Offline training data Under, compared to traditional location algorithm, precision is significantly improved.This method is mainly comprised the following steps:Step one:It is determined that positioning Region, by positioning region according to doors structure feature, spatial layout feature is divided, and obtains the stronger classification results of dependency;Step Rapid two:Offline training data is obtained, each focus (AP) the RSS signals that the reference point in different classifications region can be received are gathered Value is used as training dataset;Step 3:Data characteristicses extraction is trained using Isometric Maps algorithm, dimensionality reduction is reached and denoising is made With;Step 4:Vector machine (SVM) classification based training is supported using training data, the super ginsengs of SVM is carried out using tabu search algorithm Number search, gives up traditional grid search, improves positioning precision.Set up support vector regression model of all categories simultaneously;Step Five:Tuning on-line is carried out, the RSS signal values of each AP at target is gathered, is classified using svm classifier model, obtain target big The positioning region of cause;Step 6:According to classification results, being accurately positioned for target is carried out using support vector regression model.

Claims (5)

1. it is a kind of based on manifold learning with improve support vector machine indoor orientation method, it is characterised in that including following Step:
Step 1:Determine positioning region, root doors structure feature is classified to positioning region, obtains specification area;
Step 2:Offline training data is obtained, each focus RSS letters that the reference point in different classifications region can be received are gathered Number value as training dataset;
Step 3:Offline training data feature extraction is carried out using Isometric Maps algorithm to the training dataset in step 2, is carried out Dimensionality reduction and denoising;
Step 4:Vector machine classification based training is supported using the Offline training data in step 2, is entered using tabu search algorithm Row support vector machine hyper parameter search for, obtain the support vector cassification model for training, at the same set up support of all categories to Amount regression model;
Step 5:Carry out tuning on-line, gather the RSS signal values of each focus at target, using the support trained in step 4 to Amount machine disaggregated model is classified, and obtains target positioning region substantially;Wherein, multi-direction adopting is carried out when gathering echo signal value Collection;
Step 6:According to classification results, being accurately positioned for target is carried out using the support vector regression model set up in step 4.
2. according to claim 1 based on manifold learning and the indoor orientation method for improving support vector machine, its feature exists In carrying out classifying according to whether signal can be mutated to classify to positioning region in the step 1.
3. according to claim 1 based on manifold learning and the indoor orientation method for improving support vector machine, its feature exists In the step 2 is specially:According to the size cases of each specification area, arrange with reference to points, it is multi-party in each reference point To gathered data, it is ensured that the data multidirectional of reference point;Signal strength threshold is set simultaneously, when signal intensity is less than threshold value When, the focus is 0 for the positioning point of reference of the reference point, and its RSS value is set to 0, training dataset is obtained.
4. according to claim 1 based on manifold learning and the indoor orientation method for improving support vector machine, its feature exists In, in the step 3 to the training dataset in step 2 using Isometric Maps algorithm be trained data characteristicses extraction, carry out Dimensionality reduction and denoising;Specially:
To the RSS signal data collection X for givingi, i=1,2 ..., k, k are sample total number:
1) the field figure of each point in data sample X is set up by K nearest neighbor algorithms,
2) geodesic curve distance is calculated, sets up X*X distance relation matrixes;If sample point and each other neighborhood sample point, therebetween Geodesic curve distance be Euclidean distance, otherwise, using heap optimization dijkstra's algorithm calculate shortest path, with shortest path come table Show geodesic distance;
3) Nonlinear Dimension Reduction is carried out by multidimensional scaling algorithm, obtains the low-dimensional character representation of high dimensional data;Adjust the distance in matrix The heart simultaneously carries out singular value decomposition, calculates d eigenvalue of maximum (λ of distance matrix1,...,λd) and corresponding characteristic vector U, constructs diagonal matrix Λ=diag (λ1,...,λd), obtain the low-dimensional character representation of sample data X.
5. according to claim 1 based on manifold learning and the indoor orientation method for improving support vector machine, its feature exists In the step 4 is specially:
In the case of linear separability, an optimum Optimal Separating Hyperplane is found, be spaced geometry and reach maximum, hyperplane is:
ωTX+b=0 formulas (1)
Wherein, w represents weight vector, and x is input vector, and b represents bias vector;
If set of data samples and its being expressed as (xi,yi), i=1 ..., k, x ∈ Rd, y ∈ { -1 ,+1 }, wherein x is The input of data, y are exported for result;In order to find this hyperplane, the quadratic programming problem of following formula is solved:
Above-mentioned quadratic programming problem is solved by Lagrange duality method, abbreviation obtains following formula:
Wherein α is Lagrange multiplier, and the α in formula (3) is fixed, partial derivative is sought to w and b respectively, and the partial derivative for making them is Zero, result is rewinded into formula (3), abbreviation is obtained:
α is obtained according to formula (5)iValue, just can rewind the value that formula (4) obtains w, the value of b is solved by supporting vector, and bright lattice Bright day multiplier αiThe corresponding sample that is not zero is supporting vector;
Under nonlinear situation, a mapping phi () is introduced, nonlinear data a higher dimensional space is mapped to into so as to become Linear separability, solves the problems, such as the Nonlinear Classification of indoor positioning by introducing gaussian kernel function, introduces slack variable to solve This problem, obtains the abbreviation form of non-linear Lagrange function:
In formula:C be introduce slack variable after parameter, referred to as penalty factor, and support vector machine need optimization amount, Gauss Kernel functionNonlinear Classification function is accordingly:
Wherein, it is as follows come the Nonlinear Classification flow process for solving the problems, such as indoor positioning by introducing gaussian kernel function:
(1) Tabu search algorithm initialization;According to parameter searching step sizes, the initialization of taboo list length is carried out, generate hyper parameter The search grid of (C, σ), designs the producing method and maximum iteration time of neighborhood;Neighborhood generating mode is:It is big according to data Little setting radius, the square neighborhood centered on current solution;
(2) fitness is calculated, checks estimated performance;According to each solution design support vector machine, cross validation side is rolled over using k- The fitness of method inspection solution, k value determine according to training data size, limit minimum and must not be less than 4, if k- folding cross validations is pre- The inverse of probability is surveyed as fitness value;
(3) update the taboo element of taboo list and record current optimal value, adopt and changing value is solved as taboo element;
(4) cyclic search, checks whether and reaches termination condition or reach maximum iteration time, reach, stop search, output Optimal solution.
CN201611225272.7A 2016-12-27 2016-12-27 Indoor positioning method based on manifold learning and improved support vector machine Active CN106604229B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611225272.7A CN106604229B (en) 2016-12-27 2016-12-27 Indoor positioning method based on manifold learning and improved support vector machine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611225272.7A CN106604229B (en) 2016-12-27 2016-12-27 Indoor positioning method based on manifold learning and improved support vector machine

Publications (2)

Publication Number Publication Date
CN106604229A true CN106604229A (en) 2017-04-26
CN106604229B CN106604229B (en) 2020-02-18

Family

ID=58604323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611225272.7A Active CN106604229B (en) 2016-12-27 2016-12-27 Indoor positioning method based on manifold learning and improved support vector machine

Country Status (1)

Country Link
CN (1) CN106604229B (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107333238A (en) * 2017-07-03 2017-11-07 杭州电子科技大学 A kind of indoor fingerprint method for rapidly positioning based on support vector regression
CN107688820A (en) * 2017-07-11 2018-02-13 浙江新再灵科技股份有限公司 A kind of Elevator Fault Diagnosis method based on BCSA Support Vector Machines Optimizeds
CN109143164A (en) * 2018-10-11 2019-01-04 哈尔滨工业大学 The radio signal source localization method returned based on Gaussian process
CN109214268A (en) * 2018-07-12 2019-01-15 浙江工业大学 A kind of packed tower liquid flooding state on_line monitoring method based on integrated manifold learning
CN109239655A (en) * 2018-09-21 2019-01-18 上海中兴易联通讯股份有限公司 A kind of wireless signal acquiring for indoor positioning and processing method and system
CN109511095A (en) * 2018-11-30 2019-03-22 长江大学 A kind of visible light localization method and system based on Support vector regression
CN109508730A (en) * 2018-09-27 2019-03-22 东南大学 A kind of ionosphere phase scintillation detection method based on non-linear SVM algorithm
CN109991591A (en) * 2018-01-02 2019-07-09 中兴通讯股份有限公司 Localization method, device, computer equipment and storage medium based on deep learning
CN110221266A (en) * 2019-06-11 2019-09-10 哈尔滨工程大学 A kind of marine radar target rapid detection method based on support vector machines
CN110263826A (en) * 2019-05-31 2019-09-20 河南大学 The construction method and its detection method of Noise non-linear procedure fault detection model
CN111309850A (en) * 2020-02-10 2020-06-19 深圳云天励飞技术有限公司 Data feature extraction method and device, terminal equipment and medium
CN111757328A (en) * 2020-06-23 2020-10-09 南京林业大学 Cross-technology communication cheating attack detection method
CN111885700A (en) * 2020-06-08 2020-11-03 广州杰赛科技股份有限公司 Mobile terminal positioning method and device combined with support vector machine
CN112200282A (en) * 2020-10-14 2021-01-08 上海海事大学 RFID intelligent book positioning method based on feature weighting support vector machine
CN112394320A (en) * 2020-04-26 2021-02-23 南京邮电大学 Indoor high-precision centroid positioning method based on support vector machine
CN114254265A (en) * 2021-12-20 2022-03-29 军事科学院***工程研究院网络信息研究所 Satellite communication interference geometric analysis method based on statistical manifold distance
CN114609602A (en) * 2022-03-09 2022-06-10 电子科技大学 Feature extraction-based target detection method under sea clutter background

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2007214360A1 (en) * 2006-09-01 2008-03-20 Innovative Dairy Products Pty Ltd Whole genome based genetic evaluation and selection process
CN103164709A (en) * 2012-12-24 2013-06-19 天津工业大学 Method for optimizing support vector machine based on tabu search algorithm
CN104619014A (en) * 2015-01-09 2015-05-13 中山大学 SVM-KNN (Support Vector Machine-K Nearest Neighbor)-based indoor positioning method
CN105301558A (en) * 2015-09-22 2016-02-03 济南东朔微电子有限公司 Indoor positioning method based on bluetooth position fingerprints

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2007214360A1 (en) * 2006-09-01 2008-03-20 Innovative Dairy Products Pty Ltd Whole genome based genetic evaluation and selection process
CN103164709A (en) * 2012-12-24 2013-06-19 天津工业大学 Method for optimizing support vector machine based on tabu search algorithm
CN104619014A (en) * 2015-01-09 2015-05-13 中山大学 SVM-KNN (Support Vector Machine-K Nearest Neighbor)-based indoor positioning method
CN105301558A (en) * 2015-09-22 2016-02-03 济南东朔微电子有限公司 Indoor positioning method based on bluetooth position fingerprints

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
FENG ZIYONG: "DLANet: A manifold-learning-based discriminative feature learning network for scene classification", 《WEB OF SCIENCE》 *
桑楠等: "基于SVM分类和回归的WiFi室内定位方法", 《计算机应用研究》 *
邓志安: "基于学习算法的WLAN室内定位技术研究", 《信息科技辑》 *

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107333238A (en) * 2017-07-03 2017-11-07 杭州电子科技大学 A kind of indoor fingerprint method for rapidly positioning based on support vector regression
CN107688820B (en) * 2017-07-11 2020-09-08 浙江新再灵科技股份有限公司 Elevator fault diagnosis method based on BCSA optimized support vector machine
CN107688820A (en) * 2017-07-11 2018-02-13 浙江新再灵科技股份有限公司 A kind of Elevator Fault Diagnosis method based on BCSA Support Vector Machines Optimizeds
CN109991591B (en) * 2018-01-02 2023-08-22 中兴通讯股份有限公司 Positioning method and device based on deep learning, computer equipment and storage medium
CN109991591A (en) * 2018-01-02 2019-07-09 中兴通讯股份有限公司 Localization method, device, computer equipment and storage medium based on deep learning
CN109214268A (en) * 2018-07-12 2019-01-15 浙江工业大学 A kind of packed tower liquid flooding state on_line monitoring method based on integrated manifold learning
CN109214268B (en) * 2018-07-12 2021-08-03 浙江工业大学 Packed tower flooding state online monitoring method based on integrated manifold learning
CN109239655A (en) * 2018-09-21 2019-01-18 上海中兴易联通讯股份有限公司 A kind of wireless signal acquiring for indoor positioning and processing method and system
CN109508730A (en) * 2018-09-27 2019-03-22 东南大学 A kind of ionosphere phase scintillation detection method based on non-linear SVM algorithm
CN109508730B (en) * 2018-09-27 2021-07-27 东南大学 Ionosphere phase scintillation detection method based on nonlinear SVM algorithm
CN109143164A (en) * 2018-10-11 2019-01-04 哈尔滨工业大学 The radio signal source localization method returned based on Gaussian process
CN109511095B (en) * 2018-11-30 2021-06-04 长江大学 Visible light positioning method and system based on support vector machine regression
CN109511095A (en) * 2018-11-30 2019-03-22 长江大学 A kind of visible light localization method and system based on Support vector regression
CN110263826A (en) * 2019-05-31 2019-09-20 河南大学 The construction method and its detection method of Noise non-linear procedure fault detection model
CN110221266B (en) * 2019-06-11 2022-12-13 哈尔滨工程大学 Marine radar target rapid detection method based on support vector machine
CN110221266A (en) * 2019-06-11 2019-09-10 哈尔滨工程大学 A kind of marine radar target rapid detection method based on support vector machines
CN111309850A (en) * 2020-02-10 2020-06-19 深圳云天励飞技术有限公司 Data feature extraction method and device, terminal equipment and medium
CN111309850B (en) * 2020-02-10 2022-03-25 深圳云天励飞技术股份有限公司 Data feature extraction method and device, terminal equipment and medium
CN112394320B (en) * 2020-04-26 2023-06-23 南京邮电大学 Indoor high-precision centroid positioning method based on support vector machine
CN112394320A (en) * 2020-04-26 2021-02-23 南京邮电大学 Indoor high-precision centroid positioning method based on support vector machine
CN111885700A (en) * 2020-06-08 2020-11-03 广州杰赛科技股份有限公司 Mobile terminal positioning method and device combined with support vector machine
CN111885700B (en) * 2020-06-08 2022-04-12 广州杰赛科技股份有限公司 Mobile terminal positioning method and device combined with support vector machine
CN111757328A (en) * 2020-06-23 2020-10-09 南京林业大学 Cross-technology communication cheating attack detection method
CN112200282A (en) * 2020-10-14 2021-01-08 上海海事大学 RFID intelligent book positioning method based on feature weighting support vector machine
CN114254265A (en) * 2021-12-20 2022-03-29 军事科学院***工程研究院网络信息研究所 Satellite communication interference geometric analysis method based on statistical manifold distance
CN114609602A (en) * 2022-03-09 2022-06-10 电子科技大学 Feature extraction-based target detection method under sea clutter background

Also Published As

Publication number Publication date
CN106604229B (en) 2020-02-18

Similar Documents

Publication Publication Date Title
CN106604229A (en) Indoor positioning method based on manifold learning and improved support vector machine
CN112784881B (en) Network abnormal flow detection method, model and system
CN112101430B (en) Anchor frame generation method for image target detection processing and lightweight target detection method
CN103648106B (en) WiFi indoor positioning method of semi-supervised manifold learning based on category matching
CN104732244B (en) The Classifying Method in Remote Sensing Image integrated based on wavelet transformation, how tactful PSO and SVM
CN103336968A (en) Hyperspectral data dimension reduction method based on tensor distance patch calibration
CN104794368A (en) Rolling bearing fault classifying method based on FOA-MKSVM (fruit fly optimization algorithm-multiple kernel support vector machine)
Miranda et al. A hybrid meta-learning architecture for multi-objective optimization of SVM parameters
Martínez-Ballesteros et al. Selecting the best measures to discover quantitative association rules
Yi et al. An improved initialization center algorithm for K-means clustering
CN101833667A (en) Pattern recognition classification method expressed based on grouping sparsity
Peng et al. A hierarchical particle swarm optimizer with latin sampling based memetic algorithm for numerical optimization
CN111079074A (en) Method for constructing prediction model based on improved sine and cosine algorithm
CN116821715A (en) Artificial bee colony optimization clustering method based on semi-supervision constraint
Jia et al. TTSL: An indoor localization method based on Temporal Convolutional Network using time-series RSSI
Ding et al. Fruit fly optimization algorithm based on a novel fluctuation model and its application in band selection for hyperspectral image
Hu et al. Fuzzy clustering with knowledge extraction and granulation
Yi et al. New feature analysis-based elastic net algorithm with clustering objective function
Fiannaca et al. Simulated annealing technique for fast learning of SOM networks
Doshi et al. Graph neural networks with parallel neighborhood aggregations for graph classification
CN108549936A (en) The Enhancement Method that self organizing neural network topology based on deep learning is kept
CN104881688A (en) Two-stage clustering algorithm based on difference evolution and fuzzy C-means
Celik et al. Change detection without difference image computation based on multiobjective cost function optimization
CN111639712A (en) Positioning method and system based on density peak clustering and gradient lifting algorithm
CN107884744B (en) Passive indoor positioning method and device

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