TA的每日心情 | 开心 2022-1-29 15:04 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
我的M文件:
- Q, i- H1 f3 x) }8 v; by0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
6 ^! d; o) v0 ^/ [4 x5 Eyp0est = zeros(10,1);& V) P* t9 f1 w$ n' ~
opt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7)); ^) @9 x2 L; b2 @+ j
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)- z! J; f3 J8 `- r3 Q- A1 L E, P( m
5 b; s# s' e' K u; J
& @3 j" ~! q" G# U# o4 [
运行后报错:! f! O+ h, c( m3 u* _* z' ~
错误使用 decic>sls (line 170)
p5 b$ N' m( t! f2 ]/ S- H索引可能大于 1。
0 ~$ m" r$ q' ~ O% t2 r
# [; g' d6 Z& o! T. S出错 decic (line 77)6 g H$ z9 `4 X3 @9 s3 S
[dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);" O7 ]5 O5 c [7 q& L
5 l& Z; U& P6 { T出错 DAEsL (line 39)
, K$ b0 C* d3 t" m# h[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)
4 Q7 f; a9 | z! A) z: {, A1 P0 j4 [! \0 M5 B7 `5 e
大家帮分析一下是什么问题?- ^, s8 w0 r7 A- k7 H7 K+ `6 F
|
|