ODE45解方程出错

llbnmo

函数文件，
function dy=qipao(t,y,c,Pv,PA,sigma,rho,n,mu,w,P0)
dy(1,1)=y(2);
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))...
-4*mu*y(2)/(rho*y(1)^2))...
, o) t* ]  V7 H, c+((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);
& H' a+ g3 W; @& |1 @/ Iend
这是运行脚本文件，
" ~! j5 c* l# k, gclear all;
. S( B5 ~  s3 I: fclc;
4 o9 [/ }, O# Z- Gf=0.025*10^6;
mu=1.31*10^(-3);
R0=10;
tspan=[0,6];
* f- u; P. V7 K% A/ F. h7 lx0=[10;0];
( K+ H" r( W- K* e/ E2 ?PA=2.41*10^5;
! X1 I0 ?4 u& b0 l4 g4 \Pv=1.23*10^3;
n=1.33;
sigma=7.42*10^(-2);
" }* X$ ~* n' g# wrho=1000;
9 }2 G. Q6 D) M" S" F* d& TP0=1.013*10^5;
  {/ w, y+ z9 j' n, Hw=2*pi*f;
[t,y]=ode45(@qipao,tspan,x0);
) T- ^* \( x& v3 lplot(t,y(:,1),'r-');
" [" i; ]6 u4 o# g7 b- |) wxlabel('t=Time/T');ylabel('Solution R(t)/RO');legend('R=R(t)/RO')
4 w6 T- _. z: Z* ~: t, A7 O
! O" |8 ]6 u7 Q' d  K% ~) ]. ?/ Z运行的时候出现错误：
未定义函数或变量 'PO'。

出错 qipao (line 4)
3 |& r  d& m: zdy(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))...

/ d- O" v7 m3 Z7 e. C- A出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

出错 ode45 (line 115)
+ C6 j, C0 s) q: e    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

出错 qipao1 (line 15)
% I( n4 @4 D" j( l- i[t,y]=ode45(@qipao,tspan,x0);


8 y- B( L) _1 z; T4 W( `请问前辈们怎么解决，我查阅了实例和书本，初值都赋值了，ODE45调用也是参考书本来的，一模一样，为什么还会出错？？？？

nkkopd

帮你顶一下

nkkopd

不够仔细啊。
0 c4 J+ ]0 T2 {: _3 t& {2 i  ~& }1 ?自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
dy(2,1)=((PO...

Uifhjvv

来学习一下

llbnmo

nkkopd 发表于 2020-8-11 13:31
不够仔细啊。
自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
1 [% m0 x8 Y; p, K3 H9 A3 qdy(2,1)=((PO...

8 H' w5 C% v# ?* N错 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)
# I8 O! Z! |$ j& Vf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

' R7 B. H, t( K1 _! n出错 ode45 (line 115)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

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