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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
我想利用ode45函数求解一个方程组,方程组中一个方程的系数的参数来自函数外部,不管怎样都无法调用,请问应该如何解决?function dy=fun(t,y,a)
  K, J- y3 J$ }% qdy=zeros(2,1);' z/ M! Z6 [, }+ v& [
dy(1)=sin(y(2));. e( n  F$ A' j4 W7 R3 P$ Z
dy(2)=a(4)+a(3)*t+a(2)*t.^2+a(1)*t.^3;
; {% o5 k0 {& O. _8 i
5 I8 ~3 m1 g; ?2 ~/ ?' _. l; V! O% W6 j: N
[t,f]=ode45('fun',[t0 tf],[y0 yf],a);
: w- L) O7 z, X. e. U运行一直出错0 v9 Y  |% C, _: Z+ C; Y% D
错误提示:
3 g4 ]' M0 ^) O! N! _; E( U/ oAttempted to access a(4); index out of bounds because numel(a)=0.- i4 o" w# `- _# s
  • TA的每日心情
    开心
    2022-1-29 15:05
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2022-10-18 13:13 | 只看该作者
    试试5 l0 e. g2 E" W! p9 f
    [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: O. ^+ E  w: j; ]; g2 _! g
    试试4 j- |# W$ c. C1 i- S5 E3 j$ D) i
    [t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);
    , j% H  q5 t  Z$ o! ^, n
    我现在想用ode45求解一个  微分方程  
    1 M9 P+ \- R/ ep=1;
      }. h3 J$ L# o2 o; f% }a=0.826;: ~+ q9 ~1 L& b5 x5 x5 w& y
    s=a*cos(p*t);8 @) Y- y3 k. t9 ]4 F
    s1=0.001*cos(1*t);
      ~  w& E+ C* r3 [5 \9 T  U, ^* l. Px1=s1+sqrt(0.01)*randn(length(t));
    6 e& M) c1 o1 h9 ]6 t* Uyprime=[p*y(2);p*(-0.5*y(2)+y(1)^1-y(1)^3+s+x1)];
    6 T  s! A6 \* l& K这里面如果将s1改为一串序列(从外面调用一串数据)  而不是用公式表示  还能不能求解呢?要怎么改呢8 l" K2 O0 L' E$ g& A
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-7-10 10:23 , Processed in 0.125000 second(s), 28 queries , Gzip On.

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

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

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