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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 有限差分法 第二类边界条件
51/1返回列表
查看: 5215  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-05-09 10:48:18
已阅   同方向广播   申请计算强帖   回复此楼   编辑   查看我的主页

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
3楼2014-05-09 21:24:53
已阅   申请计算强帖   回复此楼   编辑   查看我的主页
查看全部 13 个回答

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
fegg7502: 金币+3, 应助指数+1, 专家考核, 3ks 2014-05-13 09:16:37
这个可以考虑用bvp解算器做
BVP解算器
solinit = bvpinit(x, yinit, params)
sol = bvpsolver(odefun,bcfun,solinit,options)
由于边值问题可能有多解，为了便于我们确定那个解是我们需要的，所以必须使用bvpinit函数对初值进行估计
解算器(bvpsolver)：Matlab中提供了bvp4c和bvp5c，后者误差控制更好些
输入参数：
x：需要计算的网格点，相当于ode**的tspan
yinit：猜测的值，可以是具体值，也可以是函数，类似与 ode**的 x0
params：其它未知参数，也是一个猜测值
odefun：描述边值问题微分方程的函数句柄

bcfun：边值函数，一般是双边值(x 的上下限即认为两个边界)，但也支持多边值(具体看帮助)
bc(y(a),y(b))=0
举个例子，y(0)=0,y(4)=-2 ---bcfun可写成 res =@(ya,yb) [ ya(1) ; yb(1) + 2];
y'(0)=y'(pi)=0,y(0)=1-----------bcfun可写成res = @(ya,yb) [  ya(2); yb(2);ya(1)-1 ];
solinit：由bvpinit生成的初始化网格
solinit = bvpinit(linspace(pa,pb,10),[1 0]);
bvpinit的第一个参数表示求解域及中心节点，第二个参数表示猜测初始值，基本都写成[1 0]，它的变化一般对解无影响
options：BVP解算器优化参数，可以通过bvpset设置，具体参数查看帮助

降阶法表达二阶的常微分方程，ya(1)表示函数在求解域[a,b]的左边界a上的值，ya(2)表示函数的一阶导数在a上的值，“ya(1)-5”实际上是“ya(1)-5=0”的简写，表示y(x=a)=5，即函数在左边界上的值为5。yb是右边界，其他与ya相同。
用bvp4c求得的解sol是一个结构体，其中的子项y是一个2行多列的矩阵，第一行表示函数值，第二行表示函数的导数。另一个子项yp也是一个2行多列的矩阵，第一行表示函数一阶导数值，第二行表示函数二阶导数值。用deval求出的y与sol.y相同
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼2014-05-09 19:47:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
478107371: 金币+60, ★★★★★最佳答案, 大神，请接受我的膝盖，分全是你一人儿的 2014-05-10 00:11:25
fegg7502: 金币+2, 应助指数+1, 专家考核, 3ks 2014-05-13 09:16:54
举个例子吧
y''-3y'+2y=x， y‘(0)=1,y’(1)=2
CODE:
%解析解
yy=dsolve('D2y-3*Dy+2*y=x','Dy(0)=1','Dy(1)=2','x');
x=0:0.01:1;hold on;plot(x,eval(yy),'bo')
%bvp 数值解
f1=@(x,y)[y(2);x-2*y(1)+3*y(2)];
solinit = bvpinit(linspace(0,1,20),[-1 1]);
f2=@(ya,yb)[ya(2)-1;yb(2)-2];
sol = bvp4c(f1,f2,solinit);
xint = linspace(0,1,25);
yint = deval(sol,xint);
hold on; plot(xint,yint(1,:),'rp-');
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
4楼2014-05-09 21:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)
本帖仅楼主可见
赞一下
5楼2014-05-10 01:25:56
已阅   申请计算强帖   回复此楼   编辑   查看我的主页
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 6/300 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +8 a089 2026-08-07 8/400 2026-08-14 18:46 by beefly
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭