求解二元微分方程组
题目
求解二元微分方程组
C1(P0-P1)dt=VdP1+C2(P1-P2)dt (1)
C2(P1-P2)dt=VdP2+C3P2dt (2)
求P1和P2关于t的函数,其余均为常数.
第一个式子V为V1
第一个式子V为V2
请最终分别给出P1和P2的函数式。
答案
移项一下,将得到:
dP1/dt=[c1*(p0-P1)-c2*(P1-P2)]/V1;
dP2/dt=[c2*(p1-P2)-c3*P2]/V1;
这个是常微分方程组,可以用龙格库塔法求解.
以下过程是用MATLAB实现的
新建一个m文件:equ.m
function dx = equ(t,x)
%下面是设置参数,可以根据你的实际情况自行设置
c1=1;
c2=2;
c3=3;
P0=0.5;
V1=1;
V2=2;
dx=[(c1*(P0-x(1))-c2*(x(1)-x(2)))/V1;(c2*(x(1)-x(2))-c3*x(2))/V2];
然后返回MATLAB的命令窗口,输入
t_end=10;
x0=[0;1]; %系统的初值
[t,P]=ode45('equ',[0,t_end],x0);
plot(t,P) %作图两个变量的图
figure;
plot(P(:,1),P(:,2)); %相轨迹
举一反三
已知函数f(x)=x,g(x)=alnx,a∈R.若曲线y=f(x)与曲线y=g(x)相交,且在交点处有相同的切线,求a的值和该切线方程.
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
最新试题
- Lucy,can you take my book ___the office.
- 1.Would you please ()the window?It’s too cold outside.A open B opening C not open D not to open
- 地球无私地提供矿产资源给人类.改成"把"字句是什么?
- 将与门,或门,非门,与非门,或非门的逻辑综述成一句话
- 找出课文中描写小孩子过春节的部分读一读,再联系实际,说说你是怎样过春节的
- 二元一次方程,求解答,在线等
- 我们不如搞一次户外烧烤翻译成英语
- 醇酸酯化时,是醇去OH酸去H,还是是醇去H酸去OH
- 翻译Rock is the music of teenage rebellion
- 一个简单的英语句子分析
热门考点