
1个回答
展开全部
哥们,这没有解析解。我算了一下,但是结果是x与y值保持0.5不变。你瞅瞅咋回事,我也很想知道这是为什么
追问
大哥,这个是这篇论文的结论,我不知道怎么设计程序语言,能不能麻烦把你的程序发给我,让我试试先,谢谢了
追答
function dy=myfun(t,y)
dy=zeros(4,1);
dy(3)=y(1);%y(3)表示的是x,y(1)表示的是x的导数
dy(4)=y(2);%y(4)表示的是y,y(2)表示的是y的导数
dy(1)=0.75*y(3)^3+0.75*y(3)*y(4)^2-6*y(3);
dy(2)=0.75*y(4)^3+0.75*y(4)*y(3)^2-6*y(4);
保存
然后在窗口调用[t y]=ode45(@myfun,[0 12],[0 0 0.5 0.5])%t是得到的时间序列,y是x、y及其导数的值
然后plot(T,y(:,3))%x的图像
plot(T,y(:,3))%y的图像
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询