之前我一般是把sqlite的数据导出到excel里面，然后把数据复制到matlab里面再跑算法的。

感觉有点低效，所以最近两天安装了一个sqlite驱动，用matlab打开db文件，读出数据之后直接跑算法。

%% main.m

%% 数据初始化

clear

close all

%获取数据库中的表名

len = input('单个位置采样点数为len: \ne.g.200\n');%单个位置采样点数为len

fileFolder=fullfile(pwd);

dirOutput=dir(fullfile(fileFolder,'*.db'));

fileNames={dirOutput.name}'

dbName = input('请输入db文件的代号: \ne.g.1\n');

dbName = fileNames(dbName)

dbName = strcat(strcat(pwd,'/'),dbName);

dbName = strrep(dbName,'\','/');

dbName = char(dbName);

tablename = fromdbNameGetTableName(dbName);

[m,n]=size(tablename);

len = 200

result = [];

for i=1:m

SingleAPData = [];

for j = 1:8

%先处理db文件转化为矩阵datalevel和r

datalevel = yy_db(dbName,tablename(i,:),j);

SingleAPData = [SingleAPData;datalevel];

end

SingleAPData = SingleAPData';

result = [result;ProcessSingleAP( SingleAPData,i,0.9 )];

end

%% fromdbNameGetTableName.m

%这是数据库获取表单的函数

%dbname 数据库文件名称

function tablename = fromdbNameGetTableName(dbName)

%连接数据库

conn = database(dbName,'yy','','org.sqlite.JDBC',strcat('jdbc:sqlite:',dbName));

ping(conn);

xs=exec(conn,'SELECT name FROM sqlite_master WHERE type=''table'' ORDER BY name');

xs=fetch(xs);

x=xs.Data;

tablename = cell2mat(x(3:end,1));

%% yy_db.m数据库操作都在这里了

%这是数据库数据处理函数

%将数据库中的对应数据转换为行向量datalevel,R

%dbname 数据库文件名称

%tablename 数据表的名称

%row 对应datalevel数据在第几列

%Length 对应单个位置采样点数长度

function [datalevel] = yy_db(dbName,tablename,row)

%连接数据库

conn = database(dbName,'yy','','org.sqlite.JDBC',strcat('jdbc:sqlite:',dbName));

ping(conn);

users = exec(conn,['select * from ',tablename]);

users=fetch(users);

user=users.Data;

%获取Level以及r

if row == 1

datalevel=transpose(cell2mat(user(:,row)));

else

datalevel=transpose(str2num(char(user(:,row))));

end

end

%% GetNumfromAllSameMetre

function [ num,ind ] = GetNumfromAllSameMetre( AP_mat,metre )

%UNTITLED2 此处显示有关此函数的摘要

% 此处显示详细说明

ind = find (AP_mat(:,2)==metre);

num = length(ind);

end

%% ProcessSingleAP

function Mean = ProcessSingleAP( AP,APnum,yuzhi )

%UNTITLED5 此处显示有关此函数的摘要

% 此处显示详细说明

Mean = [];

need_list = [0.5,1,3,5,7,9,11,13];

for i = 1:length(need_list)

ind = find(AP(:,2) == need_list(i));

tempdata = AP(ind,APnum+2);

if length(tempdata) == 0

Mean = [Mean,0];

else

Mean = [Mean,GetMeanfromFrequency(tempdata,yuzhi)];

end

end

end

%% GetMeanfromFrequency

function [ Mean ] = GetMeanfromFrequency( raw_data,yuzhi )

%UNTITLED4 此处显示有关此函数的摘要

% 此处显示详细说明

Mean = 0;

table = tabulate(raw_data);

temp_table = table;

gailv = 0;

yuzhi = yuzhi*100;

while gailv < yuzhi

[val,ind] = max(temp_table(:,3));

gailv = gailv + val;

Mean = temp_table(ind,1)*val/yuzhi + Mean;

temp_table(ind,:) = [];

end

end