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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 二阶常微分方程的编写思路
71/1返回列表
查看: 972  |  回复: 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的回帖

angocn

至尊木虫 (文坛精英)

【答案】应助回帖

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

依旧是则徐

铜虫 (小有名气)
引用回帖:
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的回帖

依旧是则徐

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by 月只蓝 at 2014-12-10 13:52:56
类似的思路，先降阶 再求解。...

版主，按照你的帮助，我这样写的程序，但是子程序里的i，j 怎样设置才能使上面传入的i，j 值呢。
!求解Y''(T)=-(i+j*cos(2*T))*Y
program main
  use IMSL
  implicit none
  integer, parameter :: MXPARM = 50
  integer, parameter :: N = 2
  integer IDO, ISTEP, NOUT
  real PARAM(MXPARM), T, TEND, Y(N)
  real, parameter :: TOL = 1E-8
  external FCN
  real :: i,j

  do i=0,1,0.001
     do j=0,1.5,0.1
           write(*,"(I2,I2)" i,j

       PARAM = 0
       WRITE (*,"(' Time', 9X,'Y')"
       Y(1) = 0 ! Initial condition
       Y(2) = 2.0
       T = 0.0    ! Y(T)=4.0中的时间T值
       IDO = 1
         do ISTEP=0,100,10
           TEND = ISTEP
               call IVPRK (IDO, N, FCN, T, TEND, TOL, PARAM, Y)
           write (*,'(F5.1,F12.4)') T, Y
         end do

       call IVPRK (3, N, FCN, T, TEND, TOL, PARAM, Y)
         end do
  end do
   
  stop
  end program
! 计算微分方程式
! YPRIME=Y'=f(t,y)
subroutine FCN (N, T, Y, YPRIME)
  implicit none
  integer N
  real T, Y(N), YPRIME(N)
  real i,j
  !Y''(T)=-(i+j*cos(2*T))*Y
  YPRIME(2) =-(i+j*cos(2*T))*Y(1)
  return
end subroutine
赞一下 回复此楼
6楼2014-12-14 11:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

猪猪会飞
祝福。
7楼2014-12-14 11:47:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 依旧是则徐 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[考博] 博士申请 +4 星…… 2026-05-18 5/250 2026-05-18 17:34 by 炎甲00
[基金申请] 评审有感 +3 popular289 2026-05-18 6/300 2026-05-18 15:17 by gy116024
[基金申请] 今年审到国自然15份，谈谈感受 +16 国自然国社科中 2026-05-17 16/800 2026-05-18 14:58 by gy116024
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 重磅！青年科学基金项目（C类）资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 10/500 2026-05-18 07:50 by 水和泥不是水泥
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-17 19:39 by Equinoxhua
[考博] 光量子物理方向 博士招生 1人（2026.09） +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[考博] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:11 by 11n4dfd8yn
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[考博] 现在不知道怎么办，感觉很痛苦 +5 qweww 2026-05-11 6/300 2026-05-16 19:57 by Equinoxhua
[有机交流] 求助2，4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[考研] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 11:29 by h3oerqvkv9
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
信息提示
关闭
请填处理意见
关闭