找回密码
 注册
关于网站域名变更的通知
查看: 466|回复: 4
打印 上一主题 下一主题

ODE45解方程出错

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-8-10 15:56 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
函数文件,
4 s' L# C  [! }% h8 U. Y! G/ Zfunction dy=qipao(t,y,c,Pv,PA,sigma,rho,n,mu,w,P0)6 n5 ~$ S8 l$ z3 d+ `- O7 ?
dy(1,1)=y(2);
0 z1 f7 X, _+ F& M/ B& ~4 }dy(2,1)=((PO+2*sigma/RO-Pv)*((RO/y(1))^(3*n)+Pv+PA*sin(w*t)-P0-2*sigma/y(1))/(rho*y(1))...
; x, t2 ?: C! m  w-4*mu*y(2)/(rho*y(1)^2))...; {1 d4 e" P: X3 g
+((P0+2*sigma/R0-Pv)*(R0^(3*n))*(-3)*n*(y(2))^((-3)*n-1)-w*PA*cos(w*t))/(rho*c)...( f. S9 D' K3 A$ {8 G
-3/2*y(2)^2/y(1);9 ^+ V  l4 v# {/ v
end
- o& H: H7 r6 c+ |$ @这是运行脚本文件,
' r5 \/ M3 R9 C. U% gclear all;, h  T( q% |; s& z
clc;
6 H: n0 k6 A& _% of=0.025*10^6;
, a& _/ m7 j1 F8 O) b2 {mu=1.31*10^(-3);6 E3 C3 K$ v. C8 T0 Y7 R
R0=10;$ N" L/ K# L9 E' A8 d7 ^8 J
tspan=[0,6];
& O& O( A0 p. h% [x0=[10;0];
# B* e% R. e0 g3 Z. \; P& yPA=2.41*10^5;  d0 O" ]% b  N$ N( t
Pv=1.23*10^3;( K4 x$ \! b# a6 X- u
n=1.33;7 N1 s; d6 y" n( t/ ~( [
sigma=7.42*10^(-2);
3 }- y- }7 f& j- P6 a4 nrho=1000;* P( f9 s  R+ Y5 f4 J0 w) p7 u
P0=1.013*10^5;: }( j3 i0 u5 ?( n/ c& g, s+ T- O& N
w=2*pi*f;( B0 i4 v, L; j
[t,y]=ode45(@qipao,tspan,x0);5 L% P( z0 h& `% F% h# V2 J
plot(t,y(:,1),'r-');
# n0 `( V3 S- Z- A% oxlabel('t=Time/T');ylabel('Solution R(t)/RO');legend('R=R(t)/RO')/ X3 ]& {/ k7 u: K
& Q1 m- w3 Y, K: ]% W
运行的时候出现错误:* j" \9 T% d5 x% J7 H
未定义函数或变量 'PO'。$ R1 H: p7 Y. a- m8 _

- H9 w( C" h6 z# J  E出错 qipao (line 4)0 A  b1 ^4 D6 M: }$ m! N8 h
dy(2,1)=((PO+2*sigma/RO-Pv)*((RO/y(1))^(3*n)+Pv+PA*sin(w*t)-P0-2*sigma/y(1))/(rho*y(1))...
# }2 j+ g: q  {& P. g. |  X
  _. j: D  A: \出错 odearguments (line 90)9 O1 U& W* L$ o# Q
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
7 Q" S1 e/ _% c5 n. \3 C4 J3 i
出错 ode45 (line 115): H9 t9 W7 V+ W; U! |* y
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
: G/ s+ V) e2 T+ y( f+ ?2 L
: {8 r$ u, m3 S; H, g出错 qipao1 (line 15)' ^4 r- |$ ]' I
[t,y]=ode45(@qipao,tspan,x0);- t! b. k8 j& H* D

) l& }  [( b0 z% U' v$ a; L: g( h2 G2 R7 [
请问前辈们怎么解决,我查阅了实例和书本,初值都赋值了,ODE45调用也是参考书本来的,一模一样,为什么还会出错????1 D# L' R* o  o5 ~$ @

该用户从未签到

2#
发表于 2020-8-10 17:47 | 只看该作者
帮你顶一下

该用户从未签到

3#
发表于 2020-8-11 13:31 | 只看该作者
不够仔细啊。
: t6 B8 ^* B& L% {0 C  G自定义函数里面的输入变量是P0,结果你下面就打字把数字零打成了字母欧( d. u5 _% h% C& Z
dy(2,1)=((PO...

点评

错 qipao (line 4) dy(2,1)=((P0+2*sigma/RO-Pv)*((RO/y(1))^(3*n)+Pv+PA*sin(w*t)-P0-2*sigma/y(1))/(rho*y(1))... 出错 odearguments (line 90) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to  详情 回复 发表于 2020-8-12 10:19

该用户从未签到

4#
发表于 2020-8-12 09:51 | 只看该作者
来学习一下

该用户从未签到

5#
 楼主| 发表于 2020-8-12 10:19 | 只看该作者
nkkopd 发表于 2020-8-11 13:31% x; J9 C8 H$ ]8 r: g  f9 _
不够仔细啊。7 j% Y; T. ~" U2 b' X1 a9 T4 D) a
自定义函数里面的输入变量是P0,结果你下面就打字把数字零打成了字母欧
& y' L& f4 v9 L# r  X' k6 y9 _, K, Mdy(2,1)=((PO...
/ a! w+ J+ x4 r! N0 c
错 qipao (line 4). O- G4 }. o- i. t
dy(2,1)=((P0+2*sigma/RO-Pv)*((RO/y(1))^(3*n)+Pv+PA*sin(w*t)-P0-2*sigma/y(1))/(rho*y(1))...
7 T$ L; }5 f2 O, a8 `, q: q7 `- W6 W. V) n4 R
出错 odearguments (line 90)  t, Y5 w( D) M7 q6 k
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.1 J: [5 W- l0 n) r/ `1 ?9 B

" i' Q& n0 F3 \出错 ode45 (line 115)
* u* {2 t% p4 O    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
$ D& |8 O/ p( Q, X. z- s! d1 {" J- l; q( F0 ~) Z; f/ y
出错 qipao1 (line 17)0 j, r+ m. n* Z
[t,y]=ode45(@qipao,tspan,x0);
' O$ |$ |3 B8 j$ ~1 x) t
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 18:30 , Processed in 0.171875 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表