bpzf 发表于 2009-12-4 09:21:08

求助二元四阶微分方程求解

本帖最后由 TBE_Legend 于 2009-12-11 09:38 编辑

二元四阶微分方程如附件,怎么结果中有积分和复数啊?
后面是代码
omega := f*sin(Pi*m*x/a)*(sin(Pi*y/b))^2;
eq1 := E*((diff(diff(omega, y), x))^2-(diff(diff(omega, x), x))*(diff(diff(omega, y), y)));
eq2:=combine(eq1);
eq3 := diff(diff(diff(diff(F(x, y), x), x), x), x)+2*(diff(diff(diff(diff(F(x, y), y), y), x), x))+diff(diff(diff(diff(F(x, y), y), y), y), y);
pde := eq2 = eq3;
sol := pdsolve(pde,F(x,y));

maplelab 发表于 2009-12-4 10:14:58

最好用代码或者图片说明意图.

feiyuzhen 发表于 2009-12-4 13:03:55

第二张图太小,不能说明问题

maplelab 发表于 2009-12-4 14:04:10

截图真的很小,无法看清.

bpzf 发表于 2009-12-6 21:34:00

顶起来,希望高手帮忙,谢谢

caoer 发表于 2009-12-7 06:59:44

不是所有的pde都有解的,4阶2元出现复数很正常,君不见2元ODE在delta<0的情况。

bpzf 发表于 2009-12-7 09:06:42

那后面那部分怎么是积分型式呢?

maplelab 发表于 2009-12-7 18:09:08

用build选项尝试获取显式.

bpzf 发表于 2009-12-11 09:22:31

各位,我已经把代码上传,希望大家帮帮忙

maplelab 发表于 2009-12-11 13:01:42

9# bpzf

pdetest(sol, pde)

解没有问题.

bpzf 发表于 2009-12-17 10:38:57

哦,谢谢maplelab 了,我再看看

bpzf 发表于 2010-1-20 10:48:26

本帖最后由 bpzf 于 2010-1-20 11:29 编辑

为什么用显式的就计算不出来了,结果用pdetest检验不等于0???
omega := f*sin(Pi*m*x/a)*(sin(Pi*y/b))^2;

eq1 := E*((diff(diff(omega, y), x))^2-(diff(diff(omega, x), x))*(diff(diff(omega, y), y)));

eq2:=combine(eq1);

eq3 := diff(diff(diff(diff(F(x, y), x), x), x), x)+2*(diff(diff(diff(diff(F(x, y), y), y), x), x))+diff(diff(diff(diff(F(x, y), y), y), y), y);

pde := eq2 = eq3;

sol := pdsolve(pde, HINT = `+`, INTEGRATE, build);
页: [1]
查看完整版本: 求助二元四阶微分方程求解