版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者学渣进行中将赠送您 10 个金币

学渣进行中

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 279.6
散金: 20
帖子: 95
在线: 14.3小时
虫号: 29559693
注册: 2022-05-03
专业: 凝聚态物性I:结构、力学和

[求助] 超导三元凸包图已有1人参与

请问各位老师如何画超导的三元凸包图

回复此楼

» 猜你喜欢

1楼 2024-04-17 22:17:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

专家经验: +83
1ST强帖: 2
应助: 251 (大学生)
金币: 10417.6
散金: 155
红花: 65
帖子: 3137
在线: 757.4小时
虫号: 936718
注册: 2010-01-04
专业: 凝聚态物性I:结构、力学和
管辖: 第一性原理

【答案】应助回帖

分享一个我自己写的matlab代码：
clear
clc
load tmp.txt
a=tmp;
nn=0;
for i=1:length(a)
if a(i,4)<0.01
      nn=nn+1;
      b(nn,1)=a(i,1);
      b(nn,2)=a(i,2);
      b(nn,3)=a(i,4);
end
end
DT = delaunayTriangulation(b);
[p,v]=convexHull(DT);
%trisurf(p,DT.Points(:,1),DT.Points(:,2),DT.Points(:,3),...
%    'FaceColor','cyan');
k=convhull(b);
for i=1:length(p)
for j=1:3
xx(1:2)=b(p(i,j),1:2);
x(j,1:2)=[xx(1)+xx(2)/2 xx(2)/2*sqrt(3)];
end
x(4,1:2)=x(1,1:2);
plot(x(:,1),x(:,2));
hold on
end
hold off

如果用了上述代码，请引用
HF Li et al, Ultrahigh-pressure magnesium hydrosilicates as reservoirs of water in early Earth. Phys. Rev. Lett. 128 (3), 035703 (2022)
J Hou et al, H4O2+ ion stabilized by pressure. Phys. Rev. B 109 (17), 174102 (2024)
谢谢

赞一下

回复此楼

2楼2024-05-09 00:29:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主学渣进行中的主题更新