请问如何实现使用matlab对下面的一阶微分方程组进行拉普拉斯变换?diff(P1)=-1i*kai*P+omega/sqrt(2)*P0diff(P0)=-1i*sqrt(2)*omega*P1得到P0(s)和P1(s)的表达式后如何在进行反变换?最终得到P1(t)和P0(t
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/18 12:43:17
请问如何实现使用matlab对下面的一阶微分方程组进行拉普拉斯变换?diff(P1)=-1i*kai*P+omega/sqrt(2)*P0diff(P0)=-1i*sqrt(2)*omega*P1得到P0(s)和P1(s)的表达式后如何在进行反变换?最终得到P1(t)和P0(t
请问如何实现使用matlab对下面的一阶微分方程组进行拉普拉斯变换?
diff(P1)=-1i*kai*P+omega/sqrt(2)*P0
diff(P0)=-1i*sqrt(2)*omega*P1
得到P0(s)和P1(s)的表达式后如何在进行反变换?最终得到P1(t)和P0(t)
请问如何实现使用matlab对下面的一阶微分方程组进行拉普拉斯变换?diff(P1)=-1i*kai*P+omega/sqrt(2)*P0diff(P0)=-1i*sqrt(2)*omega*P1得到P0(s)和P1(s)的表达式后如何在进行反变换?最终得到P1(t)和P0(t
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dxdt=zhidao_rk4_5(t,x)
global a b c d r
x1=x(1);
x2=x(2);
x3=x(3);
x4=x(4);
dxdt=[
a*(x2-x1)+x4;
d*x1-x1*x3+c*x2;
x1*x2-b*x3;
x2*x3+r*x4;
];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
global a b c d r
a=1;
b=2;
c=3;
d=4;
r=5;
t_end=10; %时间计算的终点,即算到这个时间为止
x0=[1;1;1;1]; %x的初始值
[t,x]=ode45('zhidao_rk4_5',[0,t_end],x0);