版块导航: 正在加载中...

客户端APP下载

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

chouwa

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 36.5
散金: 4
红花: 2
帖子: 275
在线: 33.7小时
虫号: 72525
注册: 2005-06-04
性别: GG
专业: 系统可靠性与管理

[交流] 【求助】Volterra方程的求解

本人最近遇到一个问题与第二类Volterra方程有关，经过几次仿真总是不理想，怀疑是Vollterra解存在问题。请各位高人指导一下。
方程的形式为
$f(t)=\int_{a}^{t}K(t,s)f(s)ds+g(t)$
其中已知条件f（0）＝0；当t＝s时，K(t,s)＝0。
由于对精度要求不是很高，所以我采用梯形法进行求解。
例如求f（b）
取步长
$h=\frac{b-a}{N}$
有
$t_{i}=a+i\times h$
依据
$\int_{a}^{t_{i}}K(t_{i},s)f(s)ds=h(\frac{1}{2}K_{i,0}f_{0}+\sum_{j=1}^{i-1}k_{i,j}f_{j}+\frac{1}{2}K_{i,i}f_{i})$
则有
$f_{i}=h(\frac{1}{2}K_{i,0}f_{0}+\sum_{j=1}^{i-1}k_{i,j}f_{j}+\frac{1}{2}K_{i,i}f_{i})+g_{i}$
利用条件简化得
$f_{i}=h(\frac{1}{2}k_{i,0}f_{0}+\sum_{j=1}^{i-1}k_{i,j}f_{j})+g_{i}$
这里
$h=\frac{y}{N}$
我通过Mathematic用
[img] http://latex.codecogs.com/gif.la ... e;h]*f[j*h])][/img]
求f（y)
不知道是不是求解错了，请各位高人指点

[ Last edited by chouwa on 2010-10-16 at 13:17 ]

» 猜你喜欢

面上项目申报已经有3人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有5人回复
酰胺脱乙酰基已经有9人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有5人回复
有70后还继续奋斗在职场上的吗？已经有6人回复
博士延得我，科研能力直往上蹿已经有7人回复
退学或坚持读已经有27人回复
面上基金申报没有其他的参与者成吗已经有5人回复
遇见不省心的家人很难过已经有22人回复

» 本主题相关价值贴推荐，对您同样有帮助:

奋斗今天，希冀明天

1楼 2010-10-16 00:48:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chouwa

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 36.5
散金: 4
红花: 2
帖子: 275
在线: 33.7小时
虫号: 72525
注册: 2005-06-04
性别: GG
专业: 系统可靠性与管理

抱歉

Mathematic的程序表示不出来，不知道为何

[ Last edited by chouwa on 2010-10-16 at 01:19 ]

奋斗今天，希冀明天

2楼2010-10-16 01:05:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

数学EPI: 17
应助: 256 (大学生)
贵宾: 1.36
金币: 2336.2
散金: 9224
红花: 92
帖子: 4583
在线: 1303.8小时
虫号: 394701
注册: 2007-06-07
专业: 流体力学

chouwa(金币+1):谢谢 2010-10-16 13:06:36

用迭代方法自己求解就是了。

长期招收博士生，参见http://fsl-unsw.com

3楼2010-10-16 11:17:18

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 chouwa 的主题更新