| 查看: 1767 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
ysjuuplj银虫 (小有名气)
|
[求助]
求问FFT变换是怎么整的？ 已有2人参与
|
|||
我想把一个方程的时域解变换到频域，，我不会编程，求各位编一下图中这个问题简单的FFT程序，不胜感激 1.jpg [ Last edited by ysjuuplj on 2014-3-10 at 19:15 ] |
回复此楼» 猜你喜欢
求助:我三月中下旬出站，青基依托单位怎么办？
已经有8人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 新手求问 NGX,Y,Z 和 kpoints 有关联吗？
已经有6人回复
- 利用matlab进行傅里叶变换时，fft函数不需要输入透镜焦距
已经有6人回复
- 各位大神，有个问题请教，数据类型的类似正弦波的信号，能在matlab中进行ＦＦＴ变换吗
已经有8人回复
- TEM图去背景？？
已经有11人回复
- 单个数列的傅里叶变换（FFT）
已经有6人回复
- 快速傅立叶变换（FFT）详解，例子和程序，好资源
已经有31人回复
- 敲击杯子的声音，经fft变换后有很多个频率值，应该取哪个频率值？
已经有14人回复
- FFT变换，是否能输出想要的某个频率结果？
已经有12人回复
- 求助：谁用过origin傅里叶变换图（FFT)
已经有7人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
feixiaolin: 金币+50, reward 2014-03-14 08:58:49
感谢参与，应助指数 +1
feixiaolin: 金币+50, reward 2014-03-14 08:58:49
|
按照你的思路，尝试了下， 对你的ODE，时间区间[0,60], 初值（0.0，0.0），然后OMEGA=[0.1:0.01:3]。 解ODE采用比RK方法精度更高的Bulirsch-Stoer方法，均采样1024个时域信号.，对输出的结果直接构造y2。 做实值FFT，然后求最大值。 求出的结果抖动很厉害，这个我也不知道（加大时间区间也是如此）。第一张图是原始数据，第二张是bezier曲线光滑所致。 这是C++写的，gnuplot绘制，Matlab应该更方便，不过我不熟悉Matlab，  只能祝好运呵呵。 11.png  22.png |
ysjuuplj
4楼2014-03-12 21:53:40
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2014-03-10 21:38:01
ysjuuplj
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1081.5
- 帖子: 59
- 在线: 27.4小时
- 虫号: 830038
- 注册: 2009-08-18
- 性别: MM
- 专业: 冠状动脉性心脏病
3楼2014-03-11 18:40:11
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
ysjuuplj: 金币+100, 虽然看不懂，但还是把分给你吧，我没用FFT，用另一种方法画出来了 2014-03-13 23:17:43
感谢参与，应助指数 +1
ysjuuplj: 金币+100, 虽然看不懂，但还是把分给你吧，我没用FFT，用另一种方法画出来了 2014-03-13 23:17:43
|
* 可以求出解析解 令 Y=InvRe[y]=exp(iz)= exp(iz)，其中z=x(t)+iy(t); InvRe[y] 代表复数化y Y= exp(iz)，则Y’= i*exp(iz), Y’’= - exp(iz) 原方程被复化， - exp(iz) + exp(iz)+0.1* i*exp(iz)+0.2*exp(iz)* exp(-iz)* exp(iz)=exp(i*Ohmiga*t) 0.1* i*exp(iz)+0.2*exp(-2y)* exp(iz)=exp(i*Ohmiga*t) 0.1* i*exp(ix-y)+0.2*exp(-2y)* exp(ix-y)=exp(i*Ohmiga*t) 0.1* exp(i(*x+pi/2)-y)+0.2* exp(ix-3y)=exp(i*Ohmiga*t) 0.1* exp(-y)*exp(i(*x+pi/2))+0.2*exp(-3y)* exp(ix)=exp(i*Ohmiga*t) 取 a=0.1* exp(-y), b=0.2* exp(-3y) 上式变为 a*[cos(x+pi/2)+i*sin(x+pi/2)]+b[cos(x)+i*sin(x)]= exp(i*Ohmiga*t) [-a*sin(x)+b*cos(x)] + i*[a*cos(x)+b*sin(x)]= exp(i*Ohmiga*t) 推出 1） A^2+b^2=1, [0.1* exp(-y)]^2 + [0.2* exp(-3y)]^2=1, y(t)= -0.531582 = constant a=0.170262, b=0.985415 2） [a*cos(x)+b*sin(x)]/ [-a*sin(x)+b*cos(x)] =tan(Ohmiga*t) [a+b*tan(x)]/ [-a*tan(x)+b] =tan(Ohmiga*t) 最后，y(t)= Re[Y]=exp(-y)cos(x)=1.70162*cos(x) |
» 本帖已获得的红花（最新10朵）
5楼2014-03-13 08:29:22