I have an Anonymous Function and want to integrate it

For a simple example, I have an Anonymous function A = @ (x) [x 2*x;3*x 4*x], and I want to integrate it from 0 to 1. If I use QUADV it will integrate the entire matrix A. However, it is not efficient for very large problem. I am wondering if there is a way that I can only integrate some specific terms, say only A(1,1).
The background is that I have a very large 80000 by 80000 matrix, and only a few terms need to be calculated.

回答(2 个)

Can you not just call:
A( myMatrix(1,1) );
to evaluate just the first element of your matrix? Or
A( myMatrix( xStart:xEnd, yStart:yEnd ) );
to evaluate a rectangular subset of your matrix? Or any other form of indexing (especially linear indexing) to choose more specifically.

2 个评论

Thanks Adam,
what is myMatrix?
Could u give me more details?
now I have
A = @ (x) [x 2*x;3*x 4*x]; >> quadv(A,0,1)
ans =
0.5000 1.0000
1.5000 2.0000
Could u just integrate A(1,1)?
I'm not familiar with using quadv so its usage does complicate matters a little though I would think it should be perfectly possible.
'myMatrix' is just a place holder for whatever your matrix is which I assumed was the input x though scanning through quadv it seems your function expects only a scalar input as x so I'm not sure where your matrix comes into it.

请先登录,再进行评论。

clear all
close all
clc;
syms x
A=[x 2*x;3*x 4*x];
int(A(1,1),0,1)

3 个评论

Thanks M.A.!
I think int is not as efficient as quad for very large problem, so do u know if I can use numerical integration to do it?
for numerical integration you can use quad , quadl or double(int())
syms x
A=[x 2*x;3*x 4*x];
double(int(A(1,1),0,1))
Quad, quadl, and quadv are all deprecated. The function is integral(), and for array valued problems, use the 'ArrayValued',true option.
It is easy to integrate selected terms in a way that unfortunately evaluates the entire matrix, anyway. The trick is to avoid that. I'm not sure how to do it in the general context.

请先登录,再进行评论。

类别

帮助中心File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

提问:

2014-11-18

评论:

2014-11-20

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by