ODE45解方程出错

llbnmo

函数文件，
function dy=qipao(t,y,c,Pv,PA,sigma,rho,n,mu,w,P0)
! R9 n3 I+ M) J: \# }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))...
+((P0+2*sigma/R0-Pv)*(R0^(3*n))*(-3)*n*(y(2))^((-3)*n-1)-w*PA*cos(w*t))/(rho*c)...
; r$ {! d, W$ M% K' H3 O-3/2*y(2)^2/y(1);
8 K  Y/ E5 Z3 r2 w3 X7 O, J) {5 Xend
: j- l/ H; `5 A- q% \这是运行脚本文件，
clear all;
8 N6 R/ e" |1 Zclc;
  x2 N. z  C3 Q0 y( w, vf=0.025*10^6;
- E3 B: S  e0 T6 Q/ wmu=1.31*10^(-3);
4 L1 j* t/ v* q* H  c2 A' jR0=10;
4 x+ P& X# p* q  v6 X$ j) R$ @tspan=[0,6];
x0=[10;0];
! s5 n2 r0 W0 Y# ]( M3 j- [PA=2.41*10^5;
, E9 u2 ~' d5 `. f, [2 jPv=1.23*10^3;
' S6 ?3 T$ @) ?7 T; y" @n=1.33;
& x3 k) d, O* w3 I# V) Lsigma=7.42*10^(-2);
rho=1000;
% Q9 i7 u% Q# PP0=1.013*10^5;
w=2*pi*f;
& \6 v* g5 X: ]& ]; M[t,y]=ode45(@qipao,tspan,x0);
" |" z8 A- [6 ~& ~$ r: f- J, c1 `plot(t,y(:,1),'r-');
xlabel('t=Time/T');ylabel('Solution R(t)/RO');legend('R=R(t)/RO')

运行的时候出现错误：
0 m( U2 Y% @& U3 p$ d  i/ B5 z未定义函数或变量 'PO'。
, Y+ P/ R1 i! u8 X( r& D8 _( w
3 u: G/ X2 ?0 U出错 qipao (line 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))...
, E& n, B- g; H) V
出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

出错 ode45 (line 115)
9 R. h  a; U" p% L4 s, _) Q    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
  n6 N8 e1 E+ w( r  v' m
出错 qipao1 (line 15)
) r) F( _% a  U6 b$ U  c[t,y]=ode45(@qipao,tspan,x0);


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

nkkopd

帮你顶一下

nkkopd

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

Uifhjvv

来学习一下

llbnmo

nkkopd 发表于 2020-8-11 13:31
不够仔细啊。
自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
( V7 X/ `+ n. s. edy(2,1)=((PO...

错 qipao (line 4)
9 i1 m; W. T& U7 K0 [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))...
/ w( r, R5 T8 ?. n" N0 c
出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

出错 ode45 (line 115)
7 ]! B$ ~8 w! C1 }: _7 y! d$ K    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
6 t& W. M- X9 T# @) I( e3 ^
0 G0 z4 Z' J: Z/ O3 g( \0 }2 {出错 qipao1 (line 17)
& l: f8 Y! {9 N9 I2 A" P8 k[t,y]=ode45(@qipao,tspan,x0);