Question

用matlab求解微分方程组

t的范围为[0,pi/6] x(0)=0, y(0)=x(pi/6) (y的初值等于t=pi/6时x的值) (如果采用龙格库塔法ode45，初值条件是否可以是非确定值，如:x(0)=0, y(0)=x(pi/6) ;如果堪称边值(BVP)问题，该如何编写程序? ) 请问各位大神，如何求解，并画出图形

Answer 1

首先
建立M文件:
function
dx=myfun(t,x)
dx=zeros(2,1);
dx(1)=x(1)*cos(t)+5*x(2)*cos(t+pi/6);
dx(2)=5*x(1)+5*x(2)*cos(t+pi/6);
主窗口:
(这里y(0)=x(pi/6)我取的是1,,初值必须是已知的，，所以你可以大致算下它)
[t,x]=ode45('myfun',[0
pi/6],[0
1]);
Y=0:0.01:pi/6;
plot(x(:,1),x(:,2))
ode解微分求值错，，要好好看的,呵呵