Dsb3.m
上传用户:loeagle
上传日期:2013-03-02
资源大小:1236k
文件大小:2k
- % dsb1.m
- % Matlab demonstration script for DSB-AM modulation. The message signal
- % is +1 for 0 < t < t0/3, -2 for t0/3 < t < 2t0/3 and zero otherwise.
- echo on
- t0=.15; % signal duration
- ts=0.0005; % sampling interval
- fc=250; % carrier frequency
- snr=10; % SNR in dB (logarithmic)
- fs=1/ts; % sampling frequency
- t_long=20*t0; % extended signal duration
- n0=floor(t_long/ts)+1; % length of extended signal vector
- n00=2^(ceil(log2(n0)));
- df=fs/n00; % frequency resolution
- t_long1=[0:ts:t_long]; % extended time vector
- f=[-fs/2+df:df:fs/2]; % frequency vector
- % message signal
- m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)];
- m_long=[m,zeros(1,n0-length(m))]; % extended message signal
- c_long=cos(2*pi*fc.*t_long1); % extended carrier signal
- u_long=m_long.*c_long; % extended modulated signal
- du_long=u_long.*c_long;
- DU=fft(du_long,n00)/fs;
- M=fft(m_long,n00)/fs; % spectrum of the extended message signal
- U=fft(u_long,n00)/fs; % spectrum of the extended modulated signal
- pause % Press any key to see a plots of the magnitude of the message and the
- % modulated signal in the frequency domain.
- subplot(2,1,1)
- plot(f,abs(fftshift(M)))
- xlabel('Frequency')
- title('Spectrum of the message signal')
- subplot(2,1,2)
- plot(f,abs(fftshift(U)))
- title('Spectrum of the modulated signal')
- xlabel('Frequency')
- pause % Press a key to see the modulated and modulated X carrier in freq. domain
- subplot(2,1,1)
- plot(f,abs(fftshift(U)))
- title('Modulated signal spectrum')
- xlabel('Frequency')
- subplot(2,1,2)
- plot(f,abs(fftshift(DU)))
- title('Modulated signal X Carrier spectrum')
- xlabel('Frequency')