% ch5example7prog1.m
clear;
Fs=10000;      % 仿真的采樣率
t=1/Fs:1/Fs:1; % 仿真時間點
m_t(Fs*1)=0;   % 基帶信號變量初始化
for F=150:400  % 基帶信號發生：頻率150Hz～400Hz
    m_t=m_t+0.003*sin(2*pi*F*t)*(400-F); % 幅度隨線性遞減
end
m_t90shift=imag(hilbert(m_t)); % 基帶信號的希爾伯特變換
carriercos=cos(2*pi*1000*t);   % 1000Hz載波 cos
carriersin=sin(2*pi*1000*t);   % 1000Hz正交載波 sin
s_SSB1=m_t.*carriercos-m_t90shift.*carriersin; % 上邊帶SSB
s_SSB2=m_t.*carriercos+m_t90shift.*carriersin; % 下邊帶SSB
% 下面作出各波形以及頻譜
figure(1);
subplot(4,2,1); plot(t(1:100),carriercos(1:100),...
t(1:100),carriersin(1:100),'--r'); % 載波
subplot(4,2,2); plot([0:9999],abs(fft(carriercos))); % 載波頻譜
axis([0 2000 -500 6000]);
subplot(4,2,3); plot(t(1:100),m_t(1:100));           % 基帶信號
subplot(4,2,4); plot([0:9999],abs(fft(m_t)));        % 信號頻譜
                axis([0 2000 -500 6000]);
subplot(4,2,5); plot(t(1:100),s_SSB1(1:100));    % SSB波形上邊帶
subplot(4,2,6); plot([0:9999],abs(fft(s_SSB1))); % SSB頻譜上邊帶
                axis([0 2000 -500 6000]);
subplot(4,2,7); plot(t(1:100),s_SSB2(1:100));    % SSB波形下邊帶
subplot(4,2,8); plot([0:9999],abs(fft(s_SSB2))); % SSB頻譜下邊帶
axis([0 2000 -500 6000]);