傅里叶变换matlab程序实例

下面是一个简单的 Matlab 程序实例，实现对信号的离散傅里叶变换（DFT）：

% 定义信号及采样频率
fs = 1000; % 采样频率为 1000 Hz
t = 0:1/fs:1-1/fs; % 信号持续时间为 1 秒
f1 = 50; % 第一个正弦波频率为 50Hz
f2 = 120; % 第二个正弦波频率为 120Hz
x = 0.7*sin(2*pi*f1*t) + sin(2*pi*f2*t); % 生成信号

% 计算傅里叶变换
N = length(x); % 信号长度
X = fft(x); % 对信号进行快速傅里叶变换

% 计算频率谱
f = (0:N-1)*(fs/N); % 计算频率范围
power = abs(X).^2/N; % 计算功率谱

% 绘制信号波形图和频率谱图
subplot(2,1,1);
plot(t,x);
xlabel('时间 (秒)');
ylabel('幅度');
title('信号波形图');
subplot(2,1,2);
plot(f,power);
xlabel('频率 (Hz)');
ylabel('功率');
title('频率谱图');

这个程序生成了一个包含两个正弦波的信号，然后计算了这个信号的傅里叶变换，并绘制了信号波形图和频率谱图。你可以通过修改生成信号的频率或添加其他信号来尝试不同的情况。

相关问题

傅里叶变换计算幅频、相频特性matlab程序实例

下面是一个计算傅里叶变换的幅频、相频特性的 MATLAB 程序实例：

% 定义信号
t = 0:0.001:1;
f = 10;
x = sin(2*pi*f*t);

% 计算傅里叶变换
N = length(x);
X = fft(x)/N;
f_axis = (0:N-1)*(1/N);
f_axis = f_axis - max(f_axis)/2;

% 计算幅频特性
X_amp = abs(X);
X_amp = X_amp(1:N/2+1);
X_amp(2:end-1) = 2*X_amp(2:end-1);
figure;
plot(f_axis(1:N/2+1), X_amp);
xlabel('Frequency (Hz)');
ylabel('Amplitude');

% 计算相频特性
X_phase = angle(X);
X_phase = X_phase(1:N/2+1);
figure;
plot(f_axis(1:N/2+1), X_phase);
xlabel('Frequency (Hz)');
ylabel('Phase');

这个程序会输出两个图像，第一个是幅频特性图像，第二个是相频特性图像。在幅频特性图像中，横轴表示频率，纵轴表示幅度。在相频特性图像中，横轴同样表示频率，纵轴表示相位。