g*****a 发帖数: 340 | 1 用cosine transform解偏微分方程。用periodic边界条件没的说,不用管就好。但如果
加入一定边界条件如何在变换空间中体现出来?
譬如说在边界上dC/dx=0。一般方法中直接代进去就好了。但如果用余弦变换的话整个
区域一起处理,怎么搞?有弄国的么?佛立业变换的情况应该是类似的
可能有点乱,总归有任何提示和意见都请不吝赐教 |
r****y 发帖数: 1437 | 2 I coded similar b.c. and PDE before.
There are two ways, both make use C(x)=Const, i.e., C(x) only
has DC component, wavenumber0, no AC components.
First is to make sure of eigenfunctions for this B.C., note
the different treatments for AC and DC components. And you can easily prove
d^nC/dx^n=0 for n = 2i-1
Second, find the quantity conserved in your equation, dC/dx=0 means
integral over x at B.C. is a value never change with time. Then you can make
use of this to
【在 g*****a 的大作中提到】 : 用cosine transform解偏微分方程。用periodic边界条件没的说,不用管就好。但如果 : 加入一定边界条件如何在变换空间中体现出来? : 譬如说在边界上dC/dx=0。一般方法中直接代进去就好了。但如果用余弦变换的话整个 : 区域一起处理,怎么搞?有弄国的么?佛立业变换的情况应该是类似的 : 可能有点乱,总归有任何提示和意见都请不吝赐教
|
l******n 发帖数: 9344 | 3 边界条件,像你说的这个一般是
C(1)=C(-1)
用cosine transform或者sine transform的实质是对角化所得的矩阵,解了之后,在用
反变幻回去就是解,根据你自己的情况看吧
fftw的document里边有介绍,比如2维有00,01,10,11,不同的选择,就是根据不同的
边值,需要shift一下
【在 g*****a 的大作中提到】 : 用cosine transform解偏微分方程。用periodic边界条件没的说,不用管就好。但如果 : 加入一定边界条件如何在变换空间中体现出来? : 譬如说在边界上dC/dx=0。一般方法中直接代进去就好了。但如果用余弦变换的话整个 : 区域一起处理,怎么搞?有弄国的么?佛立业变换的情况应该是类似的 : 可能有点乱,总归有任何提示和意见都请不吝赐教
|
g*****a 发帖数: 340 | 4 多谢
第二个还是不太明白,如何控制整个域内的C(x)平均值呢?而且'dC/dx=0'应该是'C(x)
在域内总和不变'的充分非必要条件吧?
有没有更进一步的资料我可一看呢?你关于这个方法有文章么?
如果可以的话可以引用的哈,呵呵
prove
means
make
如果
整个
【在 r****y 的大作中提到】 : I coded similar b.c. and PDE before. : There are two ways, both make use C(x)=Const, i.e., C(x) only : has DC component, wavenumber0, no AC components. : First is to make sure of eigenfunctions for this B.C., note : the different treatments for AC and DC components. And you can easily prove : d^nC/dx^n=0 for n = 2i-1 : Second, find the quantity conserved in your equation, dC/dx=0 means : integral over x at B.C. is a value never change with time. Then you can make : use of this to
|
g*****a 发帖数: 340 | 5 多谢
还是不太了
大致来说就是强制最外层节点的值同次外层相同,dC/dx向外层的向量,在次外层就等
于零了?
fftw的document是在matlab里的么?木有找到相关内容哈?能否给个其它链接呢?
如果
整个
【在 l******n 的大作中提到】 : 边界条件,像你说的这个一般是 : C(1)=C(-1) : 用cosine transform或者sine transform的实质是对角化所得的矩阵,解了之后,在用 : 反变幻回去就是解,根据你自己的情况看吧 : fftw的document里边有介绍,比如2维有00,01,10,11,不同的选择,就是根据不同的 : 边值,需要shift一下
|
l******n 发帖数: 9344 | 6 google fftw
it is not matlab
【在 g*****a 的大作中提到】 : 多谢 : 还是不太了 : 大致来说就是强制最外层节点的值同次外层相同,dC/dx向外层的向量,在次外层就等 : 于零了? : fftw的document是在matlab里的么?木有找到相关内容哈?能否给个其它链接呢? : : 如果 : 整个
|
g*****a 发帖数: 340 | 7 thanks!
就等
【在 l******n 的大作中提到】 : google fftw : it is not matlab
|