ODE45解方程出错

llbnmo

函数文件，
function dy=qipao(t,y,c,Pv,PA,sigma,rho,n,mu,w,P0)
" h, y+ ~: e0 b% h8 O% Z/ ndy(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)...
, E, C/ ?8 p' R" K8 y! q+ i& N3 W-3/2*y(2)^2/y(1);
end
! K* Z* q2 y) ]7 g4 q# A& e- c: c0 ?$ L这是运行脚本文件，
- q( p. J9 A# {2 e: T6 \8 \5 y: Dclear all;
7 b$ {; y$ t4 U) ^% _+ M9 n. Hclc;
f=0.025*10^6;
: Z5 \  _; T2 \" H7 g$ q# z. F7 x) k$ ?mu=1.31*10^(-3);
R0=10;
tspan=[0,6];
x0=[10;0];
5 W+ Z: G6 f- z$ hPA=2.41*10^5;
Pv=1.23*10^3;
0 T# @- G" h3 S, ?0 G5 g6 u: [n=1.33;
sigma=7.42*10^(-2);
rho=1000;
2 ]* R1 H5 y9 F) EP0=1.013*10^5;
$ q% U9 _! J; Kw=2*pi*f;
[t,y]=ode45(@qipao,tspan,x0);
plot(t,y(:,1),'r-');
6 v7 |# |* t1 t% }xlabel('t=Time/T');ylabel('Solution R(t)/RO');legend('R=R(t)/RO')

2 }0 c  {  H9 Y: R- z1 e/ K+ A6 Q运行的时候出现错误：
; n7 G% E* Y( R8 ^) P未定义函数或变量 'PO'。
3 E. Y% O! _5 m, ]0 w) H
! N/ j! p4 g0 U; J/ Q" J出错 qipao (line 4)
5 ^: H1 y7 v4 q9 A* k3 k* {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 w; f% F; h; G# t9 [2 w! S
出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

% Y% f% Y) `4 x& k5 h" ?, z出错 ode45 (line 115)
' O+ j- u" t6 q( }    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

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

; ?" `0 b' w  @: E( K' r! o
, U, c( u' F  n3 G9 P5 {' S请问前辈们怎么解决，我查阅了实例和书本，初值都赋值了，ODE45调用也是参考书本来的，一模一样，为什么还会出错？？？？
5 i4 b7 o( l6 s& B: Y- y2 E

nkkopd

帮你顶一下

nkkopd

不够仔细啊。
0 @! ^' h5 W* w. [9 U2 N自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
dy(2,1)=((PO...

Uifhjvv

来学习一下

llbnmo

nkkopd 发表于 2020-8-11 13:31
+ q' W, z( M$ Q( H- [  n  I不够仔细啊。
自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
8 d; j7 d' v8 I. L8 Zdy(2,1)=((PO...

$ r- _7 g1 [) r; Y错 qipao (line 4)
/ \) w9 G& Q9 \7 V3 C# ady(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))...
5 v6 H' c0 V( \( N6 C9 u- p" E
出错 odearguments (line 90)
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
- y; w  ]2 `6 m* t( O/ `" ?
出错 ode45 (line 115)
. i# K4 Q) w( E; \    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

出错 qipao1 (line 17)
[t,y]=ode45(@qipao,tspan,x0);
. T$ [7 f% e" Z0 r' l