ODE45解方程出错

llbnmo

函数文件，
" d- d( J" b: k$ i; pfunction dy=qipao(t,y,c,Pv,PA,sigma,rho,n,mu,w,P0)
dy(1,1)=y(2);
5 s' T  A8 h5 G6 y) u/ `$ edy(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))...
" Q# u2 }- ]" S# k-4*mu*y(2)/(rho*y(1)^2))...
, T0 \0 c0 \7 p- q+((P0+2*sigma/R0-Pv)*(R0^(3*n))*(-3)*n*(y(2))^((-3)*n-1)-w*PA*cos(w*t))/(rho*c)...
-3/2*y(2)^2/y(1);
1 `! ~! R2 i" {end
这是运行脚本文件，
& e0 O: R1 h+ q, d* A+ dclear all;
$ I( C1 i: G6 s! M- p, Iclc;
f=0.025*10^6;
. I+ [8 O( P  s3 i# Z, X  mmu=1.31*10^(-3);
; N, D, R- f) N, W" JR0=10;
6 z! p7 [( U; v$ Qtspan=[0,6];
x0=[10;0];
PA=2.41*10^5;
% F2 r  N7 c* ~' ~& [* A0 d" T3 FPv=1.23*10^3;
n=1.33;
8 g5 c! p0 u# g( @& D* Esigma=7.42*10^(-2);
rho=1000;
# C# C% J; O$ j- P9 _1 [( a/ M& `5 sP0=1.013*10^5;
% O6 z* g; `2 w" l* ^) Vw=2*pi*f;
* B  S1 A2 `. P7 V# _: j$ g[t,y]=ode45(@qipao,tspan,x0);
plot(t,y(:,1),'r-');
5 B( H0 N* r" i3 q: {5 h' txlabel('t=Time/T');ylabel('Solution R(t)/RO');legend('R=R(t)/RO')
5 z' N, ^. @1 S! [  ]
4 ^2 k8 A3 U$ m3 J' q  i0 }( k& `运行的时候出现错误：
' g& J6 Q& T9 C$ o* C. g未定义函数或变量 'PO'。
8 P3 ^+ `: x. I1 r5 v3 V
0 z  o( Q0 E) @  ^2 ~- b出错 qipao (line 4)
& R0 b; g( b& C# w* _0 |' Pdy(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))...
7 B- p  ?6 C. g; o! t2 V4 v. v
出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

出错 ode45 (line 115)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
/ j9 s0 l8 j& L
. P+ u: X! i9 e$ \出错 qipao1 (line 15)
! @& l3 M$ _( o2 o5 F# |4 _[t,y]=ode45(@qipao,tspan,x0);


请问前辈们怎么解决，我查阅了实例和书本，初值都赋值了，ODE45调用也是参考书本来的，一模一样，为什么还会出错？？？？

nkkopd

帮你顶一下

nkkopd

不够仔细啊。
. l( O" D3 _) |. s2 F自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
dy(2,1)=((PO...

Uifhjvv

来学习一下

llbnmo

nkkopd 发表于 2020-8-11 13:31
不够仔细啊。
8 j2 |4 A8 O* M3 u自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
dy(2,1)=((PO...

4 g6 s% O8 J( x* C错 qipao (line 4)
+ t7 V% l% ^8 w- d5 p2 l  `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))...
( I4 c/ l( A, d  [* `& ?& g! [
出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

0 V2 V* _/ t& \, N% \' R. k( y# y5 ]出错 ode45 (line 115)
" M( n+ i! K& I, i5 _5 o' O3 X    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

出错 qipao1 (line 17)
[t,y]=ode45(@qipao,tspan,x0);