MATLAB编程：出血性脑卒中病情分析与血肿扩张预测代码解析

第一个问题：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. 在处理后续检查数据时，要考虑到可能存在的空值情况，避免程序出错。 ######[AI写代码神器 | 1690点数解答 | 2025-06-22 13:47:40]