close all; clear all; clc
data = csvread('G:./matlabdata./data3.8.csv', 2, 1); %K-T-model yingbi
Km = data(:, 9);
Gm = data(:, 12);
Apa= data(:,33);
X0 = data(:, 17); % 0MPa por
X5 = data(:, 20); % 5MPa por
X10 = data(:, 23); % 10MPa por
X15 = data(:, 26); % 15MPa por
X20 = data(:, 29); % 20MPa por
X25 = data(:, 32); % 25MPa por
Kd0 = abs((Km.*(-4.*Apa.*X0.*Gm.^3.*Km-12.*Apa.*X0.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+24.*pi.*Gm.*Km+9.*pi.*Km.^2))./(3.*Apa.*X0.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+9.*Apa.*X0.*Gm.*Km.^3+24.*pi.*Gm.*Km+9.*pi.*Km.^2)); % 0MPa Kd
Kd5 = abs((Km.*(-4.*Apa.*X5.*Gm.^3.*Km-12.*Apa.*X5.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+24.*pi.*Gm.*Km+9.*pi.*Km.^2))./(3.*Apa.*X5.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+9.*Apa.*X5.*Gm.*Km.^3+24.*pi.*Gm.*Km+9.*pi.*Km.^2)); % 5MPa Kd
Kd10 = abs((Km.*(-4.*Apa.*X10.*Gm.^3.*Km-12.*Apa.*X10.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+24.*pi.*Gm.*Km+9.*pi.*Km.^2))./(3.*Apa.*X10.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+9.*Apa.*X10.*Gm.*Km.^3+24.*pi.*Gm.*Km+9.*pi.*Km.^2)); % 10MPa Kd
Kd15 = abs( (Km.*(-4.*Apa.*X15.*Gm.^3.*Km-12.*Apa.*X15.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+24.*pi.*Gm.*Km+9.*pi.*Km.^2))./(3.*Apa.*X15.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+9.*Apa.*X15.*Gm.*Km.^3+24.*pi.*Gm.*Km+9.*pi.*Km.^2)); % 15MPa Kd
Kd20 = abs((Km.*(-4.*Apa.*X20.*Gm.^3.*Km-12.*Apa.*X20.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+24.*pi.*Gm.*Km+9.*pi.*Km.^2))./(3.*Apa.*X20.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+9.*Apa.*X20.*Gm.*Km.^3+24.*pi.*Gm.*Km+9.*pi.*Km.^2)); % 20MPa Kd
Kd25 = abs((Km.*(-4.*Apa.*X25.*Gm.^3.*Km-12.*Apa.*X25.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+24.*pi.*Gm.*Km+9.*pi.*Km.^2))./(3.*Apa.*X25.*Gm.^2.*Km.^2+16.*pi.*Gm.^2+9.*Apa.*X25.*Gm.*Km.^3+24.*pi.*Gm.*Km+9.*pi.*Km.^2)); % 25MPa Kd
Gd0 = -(512.*Gm.^4.*X0+2592.*Gm.^2.*Km.^2.*X0-600.*pi.*Apa.*Gm.^4+972.*Gm.*Km.^3.*X0+2112.*Gm.^3.*Km.*X0-4725.*pi.*Apa.*Gm.^2.*Km.^2-2025.*pi.*Apa.*Gm.*Km.^3-3150.*pi.*Apa.*Gm.^3.*Km+48.*pi.*Apa.*Gm.^4.*X0+243.*pi.*Apa.*Gm.*Km.^3.*X0+270.*pi.*Apa.*Gm.^3.*Km.*X0+459.*pi.*Apa.*Gm.^2.*Km.^2.*X0)./(3.*(256.*Gm.^3.*X0+216.*Km.^3.*X0+200.*pi.*Apa.*Gm.^3+675.*pi.*Apa.*Km.^3+816.*Gm.*Km.^2.*X0+896.*Gm.^2.*Km.*X0+1575.*pi.*Apa.*Gm.*Km.^2+1050.*pi.*Apa.*Gm.^2.*Km+24.*pi.*Apa.*Gm.^3.*X0+54.*pi.*Apa.*Km.^3.*X0+162.*pi.*Apa.*Gm.*Km.^2.*X0+120.*pi.*Apa.*Gm.^2.*Km.*X0)); % 0MPa Gd
Gd5 = -(512.*Gm.^4.*X5+2592.*Gm.^2.*Km.^2.*X5-600.*pi.*Apa.*Gm.^4+972.*Gm.*Km.^3.*X5+2112.*Gm.^3.*Km.*X5-4725.*pi.*Apa.*Gm.^2.*Km.^2-2025.*pi.*Apa.*Gm.*Km.^3-3150.*pi.*Apa.*Gm.^3.*Km+48.*pi.*Apa.*Gm.^4.*X5+243.*pi.*Apa.*Gm.*Km.^3.*X5+270.*pi.*Apa.*Gm.^3.*Km.*X5+459.*pi.*Apa.*Gm.^2.*Km.^2.*X5)./(3.*(256.*Gm.^3.*X5+216.*Km.^3.*X5+200.*pi.*Apa.*Gm.^3+675.*pi.*Apa.*Km.^3+816.*Gm.*Km.^2.*X5+896.*Gm.^2.*Km.*X5+1575.*pi.*Apa.*Gm.*Km.^2+1050.*pi.*Apa.*Gm.^2.*Km+24.*pi.*Apa.*Gm.^3.*X5+54.*pi.*Apa.*Km.^3.*X5+162.*pi.*Apa.*Gm.*Km.^2.*X5+120.*pi.*Apa.*Gm.^2.*Km.*X5)); % 5MPa Gd
Gd10 = -(512.*Gm.^4.*X10+2592.*Gm.^2.*Km.^2.*X10-600.*pi.*Apa.*Gm.^4+972.*Gm.*Km.^3.*X10+2112.*Gm.^3.*Km.*X10-4725.*pi.*Apa.*Gm.^2.*Km.^2-2025.*pi.*Apa.*Gm.*Km.^3-3150.*pi.*Apa.*Gm.^3.*Km+48.*pi.*Apa.*Gm.^4.*X10+243.*pi.*Apa.*Gm.*Km.^3.*X10+270.*pi.*Apa.*Gm.^3.*Km.*X10+459.*pi.*Apa.*Gm.^2.*Km.^2.*X10)./(3.*(256.*Gm.^3.*X10+216.*Km.^3.*X10+200.*pi.*Apa.*Gm.^3+675.*pi.*Apa.*Km.^3+816.*Gm.*Km.^2.*X10+896.*Gm.^2.*Km.*X10+1575.*pi.*Apa.*Gm.*Km.^2+1050.*pi.*Apa.*Gm.^2.*Km+24.*pi.*Apa.*Gm.^3.*X10+54.*pi.*Apa.*Km.^3.*X10+162.*pi.*Apa.*Gm.*Km.^2.*X10+120.*pi.*Apa.*Gm.^2.*Km.*X10)); % 10MPa Gd
Gd15 = -(512.*Gm.^4.*X15+2592.*Gm.^2.*Km.^2.*X15-600.*pi.*Apa.*Gm.^4+972.*Gm.*Km.^3.*X15+2112.*Gm.^3.*Km.*X15-4725.*pi.*Apa.*Gm.^2.*Km.^2-2025.*pi.*Apa.*Gm.*Km.^3-3150.*pi.*Apa.*Gm.^3.*Km+48.*pi.*Apa.*Gm.^4.*X15+243.*pi.*Apa.*Gm.*Km.^3.*X15+270.*pi.*Apa.*Gm.^3.*Km.*X15+459.*pi.*Apa.*Gm.^2.*Km.^2.*X15)./(3.*(256.*Gm.^3.*X15+216.*Km.^3.*X15+200.*pi.*Apa.*Gm.^3+675.*pi.*Apa.*Km.^3+816.*Gm.*Km.^2.*X15+896.*Gm.^2.*Km.*X15+1575.*pi.*Apa.*Gm.*Km.^2+1050.*pi.*Apa.*Gm.^2.*Km+24.*pi.*Apa.*Gm.^3.*X15+54.*pi.*Apa.*Km.^3.*X15+162.*pi.*Apa.*Gm.*Km.^2.*X15+120.*pi.*Apa.*Gm.^2.*Km.*X15)); % 15MPa Gd
Gd20 = -(512.*Gm.^4.*X20+2592.*Gm.^2.*Km.^2.*X20-600.*pi.*Apa.*Gm.^4+972.*Gm.*Km.^3.*X20+2112.*Gm.^3.*Km.*X20-4725.*pi.*Apa.*Gm.^2.*Km.^2-2025.*pi.*Apa.*Gm.*Km.^3-3150.*pi.*Apa.*Gm.^3.*Km+48.*pi.*Apa.*Gm.^4.*X20+243.*pi.*Apa.*Gm.*Km.^3.*X20+270.*pi.*Apa.*Gm.^3.*Km.*X20+459.*pi.*Apa.*Gm.^2.*Km.^2.*X20)./(3.*(256.*Gm.^3.*X20+216.*Km.^3.*X20+200.*pi.*Apa.*Gm.^3+675.*pi.*Apa.*Km.^3+816.*Gm.*Km.^2.*X20+896.*Gm.^2.*Km.*X20+1575.*pi.*Apa.*Gm.*Km.^2+1050.*pi.*Apa.*Gm.^2.*Km+24.*pi.*Apa.*Gm.^3.*X20+54.*pi.*Apa.*Km.^3.*X20+162.*pi.*Apa.*Gm.*Km.^2.*X20+120.*pi.*Apa.*Gm.^2.*Km.*X20)); % 20MPa Gd
Gd25 = -(512.*Gm.^4.*X25+2592.*Gm.^2.*Km.^2.*X25-600.*pi.*Apa.*Gm.^4+972.*Gm.*Km.^3.*X25+2112.*Gm.^3.*Km.*X25-4725.*pi.*Apa.*Gm.^2.*Km.^2-2025.*pi.*Apa.*Gm.*Km.^3-3150.*pi.*Apa.*Gm.^3.*Km+48.*pi.*Apa.*Gm.^4.*X25+243.*pi.*Apa.*Gm.*Km.^3.*X25+270.*pi.*Apa.*Gm.^3.*Km.*X25+459.*pi.*Apa.*Gm.^2.*Km.^2.*X25)./(3.*(256.*Gm.^3.*X25+216.*Km.^3.*X25+200.*pi.*Apa.*Gm.^3+675.*pi.*Apa.*Km.^3+816.*Gm.*Km.^2.*X25+896.*Gm.^2.*Km.*X25+1575.*pi.*Apa.*Gm.*Km.^2+1050.*pi.*Apa.*Gm.^2.*Km+24.*pi.*Apa.*Gm.^3.*X25+54.*pi.*Apa.*Km.^3.*X25+162.*pi.*Apa.*Gm.*Km.^2.*X25+120.*pi.*Apa.*Gm.^2.*Km.*X25)); % 25MPa Gd
D = {'Kd0','Gd0','Kd5','Gd5','Kd10','Gd10','Kd15','Gd15','Kd20','Gd20','Kd25','Gd25'};
d = [Kd0,Gd0,Kd5,Gd5,Kd10,Gd10,Kd15,Gd15,Kd20,Gd20,Kd25,Gd25];
xlswrite('G:./matlabdata./KT-model.xls',D,'Penny','A1')
xlswrite('G:./matlabdata./KT-model.xls',d,'Penny','A2')
matlabK-T模型函数包.zip_K-T函数包_K._孔隙_孔隙介质_孔隙介质模型
版权申诉
5星 · 超过95%的资源 168 浏览量
2022-09-15
00:32:47
上传
评论
收藏 2KB ZIP 举报
alvarocfc
- 粉丝: 111
- 资源: 1万+
最新资源
- rfid的应用场景.zip
- SpringBlade是一个由商业级项目升级优化而来的微服务开发平台
- mod-ldap-2.4.6-45.el7.centos.x86-64.rpm
- vbscript是什么,他的作用
- mod-ldap-2.4.6-67.el7.centos.x86-64.rpm
- mod-ldap-2.4.6-31.el7.centos.1.x86-64.rpm
- 电子竞技战队网站设计与实现(ssm+jsp+mysql)含万字详细文档+运行截图
- mod-ldap-2.4.6-40.el7.centos.x86-64.rpm
- ResNet50模型训练猫狗数据集
- mod-ldap-2.4.6-80.el7.centos.x86-64.rpm
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈