| 查看: 2234 | 回复: 3 | ||
[求助]
二阶微分方程的数值解的傅里叶变换
|
我用NDsolve求得了二阶微分方程的数值解，然后定义并画出了函数p[t]的图形，如下图，我现在想求p[t]的傅里叶变换，没有 函数p[t]关于变量t的具体表达式，而且其值也不是离散的，该用什么语句呢? 先谢谢啦~  |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求数学高手 下面的微分方程的数值解 急急急
已经有3人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 【求助】求2010年注册化工工程师基础考试的考试大纲
已经有4人回复
- 【求助】一个二阶常微分方程的解法（重奖金币！）
已经有16人回复
- 【求助】大家谈谈对注册环保工程师的了解
已经有15人回复
- 【信息】环保工程师报考条件
已经有11人回复
2楼2012-08-28 13:47:59
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 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
