用matlab编程实现四阶龙格库塔解二元二阶微分方程组

用matlab编程实现四阶龙格库塔解二元二阶微分方程组,初始条件方程如图其中a=b=0.1,a0=b0=c0=2急求具体程序!谢谢... 用matlab编程实现四阶龙格库塔解二元二阶微分方程组,初始条件

方程如图

其中a=b=0.1,a0=b0=c0=2

急求具体程序!谢谢
展开
 我来答
tianxiawulang
推荐于2018-03-02 · TA获得超过2.7万个赞
知道大有可为答主
回答量:4732
采纳率:89%
帮助的人:2649万
展开全部

求解二阶微分方程,初始条件还需要给出y1'(0)和y2'(0)。这里暂时按照0处理。

 

function zd530003514

a=0.1;

b=0.1;

Y0 = [b-1; 0; b; 0];

 

% 解方程

[t,Y]= ode45(@ode,[0 10],Y0);

y1=Y(:,1);

y2=Y(:,3);

 

% 绘图 

subplot 211

plot(t,y1);

subplot 212

plot(t,y2);

 

% 微分方程定义 

function dY = ode(t, Y)

L1=5;

L2=0.01;

a0=2;

b0=2;

c0=2;

 

y1=Y(1);y2=Y(3);

dY = [

    Y(2);

    -(a0*y2+b0*y2^2+c0*y2^3) - L1^2*L2*y1 - L1^2*y1;

    Y(4);

    -(a0*y2+b0*y2^2+c0*y2^3) - L1^2*L2*y1;

];

 

本回答被提问者和网友采纳
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
东莞大凡
2024-11-14 广告
标定板认准大凡光学科技,专业生产研发厂家,专业从事光学影像测量仪,光学投影测量仪.光学三维测量仪,光学二维测量仪,光学二维测量仪,光学三维测量仪,光学二维测量仪.的研发生产销售。东莞市大凡光学科技有限公司创立于 2018 年,公司总部坐落于... 点击进入详情页
本回答由东莞大凡提供
scxczyyl
2013-03-06 · TA获得超过389个赞
知道小有建树答主
回答量:433
采纳率:0%
帮助的人:265万
展开全部
function dy=fun(t,y)
a0=2;
b0=2;
c0=2;
lamda1=5;
lamda2=0.01;
dy=zeros(4,1);
dy(1)=y(2);
dy(2)=-(a0*y(2)+b0*y(2)^2+c0*y(2)^3)-lamda1^2*lamda2*y(1)-lamda1^2*y(1);
dy(3)=y(4);
dy(4)=-(a0*y(2)+b0*y(2)^2+c0*y(2)^3)-lamda1^2*lamda2*y(1);
[T,Y]=ode45('fun',[0 100],[0 0 0.1 0]);
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
maliang351
2013-03-06 · TA获得超过239个赞
知道答主
回答量:196
采纳率:100%
帮助的人:120万
展开全部
matlab解微分方程,函数调用的形式是一样的,把方程转换为代码就可以了。四阶龙格库塔应该是众多解法中的一个。详细参阅help文件。
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
匿名用户
2013-03-06
展开全部
你也弄毕设呢?
追问
是的,所以很急
追答
我也用MATLAB做毕设,不过你这个比我的可难多了
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 2条折叠回答
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式