在MATLAB中集成功能

时间:2016-07-15 01:34:27

标签: matlab integration

我是Matlab的新手。我想整合一个功能。我尝试使用int()和integrate(),但它们都会给我带来问题 - 没有足够的参数或其他不同的错误,我尝试了很多文档组合。我目前的代码如下,我希望能够将数字p和q传递给res并获得数值结果:

syms x;
w = 1;
hbar = 1.054571800*10^(-34);
k = (w/(pi*hbar))^(1/4);
e = @(q) (w*q/hbar)^(1/2);
waveF = @(q) k*exp(-feval(e,q)*feval(e,q)*1/2.0)*1.0/1;
func = @(p,q) waveF(q-x/2)*waveF(q+x/2)*exp(1i*p*x/(hbar));

res = @(p,q) int(func(p,q), x = -Inf..Inf);

目前“x =”表示为错误,但根据文档似乎没有问题。

感谢。

1 个答案:

答案 0 :(得分:4)

您正在使用anonymous functionsSymbolic Toolbox一致,并错误地使用MuPAD int版本的Symbolic int,这实际上会产生错误,当您需要{ {3}}

虽然将匿名函数与符号调用混合在一起并不违法,但我认为通过一次坚持一种计算范式会更好:

  1. 使用符号函数的纯符号版本:

    syms x p q e(q) waveF(q) func(p,q) res(p,q);
    w         = sym(1);
    hbar      = sym('1.054571800E-34');
    k         = (w/(pi*hbar))^(1/4);
    e(q)      = sqrt(w*q/hbar);
    waveF(q)  = k*exp(-e(q)^2/2);
    func(p,q) = waveF(q-x/2)*waveF(q+x/2)*exp(1i*p*x/(hbar));
    res(p,q)  = int(func(p,q),x,-Inf,Inf);
    

    我将hbar的值包装在引号中以强制使用提供的值,而不是在计算过程中被强制使用的最近的理性表示。

    1. 使用匿名函数和数字integral函数的纯数字版本:

      w     = 1;
      hbar  = 1.054571800E-34;
      k     = (w/(pi*hbar)).^(1/4);
      e     = @(q) sqrt(w*q/hbar);
      waveF = @(q) k*exp(-e(q).^2/2);
      func  = @(p,q,x) waveF(q-x/2).*waveF(q+x/2).*exp(1i*p*x/hbar);
      res   = @(p,q) integral(@(x) func(p,q,x),-Inf,Inf);
      
    2. 在评估NaN时,这两个版本都会生成res,但这可能是被积函数的缺点。但是,上述脚本背后的功能形式和想法不会随着不同的集成而发生变化。