| 编程中国 | 业界新闻 | 技术文章 | 视频教程 | 下载频道 | 程序源码 | 个人空间 | 编程论坛
全能ASP/PHP/ASP.NET主机,支持月付专业 MSSQL 数据库空间,支持月付专业 MySQL 数据库空间,支持月付学习型 ASP/PHP/ASP.NET 主机 30元/年
高端软件开发 = 年薪十万不是梦赛孚耐:软件保护加密专家身份认证令牌USB KEY 
共有 991 人关注过本帖
标题:matlab关于ode45解二阶微分方程的困惑
收藏  订阅  推荐  打印 
watermoon
Rank: 1
等级:新手上路
帖子:1
积分:110
注册:2007-8-5
matlab关于ode45解二阶微分方程的困惑

matlab关于ode45解二阶微分方程的困惑

一个二阶微分方程:
y''+y'+y=sin(t)

初始条件为y(0)=5,y'(0)=6。

过程:
先降阶为一阶微分方程组
y'=z
z'=-z-y+sin(t)

编制如下函数m文件
function dy=weifen(t,x)
dy=zeros(2,1);
%y=x(1)
%z=x(2)
dy(1)=x(2);
dy(2)=sin(t)-x(2)-x(1);

*******************
然后用ode45解方程
[t,y]=ode45(@weifen,[0 20],[5 6])
plot(t,y)就可以画出y和y'的图像

但是我用同样的思路,是不过在函数文件的编制时候将函数的代表意义换一下,结果出来的图形就完全不一样了?

即我的m文件重新编制为
function dy=weifen(t,x)
dy=zeros(2,1);
%y=x(2)
%z=x(1)
dy(1)=x(1);
dy(2)=sin(t)-x(1)-x(2);

*******************
然后同样用ode45解方程,根据函数文件只需将初始条件换一下位置
[t,y]=ode45(@weifen,[0 20],[6 5])
plot(t,y)就可以画出y和y'的图像,但是结果完全不一样了,怎么回事?
难道必须有一定的顺序么?

搜索更多相关主题的帖子: 微分方程  matlab  二阶  困惑  
2007-8-5 21:26
abingchem
Rank: 12Rank: 12Rank: 12
等级:贵宾
威望:24
帖子:716
积分:7466
注册:2004-12-30

这个顺序是必须的,所谓置换降阶是采用了升序的顺序

此人已消失
2007-8-6 20:06
关于我们 | 广告合作 | 编程中国 | 清除Cookies | Archiver | WAP | TOP

编程中国 版权所有,并保留所有权利。鲁ICP备08000592号
Powered by Discuz, Processed in 0.058171 second(s), 9 queries.
Copyright©2004-2008, BCCN.NET, All Rights Reserved