有人用Matlab作過離散傅立葉轉換麼﹖

stubor

不是一般連續信號的FFT,是discrete FFT, 可以用frontal analyse的。。$考虑$ $考虑$
希望有人做過啊。。。頭疼中。。:mad:

flight

去这里看看：

http://bbs.matwav.com/index.jsp

stubor

原帖由 flight 于 2007-4-16 13:00 发表
去这里看看：

http://bbs.matwav.com/index.jsp

多謝。。。

eisenstange

or you can see this page and ask LZ sometime

https://www.dolc.de/forum/viewthr ... &extra=page%3D1

stubor

原帖由 eisenstange 于 2007-4-16 14:40 发表
or you can see this page and ask LZ sometime

https://www.dolc.de/forum/viewthr ... &extra=page%3D1

多謝多謝了。。我應該先搜索下的。。。$害羞$ $害羞$

熊猫羊

frontal analyse?还是Friour analyse?

stubor

原帖由 熊猫羊 于 2007-4-16 22:33 发表
frontal analyse?还是Friour analyse?

frontal analyse....use SHAH convolution
因為我

leiwang81

看了半天也不知道楼主到底有什么问题，不是每个人都学你这个专业的，chromatography是什么意思我都不知道。你要稍微详细一点描述你的问题，别人才能帮得上你的忙呀！
顺便指出，连续信号是没有FFT的，DFT和FFT只是效率问题，是同一个东西。
你说要分离信号，如果你不知道怎么分离的话，我下面具个例子，希望对你又帮助。代码是matlab写的。
有一个离散信号 x = [ 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 ];
用fft频谱。
fx = fft(x);
用一下语句绘出频谱的幅值 fig1.jpg
figure;
plot(abs(fx),'o-');
现在分离出两组信号，一组是低频信号，一组是高频信号。
这组是低频信号
fxnl = zeros(1,32);
fxnl([1:2 32]) = fx([1:2 32]);
这组是高频信号
fxnh = zeros(1,32);
fxnh(3:31) = fx(3:31);
用反变换还原信号
xnl = ifft(fxnl);
xnh = ifft(fxnh);
绘出信号 fig2.jpg
figure;
plot(1:32, x, 1:32, xnl, 1:32, xnh, 1:32, (xnh+xnl),'o');
grid on;
legend('x_n','x_{new} with Low Pass', 'x_{new} with High Pass', 'Superposition');
这里最后一个用小圆圈显示的是两组分离出来的信号的叠加。由于我提取的信号，包含了全部的频谱，所以在时域里两个分离出来的信号的叠加就是原信号。

stubor

lsd﹐我知道沒人跟我一個專業﹐沒人搞我這個﹐因為我這個目前我查出?淼奈墨I全世界也就一個人﹐是我教授的徒弟曾經搞過。所以我才頭疼啊。
我就是在這裡喊一喊﹐然後學習一下fft的一些基本知識﹐然後還要自己去琢磨怎麼寫程序。。。
你要是有興趣﹐我把曾經搞過這個的這人名字告訴你﹐呵呵。
Yien C Kwok﹐
沒有很具體的例子﹐因為本身就是理論研究。:(

stubor

忘記告訴大家了﹐我是學化學的。所以搞不懂這可?旱膄ft啊﹗﹗$怒吼$ 教授?G給我一本基本上什麼也沒說的論文﹐看了我一個月那些公式也啃不懂。。$郁闷$
chromatography就是色譜的意思。
謝謝ls的xd。。。正在學習中。。

[ 本帖最后由 stubor 于 2007-4-18 11:13 编辑 ]

月亮的背面

本来还以为就我一个学化学的搞不懂FFT呢, 还好有LZ陪我$害羞$

stubor

原帖由 月亮的背面 于 2007-4-18 11:36 发表
本来还以为就我一个学化学的搞不懂FFT呢, 还好有LZ陪我$害羞$

呵呵。你現在也搞這個啦﹖$握手$ $握手$