中華民國萬萬歲

close all;clear all;clc; % Prepare to sample a signal for two seconds, % at a rate of 100 samples per second. Fs = 100; % Sampling rate t = [0:2*Fs+1]/Fs; % Time points for sampling % Create the signal, a sum of sinusoids. x = sin(2*pi*t)+sin(8*pi*t); Fc = 10; % Carrier frequency in modulation phasedev = pi/2; % Phase deviation for phase modulation y = ssbmod(x,Fc,Fs,phasedev); % Modulate. % y = awgn(y,10,'measured',103); % Add noise. z = ssbdemod(y,Fc,Fs,phasedev); % Demodulate % Plot the original and recovered signals. figure; plot(t,x,'k-',t,y,'r--',t,z,'b--'); legend('x Original signal','y Recovered signal','z');