找回密码
 注册
关于网站域名变更的通知
查看: 383|回复: 2
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2022-10-18 11:06 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
我想利用ode45函数求解一个方程组,方程组中一个方程的系数的参数来自函数外部,不管怎样都无法调用,请问应该如何解决?function dy=fun(t,y,a)
+ [$ K& l9 s2 ~; s" U& Z; ody=zeros(2,1);4 N) o! L7 A2 ^8 b' R6 _
dy(1)=sin(y(2));1 z$ a4 x7 p1 v8 q
dy(2)=a(4)+a(3)*t+a(2)*t.^2+a(1)*t.^3;1 W3 z4 M: D  x9 [: W; Y
1 K/ ^, A2 W' ?& n. S1 ^5 i# u

7 g  T# ?# K* r( b& r! _# V[t,f]=ode45('fun',[t0 tf],[y0 yf],a);
5 c) t+ J* J' z) ^& S运行一直出错; S2 X) U5 u3 I
错误提示:
& V3 |" F( Z8 b+ {* w4 }" X8 ~Attempted to access a(4); index out of bounds because numel(a)=0.
5 }6 y# S& X# ^
  • TA的每日心情
    开心
    2022-1-29 15:05
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2022-10-18 13:13 | 只看该作者
    试试
    ' q6 y/ A# U' S& \[t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);

    点评

    我现在想用ode45求解一个 微分方程 p=1; a=0.826; s=a*cos(p*t); s1=0.001*cos(1*t); x1=s1+sqrt(0.01)*randn(length(t)); yprime=; 这里面如果将s1改为一串序列(从外面调用一串数据) 而不是用公式表  详情 回复 发表于 2022-10-18 13:17

    该用户从未签到

    3#
     楼主| 发表于 2022-10-18 13:17 | 只看该作者
    someone 发表于 2022-10-18 13:13
    4 [( p2 M7 O* g# N8 a试试; L  v  ]- u4 C) f
    [t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);

    - M" S* N: D7 v) @  } 我现在想用ode45求解一个  微分方程  " T! J2 E* m! K: |" P
    p=1;' ~$ ~2 x! e. `( L* k  x1 k7 s
    a=0.826;
    1 u' H% i) D0 G* O8 M6 bs=a*cos(p*t);# Q0 S) {2 P% D9 p. t
    s1=0.001*cos(1*t);
    5 M1 n; p) U0 S2 \* G* b3 P7 J, gx1=s1+sqrt(0.01)*randn(length(t));1 ]& T3 J* S9 l! u
    yprime=[p*y(2);p*(-0.5*y(2)+y(1)^1-y(1)^3+s+x1)];0 L" x* g4 ~+ C( `. ?; C  [
    这里面如果将s1改为一串序列(从外面调用一串数据)  而不是用公式表示  还能不能求解呢?要怎么改呢+ }& w/ |: G- Y' o/ L  c1 w
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-10-4 11:21 , Processed in 0.125000 second(s), 27 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表