作业帮 > 综合 > 作业

matlab解微分方程组

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/16 19:07:52
matlab解微分方程组
求微分方程组的数值解并作图.
要求是用ode45的方法来解
在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function dydx=zhidao_ws19221223(x,fg)
f=fg(1);
g=fg(2);
dydx=[3*f+4*g;-4*f+3*g];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
在Matlab下面输入:
x_start=0;
x_end=2;
fg0=[0;1]; %初值
[x,fg]=ode45('zhidao_ws19221223',[x_start,x_end],fg0);
plot(x,fg);
legend('f','g');
xlabel('x');
figure
plot(fg(:,1),fg(:,2))
xlabel('f')
ylabel('g')