功能定义为积分：trapz？

Question

我想计算x，y：p（x，y）的函数p，它由x和y的另一个函数的积分定义：

p（x，y）= v（x，y）dx的积分（不定）

现在，如果我在统一网格SELECT s.* FROM state_history s JOIN( SELECT entity_id, MAX(state_added_time) AS maxTime FROM state_history WHERE state_added_time < '2015-05-15 14:25:00' GROUP BY entity_id) tmp ON tmp.entity_id = s.entity_id AND tmp.maxTime = s.state_added_time; 和p上有一个表示x的矩阵，我该如何构造函数y（应该是一个矩阵） 。

显然，如果我使用p进行整合，我会获得一个向量：

trapz

给出x=linspace(-1,1,10); v=magic(10); p=trapz(x,v); size(p) 而不是10 10。

Answer 1

假设

1. p 的公式实际上定义如下：

p （ x ， y ）=积分[ x ₀ 。 。 x ] v （ξ， y ）d ξ;

1. 整合网格x×y是统一的，其标准为dx×dy;

2. 矩阵v的行有x常量，第一行对应最小的x;

那么积分是：

p = cumtrapz(v) * dx;