matlab 函数向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...

我正在将一些代码从MATLAB转换为python，并且正在努力获得一个函数，该函数以数组参数(作为参数)来使用Scipy进行集成。在

我已经将代码简化为一个在Scipy中产生相同错误的基本示例，而等效的MATLAB代码的功能与预期相同。在

我试图将长度为m的行向量参数参数参数和长度为N的列向量参数参数参数传递给被集成在另一个单独的积分参数上的函数，期望我的集成输出将具有MxN的形状。在

以下python代码会产生此错误：File "C:\Anaconda3\lib\site-packages\scipy\integrate\quadrature.py", line 196, in quadrature

if err < tol or err < rtol*abs(val):

ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

从

^{pr2}$

等效的MATLAB代码(输出MxN形状的结果)是c = [1, 2];

r = transpose([2, 1]);

out = integral(@(p)intgd(p, r, c), -pi/2, pi/2, 'ArrayValued', true);

function intgd = intgd(p, r, c)

c_bcr = repmat(c, length(r), 1);

r_bcc = repmat(r, 1, length(c));

A = ones(length(r), length(c));

s = sin(r_bcc) - sin(p);

A(s ~= 0) = sin(c_bcr(s ~= 0).*s(s ~= 0))./(c_bcr(s ~= 0).*s(s ~= 0));

intgd = 1./fun(p, c).^2.*(A.^2);

end

function d = fun(p, c)

p1 = zeros(1, length(c));

mask = sqrt(pi./(2*c)) < 1;

p1(mask) = acos(sqrt(pi./(2*c(mask))));

d = zeros(1, length(c));

mask = abs(p) <= p1;

d(mask) = 1./(pi./(2*c(mask).^2) + cos(p));

mask = and(abs(p) > p1, abs(p) <= pi/2);

d(mask) = 1./(pi./(2*c(mask).^2) + ((cos(p1(mask)) - cos(p))./2));

end

上面的MATLAB输出out是[6.58727018139280, 0.963083280848789;

6.78600314283299, 1.05994693990888]

我不确定scipy.integrate.quadrature是如何处理通过它的对象的维度的，但想法是它应该产生相同的MxN输出。在

我知道numpy有它自己的内置广播，通常可以避免这里所示的显式广播的需要，但我不确定它如何处理m=N的数组，如本例中所示，因此我将其保持为显式。任何关于这个次要问题的建议也将是受欢迎的。在

matlab 函数向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...相关推荐

python中quad_python scipy integrate.quad用法及代码示例
计算定积分. 使用Fortran库QUADPACK中的技术将func从a集成到b(可能是无限间隔). 参数: func:{function, scipy.LowLevelCallable}集成的Pyt ...
python中eps是什么意思_如何在中使用参数epsabsscipy.integrate.quad在Python中？
根据scipy手册,quad function要指定limit参数An upper bound on the number of subintervals used in the adaptive a ...
使用scipy.integrate模块中的quad和dblquad方法分别进行单重积分和多重积分
一.单重积分一般而言,使用求解微积分可以分为两大类:符号积分(即求出解析解)和数值积分(即求出数值解).在计算机的处理当中,数值解往往更有意义. quadquadquad函数是 sc ...
int main(int argc,char *argv[]),主函数的参数问题
主函数的参数在VC++中,我们有时候看到一些程序的主函数是带有参数的,虽然我们经常用到VC++但是对主函数的参数并不是很理解: int main(int argc,char *argv[]) 对于里 ...
Python自定义函数的创建、调用和函数的参数详解
这篇文章主要介绍了Python自定义函数的创建.调用和函数的参数.变量作用域等常见问题,需要的朋友可以参考下函数是组织好的,可重复使用的,用来实现单一,或相关联功能的代码段. 函数能提高应用的模块性 ...
Swift 烧脑体操（二） - 函数的参数
前言 \\ Swift 其实比 Objective-C 复杂很多,相对于出生于上世纪 80 年代的 Objective-C 来说,Swift 融入了大量新特性.这也使得我们学习掌握这门语言变得相对来说 ...
python中函数参数_Python函数的参数
本章将介绍函数中参数的用法,更多内容请参考:Python学习指南定义函数的时候,我们把参数的名字和位置确定下来,函数的接口定义就完成了.对于函数的调用者来说,只需要知道传递正确的参数,以及函数就返回 ...
函数局部变量和函数的参数在栈中的布局
#include <stdio.h> #include <iostream>using namespace std; void func(int p1, int p2, int ...
python函数的参数可以接收哪些类型的数据_python强势来袭-0015-函数中的参数-送礼开始...
接上一讲,这节我们专门介绍python函数的参数处理本节主要内容函数的形式参数和实际参数函数的参数定义规则函数参数的默认值函数的可变参数函数的关键字参数函数的命名关键字参数 1.3. 函 ...

matlab 函数向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...

matlab 函数向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...相关推荐

最新文章

热门文章

matlab 函数 向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...

matlab 函数 向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...相关推荐

最新文章

热门文章

matlab 函数向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...

matlab 函数向量参数,Scipy integrate（quad，quadration，nquad）不能集成向量参数化函数？等效函数（MATLAB works）...相关推荐