[x,y]=dsolve('Dx=x-y-x*(x^2+y^2)','Dy=x+y-y*(x^2+y^2)','x(0)=2','y(0)=1')
得到的结果是解析解没有找到.
用数值解.
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=zhidao_rk4_5(t,x)
%x,y变量分别用x(1),x(2)表示
y=[x(1)-x(2)-x(1)*(x(1)^2+x(2)^2);x(1)+x(2)-x(2)*(x(1)^2+x(2)^2)];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
t_end=10;
x0=[2;1];
[t,x]=ode45('zhidao_rk4_5',[0,t_end],x0);
plot(t,x);
legend('x','y');
xlabel('t');
figure;
plot(x(:,1),x(:,2));
xlabel('x');
ylabel('y');
用matlab解微分方程组
用matlab解微分方程组
dx/dt=x-y-x(x^2+y^2)
dy/dt=x+y-y(x^2+y^2)
x(0)=2
y(0)1
dx/dt=x-y-x(x^2+y^2)
dy/dt=x+y-y(x^2+y^2)
x(0)=2
y(0)1
其他人气:408 ℃时间:2019-12-13 18:17:13
优质解答
我来回答
类似推荐
猜你喜欢
- 1我想问一个问,一支合唱队里男生比女生多36,男生的人数是女生的4倍,合唱队里女生有多少人.
- 2take your advice.
- 3一个长方体,长4m.宽3m.它的占地面积最大是多少,最小是多少?
- 4The poor farmer lived (happily) in his old hut.(括号里的划线提问)
- 5she wants to have a ruler 怎样改为否定句
- 6NH4+与OH-生成NH3↑?还是氨水?
- 7甲乙两地相距530千米,客货车相向而行,客车先行了80千米,货车才出发,货车出发6小时两车相遇,已知客车
- 8I need time to understand how to love you~! 我说:Afraid even if it passas ,not leave~!的中文
- 9感恩老师的作文,书信格式,400字,快啊,回答的好给分
- 10将下列化学反应文字表达式用化学方程式写出来并配平.写的好再加100分