| 查看: 1289 | 回复: 3 | ||
[求助]
Matlab求解二阶常微分方程边值问题
|
|
我用bvp4c求解了一个简单的二阶常微分方程边值问题 原方程 边界条件 我想得到不同的k时的结果。需要从别的结果中提取k的值，数量较大，一个一个修改不太现实。 据我尝试，描述微分方程的function的参数只能是(x,y)，不能再添加一个k。 怎样能整理成一个function，把k也当作参数呢？ 附：k=1的时候，MATLAB计算过程 先转换成一阶微分方程组 求解过程： 1.magFun.m 描述微分方程 function dy = magFun(x,y) k = 1; dy = zeros(2,1); dy(1) = y(2); dy(2) = -k*cos(y(1)); end 2.magBC.m 描述边界条件 function res = magBC( ya,yb ) res = [ ya(1); yb(2) ]; end 3.magSol.m 求解 function magSol solinit = bvpinit(linspace(0,1,5),[1 0]); sol = bvp4c(@magFun,@magBC,solinit); x = linspace(0,1); y = deval(sol,x); plot(x,y(1,: )) |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有7人回复
-
27博士申请（市政/环境）
已经有3人回复
2楼2016-10-21 15:29:04
shan_yin
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1833.7
- 散金: 170
- 红花: 5
- 帖子: 590
- 在线: 169.7小时
- 虫号: 4856693
- 注册: 2016-07-20
- 性别: GG
- 专业: 动力学与控制
3楼2016-12-27 02:58:17
hgq2016
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 2284.4
- 散金: 123
- 红花: 23
- 帖子: 937
- 在线: 213.2小时
- 虫号: 2960148
- 注册: 2014-02-09
- 性别: GG
- 专业: 机械动力学
4楼2018-10-22 11:19:24
