matlab代码产生单音干扰并做图
时间: 2023-08-02 15:06:27 浏览: 143
可以使用MATLAB中的sin函数生成一个单频率的正弦波来产生单音干扰。以下是一个简单的 MATLAB 代码示例,产生一个频率为1kHz的正弦波作为单音干扰并绘制时域图和频谱图:
```
fs = 10e6; % 采样率
T = 1/fs; % 采样间隔
L = 1024; % 信号长度
t = (0:L-1)*T; % 时间向量
s = sin(2*pi*100e3*t); % 基带信号
s_interfere = s + 0.1*sin(2*pi*1e3*t); % 叠加单音干扰信号
% 绘制时域图
figure(1);
subplot(2,1,1);
plot(t, s);
xlabel('时间(秒)');
ylabel('幅度');
title('基带信号');
subplot(2,1,2);
plot(t, s_interfere);
xlabel('时间(秒)');
ylabel('幅度');
title('叠加单音干扰后的信号');
% 绘制频谱图
nfft = 2^nextpow2(L); % FFT长度
S = fft(s, nfft)/L; % 基带信号的FFT
f = fs/2*linspace(0,1,nfft/2+1); % 频率向量
S_interfere = fft(s_interfere, nfft)/L; % 叠加干扰信号的FFT
figure(2);
subplot(2,1,1);
plot(f, 2*abs(S(1:nfft/2+1)));
xlabel('频率(Hz)');
ylabel('幅度');
title('基带信号频谱');
subplot(2,1,2);
plot(f, 2*abs(S_interfere(1:nfft/2+1)));
xlabel('频率(Hz)');
ylabel('幅度');
title('叠加单音干扰后的信号频谱');
```
在上述代码中,使用 `sin` 函数产生1kHz的正弦波,并将其乘以 `0.1` 叠加到基带信号上作为干扰信号。使用 `subplot` 函数分别绘制了基带信号和叠加干扰信号的时域图和频谱图。在频谱图中,使用 `fft` 函数对信号进行FFT变换,并使用 `abs` 函数计算幅度谱。在实际应用中,需要根据具体情况进行参数设置和调整。
阅读全文