ODE45解方程出错

llbnmo

函数文件，
- M# o3 j  o& ~% l! Z: [function dy=qipao(t,y,c,Pv,PA,sigma,rho,n,mu,w,P0)
( T& n' k3 |  I, X( u& b+ s/ idy(1,1)=y(2);
' Q$ H+ D/ O- C' x1 M( Rdy(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))...
1 F) P: P0 i9 p* p-4*mu*y(2)/(rho*y(1)^2))...
6 s" E' I$ W) M" e- d5 z  g+((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);
& i* i% M% _9 L% n& Rend
这是运行脚本文件，
clear all;
clc;
f=0.025*10^6;
: t$ z4 w) Z! ]7 V& v! Hmu=1.31*10^(-3);
9 }1 t0 U( @; ~% v: u9 G& SR0=10;
- _5 _8 B; r" G  I6 G5 c2 itspan=[0,6];
2 {0 s6 S9 V9 ex0=[10;0];
PA=2.41*10^5;
* x8 D2 M( [. G# OPv=1.23*10^3;
  J, r1 E. T6 M& ~: [n=1.33;
sigma=7.42*10^(-2);
) Q7 L* R1 P1 M# X, Qrho=1000;
P0=1.013*10^5;
w=2*pi*f;
[t,y]=ode45(@qipao,tspan,x0);
plot(t,y(:,1),'r-');
3 i6 I+ ^" p: W9 m" {xlabel('t=Time/T');ylabel('Solution R(t)/RO');legend('R=R(t)/RO')
1 O! ~0 S7 @- Z0 Z
! v/ T- E' w# n2 R& r运行的时候出现错误：
未定义函数或变量 'PO'。
1 Q- f; I4 M1 C  e
# R: Q/ Y6 ~9 S出错 qipao (line 4)
$ C- _0 ?5 }( N8 [$ M% Z' }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))...
( H! L: ~3 z2 Z$ L
出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
5 [6 z. {, |! {
; @7 P: x6 s1 T+ B/ E出错 ode45 (line 115)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
0 i5 Y$ X) z! N+ N
' j* j& l$ T+ q  B8 _出错 qipao1 (line 15)
[t,y]=ode45(@qipao,tspan,x0);


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

nkkopd

帮你顶一下

nkkopd

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

Uifhjvv

来学习一下

llbnmo

nkkopd 发表于 2020-8-11 13:31
不够仔细啊。
# Y. X+ K4 U% |! I! ~" @4 ]自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
; [8 O7 D" j$ n& S4 B) U1 ydy(2,1)=((PO...

( Y' S! H3 Y2 Q8 J% Y2 ^7 q错 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 yp0.

出错 ode45 (line 115)
! K% h# p6 z( S    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
4 E# h' x" I& v4 v
出错 qipao1 (line 17)
[t,y]=ode45(@qipao,tspan,x0);