| 查看: 1882 | 回复: 9 | ||||
ysjuuplj银虫 (小有名气)
|
[求助]
求问FFT变换是怎么整的？ 已有2人参与
|
我想把一个方程的时域解变换到频域，，我不会编程，求各位编一下图中这个问题简单的FFT程序，不胜感激 1.jpg [ Last edited by ysjuuplj on 2014-3-10 at 19:15 ] |
回复此楼» 猜你喜欢
一志愿北交大材料工程总分358求调剂
已经有8人回复
-
专硕304找调剂，一线城市最好
已经有3人回复
-
一志愿南航，数一英一学硕317求调剂！！
已经有6人回复
-
295求调剂
已经有12人回复
-
285求调剂
已经有10人回复
-
材料求调剂
已经有12人回复
-
0703化学调剂325分
已经有13人回复
-
08600生物与医药-327
已经有8人回复
-
调剂
已经有10人回复
-
085600材料与化工301分求调剂院校
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 新手求问 NGX,Y,Z 和 kpoints 有关联吗？
已经有6人回复
- 利用matlab进行傅里叶变换时，fft函数不需要输入透镜焦距
已经有6人回复
- 各位大神，有个问题请教，数据类型的类似正弦波的信号，能在matlab中进行ＦＦＴ变换吗
已经有8人回复
- TEM图去背景？？
已经有11人回复
- 单个数列的傅里叶变换（FFT）
已经有6人回复
- 快速傅立叶变换（FFT）详解，例子和程序，好资源
已经有31人回复
- 敲击杯子的声音，经fft变换后有很多个频率值，应该取哪个频率值？
已经有14人回复
- FFT变换，是否能输出想要的某个频率结果？
已经有12人回复
- 求助：谁用过origin傅里叶变换图（FFT)
已经有7人回复
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
ysjuuplj2楼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
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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 |
» 本帖已获得的红花（最新10朵）
4楼2014-03-12 21:53:40
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 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
ysjuuplj
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1081.5
- 帖子: 59
- 在线: 27.4小时
- 虫号: 830038
- 注册: 2009-08-18
- 性别: MM
- 专业: 冠状动脉性心脏病
6楼2014-03-13 11:48:53
ysjuuplj
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1081.5
- 帖子: 59
- 在线: 27.4小时
- 虫号: 830038
- 注册: 2009-08-18
- 性别: MM
- 专业: 冠状动脉性心脏病
7楼2014-03-13 11:54:04
8楼2014-03-13 12:02:47
ysjuuplj
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1081.5
- 帖子: 59
- 在线: 27.4小时
- 虫号: 830038
- 注册: 2009-08-18
- 性别: MM
- 专业: 冠状动脉性心脏病
9楼2014-03-13 14:11:46
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3880
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
|
(Y‘)^3<>Re{ [iexp(iz)]*[iexp(iz)]*[iexp(iz)] } 应该是 sin(wt)=Re[-i*exp(iwt)] InvRe[y] 代表复数化y；所谓的相量分析 Y= exp(iz)，则Y’= i*exp(iz), Y’’= - exp(iz)，InvRe[sin(Ohmiga*t)]= -i*exp(i*Ohmiga*t) 原方程被复化， - exp(iz) + exp(iz)+0.1* i*exp(iz)+0.2*exp(iz)* exp(-iz)* exp(iz)=-i*exp(i*Ohmiga*t) 0.1* i*exp(iz)+0.2*exp(-2y)* exp(iz)= -i*exp(i*Ohmiga*t) 0.1* i*exp(ix-y)+0.2*exp(-2y)* exp(ix-y)= -i*exp(i*Ohmiga*t) 0.1* exp(i(*x+pi/2)-y)+0.2* exp(ix-3y)= -i*exp(i*Ohmiga*t) 0.1* exp(-y)*exp(i(*x+pi/2))+0.2*exp(-3y)* exp(ix)=exp(i*(Ohmiga*t-pi/2)) 取 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-pi/2)) [-a*sin(x)+b*cos(x)] + i*[a*cos(x)+b*sin(x)]= exp(i*(Ohmiga*t-pi/2)) 推出 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)] =-ctan(Ohmiga*t) [a+b*tan(x)]/ [-a*tan(x)+b] =-ctan(Ohmiga*t) tan(x)= [b* ctan(Ohmiga*t)+a]/ [a* ctan(Ohmiga*t)-b] 最后，y(t)= Re[Y]=exp(-y)cos(x)=1.70162*cos(x) =1.70162/sqrt{1+[b* ctan(Ohmiga*t)+a]^2/ [a* ctan(Ohmiga*t)-b]^2} [ Last edited by feixiaolin on 2014-3-14 at 09:39 ] |
10楼2014-03-13 14:27:33
5