matlab 解微分方程组 (高分求解)奇异性如何解决dx1/dt=(sin(x3)+cos(x3))/sin(x2)dx2/dt=cos(x3)-sin(x3)dx3/dt=1-(sin(x3)+cos(x3))*cot(x2)我想把x1,x2,x3的数值解求出来.作图.并要知道最后一个时刻的x1,x2,x3的数值.
来源:学生学帮网 编辑:学帮网 时间:2024/07/06 16:18:36
matlab 解微分方程组 (高分求解)奇异性如何解决
dx1/dt=(sin(x3)+cos(x3))/sin(x2)
dx2/dt=cos(x3)-sin(x3)
dx3/dt=1-(sin(x3)+cos(x3))*cot(x2)
我想把x1,x2,x3的数值解求出来.
作图.并要知道最后一个时刻的x1,x2,x3的数值.
初始值为(0,0,0) ,因此第一行的最后一项将除于一个零,导致无法算出来.若用eps代替,好像也会有很大的误差,这样如何解决.
全答出来的给100分
回答一楼的问题:首先感谢你的关注。
请问我要显示出最后一个时刻的值,是什么代码?
方法如下:先建立一个微分方程组的子函数:
function f=funx1_3(t,x)
f=zeros(3,1);
f(1)=sin(x(3))+cos(x(3))./sin(x(2));
x(2)=cos(x(3))-sin(x(3));
x(3)=1-(sin(x(3))+cos(x(3))).*cot(x(2));
return
再在主界面里调用:
[T,F] = ode45(@funx1_3,[0 12],[0 0 0]);
画图也很容易:
plot(T,F)
但是初值全是0的话解出来的是不是有问题就难说,给个非零数会不会好些,我提供了代码,你可以自己多尝试几组数据.
可不可以先人工算出来
吾是这样考虑的:
syms x1 x2 x3
f=Dx1-(sin(x3)+cos(x3))/sin(x2);
g=Dx2-cos(x3)+sin(x3);
h=Dx3+(sin(x3)+cos(x3))*cos(x2)-1;
[x1,x2,x3]=solve(f,g,h];
不过,Dx1好像不太对,各位高手参考。
matlab 微分方程组求解
matlab求解以下微分方程组!
求解微分方程组的解析解.matlab代码,谢谢
这个微分方程组杂用matlab编程求解
matlab怎么求解微分代数方程组
matlab矩阵求求解微分方程组{用matlab编程}数值求解微分方程组求程序
如何运用matlab求解非线性偏微分方程组
matlab大神,求解微分方程组代码求代码
用MATLAB求解符号微分方程组(求大虾帮忙)dx=ydy=x+y不太会用MATLAB
用matlab解微分方程组并作图
matlab微分方程组其中
用Matlab能求解二元二阶微分方程组,求Matlab程序!方程组如图.
求解常微分方程组
matlab 解微分方程组 (高分求解)奇异性如何解决dx1/dt=(sin(x3)+cos(x3))/sin(x2)dx2/dt=cos(x3)-sin(x3)dx3/dt=1-(sin(x3)+cos(x3))*cot(x2)我想把x1,x2,x3的数值解求出来.作图.并要知道最后一个时刻的x1,x2,x3的数值.
用matlab求解微分方程组打字不好打 我直接上图吧
Matlab用四阶龙格库塔法求解不可微分方程组的初值y1'=120-2*y1+2*y2 (0
用matlab编程解微分方程组编程解下边方程
用matlab解微分方程组,怎么写程序呢?