matlab编程题： 出⾎性脑卒的病情分析及其预测 出血性脑卒中指非外伤性脑实质内血管破裂引起的脑出血，占 全部脑卒中发病率的10-15%。其病因复杂，通常因脑动脉瘤破 裂、脑动脉异常等因素，导致血液从破裂的血管涌入脑组织，从而 造成脑部机械性损伤，并引发一系列复杂的生理病理反应。 出血性脑卒中后，血肿范围扩大是预后不良的重要危险因素之 一。在出血发生后的短时间内，血肿范围可能因脑组织受损、炎症 反应等因素逐渐扩大，导致颅内压迅速增加，从而引发神经功能进 一步恶化，甚至危及患者生命。因此，监测和控制血肿的扩张是临 床关注的重点之一。 目前，某医院已有160例（100例训练数据集+60例独立测试 数据集）出血性脑卒中患者的个人情况、发病到首次影像检查的时 间、首次影像学检查（CT平扫）的编号，该部分信息可在文件夹 “附录信息文件\3 综合应用题一素材”中的“表1-患者列表及临 床信息”中查询。 在随后的治疗过程中，还对每一名患者进行了多次的影像学检 查（CT平扫），该部分数据可在文件夹“附录信息文件”中的 “表2-患者影像信息血肿的体积及位置”中查询。在表2中每一 行的数据就是每

第一个问题：pr_qmks3_zhyy11.m ```matlab % 导入文件“表3-答案文件.xlsx”的“Sheet1”中的区域“A4:D103” jie_qmks3_zhyy11 = readtable('附录信息文件\3 综合应用题一素材\表3-答案文件.xlsx', 'Sheet', 'Sheet1', 'Range', 'A4:D103'); jie_qmks3_zhyy11 = table2cell(jie_qmks3_zhyy11); % 导入表1和表2 table1 = readtable('附录信息文件\3 综合应用题一素材\表1-患者列表及临床信息.xlsx'); table2 = readtable('附录信息文件\表2-患者影像信息血肿的体积及位置.xlsx'); lookup_table = readtable('附录信息文件\附表1-检索表格-流水号vs时间.xlsx'); % 处理训练数据集患者sub001至sub100 for i = 1:100 patient_id = strcat('sub', sprintf('%03d', i)); % 获取首次影像检查流水号和发病到首次影像检查时间间隔 first_image_serial = table1{table1.ID == patient_id, '入院首次影像检查流水号'}; time_interval = table1{table1.ID == patient_id, '发病到首次影像检查时间间隔'}; % 获取首次检查血肿体积 first_volume_index = find(strcmp(table2.ID, patient_id) & strcmp(table2.首次检查流水号, first_image_serial)); first_volume = table2{first_volume_index, 'HM_volume'} * 1e-3; % 获取后续检查血肿体积 follow_up_volumes = []; follow_up_serials = []; for j = 3:2:size(table2, 2) if ~isnan(table2{first_volume_index, j}) follow_up_volumes = [follow_up_volumes, table2{first_volume_index, j+1} * 1e-3]; follow_up_serials = [follow_up_serials, table2{first_volume_index, j}]; end end % 判断是否发生血肿扩张 expansion_occurred = 0; expansion_time_diff = 0; for k = 1:length(follow_up_volumes) absolute_increase = follow_up_volumes(k) - first_volume; relative_increase = (follow_up_volumes(k) - first_volume) / first_volume; if absolute_increase >= 6 || relative_increase >= 0.33 expansion_occurred = 1; % 获取影像检查时间点 image_time = lookup_table{lookup_table.流水号 == follow_up_serials(k), '时间'}; expansion_time_diff = image_time - table1{table1.ID == patient_id, '发病时间'}; expansion_time_diff = expansion_time_diff + time_interval; break; end end % 将结果赋值给jie_qmks3_zhyy11 jie_qmks3_zhyy11{i, 3} = expansion_occurred; if expansion_occurred jie_qmks3_zhyy11{i, 4} = sprintf('%.2f', expansion_time_diff); end end % 将赋值后的变量输出到文件 writetable(cell2table(jie_qmks3_zhyy11), '附录信息文件\3 综合应用题一素材\表3-答案文件.xlsx', 'Sheet', 'Sheet1', 'Range', 'A4:D103'); ``` 第二个问题：pr_qmks3_ zhyy12.m ```matlab % 导入文件“表3-答案文件.xlsx”的“Sheet1”中的区域“A104:D163” jie_qmks3_zhyy12 = readtable('附录信息文件\3 综合应用题一素材\表3-答案文件.xlsx', 'Sheet', 'Sheet1', 'Range', 'A104:D163'); jie_qmks3_zhyy12 = table2cell(jie_qmks3_zhyy12); % 导入表1和表2 table1 = readtable('附录信息文件\3 综合应用题一素材\表1-患者列表及临床信息.xlsx'); table2 = readtable('附录信息文件\表2-患者影像信息血肿的体积及位置.xlsx'); lookup_table = readtable('附录信息文件\附表1-检索表格-流水号vs时间.xlsx'); % 处理独立测试数据集患者sub101至sub160 for i = 101:160 patient_id = strcat('sub', sprintf('%03d', i)); % 获取首次影像检查流水号和发病到首次影像检查时间间隔 first_image_serial = table1{table1.ID == patient_id, '入院首次影像检查流水号'}; time_interval = table1{table1.ID == patient_id, '发病到首次影像检查时间间隔'}; % 获取首次检查血肿体积 first_volume_index = find(strcmp(table2.ID, patient_id) & strcmp(table2.首次检查流水号, first_image_serial)); first_volume = table2{first_volume_index, 'HM_volume'} * 1e-3; % 获取后续检查血肿体积 follow_up_volumes = []; follow_up_serials = []; for j = 3:2:size(table2, 2) if ~isnan(table2{first_volume_index, j}) follow_up_volumes = [follow_up_volumes, table2{first_volume_index, j+1} * 1e-3]; follow_up_serials = [follow_up_serials, table2{first_volume_index, j}]; end end % 判断是否发生血肿扩张 expansion_occurred = 0; for k = 1:length(follow_up_volumes) absolute_increase = follow_up_volumes(k) - first_volume; relative_increase = (follow_up_volumes(k) - first_volume) / first_volume; if absolute_increase >= 6 || relative_increase >= 0.33 expansion_occurred = 1; break; end end % 将结果赋值给jie_qmks3_zhyy12 jie_qmks3_zhyy12{i-100, 4} = expansion_occurred; end % 将赋值后的变量输出到文件 writetable(cell2table(jie_qmks3_zhyy12), '附录信息文件\3 综合应用题一素材\表3-答案文件.xlsx', 'Sheet', 'Sheet1', 'Range', 'A104:E163'); ``` 额外需要注意的逻辑点 1. 在处理数据时，要注意血肿体积的单位换算，表2中体积单位是 $10^{-3}$ ml，需要转换为ml。 2. 在查找流水号对应的时间点时，要确保流水号的匹配准确，避免出现查找错误。 3. 在处理后续检查数据时，要考虑到可能存在的空值情况，避免程序出错。 [2025-06-22 13:47:40 | AI写代码神器 | 1690点数解答]