4QAM信号时域图，频域图，时频域图，功率谱图通过MATlab产生

生成4QAM信号：

fc = 50;       %载频
fs = 1000;     %采样率
t = 0:1/fs:1;  %时间向量
m = cos(2*pi*10*t); %调制信号
I = m.*cos(2*pi*fc*t); %I路信号
Q = m.*sin(2*pi*fc*t); %Q路信号
s = I + Q;    %4QAM信号

时域图：

subplot(2,2,1);
plot(t,s);
xlabel('时间');
ylabel('幅度');
title('4QAM信号时域图');

频域图：

N = length(s);
fft_s = fft(s,N);
f = fs*(0:N-1)/N;
subplot(2,2,2);
plot(f,abs(fft_s));
xlabel('频率');
ylabel('幅度');
title('4QAM信号频域图');

时频域图：

subplot(2,2,3);
spectrogram(s,256,[],[],fs,'yaxis');
title('4QAM信号时频域图');

功率谱图：

subplot(2,2,4);
pwelch(s,[],[],[],fs,'centered','power');
title('4QAM信号功率谱图');

相关问题

4QAM信号通过MATlab产生时域图，频域图，时频域图，功率谱图

1. 时域图：
可以使用MATLAB中的plot函数绘制4QAM信号的时域波形。假设4QAM信号的调制方式为正交振幅调制（QAM）, 则可以使用MATLAB中的qammod函数生成4QAM信号。以下是一个简单的MATLAB代码示例：

% 生成4QAM信号
n = 1000;     % 信号点数
M = 4;        % 调制阶数
x = randi([0 M-1],n,1);  % 随机生成调制符号
s = qammod(x,M);         % 4QAM调制

% 绘制时域波形
t = linspace(0,1,n);     % 时间轴
figure
plot(t,real(s));
hold on
plot(t,imag(s));
xlabel('时间')
ylabel('幅度')
legend('实部','虚部')

2. 频域图：
可以使用MATLAB中的fft函数将4QAM信号转换到频域，然后使用MATLAB中的plot函数绘制频域波形。以下是一个简单的MATLAB代码示例：

% 生成4QAM信号
n = 1000;     % 信号点数
M = 4;        % 调制阶数
x = randi([0 M-1],n,1);  % 随机生成调制符号
s = qammod(x,M);         % 4QAM调制

% 绘制频域波形
fs = 1000;              % 采样率
f = linspace(-fs/2,fs/2,n); % 频率轴
S = fftshift(fft(s,n))/n;    % 4QAM信号的频谱
figure
plot(f,abs(S));
xlabel('频率')
ylabel('幅度')

3. 时频域图：
可以使用MATLAB中的spectrogram函数绘制4QAM信号的时频谱图。以下是一个简单的MATLAB代码示例：

% 生成4QAM信号
n = 1000;     % 信号点数
M = 4;        % 调制阶数
x = randi([0 M-1],n,1);  % 随机生成调制符号
s = qammod(x,M);         % 4QAM调制

% 绘制时频谱图
fs = 1000;          % 采样率
window = 128;       % 窗口大小
overlap = 64;       % 重叠大小
nfft = 256;         % FFT点数
spectrogram(s,window,overlap,nfft,fs,'yaxis');

4. 功率谱图：
可以使用MATLAB中的periodogram函数绘制4QAM信号的功率谱图。以下是一个简单的MATLAB代码示例：

% 生成4QAM信号
n = 1000;     % 信号点数
M = 4;        % 调制阶数
x = randi([0 M-1],n,1);  % 随机生成调制符号
s = qammod(x,M);         % 4QAM调制

% 绘制功率谱图
fs = 1000;          % 采样率
[Pxx,f] = periodogram(s,[],[],fs);
figure
plot(f,Pxx)
xlabel('频率')
ylabel('功率谱')