作者justlikelbj (00100020)
看板MATLAB
标题[问题]微分方程
时间Tue Sep 8 12:55:24 2015
我要解一个system of ODE
dydt(1)=-0.3*y(1)+0.5*y(4);
dydt(2)=-0.3*y(2)-0.5*y(3);
dydt(3)=0.5*y(2)-0.25*y(3)+2.5*y(6);
dydt(4)=-0.5*y(1)-0.25*y(4)-2.5*y(5);
dydt(5)=2.5*y(4)-y(5);
dydt(6)=-2.5*y(3)-y(6);
y(1)=a(t) y(2)=b(t)...y(6)=f(t)
initial condition a(0)=b(0)=c(0)=d(0)=f(0)=0,e(0)=1
现在我想多加一个条件 a(t)^2+b(t)^2+...+f(t)^2=1
也就是说
a(0)^2+b(0)^2+...+f(0)^2=1,
a(1)^2+b(1)^2+...+f(1)^2=1,
a(2)^2+b(2)^2+...+f(2)^2=1,
...
a(n)^2+b(2)^2+...+f(n)^2=1,
不知道要怎麽做才好？
google了很久 还是没有答案
感谢大大相助
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.114.127.67
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1441688126.A.198.html
1F：推 sunev: 这是线性微方，initial condition给定就决定了 09/08 16:03
2F：→ sunev: 没有空间加其它条件 09/08 16:03
3F：→ justlikelbj: 那如果去掉initial condition 只要我後来那个条件 要 09/08 16:53
4F：→ justlikelbj: 怎麽打呢 09/08 16:53
5F：推 sunev: 把f(t)表为其它五个函数，代回原式，就得到非线性微方 09/08 16:58
6F：→ sunev: 去解那个微方即可。一样代你想要的初始条件 09/08 16:58
7F：→ justlikelbj: 那dydt(6)怎麽办 我用y(6)=sqrt(1-y(1)*y(1)-...y(5) 09/08 19:48
8F：→ justlikelbj: *y(5))换掉前面的 09/08 19:48
9F：→ justlikelbj: 把y(6)代换之後 照理来说不能有dydt(6)吧 09/08 19:49
10F：推 sunev: 嗯，是我搞错了。就如你所说，微方一给，y(6)也就被决定了 09/08 23:47
11F：→ sunev: 所以你的条件还是不能乱给 09/08 23:47