| 查看: 1003 | 回复: 1 | ||||
| 【悬赏金币】回答本帖问题,作者学渣进行中将赠送您 10 个金币 | ||||
学渣进行中新虫 (小有名气)
|
[求助]
超导三元凸包图 已有1人参与
|
|||
| 请问各位老师如何画超导的三元凸包图 |
» 猜你喜欢
人气不行了
已经有8人回复
filecode,4个jtjc了
已经有15人回复
今天基金会出结果吗?20260819
已经有17人回复
售一区SCI文章T0P,我:8O.551.O54,科目全,可十急
已经有3人回复
售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
时间戳今天,20号变了
已经有5人回复
2026国自然函评费到账
已经有9人回复
放榜前的不淡定
已经有14人回复
建议基金发布提前给出明确的时间点
已经有13人回复
今日不放榜?网传国自然预计 8 月 27 日可查结果
已经有20人回复
dxcharlary
专家顾问 (职业作家)
-

专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10507.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 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









回复此楼