我现在正在做我的论文。我坚持使用积分。我的功能定义为
myfun =(exp(t*Q)*V*x)(j);
其中Q
和V
是一个矩阵(n * n),x
是一个元素为1的向量,然后在计算后我们得到该向量的第j个元素,然后我需要将函数与t积分。
我想quad
在matlab中使用。然而关键是它会报告内部矩阵的大小不同。因为这里的 A 不是一个数字?......
我怎样才能做到这一点。否则我只能针对 t 本身编写一个循环,这非常慢。
谢谢
您可以为此使用SUBSREF(不过,您仍然需要遍历所有j
's):
myfunOfT = @(t)(subsref(exp(t*Q)*V*x,struct('type','()','subs',j);
这将返回j
时间数组的第 th 个元素的值t
。