求怎么用matlab画出S与n的关系，最好附上matlabd的代码！！！？-盖德问答-化工人互助问答社区

共4个回答

岁月与世界，设备维修 2019-05-20回答

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

14 评论 0 举报

小九九，品质工程师 2019-05-20回答

你这是哪个大学在做热电

11 评论 0 举报

庆生，化工研发 2019-05-20回答

请问这几个式子是哪里的文献，我想引用一下

17 评论 0 举报

认栽，销售 2019-05-20回答

这段时间自己在花时间学习matlab，所以更新了我的代码，现在看上去要高级一些，共享给大家：

function = 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 ');
ylabel('S ');
legend(leg,'fontSize',14);
n = n';

展开

12 评论 0 举报