求二阶变系数振动微分方程组的Matlab程序
求二阶变系数振动微分方程组的Matlab程序求解振动微分方程的Matlab程序,二阶变系数微分方程组,怎么通过求解绘出频域和时域图,频谱曲线????
非常着急,哪位高手请指教,十分感谢!!!
可以用simulink里的S-Function解决,数据导出就可以绘图了。或者自己编程序积分。
回复 #1 fineshang 的帖子
继续关注这个问题,原来有过这样的帖子,但是没有得到一个好的结论! 就是就是没有好的结论 貌似楼主想要的是MATLAB代码?这个恐怕是见仁见智的事。
MATLAB里面有一个家族的ode求解器,供求解各种问题。大概来说,可以有这么几步
1.在主函数中设置求解器(odeset);
2.在主函数中调用求解器(例如ode45,ode15s等);
3.书写微分方程子函数,即du=f(u,t),此时需考虑变系数问题和需要输出为结果的参数的保存;
4.在主函数中读取结果,并显示。注意利用对全局变量和局部变量的设置。
另外请注意
1.如果问题为非线性问题,常常先求解线性问题的特征解。评估系统方程的数值刚度。参照数值刚度选择合理的求解器;
2.必要时在odeset中使用OutputFcn来控制输出。 今年的我也遇到这样的问题了,找了很长时间也没发现好的解决办法。。。 回复 1 # fineshang 的帖子
不知道求出来了没有啊? 貌似楼主想要的是MATLAB代码?这个恐怕是见仁见智的事。
MATLAB里面有一个家族的ode求解器,供求解各种问题。大概来说,可以有这么几步
1.在主函数中设置求解器(odeset);
2.在主函数中调用求解器(例如ode45,ode15s等);
3.书写微分方程子函数,即du=f(u,t),此时需考虑变系数问题和需要输出为结果的参数的保存;
4.在主函数中读取结果,并显示。注意利用对全局变量和局部变量的设置。
另外请注意
1.如果问题为非线性问题,常常先求解线性问题的特征解。评估系统方程的数值刚度。参照数值刚度选择合理的求解器;
2.必要时在odeset中使用OutputFcn来控制输出。
学习下,能不能详细些或者举个例子! 学习了......
你好,你的这个问题解决了吗
页:
[1]