matlab求解微分方程组

来源:百度知道 编辑:UC知道 时间:2024/05/04 12:17:23
论文中需要解一个方程组得到y变量的表达式,好像用dsolve解不出,请教各位高人,要怎样编程才能解出表达式。

Ay''+p(H-x)(y+z)=q
z'=By'''

初始条件: y(0)=0; y'(0)=0; y''(H)=0; y'''(H)=C

其中,A,B,C,p,q,H都是常系数。
需要求y关于x的表达式。

多谢各位了。
syms A B C p q H x;
[y,z]=dsolve('A*D2y+p*(H-x)*(y+z)=q','Dz=B*D3y','y(0)=0','Dy(0)=0','D2y(H)=0','D3y(H)=C','x')

这样不知道行不行,3个小时了。结果都没出来。

因为z'=By'''
所以z=By''+C
把它代入第一个方程里,再用dsolve解,从matlab给出的答案可以看到matlab无法得到这个方程的解析解,所以你还是放弃的好。

网上找的关于利用Matlab解微分方程和方程组的资料,可以参考一下
http://mcm.ustc.edu.cn/download/matlab30.pdf