matlab解偏微分方程请把matlab语句给出来~

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/28 22:13:20

matlab解偏微分方程
请把matlab语句给出来~

这只是求偏导数而已,你别搞笑哈.
给你代码
syms x y z
z =(x^2+y^2)*exp((x^2+y^2)/x/y);
z_x=diff(z,x)
z_xx=diff(z,x,2)
z_xy=diff(diff(z,x),y)
pretty(z)
pretty(z_x)
pretty(z_xx)
pretty(z_xy)
输出(这儿可能串行,在matlab中很好看的):
z_x =
2*x*exp((x^2+y^2)/x/y)+(x^2+y^2)*(2/y-(x^2+y^2)/x^2/y)*exp((x^2+y^2)/x/y)
z_xx =
2*exp((x^2+y^2)/x/y)+4*x*(2/y-(x^2+y^2)/x^2/y)*exp((x^2+y^2)/x/y)+(x^2+y^2)*(-2/x/y+2*(x^2+y^2)/x^3/y)*exp((x^2+y^2)/x/y)+(x^2+y^2)*(2/y-(x^2+y^2)/x^2/y)^2*exp((x^2+y^2)/x/y)
z_xy =
2*x*(2/x-(x^2+y^2)/x/y^2)*exp((x^2+y^2)/x/y)+2*y*(2/y-(x^2+y^2)/x^2/y)*exp((x^2+y^2)/x/y)+(x^2+y^2)*(-2/y^2-2/x^2+(x^2+y^2)/x^2/y^2)*exp((x^2+y^2)/x/y)+(x^2+y^2)*(2/y-(x^2+y^2)/x^2/y)*(2/x-(x^2+y^2)/x/y^2)*exp((x^2+y^2)/x/y)
2 2
2 2 x + y
(x + y ) exp(-------)
x y
2 2 / 2 2\ 2 2
x + y 2 2 | x + y | x + y
2 x exp(-------) + (x + y ) |2/y - -------| exp(-------)
x y | 2 | x y
\ x y /
2 2 / 2 2\ 2 2
x + y | x + y | x + y
2 exp(-------) + 4 x |2/y - -------| exp(-------)
x y | 2 | x y
\ x y /
/ 2 2\ 2 2
2 2 | 2 x + y | x + y
+ (x + y ) |- --- + 2 -------| exp(-------)
| x y 3 | x y
\ x y /
/ 2 2\2 2 2
2 2 | x + y | x + y
+ (x + y ) |2/y - -------| exp(-------)
| 2 | x y
\ x y /
/ 2 2\ 2 2 / 2 2\ 2 2
| x + y | x + y | x + y | x + y
2 x |2/x - -------| exp(-------) + 2 y |2/y - -------| exp(-------)
| 2 | x y | 2 | x y
\ x y / \ x y /
/ 2 2\ 2 2
2 2 | 2 2 x + y | x + y
+ (x + y ) |- ---- - ---- + -------| exp(-------)
| 2 2 2 2 | x y
\ y x x y /
/ 2 2\ / 2 2\ 2 2
2 2 | x + y | | x + y | x + y
+ (x + y ) |2/y - -------| |2/x - -------| exp(-------)
| 2 | | 2 | x y
\ x y / \ x y /