| 查看: 1241 | 回复: 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,: )) |
回复此楼» 猜你喜欢
让我中一个面上吧！
已经有10人回复
-
科研孤儿太难了
已经有14人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
-
今天基金会出结果吗？20260819
已经有15人回复
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
40