应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 二阶常微分方程的编写思路
51/1返回列表
查看: 910  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

依旧是则徐

铜虫 (小有名气)

[求助] 二阶常微分方程的编写思路 已有2人参与

求高手解决一下方程的编写思路 最好详细一些。或者给出相似的计算例子也行
y··（t）+2y·（t)+y=0   初始条件y（0）=4    y·（0）=2
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-12-09 12:11:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
令y1=y y2=y'
原二阶常微分初值问题 化为 一阶常微分方程组初值问题
化成方程组：
y1'=y2
y2'=-(2*y2+y1)

初值条件：
y1(0)=4, y2(0)=2

解析解见附图1。

数值解MATLAB代码：
CODE:
function solve_odes
clear all;clc
[t y]=ode45(@odes,[0:0.05:20],[4 2]);
figure(33)
plot(t,y(:,1),'b-s',t,y(:,2),'r-d'),axis([-0.5 20 -2 4.5]),legend('y1','y2')


function f=odes(t,y)
f(1)=y(2);
f(2)=-(2*y(2)+y(1));
f=f';

二阶常微分方程的编写思路
附图1.png


二阶常微分方程的编写思路-1
附图2.png
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2014-12-09 21:56:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

angocn

至尊木虫 (文坛精英)

【答案】应助回帖

感谢参与，应助指数 +1
这个随便翻开一本数值分析的教材都能找到吧！甚至有的还带源程序。
赞一下 回复此楼
其实我喜欢历史，可是却学了理科
2楼2014-12-09 13:23:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

依旧是则徐

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by 月只蓝 at 2014-12-09 21:56:02
令y1=y y2=y'
原二阶常微分初值问题 化为 一阶常微分方程组初值问题
化成方程组：
y1'=y2
y2'=-(2*y2+y1)

初值条件：
y1(0)=4, y2(0)=2

解析解见附图1。

数值解MATLAB代码：

function solve_odes
...

谢谢，能再请问一个问题吗，如果是V``(x)+(a+bcos2x）v=0这样的方程编写思路是怎样的，其中a，b 未知，可以变化。我会追加金币的。非常感谢！！！
赞一下 回复此楼
4楼2014-12-10 10:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
依旧是则徐: 金币+15, ★★★★★最佳答案 2014-12-14 16:01:42
引用回帖:
4楼: Originally posted by 依旧是则徐 at 2014-12-10 10:09:25
谢谢，能再请问一个问题吗，如果是V``(x)+(a+bcos2x）v=0这样的方程编写思路是怎样的，其中a，b 未知，可以变化。我会追加金币的。非常感谢！！！...

类似的思路，先降阶 再求解。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2014-12-10 13:52:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭