求解二元微分方程组C1(P0-P1)dt=VdP1+C2(P1-P2)dt (1)C2(P1-P2)dt=VdP2+C3
来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/10/03 05:24:58
求解二元微分方程组
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的函数式。
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)); %相轨迹
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)); %相轨迹
求解二元微分方程组C1(P0-P1)dt=VdP1+C2(P1-P2)dt (1)C2(P1-P2)dt=VdP2+C3
p2->next=p1;p2=p1;
while(*p1) p1++; while(*p2) { *p1=*p2; p1++; p2++; } *p1=0;
已知曲线c1:y=e*x与c2:y=-1/e*x,若c1,c2分别在p1,p2处得切线是同一条切线,试求出切线方程
k=*p1*(*p2)是什么意思?
p2->next=p1是什么意思
swap(int *p1,int *p2) {int temp; temp=*p1; *p1=*p2; *p2=temp
P1=¬P2 P1= P2 力P1 等于力p2意义和区别
有指针变量p1,p2,p1="abc";p2="abc".strcat(p1,p2)如何理解
p1=p2; \\p1 的指针指向 p2,意味着p2值改变,p1值也改变
p1+p2+p3+p4=1 0
int*p,*p1,*p2,*p=*p1,*p1=*p2,*p2=*p3哪里有错.另外p=p1,p1=p2,p2=p3呢