%% a
load data
Num=unique(Table(2:end,1),'stable');
T=[];%记录复查时间
ED=[];%记录水肿数据
A=[];%记录除此诊断信息
for i=1:length(Num)
a=find(Table(:,1)==Num(i));
A=[A;double(Table(a(1),4:end))];
T=[T;double(Table(a,3))];
ED=[ED;double(Table(a,34))];
end
figure
plot(T,ED,'*')
xlim([1,2000])
%高斯模型
hold on
gaussModel = fit(T, ED, 'gauss1')
plot(gaussModel)
xlabel('时间')
ylabel('水肿/10^-3ml')
%计算残差
ED_fit=gaussModel(T);
Error=[];
for i=1:length(Num)
a=find(Table(:,1)==Num(i))-1;
Error(i,1)=mean(abs(ED_fit(a)-ED(a)));
end
%% b
%个体差异用Tabel第4-15列的指标来进行聚类
cluster_n=5;%聚类中心
[center, U, obj_fcn] = fcm(double(A(:,1:12)), cluster_n);
figure%目标函数变化值
plot(obj_fcn)
xlabel('iteration')
ylabel('obj.fcn_value')
title('FCM聚类')
[~,u]=max(U);
ED_fit2=[];
for i=1:cluster_n
a=find(u==i);
b=find(ismember(Table(2:end,1),Num(a))==1);
figure
plot(T(b),ED(b),'*')
xlim([1,2000])
%高斯模型
hold on
disp('亚类1:')
gaussModel = fit(T(b), ED(b), 'gauss1')
plot(gaussModel)
xlabel('时间')
ylabel('水肿/10^-3ml')
title(['亚类',num2str(i)])
ED_fit2(b,1)=gaussModel(T(b));
end
%计算残差
Error2=[];
for i=1:length(Num)
a=find(Table(:,1)==Num(i))-1;
Error2(i,1)=mean(abs(ED_fit2(a)-ED(a)));
end
result2=[["患者","残差(全体)","残差(亚类)","亚类"];[Num,Error,Error2,u']];
%% c
%计算水肿指标的变化率,在不同治疗方法下,改善效率
K=[];
for i=1:length(Num)
a=find(Table(:,1)==Num(i));
if length(a)>=5%至少五次检查数据才进行计算
k=(double(Table(a(2:end),34))-double(Table(a(1:end-1),34)))./(double(Table(a(2:end),3))-double(Table(a(1:end-1),3)));
kk=find(k<0);
if length(kk)==0
K(i,1)=0;
else
K(i,1)=mean(k(kk));
end
else
K(i,1)=NaN;
end
end
c=setdiff([1:length(K)],find(isnan(K)==1));
G=double(Table(2:end,16:22));
Z=Table(1,16:22);
for i=1:length(Z)
[p(i),anovatab{i},stats{i}]=anova1(K(c),G(c,i)+1,'off');%单因素方差分析
fa=finv(0.95,anovatab{i}{2,3},anovatab{i}{3,3});%计算fa
F=anovatab{i}{2,5};%F值
if p(i)<=0.01 && F>fa
disp([Z(i)+"对水肿体积进展模式影响非常显著"])
fprintf('p值为%.4f<0.01,F值为%.2f>%.2f\n',p(i),F,fa)
elseif p(i)<=0.05 && F>fa
disp([Z(i)+"对水肿体积进展模式影响显著"])
fprintf('p值为%.4f<0.05,F值为%.2f>%.2f\n',p(i),F,fa)
else
disp([Z(i)+"对水肿体积进展模式影响不显著"])
fprintf('p值为%.4f,F值为%.2f\n',p(i),F)
end
end
disp('不同治疗对水肿进展模式的影响大小为:')
[~,q]=sort(p);
str=Z(q)+">";
disp(str)
%% d
%同上述步骤求血肿体积的
%计算水肿指标的变化率,在不同治疗方法下,改善效率
K2=[];
for i=1:length(Num)
a=find(Table(:,1)==Num(i));
if length(a)>=5%至少五次检查数据才进行计算
k=(double(Table(a(2:end),23))-double(Table(a(1:end-1),23)))./(double(Table(a(2:end),3))-double(Table(a(1:end-1),3)));
kk=find(k<0);
if length(kk)==0
K2(i,1)=0;
else
K2(i,1)=mean(k(kk));
end
else
K2(i,1)=NaN;
end
end
c=setdiff([1:length(K2)],find(isnan(K2)==1));
for i=1:length(Z)
[p2(i),anovatab2{i},stats2{i}]=anova1(K2(c),G(c,i)+1,'off');%单因素方差分析
fa=finv(0.95,anovatab2{i}{2,3},anovatab2{i}{3,3});%计算fa
F=anovatab2{i}{2,5};%F值
if p2(i)<=0.01 && F>fa
disp([Z(i)+"对血肿肿体积进展模式影响非常显著"])
fprintf('p值为%.4f<0.01,F值为%.2f>%.2f\n',p2(i),F,fa)
elseif p2(i)<=0.05 && F>fa
disp([Z(i)+"对血肿体积进展模式影响显著"])
fprintf('p值为%.4f<0.05,F值为%.2f>%.2f\n',p2(i),F,fa)
else
disp([Z(i)+"对血肿体积进展模式影响不显著"])
fprintf('p值为%.4f,F值为%.2f\n',p2(i),F)
end
end
disp('不同治疗对血肿进展模式的影响大小为:')
[~,q2]=sort(p2);
str2=Z(q2)+">";
disp(str2)
%在算下血肿指标与水肿指标的相关性
x0=double(Table(2:end,23));
y0=double(Table(2:end,34));
theta=x0'*y0/(norm(x0)*norm(y0));%余弦相似度
fprintf('血肿指标与水肿指标的相关度为:%.4f\n',theta)
没有合适的资源?快使用搜索试试~ 我知道了~
基于Matlab求解2023年华为杯E题程序(源码+数据).rar
共10个文件
xlsx:5个
m:4个
mat:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 6 下载量 155 浏览量
2023-09-26
10:14:52
上传
评论 7
收藏 846KB RAR 举报
温馨提示
1、资源内容:基于Matlab求解2023年华为杯E题程序(源码+数据).rar 2、适用人群:计算机,电子信息工程、数学等专业的学习者,作为自己程序部分功能的“参考资料”使用。 3、解压说明:本资源需要电脑端使用WinRAR、7zip等解压工具进行解压,没有解压工具的自行百度下载即可。 4、免责声明:本资源作为“参考资料”而不是“定制需求”,代码只能作为参考,不能完全复制照搬。不一定能够满足所有人的需求,需要有一定的基础能够看懂代码,能够自行调试代码并解决报错,能够自行添加功能修改代码。由于作者大厂工作较忙,不提供答疑服务,如不存在资源缺失问题概不负责,谢谢理解。
资源推荐
资源详情
资源评论
收起资源包目录
基于Matlab求解2023年华为杯E题程序(源码+数据).rar (10个子文件)
基于Matlab求解2023年华为杯E题程序(源码+数据)
表1-患者列表及临床信息.xlsx 27KB
附表1-检索表格-流水号vs时间.xlsx 30KB
表4-答案文件.xlsx 13KB
yuchuli.m 3KB
wen3.m 2KB
wen1.m 2KB
表2-患者影像信息血肿及水肿的体积及位置.xlsx 183KB
wen2.m 4KB
表3-患者影像信息血肿及水肿的形状及灰度分布.xlsx 503KB
data.mat 130KB
共 10 条
- 1
Matlab仿真实验室
- 粉丝: 2w+
- 资源: 2181
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
- 1
- 2
前往页