使用Matlab产生一个需要传输的数字输入信号，计算输入信号功率谱密度，并利用Matlab显示； （2）生成ASK数字调制信号，计算ASK调制信号功率谱密度，并利用Matlab显示； （3）生成2PSK数字调制信号，计算2PSK调制信号功率谱密度，并利用Matlab显示； （4）生成2FSK数字调制信号，计算2FSK调制信号功率谱密度，并利用Matlab显示；

以下是使用Matlab实现的代码，其中输入信号为随机生成的正弦波，采样率为1000Hz，时长为1秒：

% 生成输入信号
fs = 1000; % 采样率
t = 0:1/fs:1-1/fs; % 时间序列
x = sin(2*pi*50*t) + sin(2*pi*120*t); % 生成两个正弦波并相加
figure;
plot(t,x)
xlabel('Time (s)')
ylabel('Amplitude')
title('Input signal')

% 计算输入信号功率谱密度
N = length(x);
xdft = fft(x);
xdft = xdft(1:N/2+1);
psdx = (1/(fs*N)) * abs(xdft).^2;
psdx(2:end-1) = 2*psdx(2:end-1);
freq = 0:fs/length(x):fs/2;
figure;
plot(freq,10*log10(psdx))
grid on
xlabel('Frequency (Hz)')
ylabel('Power/Frequency (dB/Hz)')
title('Input signal power spectral density')

生成ASK数字调制信号，可以使用以下代码，其中调制信号为输入信号乘以二进制序列，并对于二进制序列中的0和1，分别使用不同的载波频率：

% 生成ASK数字调制信号
fc0 = 100; % 载波频率1
fc1 = 200; % 载波频率2
m = randi([0 1],1,N); % 生成二进制序列
s = x.*(m==1)*fc1 + x.*(m==0)*fc0; % 生成ASK数字调制信号
figure;
plot(t,s)
xlabel('Time (s)')
ylabel('Amplitude')
title('ASK modulated signal')

% 计算ASK数字调制信号功率谱密度
N = length(s);
sdft = fft(s);
sdft = sdft(1:N/2+1);
psds = (1/(fs*N)) * abs(sdft).^2;
psds(2:end-1) = 2*psds(2:end-1);
figure;
plot(freq,10*log10(psds))
grid on
xlabel('Frequency (Hz)')
ylabel('Power/Frequency (dB/Hz)')
title('ASK modulated signal power spectral density')

生成2PSK数字调制信号，可以使用以下代码，其中调制信号为输入信号乘以二进制序列，并对于二进制序列中的0和1，分别使用不同的相位：

% 生成2PSK数字调制信号
fc = 100; % 载波频率
m = randi([0 1],1,N); % 生成二进制序列
s = x.*cos(pi*m)*fc; % 生成2PSK数字调制信号
figure;
plot(t,s)
xlabel('Time (s)')
ylabel('Amplitude')
title('2PSK modulated signal')

% 计算2PSK数字调制信号功率谱密度
N = length(s);
sdft = fft(s);
sdft = sdft(1:N/2+1);
psds = (1/(fs*N)) * abs(sdft).^2;
psds(2:end-1) = 2*psds(2:end-1);
figure;
plot(freq,10*log10(psds))
grid on
xlabel('Frequency (Hz)')
ylabel('Power/Frequency (dB/Hz)')
title('2PSK modulated signal power spectral density')

生成2FSK数字调制信号，可以使用以下代码，其中调制信号为输入信号乘以二进制序列，并对于二进制序列中的0和1，分别使用不同的频率：

% 生成2FSK数字调制信号
f0 = 100; % 频率1
f1 = 200; % 频率2
m = randi([0 1],1,N); % 生成二进制序列
s = x.*cos(2*pi*f0*t + pi*m)*fc + x.*cos(2*pi*f1*t + pi*(m-1))*fc; % 生成2FSK数字调制信号
figure;
plot(t,s)
xlabel('Time (s)')
ylabel('Amplitude')
title('2FSK modulated signal')

% 计算2FSK数字调制信号功率谱密度
N = length(s);
sdft = fft(s);
sdft = sdft(1:N/2+1);
psds = (1/(fs*N)) * abs(sdft).^2;
psds(2:end-1) = 2*psds(2:end-1);
figure;
plot(freq,10*log10(psds))
grid on
xlabel('Frequency (Hz)')
ylabel('Power/Frequency (dB/Hz)')
title('2FSK modulated signal power spectral density')