在六角形区域中,随机采样2W个点,根据公式,计算费米面nesting,使用Orgin画图公式对多边形插值处理

主函数

delta = 0.01;
T = 0.0026;
N = 20000;
k=GetKPoints(N);
sus = zeros(N,1);
for i=1:N                    for j=1:NHk1=Hamiltonian_TaS2_k(k(j,1),k(j,2));[~,Ek1]=eig(Hk1);E1 = Ek1(1,1);f1 = (exp(E1/T)+1)^-1;Hk2=Hamiltonian_TaS2_k(k(j,1)+k(i,1),k(j,2)+k(i,2));[~,Ek2]=eig(Hk2);E2 = Ek2(1,1);f2 = (exp(E2/T)+1)^-1;sus(i,1) = sus(i,1) + f1*(1-f2)/(E2 - E1 + 1i*delta);end
end
plot3(k(:,1),k(:,2),real(sus),'.');
data = [k,real(sus)];

优化升级方案:

取更多的点N = 4W

计算下面平行四边形区域的布里渊区,均匀取点

二维费米面的nesting计算相关推荐

  1. Vasp二维材料单胞基本计算范例

    Vasp二维材料单胞基本计算范例,以InSe为例子. 总的来说,对于非磁性的简单二维材料体系,计算分为以下几步 1. 结构优化 2.静态计算 3.能带和态密度的计算 上面几个过程,实际需要改动的参数只 ...

  2. 基于二维矩阵的FFT计算原理

    FFT算法是信号处理领域最基本.最经典的算法,在工程实践中用处十分广泛,但是在一些对FFT点数要求较大或者计算FFT实时性要求较高的场合,按常规方式直接计算FFT难以满足工程实际的要求.本文针对长点数 ...

  3. 【C 语言】二级指针作为输入 ( 二维数组 | 二维数组内存大小计算 | 指针跳转步长问题 )

    文章目录 一.二维数组内存大小计算 二.二维数组内存大小意义 一.二维数组内存大小计算 给定一个二维数组 : 该 二维数组 中有 444 个 一维数组 , 每个一维数组有 101010 个 char ...

  4. 二维有限体积 matlab,二维有限体积法计算热传导及源码.pdf

    二维有限体积法计算热传导及源码 //#include "stdafx.h" #include #include #include #include #include using n ...

  5. 二维联合分布的密度函数计算

    二维联合分布的密度函数计算 @(概率论) 已知随机向量(X1,X2)(X_1,X_2)的概率密度为f(x1,x2)f(x_1,x_2),设Y1=2X1,Y2=13x2Y_1 = 2X_1,Y_2 = ...

  6. 按如下函数原型编程从键盘输入一个m行n列的二维数组,然后计算数组中元素的最大值及其所在的行列下标值。其中,m和n的值由用户键盘输入。已知m和n的值都不超过10。

    题目:按如下函数原型编程从键盘输入一个m行n列的二维数组,然后计算数组中元素的最大值及其所在的行列下标值.其中,m和n的值由用户键盘输入.已知m和n的值都不超过10. void InputArray( ...

  7. 练习题记录:求解距离矩阵,首先生成一百个二维坐标点,计算任意两个坐标点的距离

    碰到一道练习题:求解距离矩阵,首先生成一百个二维坐标点,计算任意两个坐标点的距离.我一开始看到这题目是一脸懵逼,然后上网找了一下,更是二脸懵逼.生成一百个二维坐标点,我首先就是想去生成一个1010的二 ...

  8. (C++)二维坐标系下,计算两点之间方位角

    (C++)二维坐标系下,计算两点之间方位角 思路 代码 思路 将其中点ori视为坐标原点,求点p相对其方位角 代码 double GlobalSettings::GetPointHAngle(Poin ...

  9. 一篇文章带你搞定二维插值的 MATLAB 计算

    前面已经学习了二维插值的基本概念:一篇文章带你认识数学建模中的二维插值 本篇文章主要实现使用MATLAB进行二维插值计算 文章目录 一.网格节点的插值计算 二.散点数据的插值计算 1. 示例 1 2. ...

最新文章

  1. 如何选择合适的损失函数,请看......
  2. C#正则表达式匹配字符串中的数字
  3. Chrome 打印PDF技巧
  4. python绘制曲线图-python绘制多个曲线的折线图
  5. Bootstrap系列 -- 34. 按钮下拉菜单
  6. CodeChef Cards, bags and coins [DP 泛型背包]
  7. CTFshow php特性 web127
  8. Java生成pgp密钥对_在Javascript中生成PGP密钥对,并使用加密的PGP私钥对文本进行签名...
  9. java线程死锁_Java线程死锁–案例研究
  10. android oat如何提取dex文件字节码,Android: 使用oatdump反编译oat文件
  11. Swift中文教程(十二) 下标
  12. sql怎么select中位数_怎么能避免写出慢SQL?
  13. 敏捷计划是否基于功能,是否更好?
  14. iTunes 10 选择自定义IPSW
  15. M1 Pro 和 M1 Max MacBook Pro,我们该选择哪款MacBook
  16. matlab编写扫雷,MATLAB版本的扫雷小游戏
  17. 天池比赛——用户情感可视化分析
  18. 开发初探-前端日记 — html、css的使用,模拟制作百度登陆界面
  19. requests使用splash
  20. 生成对抗网络GAN原理——浅显易懂,个人收藏

热门文章

  1. 男人的这十大幼稚行为 女人最讨厌
  2. 上月技术关注:Google大表
  3. Mac OS 系统 latex 转 word
  4. 【音视频开发之起始】
  5. 大鱼海棠国漫主题HTML前端网页设计成品预览与讲解DIV CSS JS
  6. 金额达 500 亿美元,AMD 完成收购赛灵思
  7. 这些酷酷的文字特效图片,制作真的很简单!
  8. 深圳公司法人变更的条件是什么?
  9. 科大讯飞2014届实习生招聘笔试题
  10. [云炬创业基础笔记]识别创业机会的方法之一头脑风暴法