（1）将BPSK换成16QAM 要求：画出频谱图，星座图，以及解码后的时域波形图。要求使用matlab

我理解了您的问题，下面是对您问题的回答：

在Matlab中，将BPSK调制方式换成16QAM的操作可以通过以下步骤实现：

1. 定义调制符号点

16QAM调制是一种基于正交幅度调制（QAM）的数字调制方式，使用16个符号进行数据传输。我们可以定义一个包含所有16个符号点的QAM调制表。在matlab中，可以使用“qammod”函数来实现。这个函数有两个输入参数，其中一个是调制阶数（例如，我们这里使用16），另一个是需要调制的数据向量。

例如，如果我们有一个长度为1000的数据向量“data”，可以使用以下代码进行调制：

symbols = qammod(data, 16);

2. 绘制星座图和频谱图

接下来，我们可以使用“scatterplot”函数在星座图中绘制调制之后的符号点。

scatterplot(symbols);

对于频谱图，我们可以通过使用“fft”函数和“plot”函数来实现。首先，我们需要将调制符号点转换成基带信号。我们可以使用“modulate”函数来执行此操作。此函数需要三个输入参数：调制符号点，载波频率和采样率。

Fs = 10000; %采样率
carrier_freq = 2000; %载波频率
modulated_signal = modulate(symbols, carrier_freq, Fs);

然后我们可以使用FFT函数来获取信号的频谱，在使用plot函数将频谱可视化。

N = length(modulated_signal);
freq = linspace(-Fs/2,Fs/2,N);
f_signal = fftshift(fft(modulated_signal));
plot(freq,abs(f_signal));

3. 解码后的时域波形图

为了解调信号，我们需要首先将信号乘以载波信号，并在经过低通滤波器之后进行采样。在这里，我们使用fir1（）函数来设计一个零相位带通滤波器，并使用filter函数进行滤波。最后，我们可以使用“demodulate”函数来解调信号

carrier = cos(2*pi*carrier_freq*(0:N-1)/Fs);
received_signal = modulated_signal .* carrier;
filtered_signal = filter(fir1(100,carrier_freq/Fs), 1, received_signal);
demodulated_signal = demod(received_signal, carrier_freq, Fs, 'qam');

为了绘制解码后的时域波形，我们可以使用下面的代码：

plot(demodulated_signal);

希望这个回答能够帮助您解决问题！如果您有任何其他问题，请随时问我。