讨论一下两种quad积分方式

vansen

本人最近刚学MATLAB，在学习ik760218转帖的MATLAB 小技巧的时候发现了一下问题，还请大虾们指点一二，原帖地址：http://forum.vibunion.com/thread-2588-1-3.html

方法一：

1. 
   
2. m=5;
   
3. a = m;
   
4. b = m + 1;
   
5. K = m.^2;
   
6. y = subs('sin(K*x)/x', 'K', sym(K, 'd') )
   
7. quad( inline(vectorize(char(y))), a, b)
   
8. 
   
9. y =
   
10. 
    
11. sin(25.*x)/x
    
12. 
    
13. 
    
14. ans =
    
15.     0.0016
    

复制代码

方法二：

1. 
   
2. m=5;
   
3. a = m;
   
4. b = m + 1;
   
5. K = m.^2;
   
6. sym x;
   
7. l=sin(K*x)/x;
   
8. y = subs('l', 'K', sym(K, 'd') );
   
9. quad( inline(vectorize(char(y))), a, b)
   
10. ans =
    
11.     5.5000
    

复制代码

两个的计算结果为什么相差这么大呢？百思不得其解，还请各位大大指点！

sigma665

第2种是错的
y = subs('l', 'K', sym(K, 'd') );
无法自动识别l=sin(K*x)/x;
因为相当于重新定义了一个字符型l，把原来的覆盖了

vansen

谢谢sigma665！
要是我被积的函数很复杂，比如不止3个参数（ a,b,k )是在运行时确定的，可能有10+个，难道对于每一个不确定的参数都要用y = subs('sin(K*x)/x', 'K', sym(K, 'd') )运行一下，再求积分吗？

编了个有很多参数的程序，运行时报Not enough inputs to inline function.

1. 
   
2. a=1;
   
3. b=2;
   
4. c=3;
   
5. x=a+b;
   
6. y=a+c;
   
7. z=a*b;
   
8. m = subs('x*y*z*k', 'x', sym(x, 'd'));
   
9. m = subs('x*y*z*k', 'y', sym(y, 'd'));
   
10. m = subs('x*y*z*k', 'z', sym(z, 'd'));
    
11. quad( inline(vectorize(char(m))), a, b)
    

复制代码

??? Error using ==> inline.subsref at 14
Not enough inputs to inline function.
Error in ==> quad at 71
y = f(x, varargin{:});

[ 本帖最后由 vansen 于 2008-7-24 19:39 编辑 ]

rocwoods

1. 举个例子：
   
2. 计算a*sin(x)+b*sqrt(log(x))+c^2*d这个表达式在x=10,20的积分。
   
3. %经过大量计算得到a,b,c,d的值后
   
4. 譬如a=1,b=2,c=3,d=4
   
5. 
   
6. 定义f=@(x) a.*sin(x)+b.*sqrt(log(x))+c.^2.*d
   
7. quad(f,10,20)
   
8. 即可
   
9. 注：@是匿名函数声明符号，7.0以上版本支持

复制代码

vansen

问题解决了，谢谢各位热心的前辈们！
刚来这里不久，发现这里的学术氛围非常好，给新手们总结的索引贴真的非常好，而且对于每一个问题总有很多热心的人来帮忙解决，很有一种家的温馨的感觉。希望赶快成长，希望能有为大家服务的时候:@)