找回密码
 加入华同
搜索
黄金广告位联系EMAIL:[email protected] 黄金广告[email protected]
查看: 2050|回复: 3

谁会编程解二阶微分方程

[复制链接]
发表于 2013-10-11 23:06:08 | 显示全部楼层 |阅读模式
vb能做吗
发表于 2013-10-12 06:56:57 | 显示全部楼层
matlab做
回复

使用道具 举报

发表于 2014-7-14 00:28:29 | 显示全部楼层
matlab 二阶微分方程8 @7 T3 k$ @2 n1 T7 M6 l" e
/ v4 C- Z# H5 S
>> clear2 z  v+ c& l; n0 `5 `
>> syms a b c d e;
  S9 `8 a2 u' R6 M. ]# D>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
" @) H' w: _" c/ C: P
. g# p9 h$ c# Cy =
+ a8 E, U8 ]5 A- \; H
( o  }" w* q8 y0 {8 ]% W(2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b + (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2))1 }# Q0 y5 e( J3 X3 \: U3 u

" o+ B  o; B6 e3 A& ~; l>> a=1;b=1;c=1;d=1;e=1;%若常数已知
4 q* [! {' L; x
6 v. x  {5 i" K: D>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
  D, O* g; o4 D>> y=eval(y)
; L4 g# N1 N) `( s0 a4 w: k4 _. E* K1 ]5 f) Q/ A1 P# V, L+ ^
y =( _1 ]7 F) W8 }3 u

, O/ d$ e: ]4 o  b   1.1932 + 0.0000i
  _1 [! e3 s9 p3 |这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的。
% w5 z7 u8 t3 I5 o9 j/ |/ T4 Dsyms a b c d e;+ G& m8 e( n. f" l: {! f
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
9 n) H# l- F- B& c# P( jx=200:0.1:400;%取值范围,步长
8 i( F- k, _! _! K8 }# R/ pa=1;b=1;c=1;d=1;e=1;%带入参数
2 d2 g5 `5 {: o2 @3 y7 P6 }5 fy=eval(y);
! _4 @: j  |& D* N6 d4 a! uplot(x,y)
回复

使用道具 举报

发表于 2014-7-31 23:39:02 | 显示全部楼层
大神跑这来问问题。。厉害~!~( O# W( o2 c9 M+ W1 Z6 Z3 O
vb做界面了,matlab做比较轻松愉快~!~
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 加入华同

本版积分规则

Archiver|手机版|小黑屋|华人同志

GMT+8, 2025-12-13 21:05 , Processed in 0.097605 second(s), 5 queries , Redis On.

Powered by Discuz! X3.5

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表