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

RWS_123

我想利用ode45函数求解一个方程组，方程组中一个方程的系数的参数来自函数外部，不管怎样都无法调用，请问应该如何解决？function dy=fun(t,y,a)
7 V7 Z6 w1 a: D! s% {% }dy=zeros(2,1);
/ {- c3 H; b% e6 ^dy(1)=sin(y(2));
- m/ f; Z6 {4 S% z+ E% e' `& Kdy(2)=a(4)+a(3)*t+a(2)*t.^2+a(1)*t.^3;
. J1 D  m% Y  K- L5 G$ G0 }

% H; H' u9 m  ^% f$ m[t,f]=ode45('fun',[t0 tf],[y0 yf],a);
! S  ~' [2 S$ L3 t2 Q8 b运行一直出错
错误提示：
9 C2 x- @# z3 pAttempted to access a(4); index out of bounds because numel(a)=0.

someone

试试
8 K* o4 e+ V6 G; l9 U[t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);

RWS_123

someone 发表于 2022-10-18 13:13
. c/ c3 A. k) B& C5 D试试
- [7 F$ \7 i- g" k: k0 l3 F, n% l[t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);

7 O% i4 o' ?3 _7 {2 z 我现在想用ode45求解一个  微分方程  
p=1;
a=0.826;
, C( s1 Y. o5 xs=a*cos(p*t);
s1=0.001*cos(1*t);
x1=s1+sqrt(0.01)*randn(length(t));
$ c! n4 q1 X/ {0 b1 m. Uyprime=[p*y(2);p*(-0.5*y(2)+y(1)^1-y(1)^3+s+x1)];
1 X0 p3 p: u  h这里面如果将s1改为一串序列（从外面调用一串数据）  而不是用公式表示  还能不能求解呢？要怎么改呢
/ N' O; Y$ e$ L6 J