| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1027 人关注过本帖
标题:[求助]用ode45解微分方程的程序,错哪儿了????
只看楼主 加入收藏
jianminliu
Rank: 1
等 级:新手上路
帖 子:5
专家分:0
注 册:2007-7-20
收藏
 问题点数:0 回复次数:1 
[求助]用ode45解微分方程的程序,错哪儿了????
微分方程为:在附件中。
我编的程序如下:
x(1)=x;x(2)=y;z=[x(1),x(2)]';
function xdot=shier(t,z)
r=1;d=0.5;a=0.1;b=0.02;
xdot=[r-a*x(2),0;0,-d+b*x(1)]*z;
ts=0:0.1:15;
x0=[25,2]';
[t,z]=ode45(@shier,ts,x0);
[t,z]
plot(t,z),grid,
gtext('\fontsize{12}x(t)'),gtext('\fontsize{12}y(t)'),
pause,plot(z(:,1),z(:,2)),grid,
xlabel('x'),ylabel('y')

运行后出现:
??? Error: File: shier.m Line: 2 Column: 1
Function definitions are not permitted at the prompt or in scripts.
这是为什么?
raP7jezv.rar (5 KB) [求助]用ode45解微分方程的程序,错哪儿了????


[此贴子已经被作者于2007-7-20 23:46:41编辑过]

搜索更多相关主题的帖子: 方程 解微分 
2007-07-20 23:44
abingchem
Rank: 6Rank: 6
等 级:贵宾
威 望:24
帖 子:716
专家分:0
注 册:2004-12-30
收藏
得分:0 
ode求解器所求解的微分方程有个固定的格式,这个格式不能改的,你把r等的值直接使用吧

此人已消失
2007-07-21 20:37
快速回复:[求助]用ode45解微分方程的程序,错哪儿了????
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.042182 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved