常微分方程组边值问题解法.docVIP

下载本文档

14
0
约8.46千字
约 21页
2019-10-07 发布于天津
举报
版权申诉

常微分方程组边值问题解法.doc

1、本文档内容版权归属内容提供方，所产生的收益全部归内容提供方所有。如果您对本文有版权争议，可选择认领，认领后既往收益都归您。。
2、本文档由用户上传，本站不保证质量和数量令人满意，可能有诸多瑕疵，付费之前，请仔细先通过免费阅读内容等途径辨别内容交易风险。如存在严重挂羊头卖狗肉之情形，可联系本站下载客服投诉处理。
3、文档侵权举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

常微分方程组边值问题解法打靶法打靶法求常微分方程的边值问题应用打靶法求解下列边值问题解将其转化为常微分方程组的初值问题命令真实解函数边值常微分方程组函数命令计算结果有限差分法同上例若划分为个区间则函数有限差分法求常微分方程的边值问题命令计算结果真实解正交配置法构造正交矩阵函数正交配置矩阵均用矩阵法求对称性与非对称性正交配置矩阵为形状因子为零点数为对称的权函数为权函数非为权函数为零点数为非对称的权函数为权函数非为权函数正交多项式求根适用于对称问题为形状因子为配置点数为权函数权函数为权函数为为多项式

PAGE 常微分方程组边值问题解法打靶法Shooting Method（shooting.m） %打靶法求常微分方程的边值问题 function [x,a,b,n]=shooting(fun,x0,xn,eps) if nargin3 eps=1e-3; end x1=x0+rand; [a,b]=ode45(fun,[0,10],[0,x0]); c0=b(length(b),1); [a,b]=ode45(fun,[0,10],[0,x1]); c1=b(length(b),1); x2=x1-(c1-xn)*(x1-x0)/(c1-c0); n=1; while (norm