TA的每日心情 | 开心
2022-1-29 15:04 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我的M文件：
9 x1 x! Y/ G# c- zy0est = [0.3133283972712190; 0.001104257865152230; 0.55; 0; 0; 0; 0; 0; 0; 0;];
/ [$ E2 i( B1 z/ u( O1 Syp0est = zeros(10,1);
6 Z4 C4 e. `# I4 k; f5 K! x+ Zopt = odeset('RelTol', 10.0^(-7), 'AbsTol' , 10.0^(-7));2 S# C4 V* D2 G# O2 x9 u8 ?9 s- U
[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)! Y+ c3 z2 w) h/ ^
5 P- n; W! V1 {$ n. D# C
5 s5 W3 [4 W. e! ^! ~( {: u运行后报错：
9 u; z, g6 z9 m. X* F错误使用 decic>sls (line 170)
" @9 @% X. h; n* b) d索引可能大于 1。7 M( p; L/ E7 R/ i! N% Z
H# i# p8 D; f! a, E" D% J- v' Q" w
出错 decic (line 77)" {2 E; s: b7 k) {
[dy,dyp] = sls(res,dfdy,dfdyp,neq,free_y,free_yp);9 }( x0 n, Q7 j3 {6 M3 e0 m
# \$ F5 W( l1 ~2 t出错 DAEsL (line 39)
. U# b" S+ f$ U" V# A) I[y0, yp0] = decic(F, 2.533037911852984, y0est, [], yp0est, [], opt)' D8 C8 }$ P; D: o( o. r- L
3 X1 W, ^! R/ i' L6 I6 c大家帮分析一下是什么问题？6 {$ a0 u* g: F
|
|
/1 
发表于 2022-9-15 13:15
收藏
淘帖
支持！
反对！
楼主