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))...
6 \: L2 }) T& M# _) y+((P0+2*sigma/R0-Pv)*(R0^(3*n))*(-3)*n*(y(2))^((-3)*n-1)-w*PA*cos(w*t))/(rho*c)...
: A2 {+ p. R+ Y6 x0 F& G% s1 S-3/2*y(2)^2/y(1);
end
& z( O$ P4 s& F) b这是运行脚本文件，
clear all;
clc;
' G+ t/ J1 n2 ]& t0 jf=0.025*10^6;
- F: {' R0 |8 }" i9 \/ O; Rmu=1.31*10^(-3);
; R1 D# x% |5 @/ U! j. m$ a5 hR0=10;
" ~+ D4 M/ P7 V3 `* c/ qtspan=[0,6];
x0=[10;0];
8 e1 f4 }8 z  R0 }. U) O- @) B3 BPA=2.41*10^5;
* k# a- n" w, b5 z# O) [- EPv=1.23*10^3;
4 A4 f- S9 b, M) r- An=1.33;
4 G1 \" q3 _$ H4 N, C$ usigma=7.42*10^(-2);
! ^7 n# ]3 C! q) arho=1000;
5 E, x4 X, S, E' z9 b1 eP0=1.013*10^5;
+ n6 f4 e: X5 b* A! ]' Fw=2*pi*f;
0 B2 `& p9 j7 t0 K1 h! [# y4 L3 v3 E[t,y]=ode45(@qipao,tspan,x0);
plot(t,y(:,1),'r-');
xlabel('t=Time/T');ylabel('Solution R(t)/RO');legend('R=R(t)/RO')

8 p2 f5 f/ M  K运行的时候出现错误：
/ w+ ~4 L& {6 c% r6 g" ], H# r未定义函数或变量 'PO'。
( V" @6 e% I+ Y  w  L: C
出错 qipao (line 4)
& p& `; E( G4 Sdy(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))...

出错 odearguments (line 90)
* d: R, Z' q6 V4 Sf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
$ D( T9 v) ?7 _. c2 s& V
出错 ode45 (line 115)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

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


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

nkkopd

帮你顶一下

nkkopd

不够仔细啊。
/ v! v. k* }7 v自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
* K3 }- D% H9 A+ o5 @dy(2,1)=((PO...

Uifhjvv

来学习一下

llbnmo

nkkopd 发表于 2020-8-11 13:31
不够仔细啊。
自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
: A3 q- T% O7 ~% Pdy(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 yp0.
1 M8 \( u. Q4 F1 O9 Z
出错 ode45 (line 115)
+ o5 k: B0 U- H* \6 u    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

出错 qipao1 (line 17)
. A4 U* r7 V! o+ X4 o: [' k( |[t,y]=ode45(@qipao,tspan,x0);
" M" }: X' B+ P$ o5 w+ L