| 查看: 2471 | 回复: 3 | |||
[求助]
二阶微分方程的数值解的傅里叶变换
|
我用NDsolve求得了二阶微分方程的数值解，然后定义并画出了函数p[t]的图形，如下图，我现在想求p[t]的傅里叶变换，没有 函数p[t]关于变量t的具体表达式，而且其值也不是离散的，该用什么语句呢? 先谢谢啦~  |
回复此楼» 猜你喜欢
科研孤儿太难了
已经有19人回复
-
能否退出参与的面上项目解除限项
已经有9人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
-
范进中举一文的中心思想
已经有6人回复
-
建议基金发布提前给出明确的时间点
已经有14人回复
-
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求数学高手 下面的微分方程的数值解 急急急
已经有3人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 【求助】求2010年注册化工工程师基础考试的考试大纲
已经有4人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
- 【求助】大家谈谈对注册环保工程师的了解
已经有15人回复
- 【信息】环保工程师报考条件
已经有11人回复
2楼2012-08-28 13:47:59
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4677.7
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★
xzhdty: 金币+1, 谢谢参与 2012-08-31 22:36:42
xzhdty: 金币+1, 谢谢参与 2012-08-31 22:36:42
|
傅立叶变换要求知道函数在整个实轴上都有定义的了 你NDSolve给出的解 x 范围到不了无穷 不够做变换吧 1. 应该可以用离散的变换 或者 人工的延拓 2. 知道在边界收敛的话。。。。可以把求解区间取大 其他取为0 延拓 3. 或者做傅立叶级数展开 4. 或者你对方程先做傅立叶变换 然后求解变换后的方程 一个例子 这个例子里面的延拓不满足原来的微分方程 Clear["Global`*"]; rule1 = NDSolve[y'[x] + 2 y[x] == 0 && y[0] == 1, y[x], {x, 0, 10}] Plot[y[x] /. rule1[[1]], {x, 0, 1}] yy[x_ /; 0 <= x <= 1] = y[x] /. rule1[[1]]; yy[x_ /; x < 0] = 0; yy[x_ /; x > 1] = 0; yybar2[w_?NumericQ] := FourierTransform[yy[x], x, w] yybar2[1.0] // N Plot[Re[yybar2[xx]], {xx, 1, 2}] 这程序里面数值积分精度有点问题 可能需要重新设置下参数 或者自己代入到傅立叶变换表达式中可能更方便 |
» 本帖已获得的红花（最新10朵）
Lucy3楼2012-08-30 11:32:24
4楼2012-08-31 01:47:12
40