使用matlab代码，用相移法生产单边带调制，并显示出单边带信号的频谱与调制信号的频谱，给出程序代码

以下是使用MATLAB代码实现移相法产生单边带调制（LSSB）并绘制频谱的示例代码：

% 设置信号参数
Fc = 100;           % 载波频率
Fs = 1000;          % 采样频率
T = 1/Fs;           % 采样时间间隔
t = 0:T:1-T;        % 采样时间序列
Am = 1;             % 调制信号幅度
Fm = 10;            % 调制信号频率
phi = pi/2;         % 相位偏移
Ac = 2;             % 载波幅度

% 生成调制信号
m_signal = Am*cos(2*pi*Fm*t); 

% 生成载波信号
c_signal = Ac*cos(2*pi*Fc*t);

% 移相法产生LSSB信号
LSSB_signal = m_signal.*cos(2*pi*Fc*t) - imag(hilbert(m_signal)).*sin(2*pi*Fc*t+phi);

% 绘制信号频谱
N = length(t);
f = Fs*(0:(N/2))/N;
M_signal = fft(m_signal)/N;
C_signal = fft(c_signal)/N;
LSSB_signal_f = fft(LSSB_signal)/N;
LSSB_signal_f_shift = fftshift(LSSB_signal_f);

figure;
subplot(3,2,1); plot(t,m_signal); xlabel('时间（s）'); ylabel('幅度'); title('调制信号');
subplot(3,2,2); plot(f,2*abs(M_signal(1:N/2+1))); xlabel('频率（Hz）'); ylabel('幅度'); title('调制信号频谱');
subplot(3,2,3); plot(t,c_signal); xlabel('时间（s）'); ylabel('幅度'); title('载波信号');
subplot(3,2,4); plot(f,2*abs(C_signal(1:N/2+1))); xlabel('频率（Hz）'); ylabel('幅度'); title('载波信号频谱');
subplot(3,2,[5,6]); plot(f,2*abs(LSSB_signal_f_shift(1:N/2+1))); xlabel('频率（Hz）'); ylabel('幅度'); title('LSSB信号频谱');

这段代码首先设置了一些信号参数，然后生成了调制信号和载波信号。接着，使用移相法将调制信号产生LSSB信号。最后，使用MATLAB的`fft`函数计算信号的频谱，并使用`subplot`将三个信号的时间域波形和频谱绘制在同一张图中。其中，LSSB信号的频谱使用了`fftshift`函数进行了中心化处理。