作者zxvc (zxvc)
看板MATLAB
标题Re: [问题] 求解微分方程式的问题
时间Wed Mar 22 16:59:06 2006
※ 引述《brucel3 (bruce N N  )》之铭言:
: 请问板上的各位高手,
: 小弟在求解微分方程式的时候
: 假设我有一个a向量 a=[1,2,3]
: 一个b向量 b=[3,4,5]
: 然後解一个微分方程式 D2y+a*y=sin(b*t)
: 初始条件为 y(0)=0 and Dy(0)=0
: 我的程式为
: a=[1,2,3];
: b=[3,4,5];
: for ii=1:3
: answer=dsolve('D2y+a(ii)*y=sin(b(ii)*t)','y(0)=0','Dy(0)=0');
这一行请修改如下:
~~~~~~~~~~~~~~~~
answer(ii)=
dsolve(strcat('D2y+',int2str(a(ii)),'*y=sin(',int2str(b(ii)),'*t)'),
'y(0)=0','Dy(0)=0');
~~~~~~~~~~~~~~~~
: end
: 理想中的解应该是由'数值'来表达的解
: 但matlab的解却是由 a(ii),b(ii) 来表示
: 例如 他的解为 sin(a(ii)*t)*b(ii)/a(ii)
: 感谢板上大大的回答 ~~
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.115.203.56