一道二阶常系数线性微分方程求解，求助

雨的印记

求朋友们帮帮忙，我不是数学专业的，但是现在需要用到这个结果，求帮助

liuyalong008

雨的印记 发表于 2014-10-27 17:44
这样可以有人求解么？帮帮忙吧各位大神

1. 可以尝试符号解和数值解法
   
2. % solution1
   
3. f=dsolve('D2y=0.05*sin(3.28859*x)-0.65836*y','y(0)=1','Dy(0)=2','x')
   
4. simplify(f)
   
5. % s2
   
6. F = @(x,y)[y(2);
   
7.            -0.05*sin(3.28859*x) - 0.65836*y(1)]
   
8. [x,y]=ode23(F,[0 10],[1 2]);
   
9. plot(x,y,'ro')

复制代码

雨的印记

我没有学过MATLAB编程，所以，如果有朋友可以帮忙编程的话，一定要帮帮忙，做论文，在这里卡住了，求帮助，谢谢大家了

雨的印记

求大神帮帮忙啊

雨的印记

有路过的朋友会做的帮帮忙啊

雨的印记

这样可以有人求解么？帮帮忙吧各位大神

TBE_Legend

雨的印记 发表于 2014-10-27 17:44
这样可以有人求解么？帮帮忙吧各位大神

算法的发生的发生的费

雨的印记

liuyalong008 发表于 2014-10-29 20:27

谢谢版主了，我之前没学过matlab，所以看起来有点困难，我用的是dsolve这个函数来解决问题，可以搞定，相对容易接受

雨的印记

liuyalong008 发表于 2014-10-29 20:27

谢谢版主了，我之前没学过matlab，所以看起来有点困难，我用的是dsolve这个函数来解决问题，可以搞定，相对容易接受

独孤神宇

雨的印记 发表于 2014-10-27 17:44
这样可以有人求解么？帮帮忙吧各位大神

这不就是一个无阻尼受迫振动么  q(x)=【-0.05/（C-3.28859^2）】sin(3.28859t)

雨的印记

独孤神宇 发表于 2014-10-30 15:54
这不就是一个无阻尼受迫振动么  q(x)=【-0.05/（C-3.28859^2）】sin(3.28859t)

嗯，是的，是一个无阻尼受迫振动