| 查看: 1114 | 回复: 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,: )) |
回复此楼» 猜你喜欢
材料284求调剂，一志愿郑州大学英一数二
已经有5人回复
-
0856材料专业298分有科研经历 硕士研究生调剂自荐信
已经有3人回复
-
317一志愿华南理工电气工程求调剂
已经有11人回复
-
高分子化学与物理调剂
已经有9人回复
-
寻找调剂
已经有4人回复
-
311求调剂
已经有10人回复
-
材料学调剂
已经有8人回复
-
材料化工调剂
已经有5人回复
-
291分工科求调剂
已经有8人回复
-
0856化工专硕求调剂
已经有5人回复
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 (幼儿园)
- 金币: 2283.6
- 散金: 123
- 红花: 23
- 帖子: 937
- 在线: 213.2小时
- 虫号: 2960148
- 注册: 2014-02-09
- 性别: GG
- 专业: 机械动力学
4楼2018-10-22 11:19:24
10