song113
2007-12-04, 23:09
各位MATLAB高手们，
你们好。小弟正在用MATLAB 7.4 计算一个数值积分。试用了3个函数, trapz, quad, quadl. Trapz的结果不太理想，因为不够精确。 quad, quadl用了之后发现一个很奇怪的问题。下面是我的函数
%-------------------%
function y = myfun(x,b)
y = ((sin(2 .* x)).^2 .* (exp((-x .^2)/(2 .* b .^2)))) ./ (b .* sqrt(2 .*pi));
%-------------------%
得到的结果
clear all
b = 1/125;
>> quadl (@(x)myfun(x,b), 0, pi/4)
ans =
2.8540e-019
clear all
>> b = 1/125;
>> quadl (@(x)myfun(x,b), 0, 0.08)
ans =
1.2797e-004
因为这个函数是一直都大于0的。我从0积分到pi/4，理论上应该比0积分到0.08要大很多。但是结果确实相反的。小弟思考了很久都不知道怎么回事？恳请各位赐教。多谢！
或者能介绍一个更好的积分方法适合我的函数的。函数就是上面给出的那个，积分区间是0到pi/4。先谢了:)
你们好。小弟正在用MATLAB 7.4 计算一个数值积分。试用了3个函数, trapz, quad, quadl. Trapz的结果不太理想，因为不够精确。 quad, quadl用了之后发现一个很奇怪的问题。下面是我的函数
%-------------------%
function y = myfun(x,b)
y = ((sin(2 .* x)).^2 .* (exp((-x .^2)/(2 .* b .^2)))) ./ (b .* sqrt(2 .*pi));
%-------------------%
得到的结果
clear all
b = 1/125;
>> quadl (@(x)myfun(x,b), 0, pi/4)
ans =
2.8540e-019
clear all
>> b = 1/125;
>> quadl (@(x)myfun(x,b), 0, 0.08)
ans =
1.2797e-004
因为这个函数是一直都大于0的。我从0积分到pi/4，理论上应该比0积分到0.08要大很多。但是结果确实相反的。小弟思考了很久都不知道怎么回事？恳请各位赐教。多谢！
或者能介绍一个更好的积分方法适合我的函数的。函数就是上面给出的那个，积分区间是0到pi/4。先谢了:)