求助：电势与量子电容

RSS

各位老师，请问下在量子电 容与电势的关系图中，为什么我拟合出来的图所有的体系都是一样的趋势呢？能否帮我看看是代码哪里出了问题嘛？以下是我的代码：cd ('F:\桌面')
% 导入态密度数据文件
data = readtable('DOSCAR 3.dat');
energy = data(:, 1).Variables; % 能量数据
density = data(:, 2).Variables; % 态密度数据
phi = linspace(-1, 1, 100); % 生成100个从-1到1的电势值

% 定义量子电容函数
quantum_capacitance = @(x, a, b, phiG, T) a * exp(-b*(x-phiG)/(T+1e-10));

% 定义拟合模型
model = fittype(@(a, b, phiG, T, x) quantum_capacitance(x, a, b, phiG, T));
% 进行拟合
fitresult = fit(energy, density, model, 'StartPoint', [-0.5, 0.5, 0, 300]);
a = fitresult.a; % 拟合参数a
b = fitresult.b; % 拟合参数b
phiG = fitresult.phiG; % 拟合参数phiG
T = 300;

% 计算电势-量子电容
C = quantum_capacitance(phi, a, b, phiG, T);

% 绘制电势-量子电容关系图
figure;
plot(phi, C);
xlabel('phi (V)')
ylabel('C (μF/cm^2)');

幻想银河

请问楼主问题解决了吗  我也在计算这个量子电容 也遇到一些问题 是否可以交流一下

W1312319875

求楼主赐教，最近也在做量子电容的计算，方便加个联系方式学习一下吗