dsolve问题解决的系统的常微分方程

2视图(30天)

显示旧的评论

詹尼斯Cakstins 2017年5月11日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/339880-dsolve-issue-with-solving-nasty-system-of-odes

评论道: 詹尼斯Cakstins2017年5月16日

答:接受 Torsten

亲爱的所有,

我一直有这个问题一段时间,我不能找到一个解决方案系统的常微分方程如下:

                         %刚度矩阵的输入
                        
                         X1 = -138.7840 e - 024;
                        
                         X2 = 123.0769 e + 003;
                        
                         X3 = -61.5385 e + 003;
                        
                         %的刚度矩阵:
                        
                         (X1 = 0 0 0;
                        
                         0 X2 X3 X3;
                        
                         0 X3 X2 X3;
                        
                         0 X3 X3 X2);
                        
                         %定义常微分方程
                        
                         信谊x (t) y (t) z (t) w (t)
                        
                         Y = [x;y;z;w];
                        
                         Q = [x;x * y;x * z;x * w];%乘以刚度矩阵
                        
                         我= Y (0) = = (16.0746;12.9888;5;0);%初始条件
                        
                         常微分方程= diff (Y) = = * Q;
                        
                         [xSol (t) ySol (t) zSol (t) wSol (t)) = dsolve(常微分方程,我);
                        
                         在10 ^ %显示解决方案
                        
                         显示(双(xSol (10 ^)));
                        
                         显示(双(ySol (10 ^)));
                        
                         显示(双(zSol (10 ^)));
                        
                         显示(双(wSol (10 ^)));

这是我得到的输出消息:

                         警告:明确的解决方案不能被发现。
                        
                         在dsolve >(第201行)
                        
                         在question_for_forum(16行)
                        
                         错误使用符号/ subsindex(第732行)
                        
                         无效的索引或函数定义。当定义一个函数,确保这一
                        
                         参数是象征性的变量和函数的主体是一个符号表达式。
                        
                         当索引,输入必须为数字、逻辑或“:”。
                        
                         错误在question_for_forum(16行)
                        
                         [xSol (t) ySol (t) zSol (t) wSol (t)) = dsolve(常微分方程,我);