matlab求解微分方程组变系数微分方程

来源：蜘蛛抓取(WebSpider) 时间：2015-03-24 09:27 标签： matlab求解方程组
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
实验七__用matlab求解常微分方程 doc
下载积分：400
内容提示：
文档格式：PDF|
浏览次数：276|
上传日期： 19:45:49|
文档星级：
该用户还上传了这些文档
实验七__用matlab求解常微分方程 doc.PDF
官方公共微信matlab求解微分方程组dr/dt=vdv/dt=(.*t)*sin(a)-4.9*10^12/(r^2)+r*w*wdw/dt=-(-2.*t))*cos(a)+2*w*w/r初始条件是r(0)=1749013,v(0)=0,w(0)=0.,t取值0到500然后a是一个变量,我想在运行时再赋值网上查了说用ode_百度作业帮
matlab求解微分方程组dr/dt=vdv/dt=(.*t)*sin(a)-4.9*10^12/(r^2)+r*w*wdw/dt=-(-2.*t))*cos(a)+2*w*w/r初始条件是r(0)=1749013,v(0)=0,w(0)=0.,t取值0到500然后a是一个变量,我想在运行时再赋值网上查了说用ode45,但是本人初学matlab实在是写不出来.求大神帮我写一下程序.
y(3)方程可能有问题,r急剧增加!dy = @(t,y,a) [y(2);(.*t)*sin(a)-4.9*10^12./(y(1).^2)+y(1).*y(3).*y(3);-(-2.*t))*cos(a)+2*y(3).*y(3)./y(1) ];
% 微分方程a =
在运行时再赋值y0 = [,0.];
% 初始条件Tend = 0.3; %
t取值0到0.3[T,Y] = ode45(@(t,y) dy(t,y,a),[0 5],y0);x1=Y(:,1);x2=Y(:,2);x3=Y(:,3);figureplot(T,x1)title(' r ')xlabel('t');ylabel('r');查看: 1198|回复: 6|关注: 0
跪求：四阶变系数微分方程数值求法
各位大侠：
& && & 小弟跪求高手能够帮我解决一个微分方程求解：四阶变系数微分方程。方程具体如下：
方程.jpg (25.51 KB, 下载次数: 0)
四阶变系数微分方程
21:12 上传
如果能够有解析解最好，没有解析解的话，用Matlab等数学工具求解的数值解也可以（求保留代码）。
小弟在此拜谢，并承诺所出论文以第二作者身份作为酬谢。（此方程涉及一个*家级科研项目。）
可以发邮件：，并告诉我您真实身份与单位。
本帖最后由 kNi0rvKOgO 于
23:57 编辑
请给出具体数据，你不会是搞岩土、隧道的吧
用1stOpt，非常容易的事。
youyouyou 发表于
用1stOpt，非常容易的事。
老大啊，啊能够要下你的QQ号啊
请问你的问题解决了没？我也遇到个问题，就是其中的一些系数需要根据变量来判断，然后得出微分方程的系数！
zhuzhubaoibei 发表于
请问你的问题解决了没？我也遇到个问题，就是其中的一些系数需要根据变量来判断，然后得出微分方程的系数！ ...
您好，您的问题解决了吗？我也有这样的问题。
罗同学 发表于
您好，您的问题解决了吗？我也有这样的问题。
我后来用矩阵传递法算了，想直接解出来还是比较麻烦的
站长推荐 /1
Powered byMATLAB&Simulink&求解微分方程&等幅震荡变成数值解衰减
在学习电力系统低频振荡时候，用MATLAB
SIMULIN搞了一下，发现当用ODE45解常微分方程时候，如果步长取如下参数时候，原本应该是等幅震荡的解变成了衰减解。
而如果参数改成如下时候：
以上都是在电力系统低频振荡模型下的仿真：
整个转子运动方程式在受到小干扰时候，在原来工作点附近线性展开，推导出的一个微分方程：
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。

我要回帖

更多关于 matlab求解方程组 的文章

 

随机推荐