关于半导体费米积分的积分,求助matlab大神 本人matbal很弱，最近做一个关于费米积分的数值计算

可以这样直接写。

u=1E=0.8n=(u-E)/0.026

F0=quadgk(@(x) x.^(-0.5)./(1+exp(x-n)),0,inf)

运行结果

如还有问题，可以进一步交流。

半导体电阻率的多种测量方法应用与注意事项依据掺杂水平的不同，半导体材料可能有很高的电阻率。有几种因素可能会使测量这些材料电阻率的工作复杂化，其中包括与材料实现良好接触的问题。已经设计出专门的探头来测量半导体晶圆片和半导体棒的电阻率。这些探头通常使用硬金属，如钨来制作，并将其磨成一个探针。在这种情况下接触电阻非常高，所以应当使用四点同线（collinear）探针或者四线隔离探针。其中两个探针提供恒定的电流，而另外两个探针测量一部分样品上的电压降。利用被测电阻的几何尺寸因素，就可以计算出电阻率。 看起来这种测量可能是直截了当的，但还是有一些问题需要加以注意。对探针和测量引线进行良好的屏蔽是非常重要的，其理由有三点： 1 电路涉及高阻抗，所以容易受到静电干扰。 2 半导体材料上的接触点能够产生二极管效应，从而对吸收的信号进行整流，并将其作为直流偏置显示出来。 3 材料通常对光敏感。 四探针技术 四点同线探针电阻率测量技术用四个等距离的探针和未知电阻的材料接触。此探针阵列放在材料的中央。图4-25是这种技术的图示。

已知的电流流过两个外部的探针，而用两个内部的探针测量电压。电阻率计算如下： 其中：V = 测量出的电压（伏特） I = 所加的电流（安培） t = 晶圆片的厚度（厘米） k = 由探头与晶圆片直径之比和晶圆片厚度与探头分开距离之比决定的修正因数。

如图4-26所示，更实际的电路还包括每个探针的接触电阻和分布电阻（r1到r4）、电流源和电压表从其LO端到大地的有限的电阻（RC和RV）和电压表的输入电阻（RIN）。依据材料的不同，接触电阻（r）可能会比被测电阻（R2）高300倍或更高。这就要求电流源具有比通常期望数值高得多的钳位电压，而电压表则必须具有高得多的输入电阻。

电流源不是与大地完全隔离的，所以当样品的电阻增加时，就更需要使用差分式静电计。存在问题的原因是样品可能具有非常高的电阻（108Ω或更高），此数值和静电计电压表的绝缘电阻（输入LO端到壳地，RV）具有相同的数量级。如图4-26所示，这样就会有交流电流从电流源的LO端，经过样品，流到电压表的LO端，再流回地。当电压表测量探头2和探头3之间的电压降时，该交流电流在r3上产生的电压降就会引起错误的结果。

使用两台静电计就解决了这个问题，如图4-27所示。电压表将读出两个静电计的缓冲输出之间的差值，该值等于R2上的电压。数值r1、r2、r3 和r4代表探头与样品材料接触的电阻。单位增益缓冲器具有很高的输入阻抗，所以几乎没有共模电流流过r3，于是可以很容易地计算出R2的数值。该缓冲器可以是一对JFET运算放大器或者是两个具有单位增益输出的静电计。

为了避免泄漏电流，使用隔离的或者带保护的探头与样品接触。电流源应当处于保护模式。

根据半导体理论，一般半导体材料的电阻率 和绝对温度 之间的关系为

（1—1）

式中a与b对于同一种半导体材料为常量，其数值与材料的物理性质有关。因而热敏电阻的电阻值 可以根据电阻定律写为

（1—2）

式中 为两电极间距离， 为热敏电阻的横截面， 。

对某一特定电阻而言， 与b均为常数，用实验方法可以测定。为了便于数据处理，将上式两边取对数，则有

（1—3）

上式表明 与 呈线性关系，在实验中只要测得各个温度 以及对应的电阻 的值，

以 为横坐标， 为纵坐标作图，则得到的图线应为直线，可用图解法、计算法或最小二乘法求出参数 a、b的值。

热敏电阻的电阻温度系数 下式给出

（1—4）

从上述方法求得的b值和室温代入式（1—4），就可以算出室温时的电阻温度系数。

热敏电阻 在不同温度时的电阻值，可由非平衡直流电桥测得。非平衡直流电桥原理图如右图所示，B、D之间为一负载电阻 ，只要测出 ，就可以得到 值。

当负载电阻 → ，即电桥输出处于开

路状态时， =0，仅有电压输出，用 表示，当 时，电桥输出 =0，即电桥处于平衡状态。为了测量的准确性，在测量之前，电桥必须预调平衡，这样可使输出电压只与某一臂的电阻变化有关。

若R1、R2、R3固定，R4为待测电阻，R4 = RX，则当R4→R4+△R时，因电桥不平衡而产生的电压输出为：

（1—5）

在测量MF51型热敏电阻时，非平衡直流电桥所采用的是立式电桥 ， ，且 ，则

（1—6）

式中R和 均为预调平衡后的电阻值，测得电压输出后，通过式（1—6）运算可得△R，从而求的 =R4+△R。

3、热敏电阻的电阻温度特性研究

根据表一中MF51型半导体热敏电阻（2.7kΩ）之电阻~温度特性研究桥式电路，并设计各臂电阻R和 的值，以确保电压输出不会溢出（本实验 =1000.0Ω， =4323.0Ω）。

根据桥式，预调平衡，将“功能转换”开关旋至“电压“位置，按下G、B开关，打开实验加热装置升温，每隔2℃测1个值，并将测量数据列表（表二）。

表一 MF51型半导体热敏电阻（2.7kΩ）之电阻～温度特性

温度℃ 25 30 35 40 45 50 55 60 65

电阻Ω 2700 2225 1870 1573 1341 1160 1000 868 748

表二 非平衡电桥电压输出形式（立式）测量MF51型热敏电阻的数据

i 1 2 3 4 5 6 7 8 9 10

温度t℃ 10.4 12.4 14.4 16.4 18.4 20.4 22.4 24.4 26.4 28.4

热力学T K 283.4 285.4 287.4 289.4 291.4 293.4 295.4 297.4 299.4 301.4

0.0 -12.5 -27.0 -42.5 -58.4 -74.8 -91.6 -107.8 -126.4 -144.4

0.0 -259.2 -529.9 -789 -1027.2 -124.8 -1451.9 -1630.1 -1815.4 -1977.9

4323.0 4063.8 3793.1 3534.0 3295.8 3074.9 2871.1 2692.9 2507.6 2345.1

根据表二所得的数据作出 ～ 图，如右图所示。运用最小二乘法计算所得的线性方程为 ，即MF51型半导体热敏电阻（2.7kΩ）的电阻～温度特性的数学表达式为 。

4、实验结果误差

通过实验所得的MF51型半导体热敏电阻的电阻—温度特性的数学表达式为 。根据所得表达式计算出热敏电阻的电阻～温度特性的测量值，与表一所给出的参考值有较好的一致性，如下表所示：

表三 实验结果比较

温度℃ 25 30 35 40 45 50 55 60 65

参考值RT Ω 2700 2225 1870 1573 1341 1160 1000 868 748

测量值RT Ω 2720 2238 1900 1587 1408 1232 1074 939 823

相对误差 % 0.74 0.58 1.60 0.89 4.99 6.20 7.40 8.18 10.00

从上述结果来看，基本在实验误差范围之内。但我们可以清楚的发现，随着温度的升高，电阻值变小，但是相对误差却在变大，这主要是由内热效应而引起的。

---