ode45('fun',[],[]) 中的函数如何调用外部参数

RWS_123

我想利用ode45函数求解一个方程组，方程组中一个方程的系数的参数来自函数外部，不管怎样都无法调用，请问应该如何解决？function dy=fun(t,y,a)
dy=zeros(2,1);
$ e, G. h& O6 Z4 N& R$ udy(1)=sin(y(2));
dy(2)=a(4)+a(3)*t+a(2)*t.^2+a(1)*t.^3;
/ p8 W. ?( g( h. e( o

9 f; D8 v" D3 o6 E[t,f]=ode45('fun',[t0 tf],[y0 yf],a);
运行一直出错
- h0 G) q- @/ ?+ d: O0 Q& G错误提示：
; d. }; M# o% y; R. WAttempted to access a(4); index out of bounds because numel(a)=0.
- M8 t$ R; w4 F* U2 k7 z( m

someone

试试
: R* C" n* `* I! \/ p[t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);

RWS_123

someone 发表于 2022-10-18 13:13
6 r' ~& o! I0 x4 Z6 n; W试试
4 C, {( x& u5 ?1 ~& I! b$ [[t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);

  t- R# a5 J, K 我现在想用ode45求解一个  微分方程  
2 ~! i# I( }; B0 s4 X; Wp=1;
8 E3 L' d: {" Q4 a/ e  B4 Ka=0.826;
; c( l9 [9 t# `& w8 h6 ns=a*cos(p*t);
% I6 e7 X8 q, Bs1=0.001*cos(1*t);
x1=s1+sqrt(0.01)*randn(length(t));
yprime=[p*y(2);p*(-0.5*y(2)+y(1)^1-y(1)^3+s+x1)];
这里面如果将s1改为一串序列（从外面调用一串数据）  而不是用公式表示  还能不能求解呢？要怎么改呢