matlab数据处理、傅里叶变换

最新推荐文章于 2024-10-09 10:15:14 发布

幻夜的梦

最新推荐文章于 2024-10-09 10:15:14 发布

阅读量712

点赞数 10

文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/2302_79519348/article/details/136317885

版权

matlab

matlab数据处理

傅里叶变换

在这里插入图片描述

% 生成正弦波信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间向量
f = 50; % 正弦波频率
x = sin(2*pi*f*t); % 正弦波信号

% 进行傅里叶变换
X = fft(x);

% 绘制时域波形图像
subplot(2, 1, 1);
plot(t, x);
title('正弦波信号 - 时域');
xlabel('时间 (s)');
ylabel('幅值');

% 绘制频域幅度谱图像
L = length(X);
frequencies = fs*(0:(L/2))/L;
amplitude = abs(X/L);
subplot(2, 1, 2);
stem(frequencies, amplitude(1:L/2+1));
title('正弦波信号 - 频域');
xlabel('频率 (Hz)');
ylabel('振幅');

在这里插入图片描述

% 生成方波信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间向量
f = 50; % 方波频率
dutyCycle = 0.5; % 占空比
x = square(2*pi*f*t, dutyCycle); % 方波信号

% 进行傅里叶变换
X = fft(x);

% 绘制时域波形图像
subplot(2, 1, 1);
plot(t, x);
title('方波信号 - 时域');
xlabel('时间 (s)');
ylabel('幅值');

% 绘制频域幅度谱图像
L = length(X);
frequencies = fs*(0:(L/2))/L;
amplitude = abs(X/L);
subplot(2, 1, 2);
stem(frequencies, amplitude(1:L/2+1));
title('方波信号 - 频域');
xlabel('频率 (Hz)');
ylabel('振幅');

在这里插入图片描述

% 生成三角波信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间向量
f = 50; % 三角波频率
x = sawtooth(2*pi*f*t, 0.5); % 三角波信号

% 进行傅里叶变换
X = fft(x);

% 绘制时域波形图像
subplot(2, 1, 1);
plot(t, x);
title('三角波信号 - 时域');
xlabel('时间 (s)');
ylabel('幅值');

% 绘制频域幅度谱图像
L = length(X);
frequencies = fs*(0:(L/2))/L;
amplitude = abs(X/L);
subplot(2, 1, 2);
stem(frequencies, amplitude(1:L/2+1));
title('三角波信号 - 频域');
xlabel('频率 (Hz)');
ylabel('振幅');

在这里插入图片描述

% 傅里叶逆变换
% 生成一个简单的信号
fs = 1000; % 采样率
t = 0:1/fs:1-1/fs; % 时间向量
f1 = 10; % 第一个频率成分
f2 = 50; % 第二个频率成分
x = sin(2*pi*f1*t) + 0.5*sin(2*pi*f2*t); % 组合两个正弦波信号

% 进行傅里叶变换
X = fft(x);

% 进行傅里叶逆变换
y = ifft(X);

% 绘制原始信号与逆变换后的信号
subplot(2,1,1);
plot(t, x);
title('原始信号 - 时域');
xlabel('时间 (s)');
ylabel('幅值');
subplot(2,1,2);
plot(t, y);
title('逆变换后的信号 - 时域');
xlabel('时间 (s)');
ylabel('幅值');

% 显示频谱图
figure;
frequencies = fs*(0:length(X)-1)/length(X); % 频率向量
plot(frequencies, abs(X));
title('信号的频谱');
xlabel('频率 (Hz)');
ylabel('振幅');