matlab二阶微分方程组解functionDy=dfun(t,y)D2y(1,1)=(0.2023*y(2)-960-23.85/(0.2023-0.2023*y(1)).^1.4-0.1*Dy(1,1))/96;Dy(2,1)=(-0.2023*Dy(1,1)-0.1*y(2)+2)/(2.4e7*y(1));t=linspace(0,0.15,100);ode113(#39;dfun#39;
<p>问题：matlab二阶微分方程组解functionDy=dfun(t,y)D2y(1,1)=(0.2023*y(2)-960-23.85/(0.2023-0.2023*y(1)).^1.4-0.1*Dy(1,1))/96;Dy(2,1)=(-0.2023*Dy(1,1)-0.1*y(2)+2)/(2.4e7*y(1));t=linspace(0,0.15,100);ode113(#39;dfun#39;,t,);这有什么问题<p>答案：↓↓↓<p class="nav-title mt10" style="border-top:1px solid #ccc;padding-top: 10px;">孙作玉的回答：<div class="content-b">网友采纳 把原方程写出来看看,你这种写法比较非主流<p class="nav-title mt10" style="border-top:1px solid #ccc;padding-top: 10px;">孙作玉的回答：<div class="content-b">网友采纳 你的初值有问题，再检查检查 clc f=@(t,x)(); =ode45(f,,); subplot(1,2,1),plot(t,x(:,3)) xlabel('t'),ylabel('p1') subplot(1,2,2),plot(t,x(:,1)) xlabel('t'),ylabel('x')
页:
[1]