matlab腳本如下：
%下面定義參數(shù)
e0=8.85e-12;
eb=5;
w=10;
kB=1.38e-23;
a1=-8.34e7;
g11=2e-10;
s11=8.6e-12;
s12=-2.8e-12;
a12=3.233e8;
a112=8.503e8;
a=2*(a1+a12+a112);
Q11=0.0812;
Q12=-0.0295;
h_=6.63*-34/(2*pi);
F12=0.5e-10;
g=g11+2*F12*F12/(s11+s12);
L=sqrt(eb*e0*g/(eb*e0*a)+1);
x=-50:0.01:50;
A=-F12*eb*e0*(Q11+Q12)/((eb*a+1)*(s11+s12)*w);
E=1/(e0*eb)*2*((+L/w)*sinh(w/L)-cosh(w/L))*L*exp(x/L).*(x>=-50&x<=-w);
F=(1/(e0*eb)*2*(x.*x/(2*w)-exp(-w/L)*(L/w+1)*L*cosh(x/L))+1/(e0*eb)*2*((+L/w)*sinh(w/L)-cosh(w/L))*L*exp(-w/L)-1/(e0*eb)*2*(w.*w/(2*w)-exp(-w/L)*(L/w+1)*L*cosh(w/L))).*(x>-w&x<w);
G=1/(e0*eb)*2*((-L/w)*sinh(w/L)+cosh(w/L))*(-L)*exp(-x/L).*(x>=w&x<=50);
V=A*(E+F+G);
%plot(x,V,'r','linewidth',2)
%axis([-50 50 -0.5e-10 0.5e-10])
%grid on
%上面的定義都沒有問題，這個函數(shù)V(x)的圖像可以做出來的，和文獻(xiàn)上的吻合。
%接下來定義另外一個函數(shù)N。文章中介紹這個函數(shù)的時候，引入了另一個變量y，N是對y從0到正無窮的積分，此外N還是V的函數(shù)，即N(y,V(x))，見下式。在對y定積分過程中，V可以看做常數(shù)。這個定積分用mathematica積不出來，所以我想著用matlab直接作出N(x)圖像來。
syms y
N=sqrt(2*0.05^3*y)/((2*pi^2*h_^3)*(1+exp((y-e0*V)/(kB*300))));
plot(x,[f(y),0,infinity],'r','linewidth',2)
axis([-50 50 -0.5e-10 0.5e-10])
grid on
%結(jié)果報錯了，這個我不會解決，請幫忙解決一下，謝謝！ 返回小木蟲查看更多