请教：数值积分的选择

beyondhxf

我现在需要做一组积分：

1. R=[f1(x);f2(x);...........];
   
2. for i=1:...
   
3. Ra(i)=double(int(R(i)*cos(x),-160*pi/180,0));
   
4. end
   

复制代码

上面用的是int函数，然而当R中元素过于复杂的时候，我想用数值积分。
比如quad就很好。
然而我遇到的问题是：调用quad时候，必须写出函数，如quad(@sin,-pi,0);
而我目前只能写出R(1,,j)这样的元素
因此请教如何变通可以使用quad这样的数值积分函数？
或者使用其他的数值积分函数。

[ 本帖最后由 beyondhxf 于 2010-6-8 21:29 编辑 ]

beyondhxf

虽然各位大侠不屑于回答我这个愚蠢的问题，但现在问题解决了，我还是回来反馈一下，以帮助后面有需要的人。
解决格式为：
quad(inline(vectorize(char(R(i)*R(j)))),-10,100000)

rocwoods

7.0以后尽量不要用inline，inline效率不高。楼主MATLAB版本如果比较新的话可以用str2func构造被积函数句柄。
不过可以在得到R的时候用cell数组存储各个积分表达式的函数句柄，譬如：

1. 
   
2. >> R = {@(x) x,@(x) x.^2};
   
3. quadl(R{2},0,2)
   
4. 
   
5. ans =
   
6. 
   
7.     2.6667
   

复制代码

[ 本帖最后由 rocwoods 于 2010-6-10 10:03 编辑 ]

beyondhxf

非常感谢 热心的回答。