matlab - 如何使用 MATLAB 对复积分进行数值求解

Question

我正在尝试解决这个等式：

K=sqrt((R*T)/(4*pi*lambda))*从-inf到exp的inf的积分(-((lambda+F*neta)/R*Tx)^2*R*T/4 *lambda)/exp(x)+1 关于 x

其中，neta 是从 0 到 1 的区间，其他符号（R、T、F、lambda 和 pi）具有恒定值。

我尝试使用这些代码：

代码 1

clear all;
close all;
clc;
F = 96485.34;
R = 8.3145;
T = 298.15;
lambda = 0.2;
neta=0:0.1:1;
pi=3.1415;
f=@(x) exp(-((lambda+F*neta)/R*T-x).^2*R*T/4*lambda)/(exp(x)+1);
Q=integral(f,-inf,inf);
k= sqrt((R*T)/(4*pi*lambda)).*Q

代码 2

clear all;
close all;
clc;
F = 96485.34;
R = 8.3145;
T = 298.15;
lambda = 0.2;
neta=0:0.1:1;
pi=3.1415;
x= 0:100;
f(x)=exp(-((lambda+F*neta)/R*T-x).^2*R*T/4*lambda)/(exp(x)+1);
q=quadl('f', 0, 100);
k= sqrt((R*T)/(4*pi*lambda)).*q

但这些代码返回我不知道要解决的错误。有人能帮助我吗？

谢谢

Answer 1

问题在于neta。如果它是一个标量值，那么代码就可以了。但是您如何建议在 neta 区间上针对不同的变量进行积分？这意味着什么？我尝试将 neta 的各个值从 0 到 1 并为每个值计算一个 Q 值，但每个 Q 返回 0。