翻译下列MATLAB代码%10hz采样 fs = 10; Ts = 1 / fs; n = -0.2 : 1/fs : 0.2;%采样点 f_10Hz = 0.5* sin(2 * pi * f1 * n) - cos(2 * pi * f2 * n); F_10Hz = f_10Hz * exp(-1i * n' * W) * Ts; F_10Hz = abs(F_10Hz); figure(6); subplot(4,3,4);stem(n,f_10Hz);title('10Hz采样信号');grid on; subplot(4,3,5);plot(W,F_10Hz);title('10Hz采样信号频谱');grid on;
时间: 2023-06-20 14:07:49 浏览: 143
数字信号处理大作业之50Hz的干扰噪声的采样的MATLAB完整代码
% 10Hz采样
fs = 10; % 采样频率为10Hz
Ts = 1 / fs; % 采样时间间隔
n = -0.2 : 1/fs : 0.2; % 采样点
% 生成信号
f1 = 1; % 第一个正弦波频率为1Hz
f2 = 2; % 第二个余弦波频率为2Hz
f_10Hz = 0.5 * sin(2 * pi * f1 * n) - cos(2 * pi * f2 * n);
% 计算频谱
W = 2 * pi * (-fs/2 : fs/length(n) : fs/2 - fs/length(n)); % 频域范围
F_10Hz = f_10Hz * exp(-1i * n' * W) * Ts; % 傅里叶变换
F_10Hz = abs(F_10Hz); % 取模
% 绘图
figure(6);
subplot(4,3,4);
stem(n, f_10Hz);
title('10Hz采样信号');
grid on;
subplot(4,3,5);
plot(W, F_10Hz);
title('10Hz采样信号频谱');
grid on;
阅读全文