Question

关于如何用matlab做傅里叶变换的问题

手头要处理数据，每个数据可看作是一段纵波。波的图像由origin处理得到。现在想把图像导入matlab中，进行傅里叶变换，从而得到傅里叶系数。这样就可以将其中不同振幅的波分离出来，从而求出其占整体的百分比。

目前这些波的图像都没有函数式，只是图像，是否需要先拟合，然后根据解析式进行变换。

我想了很久，查了很多书，可书上的东西有些太难了，我只是需要一个很简单的频谱分析就可以了。大侠帮帮忙吧。
谢谢你的回答，你给出的程序我在计算机上试一下。但是我没学过数字信号处理。看来要不不少课啊。

Answer 1

不需要解析式，你学过数字信号处理的话看看matlab的fft函数就应该会用

fs=10000;%采样频率
N=1000; %采样点数
n=0:N-1;
t=n/fs;
y=sin(2*pi*t)+0.1*sin(4*pi*t)+0.01*sin(8*pi*t);
%给定信号
figure(1);
subplot(211);
plot(t,x);%绘出信号的时域波形
xlabel('时间/t');
ylabel('幅值');
title('信号的时域波形');
grid on;
y=fft(x,N);%进行fft变换
m=abs(y(1:N/2))*2/N;%求信号的真实幅值
f=n*fs/N; %进行对应的频率转换
subplot(212);
stem(f(1:N/2),m(1:N/2));%绘出频谱图
xlabel('频率/Hz');
ylabel('幅值');
title('信号频谱图');
grid on;

Answer 2

不需要拟合 ：
设抽样频率为Fs（Hz）,信号点数为N，信号序列为x。
x=imread（‘’）
f = fftshift(fft(x));
w = linspace(-Fs/2, Fs/2, N);%频率坐标，单位Hz
plot(w,abs(f));
title('信号的频谱');
xlabel('频率(Hz)');

Answer 3

fs=10000;%采样频率
N=1000; %采样点数
n=0:N-1;
t=n/fs;
y=sin(2*pi*t)+0.1*sin(4*pi*t)+0.01*sin(8*pi*t);
%给定信号
figure(1);
subplot(211);
plot(t,x);%绘出信号的时域波形
xlabel('时间/t');
ylabel('幅值');
title('信号的时域波形');
grid on;
y=fft(x,N);%进行fft变换
m=abs(y(1:N/2))*2/N;%求信号的真实幅值
f=n*fs/N; %进行对应的频率转换
subplot(212);
stem(f(1:N/2),m(1:N/2));%绘出频谱图
xlabel('频率/Hz');
ylabel('幅值');
title('信号频谱图');
grid on;