matlab脚本如下：
%下面定义参数
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
%上面的定义都没有问题，这个函数V(x)的图像可以做出来的，和文献上的吻合。
%接下来定义另外一个函数N。文章中介绍这个函数的时候，引入了另一个变量y，N是对y从0到正无穷的积分，此外N还是V的函数，即N(y,V(x))，见下式。在对y定积分过程中，V可以看做常数。这个定积分用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
%结果报错了，这个我不会解决，请帮忙解决一下，谢谢！ 返回小木虫查看更多