【答案】应助回帖

» 收录本帖的淘帖专辑推荐

博士发论文的思考体会

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

二阶偏微分方程组matlab求解已经有4人回复
求大神解决matlab算积分方程组！急已经有6人回复
利用matlab求一阶变系数微分方程组的数值解。已经有7人回复
偏微分方程Matlab求解已经有7人回复
Matlab求解微分方程(组)及偏微分方程(组) 已经有224人回复
matlab 非线性微分方程求解已经有3人回复
求助非线性常微分方程和常系数双曲偏微分方程求解已经有7人回复
寻求方程组的解法，急急急！！已经有7人回复
用matlab 求解参数方程组已经有8人回复
偏微分方程求解！已经有4人回复
二阶非常系数非线性微分方程求解已经有4人回复
三元二次的三角方程如何用MATLAB求解？？已经有10人回复
matlab如何求解一个非线性微分方程组已经有8人回复
matlab求解高阶方程组已经有3人回复
求解一个二阶偏微分方程~~有结果，但是不知道中间过程（完美解决者答谢100￥）~ 已经有4人回复
matlab处理常微分方程作图问题已经有7人回复
复系数微分方程求解已经有3人回复
分数阶微分方程的求解（用差分法，matlab编程求解）已经有9人回复
二阶非线性微分方程求解已经有14人回复
如何用matlab求解矩阵系数的二阶微分方程已经有13人回复
求助一阶非线性微分方程的解已经有7人回复
如何求解描述振动的二阶微分方程已经有21人回复
Matlab求解偏微分方程组已经有22人回复
matlab求解非线性方程组已经有16人回复
请问如何求解二元一阶微分方程组已经有24人回复
【求助】常微分方程组的解析解已经有3人回复
【分享】PDECOL-经典数值求解非线性偏微分方程组源程序已经有142人回复

把一切交给时间！

1楼2014-03-25 15:58:00

回帖支持 ( 显示支持度最高的前 50 名 )

george1990

金虫 (小有名气)

应助: 11 (小学生)
金币: 1155.3
红花: 2
帖子: 130
在线: 75.3小时
虫号: 2142478
注册: 2012-11-22
性别: GG
专业: 机械动力学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lylpiao: 金币+10, ★有帮助, 谢谢 2014-03-30 11:47:01

这种问题你可以参考机械动力学里面的4阶龙格库塔解法它对于解这种方程很ok的

赞一下(1人)

天天多学一点

5楼2014-03-30 01:11:47

george1990

金虫 (小有名气)

应助: 11 (小学生)
金币: 1155.3
红花: 2
帖子: 130
在线: 75.3小时
虫号: 2142478
注册: 2012-11-22
性别: GG
专业: 机械动力学

【答案】应助回帖

★ ★ ★ ★ ★
lylpiao: 金币+5, 哈，那我明天也去试一下，有结果再告诉你 2014-04-01 23:11:10

这个我当时试了的消不消去都没什么问题两种结果是一样的至于多方程的情况我没试过

赞一下(1人)

天天多学一点

11楼2014-04-01 17:37:47

普通回帖

WWJ1991

木虫 (正式写手)

应助: 34 (小学生)
金币: 2417.5
散金: 66
红花: 10
帖子: 337
在线: 195.5小时
虫号: 2423760
注册: 2013-04-18
专业: 机构学与机器人

【答案】应助回帖

感谢参与，应助指数 +1

参考：http://zhidao.baidu.com/link?url ... le8kHjgTQAnHcUx-EMa。
百度知道给的是一阶的求法，你换成二阶就行了。可以将每一阶微分的分量写在一起，组成矩阵形式。

2楼2014-03-25 16:48:07

daiyouming

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 1.5
帖子: 3
在线: 3.2小时
虫号: 2716773
注册: 2013-10-11

【答案】应助回帖

感谢参与，应助指数 +1

引用回帖:

2楼: Originally posted by WWJ1991 at 2014-03-25 16:48:07
参考：http://zhidao.baidu.com/link?url=n-6qzdYoFHcIxmVsFUFbKrknOkbLdZstCBx1IgbFeDhTH3rdU8gqbP1xQ5hPe-FAZ6Wle8kHjgTQAnHcUx-EMa。
百度知道给的是一阶的求法，你换成二阶就行了。可以将每一阶微分的分量写在 ...

您好，请问你的意思是变成24X24矩阵吗,然后一个个一起求解，有这种参考案例吗:

3楼2014-03-25 18:39:05

WWJ1991

木虫 (正式写手)

应助: 34 (小学生)
金币: 2417.5
散金: 66
红花: 10
帖子: 337
在线: 195.5小时
虫号: 2423760
注册: 2013-04-18
专业: 机构学与机器人

【答案】应助回帖

★ ★ ★ ★ ★
lylpiao: 金币+5, ★有帮助, 谢谢 2014-03-30 11:47:20

引用回帖:

3楼: Originally posted by daiyouming at 2014-03-25 18:39:05
您好，请问你的意思是变成24X24矩阵吗,然后一个个一起求解，有这种参考案例吗:

:...

不是，还是12X12的。二阶微分用d2y表示啊。

4楼2014-03-25 20:44:17

lylpiao

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6761
红花: 1
帖子: 391
在线: 233.6小时
虫号: 1758002
注册: 2012-04-16
性别: MM
专业: 微/纳机械系统

ben_ladeng: 求助不能申请程序强帖，应助和精彩讲座可以申请程序强帖。欢迎发表高质量作品 2014-03-31 06:57:26

引用回帖:

5楼: Originally posted by george1990 at 2014-03-30 01:11:47
这种问题你可以参考机械动力学里面的4阶龙格库塔解法它对于解这种方程很ok的

请问下，该如何求解啊，送你一朵鲜花哈

H{W}]M5NY@[C48~~W@42B}6.jpg

把一切交给时间！

6楼2014-03-30 11:43:44

george1990

金虫 (小有名气)

应助: 11 (小学生)
金币: 1155.3
红花: 2
帖子: 130
在线: 75.3小时
虫号: 2142478
注册: 2012-11-22
性别: GG
专业: 机械动力学

【答案】应助回帖

function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(30*(300*sin(t)-(0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3))))/2999;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;

close all;clc;
y0=[0 0 0 0]; %x初值和x导数的初值
tspan=[0 3];
options = odeset('RelTol',1e-8,'AbsTol',1e-2);
[t y]=ode23s(@lyl,tspan,y0,options);
%plot(t,y(:,2),'r');
我这里设置的初值是随意的导致计算到2.7s就不能计算了你看看你那边的初值是多少因为没给这个我是随便设置的

天天多学一点

7楼2014-03-31 18:46:58

george1990

金虫 (小有名气)

应助: 11 (小学生)
金币: 1155.3
红花: 2
帖子: 130
在线: 75.3小时
虫号: 2142478
注册: 2012-11-22
性别: GG
专业: 机械动力学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lylpiao: 金币+5, ★有帮助, 谢谢~ 2014-04-01 10:37:09
lylpiao: 金币+5, ★有帮助, 谢谢~ 2014-04-01 10:40:34

额第一段程序发错了应该是这样的：
function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;

天天多学一点

8楼2014-03-31 18:53:11

lylpiao

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6761
红花: 1
帖子: 391
在线: 233.6小时
虫号: 1758002
注册: 2012-04-16
性别: MM
专业: 微/纳机械系统

引用回帖:

8楼: Originally posted by george1990 at 2014-03-31 18:53:11
额第一段程序发错了应该是这样的：
function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(10 ...

你好，我有个疑惑，按照你的思路，一个方程里面同时含有二阶量，dy(2)和dy(4)，这可以吗
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;

把一切交给时间！

9楼2014-04-01 10:32:53

lylpiao

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 6761
红花: 1
帖子: 391
在线: 233.6小时
虫号: 1758002
注册: 2012-04-16
性别: MM
专业: 微/纳机械系统

引用回帖:

9楼: Originally posted by lylpiao at 2014-04-01 10:32:53
你好，我有个疑惑，按照你的思路，一个方程里面同时含有二阶量，dy(2)和dy(4)，这可以吗
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2) ...

查阅相关资料：若两个高阶微分方程同时含有最高价导数项，需要先进行相应处理，先消去某个高阶导数。但是方程组少的话计算还行，多的话消去就太麻烦了。所以还是不知道求解。不知道我理解的对不对呢

把一切交给时间！

10楼2014-04-01 11:13:33