B 题:大学生平衡膳食食谱的优化设计及评价
B题完整解析分享给大家,包括具体的求解过程,以及全部的代码与技术文档,需要完整代码与技术文档,关注【云顶数模】公众号,即可获得哦
1 数据预处理
首先,对“食物编码”进行处理。通过对“中国食物编码表”进行查询,发现编号中的“x”属于冗余信息。故对该信息进行删除。
同时,存在信息错误。附件1和附件3中“土豆”的食物编码出现重复。更正为“021101”。
其次,为了提升运算速度,对数据中的文本信息进行数字编码。“食物名称”替换信息如下表所示:
对附件3中的“是”统一替换成1,“否”统一替换成“0”。
对于缺失的氨基酸含量,由于同类别的食品的含量类似,故利用同类的氨基酸量替代缺失的量。
2 问题一
2.1 问题1.1
2.1.1 评价体系的构建
为了对两位同学的食物摄入情况进行评价,构建评价体系如下:
图1 评价体系
在这一节中,本文设计了一个全面的食物摄入情况评价体系,涵盖了膳食营养评价的关键方面。这个评价体系主要包括以下四个部分:
首先是对食物结构的分析。本文按照附件4中提到的五大类别,包括谷薯类、蔬菜水果类、畜禽蛋类、奶类干豆类、以及植物油类,对给定的食谱中的各种食物进行归类统计。这样可以检查食物种类是否达到每天12种以上的标准,从而评估膳食的多样性和均衡性。
其次是对能量来源结构的评价。本文计算膳食中碳水化合物、蛋白质、脂肪三大宏量营养素的供能占比,并将结果与附件4中给出的参考范围进行对比分析。这有助于判断膳食的能量结构是否合理,是评估膳食质量的重要指标。
同时,本文也关注膳食的总能量摄入量以及在三餐中的分配情况。本文将总能量摄入与男生2400千卡/天或女生1900千卡/天的目标值进行比较,并检查早中晚三餐的能量占比是否满足30%-40%的参考区间。这有助于评估膳食是否能够满足身体的能量需求。
最后,本文还单独关注膳食中蛋白质的摄入情况。本文计算每餐混合食物的蛋白质氨基酸评分,并根据附件4中给出的标准进行评价。这可以反映膳食中蛋白质的质量和利用率。
2.1.2 指标计算
《指南》中明确了平衡膳食的基本准则,其中定义了食物的五大类别:
(1)谷、薯类;
(2)蔬菜、菌藻、水果类;
(3)畜、禽、鱼、蛋类及制品;
(4)奶、干豆、坚果、种子类及制品;
(5)植物油类。
为了进行类别统计,需要根据食物编号对每个食物的类别进行区分。
2.1.3 指标计算结果
2.1.4 基于层次分析法的膳食营养评价模型
2.1.5 评价模型的求解
2.2 问题1.2
2.2.1 食物与成分间拓扑关系的构建
2.2.2 微调模型的建立
2.2.3 基于启发式算法的模型求解
部分代码
clc
clear
%% 读入数据
boy = xlsread('附件1:1名男大学生的一日食谱.xlsx');
boy(:,2) = [];
bound_boy = find(isnan(boy(:,1)));% 找到下午和晚上的开始索引
bound_boy(2) = bound_boy(2)-1;
boy = boy(~isnan(boy(:,1)),:);
girl = xlsread('附件2:1名女大学生的一日食谱.xlsx');
girl(:,2) = [];
bound_girl = find(isnan(girl(:,1)));
bound_girl(2) = bound_girl(2)-1;
girl = girl(~isnan(girl(:,1)),:);
element = xlsread('成分元素查询(无缺失,每克).xlsx');
%% 记录矩阵
boy_set = [];
girl_set = [];
%% 计算大类数
un = unique(boy(:,1));
for i=1:size(un,1)
t = un(i,1);
if t==1
un(i,2) = 1;
elseif t==2||t==4||t==5||t==6
un(i,2) = 2;
elseif t==8||t==9||t==11||t==12
un(i,2) = 3;
elseif t==3||t==7||t==10||t==18
un(i,2) = 4;
else
un(i,2) = 5;
end
end
boy_set(1,1) = size(unique(un(:,1)),1);
un = unique(girl(:,1));
for i=1:size(un,1)
t = un(i,1);
if t==1
un(i,2) = 1;
elseif t==2||t==4||t==5||t==6
un(i,2) = 2;
elseif t==8||t==9||t==11||t==12
un(i,2) = 3;
elseif t==3||t==7||t==10||t==18
un(i,2) = 4;
else
un(i,2) = 5;
end
end
girl_set(1,1) = size(unique(un(:,1)),1);
%% 计算小类数
boy_set(1,2) = size(unique(boy(:,2)),1);
girl_set(1,2) = size(unique(girl(:,2)),1);
%% 计算脂肪营养素
boy_set(1,3) = 0;
for i=1:size(boy,1)
g = boy(i,3)*boy(i,4);
for ii = 1:size(element,1)
if element(ii,1) == boy(i,2) % 找到了
fat = element(ii,5)*g;
boy_set(1,3) = boy_set(1,3)+fat;
break
end
end
end
boy_set(1,3) = boy_set(1,3)*9; %乘转换系数
数据与代码
完整资料
另外学姐也完成了电工杯A题的代码与技术文档,完整电工杯助攻资料,关注【云顶数模】公众号,即可获得哦~