研究生论文 一类变上限积分函数图像处理！

bobbin2002

:loveliness: 欲做下积分式子函数图像进行研究：
syms x n t;
t=0:0.05:1;
a1=1/gamma(0.2);
f1=2.^(-0.5*n);
f2=x.^(-0.8);
f3=sin((2.^n).*(t-x));
f4=int(f2.*f3,x,0,t);
f5=symsum(f1.*f4,n,1,30);
y=a1.*f5;
plot(t,subs(y));
函数可以运行，但在运算f4的时候结果出现虚数，以致予函数图像和想象中的有一定差距，应该不对；不知道我这个积分表达式写的正确不，f4=int(f2.*f3,x,0,t),其中的X为 f2.*f3中的未知数，其取值从0到t；
小弟在此恭候各位大虾，如能解开此谜团，小弟不甚感激！！

maigicku

这个偶用quadl试了试，出现了搜索步长为0的情况。。期待高人来解答。。

[ 本帖最后由 maigicku 于 2009-9-21 17:29 编辑 ]

bobbin2002

高人高人你在何方 不知道传说中的高手是什么模样 万分期待您的到来 我等。。。。。。。在此先感谢楼上朋友的：）

maigicku

感觉LZ的被积函数比较差，总是出现“Warning: Maximum function count exceeded; singularity likely.”这样的警告。。

bobbin2002

出问题应该的地方我想应该就是这个积分函数，但是又不知道错在哪儿?还有我运算的时候并没有出现楼上的情况，只是说“警告，出现了虚数部分被忽略的情况”；起初我认为程序没有问题，然后把函数中部分数值改动了，结果发现out of memory的情况，应该是因本程序积分部分出现问题导致运算产生了大量的中间结果，致使图像偏差了。不知道有哪位仁兄能帮帮小弟把程序合理化一下 继续等候：）

[ 本帖最后由 bobbin2002 于 2009-9-22 18:13 编辑 ]

bobbin2002

程序已解决 多谢各位参与  现在看来问题不是出在积分表达式上，表达式没有问题，问题在于运算的方式，由于太多中间结果需要保留和处理，导致函数精度不够而与理论图像偏差太大，所以改成数值计算后，精度得到提升，图形也逼近理论了，所以 多谢大家了：）:lol