|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
最近看关于求解隐式微分方程的求解的内容,遇到了这样的问题,对于求解问题/ U. L& n: ? z, j4 o
' B' I' u3 n: h* ]7 b nx0= [1;0;0;1]
2 b& G: ?* p. ^% u/ n t. {1 B2 Z使用ode15i进行求解,具体的命令如下:) r3 W3 U, O6 O7 s0 b1 [
6 g+ Z9 Z) I/ S W. L0 i
+ O; C0 c/ `& Q* L `
odefun=@(t,x,dx)[dx(1)-x(2)
1 Z( ]' ^5 Y7 ydx(2)*sin(x(4))+dx(4)^2+2*x(1)*x(3)-x(1)*dx(2)*x(4)
: w8 c( o2 z3 ddx(3)-x(4), c1 n6 a d8 d/ z1 ^& @
x(1)*dx(2)*dx(4)+cos(dx(4))-3*x(3)*x(2)];9 l2 S/ p" ]" _
+ E1 Y$ Q4 G: H* j- O4 z
$ [% I M2 B) I: u7 E
t0=0* i* |2 v7 c# b- Y
x0=[1 0 0 1]';%本题初值 x0 的都给出了0 ?' K D! I, J: Q
fix_x0=ones(4,1);%本题中 x0 都给出了" a% H. d* T! s( r6 e+ u' l
* U9 e) q+ @! l f1 k! c4 s( X6 \. Q" |dx0=[0 0 1 1]';%本题中初值 dx0 一个都没有给出,那么全部任意写
0 ]5 R* D G4 Nfix_dx0=zeros(4,1);%本题中 dx0 一个没有给出,故全部为 0[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);/ B: m A) Y Q$ x" J
+ j, V# U6 p! Q5 Y
2 k L* |6 f- n' Q9 d8 F: z @
但是执行[x02,dx02]=decic(odefun,t0,x0,fix_x0,dx0,fix_dx0);后总是提示* a8 ~2 b8 O$ U+ P( }5 ~; ?' j! ?
Error using decic (line 109)
7 a0 {2 r( W% _' ^( T- h; z- gConvergence failure in DECIC.
# W( t3 G* J* ?/ j9 A$ v( Y* ~请问出错是为什么呢? P5 L, G. w9 k9 D+ f1 ]
|
|