| 查看: 2420 | 回复: 4 | |||
| 【悬赏金币】回答本帖问题，作者反时差将赠送您 5 个金币 | |||
[求助]
求怎么用matlab画出S与n的关系，最好附上matlabd的代码！！！已有1人参与
|
|||
该怎么输入函数，费米积分呢？求各位大佬帮助，感激不尽！！！ s与n的关系-SPB模型.png  已知参数.png |
回复此楼» 收录本帖的淘贴专辑推荐
 热电材料知识 |
» 猜你喜欢
急求博导收留
已经有15人回复
-
蹲一个男朋友
已经有35人回复
-
去年申请基金的评审意见是ChatGPT在国内是禁止的，研究方案中有使用ChatGPT不合理
已经有4人回复
-
申请基金代表性成果
已经有20人回复
-
word转成pdf之后公式里面的字体变了，正文字体没变。
已经有7人回复
-
基金委也是打工的
已经有33人回复
-
被导师抓着一个褥羊毛该咋办
已经有25人回复
-
青椒的未来在哪里？
已经有18人回复
-
审不上青基又非升即走的青椒 和 牢里踩缝纫机的犯人哪个活的更舒服一点？
已经有8人回复
-
绍兴文理学院怎么样？有没有坑啊
已经有23人回复
穿越千年
版主 (著名写手)
寂寞读书的小和尚
- FPI: 2
- 应助: 238 (大学生)
- 贵宾: 0.515
- 金币: 9994.5
- 散金: 8750
- 红花: 119
- 沙发: 10
- 帖子: 2398
- 在线: 1369.1小时
- 虫号: 1393528
- 注册: 2011-09-07
- 性别: GG
- 专业: 半导体材料
- 管辖: 功能材料
【答案】应助回帖
|
Kb=138.06488; e=1.6021892; K=Kb/e; n=59;%n 10^18/cm^3 Seebeck=139;% uV/K h=6.6260755; me=9.10938215; Fn=@(n,yita)quadgk(@(x)x.^n./(1+exp(x-yita)),0,inf); f1=@(yita)(2*Fn(1,yita)/Fn(0,yita)-yita)-Seebeck/K;%Seebeck uV/K yita=fsolve(f1,0); m=(n^(2/3)*h^2)/(4*pi*Fn(1/2,yita))^(2/3)/(2*Kb*300); m2=(m/me)*10^4; %% normalized result 你可以试试用这段代码计算有效质量，当然你自己给一个有效质量，然后计算不同浓度下理论的Seebeck可以用下面这个。我的代码也很初级，不能自动计算曲线，得手动一个一个输入。 Kb=138.06488; e=1.6021892; K=Kb/e; n=35;%n input 10^18/cm^3 h=6.6260755; me=9.10938215; m2=1.6;% input m=(m2*me)*10^-4; Fn=@(n,yita)quadgk(@(x)x.^n./(1+exp(x-yita)),0,inf); f1=@(yita)4*pi*(2*Kb*300*m)^(3/2)/h^3*Fn(1/2,yita)-n; yita=fsolve(f1,0); S=K*(2*Fn(1,yita)/Fn(0,yita)-yita);%Seebeck uV/K |
2楼2018-05-25 21:10:30
旭阳砜
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 2429.1
- 散金: 323
- 红花: 1
- 沙发: 1
- 帖子: 1129
- 在线: 135.7小时
- 虫号: 4114227
- 注册: 2015-10-01
- 专业: 半导体材料
3楼2018-09-17 01:26:48
穿越千年
版主 (著名写手)
寂寞读书的小和尚
- FPI: 2
- 应助: 238 (大学生)
- 贵宾: 0.515
- 金币: 9994.5
- 散金: 8750
- 红花: 119
- 沙发: 10
- 帖子: 2398
- 在线: 1369.1小时
- 虫号: 1393528
- 注册: 2011-09-07
- 性别: GG
- 专业: 半导体材料
- 管辖: 功能材料
【答案】应助回帖
|
这段时间自己在花时间学习matlab，所以更新了我的代码，现在看上去要高级一些，共享给大家： function [n,S] = PisarenkoPlot(n1,n2,n_step,m2) Kb=138.06488; e=1.6021892; K=Kb/e; n = n1:n_step:n2; h=6.6260755; me=9.10938215; T=300;% unit in K leg = cell (length(m2),1); for j=1:length(m2) mm2 = m2(j); m=(mm2*me)*10^-4; Fn=@(no,yita)quadgk(@(x)x.^no./(1+exp(x-yita)),0,inf); for i=1:length(n) f1=@(yita)4*pi*(2*Kb*T*m)^(3/2)/h^3*Fn(1/2,yita)-n(i); yita(i)=fsolve(f1,0); S(i,j)=K*(2*Fn(1,yita(i))/Fn(0,yita(i))-yita(i)); end plot(n,S(:,j),'LineWidth',1) hold on leg{j} = num2str(m2(j)); end hold off xlabel('n [10^{18} cm^{-3}]'); ylabel('S [\mu V \cdot K^{-1}]'); legend(leg,'fontSize',14); n = n'; |
4楼2018-09-26 15:58:24
吴一i
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 308.1
- 散金: 40
- 帖子: 182
- 在线: 15.5小时
- 虫号: 6362357
- 注册: 2017-04-21
- 性别: GG
- 专业: 无机非金属材料
5楼2019-02-16 09:34:58
50