用matlab求解微分方程!dy/dx=((1-0.786)^1.33)*((3*p1*x^2+2*p2*x+p3)/(p1*x^3+p2*x^2+p3*x+p4))*y其中 x=2008时,y=20.83*10^6 p1 = -57.924 p2 = 3.4558e+005 p3 = -6.8719e+008 p4 = 4.5547e+011求出y与x的方程,和画图 x的范围1990-2025
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/23 16:51:36
![用matlab求解微分方程!dy/dx=((1-0.786)^1.33)*((3*p1*x^2+2*p2*x+p3)/(p1*x^3+p2*x^2+p3*x+p4))*y其中 x=2008时,y=20.83*10^6 p1 = -57.924 p2 = 3.4558e+005 p3 = -6.8719e+008 p4 = 4.5547e+011求出y与x的方程,和画图 x的范围1990-2025](/uploads/image/z/8813371-67-1.jpg?t=%E7%94%A8matlab%E6%B1%82%E8%A7%A3%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%21dy%2Fdx%3D%28%281-0.786%29%5E1.33%29%2A%28%283%2Ap1%2Ax%5E2%2B2%2Ap2%2Ax%2Bp3%29%2F%28p1%2Ax%5E3%2Bp2%2Ax%5E2%2Bp3%2Ax%2Bp4%29%29%2Ay%E5%85%B6%E4%B8%AD+x%3D2008%E6%97%B6%2Cy%3D20.83%2A10%5E6++p1+%3D+-57.924++p2+%3D+3.4558e%2B005++p3+%3D+-6.8719e%2B008++p4+%3D+4.5547e%2B011%E6%B1%82%E5%87%BAy%E4%B8%8Ex%E7%9A%84%E6%96%B9%E7%A8%8B%2C%E5%92%8C%E7%94%BB%E5%9B%BE+x%E7%9A%84%E8%8C%83%E5%9B%B41990-2025)
用matlab求解微分方程!dy/dx=((1-0.786)^1.33)*((3*p1*x^2+2*p2*x+p3)/(p1*x^3+p2*x^2+p3*x+p4))*y其中 x=2008时,y=20.83*10^6 p1 = -57.924 p2 = 3.4558e+005 p3 = -6.8719e+008 p4 = 4.5547e+011求出y与x的方程,和画图 x的范围1990-2025
用matlab求解微分方程!dy/dx=((1-0.786)^1.33)*((3*p1*x^2+2*p2*x+p3)/(p1*x^3+p2*x^2+p3*x+p4))*y
其中 x=2008时,y=20.83*10^6
p1 = -57.924
p2 = 3.4558e+005
p3 = -6.8719e+008
p4 = 4.5547e+011
求出y与x的方程,和画图 x的范围1990-2025
用matlab求解微分方程!dy/dx=((1-0.786)^1.33)*((3*p1*x^2+2*p2*x+p3)/(p1*x^3+p2*x^2+p3*x+p4))*y其中 x=2008时,y=20.83*10^6 p1 = -57.924 p2 = 3.4558e+005 p3 = -6.8719e+008 p4 = 4.5547e+011求出y与x的方程,和画图 x的范围1990-2025
%%%%%%%%%%%%%%%%%%% 求解主体代码 %%%%%%%%%%%%%%%%%
clear;clc;
close all;
options = odeset('RelTol',1e-4,'AbsTol',1e-5);
[X1,Y1] = ode45(@aatest,[2008 2025],20.83E6,options);
[X2,Y2] = ode45(@aatest,[2008 1990],20.83E6,options);
figure;
hold on;
grid on;
plot(X1,Y1/1E6)
plot(X2,Y2/1E6)
%%%%%%%%%%%%%%%%%%% 微分方程对应函数 %%%%%%%%%%%%%%%%%
function dy=aatest(x,y)
p1 = -57.924;
p2 = 3.4558e+005;
p3 = -6.8719e+008;
p4 = 4.5547e+011;
dy=((1-0.786)^1.33)*((3*p1*x^2+2*p2*x+p3)/(p1*x^3+p2*x^2+p3*x+p4))*y;
结果图