ODE45解方程出错

llbnmo

函数文件，
# o4 P5 K* Q. \) V  O4 v( {function dy=qipao(t,y,c,Pv,PA,sigma,rho,n,mu,w,P0)
dy(1,1)=y(2);
2 B( l1 S6 ]5 m. r8 A% V3 mdy(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))...
& j- g0 Y: |2 _7 J1 f0 d$ g  C( }7 y-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+ G/ F& k. h+ P# R, X- ~8 l-3/2*y(2)^2/y(1);
end
% `  ~* u3 p3 i' p+ a" W这是运行脚本文件，
! A  \/ e7 u, G) `; k) J  \clear all;
; f6 ~4 j/ u5 H  m1 k# N7 `* y0 |clc;
9 C4 w) `. C/ Z& V  @- {f=0.025*10^6;
+ h( A3 v$ ?" d7 Q6 K  Amu=1.31*10^(-3);
R0=10;
; s: m. b# P9 R$ h5 G7 t) p# ?tspan=[0,6];
x0=[10;0];
8 Y) y9 n& i$ R2 `* zPA=2.41*10^5;
Pv=1.23*10^3;
5 ^" L$ m" C& xn=1.33;
sigma=7.42*10^(-2);
) C- w; Y8 N7 t/ Nrho=1000;
: P# o4 r& C/ fP0=1.013*10^5;
# T( \2 e; _' Q1 ?) Ew=2*pi*f;
" k) ?- [2 N% i[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')
/ O0 d) t  c" Y9 l
运行的时候出现错误：
; Y! v) Q: y# V/ k未定义函数或变量 'PO'。
; \1 s8 i0 y; V! g( y
出错 qipao (line 4)
8 o* ?( T& u" [5 n9 k. {  Jdy(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))...
: ~' M2 f! S, v
出错 odearguments (line 90)
: f- o3 p6 t4 Z% z2 J( W* v: Rf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
/ I4 ~( x, b' a* d) v5 g
出错 ode45 (line 115)
; ^0 M0 x1 s5 M3 K8 u8 `    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

; F9 I5 T7 |; I出错 qipao1 (line 15)
7 `8 L+ H% n' V9 p# w1 D[t,y]=ode45(@qipao,tspan,x0);
- ^9 W& C! s! G3 V

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

nkkopd

帮你顶一下

nkkopd

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

Uifhjvv

来学习一下

llbnmo

nkkopd 发表于 2020-8-11 13:31
, ]. r6 ]) I+ t4 A: q0 U0 _2 v不够仔细啊。
  s1 S8 Q% C# x; `# }) \& i* i自定义函数里面的输入变量是P0，结果你下面就打字把数字零打成了字母欧
) ~5 P/ d& ]5 D; O" v4 mdy(2,1)=((PO...

( B" j7 H% K, R) ~# R) _错 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))...
& h3 A3 n. c$ A1 Y
出错 odearguments (line 90)
& W9 F  J7 t- ~. [3 _2 N1 }5 `f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
6 O9 p' e9 Q( ~
出错 ode45 (line 115)
    odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);

) Z/ q% N% k8 Q4 |) M出错 qipao1 (line 17)
% [3 |& j, G9 a' e8 S/ y1 y[t,y]=ode45(@qipao,tspan,x0);
7 Q* W* G/ w4 `, m, a