zw123 发表于 2006-3-24 19:29

[求助]积分方程求解

我还有个问题请教<br>y2=(a1*(x2-x)*a+a2*(x-x1)*b)/(a1*(x2-x1)+a2*(x-x1));<br>y=sqrt(1+(diff(y2))^2);<br>现在需要对x求积分,积分限是x1~x2,<br>我使用j=int(y,x1,x2)<br>为什么不行,<br>应该怎么处理?<br>谢谢各位!<br>
[此贴子已经被aspen于2006-3-25 19:11:56编辑过]

yqchenlolo 发表于 2006-3-24 21:31

试了一下可以的<BR>你具体的函数是怎么定义的？a1,a2,a,b都是常数么？<BR>检查一下是不是哪里不小心搞错了<BR>或者把你之前各个变量的具体定义贴出来看看

happy 发表于 2006-3-24 21:38

回复：（zw123）[求助]积分方程

系数太多了，是解出来，不过估计有几千项

freelive 发表于 2006-3-25 10:32

对于^1/2的积分计算机很慢,而且产生虚数<BR>解决方法对 y 在(x1+x2)/2处8阶用taylor展开,在用int <BR>注意taylor只能对符号展开,要先将y定义为符号<BR>sym x;<BR>taylor( y(x),(x1+x2)/2,8);
页: [1]
查看完整版本: [求助]积分方程求解