Question

MATLAB 自定义函数有两个参数,固定一个,对另外一个进行数值积分

函数较复杂

Answer 1

举个这种情况的例子：

f = @(x,y)x+y;
quad(@(x)f(x,1), 0, 1)

追问

函数里面有BesselK函数，其参数包含变量

追答

能把具体要积分的函数给出来吗？有Bessel函数怎么了？

追问

Ee=0.8; I=0.3; B=1.2;
gama=1957*Ee;
Ec=0.665*Ee^2*B;
y=Energy/Ec;
a=gama^2*phy.^2;
b=1+a;
kexi=1/2*y.*b.^(3/2);
FluxDensity=1.327E13*Ee^2*I*b.^2.*(BesselK(2/3,kexi).^2+a./b.*BesselK(1/3,kexi).^2);

追答

要积分的函数是FluxDensity？？ 变量是Energy，phy？？如果是的话，看如下示例：

1）新建个函数文件，

function FluxDensity = f( Energy, phy )
    Ee=0.8; I=0.3; B=1.2;
    gama=1957*Ee;
    Ec=0.665*Ee^2*B;
    y=Energy/Ec;
    a=gama^2*phy.^2;
    b=1+a;
    kexi=1/2*y.*b.^(3/2);
    FluxDensity=1.327E13*Ee^2*I*b.^2.*(BesselK(2/3,kexi).^2+a./b.*BesselK(1/3,kexi).^2);
end

2）固定phy为2，求Energy积分（0到1）：

quad(@(E)f(E,2), 0, 1)

追问

clear
clc
Energy=1:10:1000;
for i=1:length(Energy)
F(i)=quadgk(@(phy)BMFluxDensity(phy,Energy(i)),-inf,inf);
end
被积分函数没问题吧

追答

你又没给出你被积函数怎么写的？你问问题太没头没尾啦，还要别人去猜你要问什么？不清楚的话发qq邮件给我吧。

追问

你的QQ是多少

追答

我的百度资料里面有，8 5 2 6 1 0 1 8 3