作业帮单摆微分方程求解：x''+(g/l)sin(x)=0,用Matlab求解,function dphi=Pendel_DGL(t,x)g= 9.81;l= 1;dphi_1=x(2);dphi_2=-g/l*sin(x(1));dphi= [dphi_1;dphi_2];[t,x]=ode45(Pendel_DGL,[0,4],[pi/2,0])我觉得思路没有错误呀
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/06 05:51:08
![单摆微分方程求解：x''+(g/l)sin(x)=0,用Matlab求解,function dphi=Pendel_DGL(t,x)g= 9.81;l= 1;dphi_1=x(2);dphi_2=-g/l*sin(x(1));dphi= [dphi_1;dphi_2];[t,x]=ode45(Pendel_DGL,[0,4],[pi/2,0])我觉得思路没有错误呀](/uploads/image/z/6925959-63-9.jpg?t=%E5%8D%95%E6%91%86%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E6%B1%82%E8%A7%A3%EF%BC%9Ax%27%27%2B%28g%2Fl%29sin%28x%29%3D0%2C%E7%94%A8Matlab%E6%B1%82%E8%A7%A3%2Cfunction+dphi%3DPendel_DGL%28t%2Cx%29g%3D+9.81%3Bl%3D+1%3Bdphi_1%3Dx%282%29%3Bdphi_2%3D-g%2Fl%2Asin%28x%281%29%29%3Bdphi%3D+%5Bdphi_1%3Bdphi_2%5D%3B%5Bt%2Cx%5D%3Dode45%28Pendel_DGL%2C%5B0%2C4%5D%2C%5Bpi%2F2%2C0%5D%29%E6%88%91%E8%A7%89%E5%BE%97%E6%80%9D%E8%B7%AF%E6%B2%A1%E6%9C%89%E9%94%99%E8%AF%AF%E5%91%80)
单摆微分方程求解：x''+(g/l)sin(x)=0,用Matlab求解,function dphi=Pendel_DGL(t,x)g= 9.81;l= 1;dphi_1=x(2);dphi_2=-g/l*sin(x(1));dphi= [dphi_1;dphi_2];[t,x]=ode45(Pendel_DGL,[0,4],[pi/2,0])我觉得思路没有错误呀
单摆微分方程求解：x''+(g/l)sin(x)=0,用Matlab求解,
function dphi=Pendel_DGL(t,x)
g= 9.81;
l= 1;
dphi_1=x(2);
dphi_2=-g/l*sin(x(1));
dphi= [dphi_1;dphi_2];
[t,x]=ode45(Pendel_DGL,[0,4],[pi/2,0])
我觉得思路没有错误呀
单摆微分方程求解：x''+(g/l)sin(x)=0,用Matlab求解,function dphi=Pendel_DGL(t,x)g= 9.81;l= 1;dphi_1=x(2);dphi_2=-g/l*sin(x(1));dphi= [dphi_1;dphi_2];[t,x]=ode45(Pendel_DGL,[0,4],[pi/2,0])我觉得思路没有错误呀
1.这段程序基本没有什么错误,只是在最后调用ode45求解时候,格式有点错误,修改一下就能运行了：[t,x]=ode45(@Pendel_DGL,[0,4],[pi/2,0])
2. 在编程时候,在help里面可以找到很多例子,一般来说照着例子来编程,基本不会发生错误,如下：
子程序：
function dx=myfun(t,x)
g=9.81;
l=1;
dx=zeros(2,1);
dx(1)=x(2);
dx(2)=-g/l*sin(x(1));
计算程序：
>> tspan=[0 4];
>> x0=[pi/2 0];
>> [t,x]=ode45(@myfun,tspan,x0);
>>
3. 目前使用匿名函数比较多,而且比较方便：
>> clear
>> g=9.81;
>> l=1;
>> fun=@(t,x)[x(2);-g/l*sin(x(1))];
>> tspan=[0 4];
>> x0=[pi/2 0];
>> [t,x]=ode45(fun,tspan,x0);
>>