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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
我想利用ode45函数求解一个方程组,方程组中一个方程的系数的参数来自函数外部,不管怎样都无法调用,请问应该如何解决?function dy=fun(t,y,a)0 b# R& B4 K7 A: E7 S" ?; M- @
dy=zeros(2,1);
- h: g- Z, K/ w, vdy(1)=sin(y(2));
/ I  s1 e. s- l* ]& X7 }2 bdy(2)=a(4)+a(3)*t+a(2)*t.^2+a(1)*t.^3;
. d/ g/ J9 R" ~
! p' k5 p; t, W3 `7 G& o; Z
* L& V' L$ z* Y' T) q' @  ~' y[t,f]=ode45('fun',[t0 tf],[y0 yf],a);
! g3 O$ H! n) N. l运行一直出错/ ~0 L: f; j- X* r. w" g% J" [
错误提示:6 }! ~3 K! c8 \2 C( D5 \" g
Attempted to access a(4); index out of bounds because numel(a)=0.
7 A. l2 @2 F5 D3 h# U# T
  • TA的每日心情
    开心
    2022-1-29 15:05
  • 签到天数: 2 天

    [LV.1]初来乍到

    2#
    发表于 2022-10-18 13:13 | 只看该作者
    试试
    9 ?! W1 C/ f2 h7 h2 w- ][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:138 c% l( L0 l& ]
    试试
    4 f) p! P$ |  i8 i# F+ Z[t,f]=ode45(@(t,y)fun(t,y,a),[t0 tf],[y0 yf]);
    ! l+ G2 Z  ~9 n. Z
    我现在想用ode45求解一个  微分方程  
    4 J$ ]/ g; U7 \p=1;  C$ n# r4 \$ ~4 U3 X
    a=0.826;4 ]* g7 E' R6 G, Z' ]: ?' \6 B
    s=a*cos(p*t);# N" e% N: F9 ]/ [# F% m! Q
    s1=0.001*cos(1*t);
    " S; r% _6 j& E  [7 l8 R, gx1=s1+sqrt(0.01)*randn(length(t));% l" D3 _7 R, q+ f% e8 ~: x
    yprime=[p*y(2);p*(-0.5*y(2)+y(1)^1-y(1)^3+s+x1)];
    ; v1 a6 d# x% S6 Z. e9 J: ?# c) o这里面如果将s1改为一串序列(从外面调用一串数据)  而不是用公式表示  还能不能求解呢?要怎么改呢
    # \* \' e- z. c
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

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

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

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

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