我现在正在写论文。我坚持使用积分。我的函数定义为
myfun =(exp(t*Q)*V*x)(j);
其中Q
和V
是一个矩阵(n*n),x
是一个元素为1的向量,经过计算我们得到了这个向量的j_th元素,然后我需要对t进行函数积分。
我想在matlab中使用quad
。然而,重点是它将报告内部矩阵不是相同的大小。因为这里的A不是一个数字?
我该怎么做呢。否则,我只能针对t本身编写循环,这是非常慢的。
谢谢
发布于 2011-08-19 11:56:34
为此,您可以使用SUBSREF (不过,您仍然需要遍历所有的j
):
myfunOfT = @(t)(subsref(exp(t*Q)*V*x,struct('type','()','subs',j);
它返回在时间t
时数组的第j
个元素的值。
https://stackoverflow.com/questions/7108090
复制