搜狗做题网matlab 二阶微分方程组解
来源：学生作业帮 编辑：搜狗做题网作业帮 分类：数学作业 时间：2024/06/16 15:20:46
matlab 二阶微分方程组解
function Dy=dfun(t,y)
D2y(1,1)=(0.0176*y(2)-960-23.85/(0.0088-0.0072*y(1)).^1.4-0.1*Dy(1,1))/96;
Dy(2,1)=(-0.0176*Dy(1,1)-0.1*y(2)+2)/(2.4e7*y(1));
t=linspace(0,0.15,100);
ode113('dfun',t,[0,16e6]);
这有什么问题
function Dy=dfun(t,y)
D2y(1,1)=(0.0176*y(2)-960-23.85/(0.0088-0.0072*y(1)).^1.4-0.1*Dy(1,1))/96;
Dy(2,1)=(-0.0176*Dy(1,1)-0.1*y(2)+2)/(2.4e7*y(1));
t=linspace(0,0.15,100);
ode113('dfun',t,[0,16e6]);
这有什么问题
把原方程写出来看看,你这种写法比较非主流
再问：
t=0时 x=0 p1=16e6 x的一阶导数也等于0 t取0到0.15 求p1与t的关系x与t的关系
再答： 你的初值有问题，再检查检查clcf=@(t,x)([x(2);-.0027/96*x(3)+10-9.5e-6/96/(.0088-.0072*x(1))^1.4-.1/96*x(2);(2-.0176*x(2)-.1*x(3))/x(1)/2.4e7]);[t,x]=ode45(f,[0 .15],[1 1 16e6]);subplot(1,2,1),plot(t,x(:,3))
xlabel('t'),ylabel('p1')
subplot(1,2,2),plot(t,x(:,1))
xlabel('t'),ylabel('x')
再问：
t=0时 x=0 p1=16e6 x的一阶导数也等于0 t取0到0.15 求p1与t的关系x与t的关系再答： 你的初值有问题，再检查检查clcf=@(t,x)([x(2);-.0027/96*x(3)+10-9.5e-6/96/(.0088-.0072*x(1))^1.4-.1/96*x(2);(2-.0176*x(2)-.1*x(3))/x(1)/2.4e7]);[t,x]=ode45(f,[0 .15],[1 1 16e6]);subplot(1,2,1),plot(t,x(:,3))
xlabel('t'),ylabel('p1')
subplot(1,2,2),plot(t,x(:,1))
xlabel('t'),ylabel('x')