小弟現在有6組二元7階PDE組(其中4個由GR來的所以是雙曲PDE,2個是Laplace)
擁有的DOF是Qx,Qz,Qt,Qc,F,K
變數是x,z
由於D[F[x,z],z]/.z->0是sin func.
所以想說可以對x做discrete F.T.來簡化問題
可是PDE的order太高了,做FT會有很多convolution
偏偏Mathematica並不會自動幫你把FT[f[x]g[x]]換成Convolution[FT[f[x]],FT[g[x]]]
每次做FT我都得要手動把A*B換成Sum[A[k],B[m-k],k]實在很麻煩
不知道有沒有什麼辦法能夠自動地把相乘換成discrete convolution