[求助][原创]菜鸟请教各位ＭＡＴＬＡＢ高手

yliu593 · 发表于 2005-9-12 05:40:33

<

>菜鸟请教各位ＭＡＴＬＡＢ高手：
<

>求下列高阶微分方程的数值解：
<

>h''(t)*1050=1/2*(0.00025*h^2-0.0107045*h+0. 0124)*h'(t)^2*26.67; h(t)=20000;h'(t)=100;
我是这样解的：
定义了一个函数： 
function dy=ly2(t,y) dy=zeros(2, 1); dy(1)=y(2); dy(2)*１０５０=0.5*0.645*(0.00025*y(1)^2-0.0107045*y(1)+0.0124)*y(2)^2*26.67;
然后运行脚本命令：
[T Y]=ode23t('ly',[0 1000],[20000 100]); plot(T,Y(:,1),'-',T,Y(:,2),'+');
但显示的错误信息是：
??? Error: File: C:\MATLAB6p5\work\ly.m Line: 4 Column: 11 Missing operator, comma, or semicolon.
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\private\odearguments.m On line 76 ==> if (nargin(ode) == 2) 
Error in ==> C:\MATLAB6p5\toolbox\matlab\funfun\ode23t.m On line 202 ==> [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, args, ...
Error in ==> C:\MATLAB6p5\work\LY0.m On line 1 ==> [T Y]=ode23t('ly',[0 1000],[20000 100]); 
这是什么意思呀，如何修改呀？或者怎样才能更好的实现呀！
特别是第一个错误，搞不懂，仔细检查下好象没有少什么东西．

cx19860824 · 发表于 2005-9-12 15:29:08

<

>一个简单又隐秘的错误，你的数字1050输入格式用的中文，改回来！

xhjj8410 · 发表于 2005-9-12 18:14:24

<

>不是呀把1050除到右边去！

		自动登录	找回密码
密码			注-册-帐-号

[求助][原创]菜鸟请教各位ＭＡＴＬＡＢ高手

dd