辽宁朝阳网站建设,wordpress主页广告,朝阳区互联网大厂,化妆品行业网站建设2018年亚太杯APMCM数学建模大赛
A题 老年人平衡能力的实时训练模型
原题再现 跌倒在老年人中很常见。跌倒可能会导致老年人出现许多并发症#xff0c;因为他们的康复能力通常较差#xff0c;因此副作用可能会使人衰弱#xff0c;从而加速身体衰竭。此外#xff0c;对跌倒…2018年亚太杯APMCM数学建模大赛
A题 老年人平衡能力的实时训练模型
原题再现 跌倒在老年人中很常见。跌倒可能会导致老年人出现许多并发症因为他们的康复能力通常较差因此副作用可能会使人衰弱从而加速身体衰竭。此外对跌倒的恐惧可能会削弱行动能力限制行动范围从而显著恶化生活质量。因此对老年人进行平衡能力评估以帮助他们改善行动状态、纠正姿势和防止意外跌倒具有重要的现实意义。 目前还没有一个包罗万象的平衡定义。在医学中平衡有两层含义。一种是当人体保持稳定姿势时的静态平衡状态。另一种是当身体在运动或受到外力时自动调整自己以保持姿势时的动态平衡状态。在力学中当物体上的合力为零时就会产生平衡。身体的平衡或稳定性与重心的位置和支承面的面积有关。如果重心线落在轴承表面内或者发生不平衡则保持机体平衡。 一家研究机构通过在老年受试者身上部署42个监测点进行了随机抽样测试。见下图所示各点的布局。 根据附件中的实验数据完成以下任务 问题1根据附件2中的数据分析老年人的平衡特征。基于对步数、重心和运动的分析建立特征提取模型。应用由42个监测点组成的系统提取25个身体平衡特征对老年人的身体平衡进行综合评估。 问题2建立基于25个指标的平衡风险评估体系评估老年人的平衡能力。给出相应的建议。 问题3根据提供的实际数据对身体平衡力进行模拟计算和比较分析。为平衡能力较弱的老年人提供有效的建议。 附件说明 附件1《老年人基本数据》格式.xlsx包含所有老年受试者的基本数据。 附件2附件2包含每个受试者在自由行走状态下的校准原始数据。有三个完整的gait文件后缀.trc可以用Excel打开。第一列是帧序列第二列是时间。从第三列开始每三列代表监测点的运动坐标xyz。共有42个监测点。
整体求解过程概述(摘要) 跌倒会造成成吨的伤害这将影响老年人的生活水平。因此结合步行时的身体姿势等因素为他们建立一个平衡能力评估体系具有重要意义。 我们的模型首先采用R型聚类方法和马氏距离从42个监测点中提取出25个具有医学意义的指标。我们的25项指标包含90.86%的全身特征信息可靠性高。这一结果使测量每个受试者的基本身体状况成为可能。 然后基于这25个指标采用改进的多项式曲线拟合方法建立了一个平衡能力体系并将稳定配速曲线与实验配速曲线进行了比较以塑造一个人的平衡能力。然后我们的研究将我们的结果与实际跌倒时间进行了比较对于所有有跌倒风险的老年人来说准确率达到了近80%。然而当应用于中等风险的老年人时准确性会迅速下降。由于这种现象我们的模型做了更多的改进。 在我们改进的评估系统中我们的模型考虑了TOPSIS的实际数据。我们使用来自实际数据的7个参数和附录2的追踪数据来衡量老年人的平衡能力。与实际跌倒时间相比我们在几乎所有老年人中获得了80%的准确率。此外我们改进的评估可以从整体上划分出平衡能力高和平衡能力低的老年人。 改进模型的敏感性分析表明模型中年龄和BMI参数的微小变化对结果的影响较小。3%的干扰只会影响2.5%的排名5%的干扰会影响7.5%的排名。最后我们从医学角度对不同类型的跌倒或平衡能力低下的老年人提出了不同的建议。
模型假设 1.记录跟踪数据时没有时间差。 2.25个指标具有医学意义和物理意义能够预测重心等物理参数。 3.跟踪数据能够预测平衡能力。由于我们的R-TYPE CLUSTER算法、PCA算法都是基于实际数据的。 4.追踪数据是每个受试者的正常步骤即每个人在实验室里都像平时一样行走尤其是在摔倒前。 5.在长者跌倒和被追踪的这些年里身体没有发生大的变化。
问题重述 背景 在过去的几十年里我国的老龄化程度逐渐加快。预计到2050年60岁以上人口将占总人口的33%。这些年来老年人摔倒事故也在大规模增加。当老年人跌倒时他们会不稳定失去平衡对老年人来说下半身肌肉衰退会导致平衡能力下降不足以支撑老年人在跌倒时下半身同时由于老年人神经系统控制能力下降神经传导减慢运动反应时间延长也会导致老年人在跌倒时不能及时调整身体保持平衡。这样就有必要衡量老年人的平衡能力为他们提供适当的建议让他们保持平衡更加注意他们的平衡。 此外对老年人平衡能力的预测研究较少。此外大多使用单个传感器检测精度有限容易误判。 问题重述 建立一个模型其中包含来自给定42个监测点的25个主要指标。该模型需要用于分析每个受试者的步数、重心、速度和加速度等。 根据实际实验数据建立一个能够评估老年人平衡能力的模型。 模拟那些几年前摔倒的老年人并向我们的模型证明我们的模型是否能够监测他们的平衡能力。 为那些平衡能力较弱的老年人提供建议。
模型的建立与求解整体论文缩略图 全部论文及程序请见下方“ 只会建模 QQ名片” 点击QQ名片即可
程序代码
部分程序如下:
%% import dataresultvar_yy[];resultvar_zz[];A dir(fullfile(*.trc));
for iiii1:1:76data1load(A(iiii).name);result_yy[];result_zz[]; resultfit_yy[];resultfit_zz[]; data2data1(:,3:end);%matrix1*126 for i1:1:25%rowyy[];zz[];for ii2:1:5%columnif data3(i,ii)~0yy[yy,data2(:,3*data3(i,ii)-1)];zz[zz,data2(:,3*data3(i,ii))];endendresult_yy[result_yy,mean(yy,2)];result_zz[result_zz,mean(zz,2)];yytofitresult_yy(:,i);zztofitresult_zz(:,i);% Create the model.fun (x,xdata)x(1)x(2)*xdatax(3)*(xdata).^2x(4)*(xdata).^3 x(5)*
(xdata).^4 x(6)*(xdata).^5x(7)*(xdata).^6x(8)*(xdata).^7;tdata data1(:,2);m0 1.0e04 *[0.0280 0.1155 -0.6000 1.1672 -1.1046
0.5433 -0.1334 0.0129];% Solve the bounded fitting problem.my1 lsqcurvefit(fun,m0,tdata,yytofit);my2 lsqcurvefit(fun,my1,tdata,yytofit);yyfittedfun(my2,tdata);mz1 lsqcurvefit(fun,m0,tdata,zztofit);mz2 lsqcurvefit(fun,mz1,tdata,zztofit);
zzfittedfun(mz2,tdata);resultfit_yy[resultfit_yy,yyfitted];resultfit_zz[resultfit_zz,zzfitted]; difference_yyabs(resultfit_yy-result_yy);difference_zzabs(resultfit_zz-result_zz);var_yyvar(difference_yy,0,1);var_zzvar(difference_zz,0,1);endresultvar_yy[resultvar_yy;var_yy];resultvar_zz[resultvar_zz;var_zz];
end
%%CLUSTER
clc, clear all
%read in circle
A dir(fullfile(*.trc));
resultfinal[];
for i1:1:2data1load(A(i).name);[m,n]size(data1);a[]; for ii1:1:42X[data1(:,3*ii1),data1(:,3*ii2)];a(:,ii)pdist(X,mahal);end%%bzscore(a); %standard the datarcorrcoef(b); %coefficient matrixdpdist(b,correlation); %计算相关系数导出的距离zlinkage(d,average); %按类平均法聚类figurehdendrogram(z); %画聚类图set(h,Color,k,LineWidth,1.3) %把聚类图线的颜色改成黑色线宽加粗
Tcluster(z,maxclust,25); %cluster into 25 catagoryresults [];for iii1:25tmfind(Tiii); %find the ith coatogory’s subjetcttmreshape(tm,1,length(tm)); %in row
results[results;tm,zeros(1,15-length(tm))];end resultfinal[resultfinal;results];end
%% topsis
clc, clear
dataload(datatopsis.txt);
[m,n]size(data);
fun(qujian,lb,ub,x)(1-(qujian(1)-x)./(qujian(1)-lb)).*(xlb
xqujian(1))(xqujian(1) xqujian(2))(1-(x-qujian(2))./(ubqujian(2))).*(xqujian(2) xub);
%properties trans
qj2[0.41,0.48]; lb20.35; ub20.52;
data(:,2)fun(qj2,lb2,ub2,data(:,2));
qj3[2.5,5]; lb32; ub35.8;
data(:,3)fun(qj3,lb3,ub3,data(:,3));
qj6[19,28]; lb617; ub632;
data(:,6)fun(qj6,lb6,ub6,data(:,6));
for j1:nb(:,j)data(:,j)/norm(data(:,j)); % normalize the matrix
end
%% weight
% data analysis in 1
maxdatarepmat(max(data),m,1);
mindatarepmat(min(data),m,1);
max_minmaxdata-mindata;
stddata(data-mindata)./max_min;
%calculatee the weight
f(1stddata)./repmat(sum(1stddata),m,1);
e-1/log(m)*sum(f.*log(f));
d1-e;
wd/sum(d); % 权重向量
%%
cb.*repmat(w,m,1); %求加权矩阵
Cstarmax(c) %求正理想解
Cstar(1)min(c(:,1)); Cstar(4)min(c(:,4));
Cstar(5)min(c(:,5));
Cstar(7)max(c(:,7)); %属性 1,3,5,6 为成本型, 属性 2 为效益型
C0min(c) %q 求负理想解
C0(1)max(c(:,1)); C0(4)max(c(:,4));
C0(5)max(c(:,5));
C0(7)min(c(:,7)); %属性 1,3,5,6 为成本型, 属性 8 为效益型
for i1:mSstar(i)norm(c(i,:)-Cstar); %求到正理想解的距离S0(i)norm(c(i,:)-C0); %求到负理想的距离
end
fS0./(SstarS0);
[sf,ind]sort(f,descend); %求排序结果
Y[ind,sf];全部论文及程序请见下方“ 只会建模 QQ名片” 点击QQ名片即可
