(MATLAB]매트랩을 이용한 DSB-SC 변복조시스템 구현

소개글

(MATLAB]매트랩을 이용한 DSB-SC 변복조시스템 구현에 대한 보고서 자료입니다.

목차

1. message 신호에 대하여 m(t), |M(f)|의 출력을 보이시오. 신호의 bandwidth 및 평균파워를 계산해보시오.
2. s(t), |S(f)|의 출력을 보이시오. s(t)의 bandwidth를 계산해보시오.
3. 채널 통과 이후 신호 r(t)의 출력을 보이시오.
4. y(t), |Y(f)|의 출력을 보이시오.
5. (t), |(f)|의 출력을 보이시오.
6. Receiver를 포락선 검파기로 교체 했을 경우의 (t) 출력을 보이시오. 문제 5의 결과와 비교해 보고 차이가 있다면 이류를 설명하시오.

7. Conventional AM (amplitude modulation) 변조를 할 경우 진폭 변조된 파형은 다음과 같다.

여기서 는 amplitude sensitivity constant.
그림 2의 message 신호에 대하여 conventional AM 변조 후, 수신부에서 포락선 검파기를 이용하여 복조시 envelope distortion이 일어나지 않게 하기 위한 송신단의 의 범위를 설정하고 그 이유를 설명하시오.
MATLAB 코드

본문내용

ata object.
pwr = avgpower(hpsd)
average_power=mean(m.*m)
%위의 두 공식을 통해 평균파워를 계산할 수 있다.
% m(t)를 plot
figure('Position',[0 300 650 650]);
subplot(2,1,1), plot(t,m,'b');
hold on;
title('Message Signal(t)');
xlabel('time domain(s), Tb=25 msec');
ylabel('m(t)');
xlim([0 8.5*tb]);
set(gca,'XTick',0:tb:8*tb);
set(gca,'XTickLabel',{'0','Tb','2Tb','3Tb','4Tb','5Tb','6Tb','7Tb','8Tb'});
text(0.15,2,{'average power of m(t)=',num2str(average_power)});
ylim([-3.5 3.5]);
%FT
nfft=length(t);
Mft=fftshift(fft(m,nfft+1));
mag_Mft=abs(Mft)/nfft;
freq=-fs/2:fs/nfft:fs/2; % 주파수 범위 설정
subplot(2,1,2), plot(freq,mag_Mft,'r');
axis([-fs/2,fs/2,0,1]);
title('FT of m(t)');
xlabel('frequency domain(Hz)');
ylabel('Magnitude');
%s(t)
s=m.*cos(2*pi*fc*t);
figure('Position',[0 300 650 650]);
subplot(2,1,1), plot(t,s,'b');
hold on;
title('s(t)');
xlabel('time domain(s), Tb=25 msec');
ylabel('s(t)');
xlim([0 8.5*tb]);
set(gca,'XTick',0:tb:8*tb);
set(gca,'XTickLabel',{'0','Tb','2Tb','3Tb','4Tb','5Tb','6Tb','7Tb','8Tb'})
ylim([-3.5 3.5]);
%FT of s(t)
Sft=fftshift(fft(s,nfft+1));
mag_Sft=abs(Sft)/nfft;
subplot(2,1,2), plot(freq,mag_Sft,'r');
axis([-fs/2,fs/2,0,1]);
title('FT of s(t)');
xlabel('frequency domain(Hz)');
ylabel('Magnitude');
%r(t)
r=awgn(s,snr);
figure('Position',[0 300 650 650]);
plot(t,r);
hold on;
title('r(t)');
xlabel('time domain(s)');
ylabel('r(t)');
xlim([0 8.5*tb]);
set(gca,'XTick',0:tb:8*tb);
set(gca,'XTickLabel',{'0','Tb','2Tb','3Tb','4Tb','5Tb','6Tb','7Tb','8Tb'})
ylim([-3.5 3.5]);
%y(t)
y=r.*2.*cos(2*pi*fc*t);
figure('Position',[0 300 650 650]);
subplot(2,1,1), plot(t,y,'b');
hold on;
title('y(t)');
xlabel('time domain(s), Tb=25 msec');
ylabel('y(t)');
xlim([0 8.5*tb]);
set(gca,'XTick',0:tb:8*tb);
set(gca,'XTickLabel',{'0','Tb','2Tb','3Tb','4Tb','5Tb','6Tb','7Tb','8Tb'})
ylim([-7 7]);
%FT of y(t)
Yft=fftshift(fft(y,nfft+1));
mag_Yft=abs(Yft)/nfft;
subplot(2,1,2), plot(freq,mag_Yft,'r');
axis([-fs/2,fs/2,0,1]);
title('FT of y(t)');
xlabel('frequency domain(Hz)');
ylabel('Magnitude');
% after LPF
H_t=fir1(48,0.02);
filtered_m=filter(H_t,1,y);
figure('Position',[0 300 650 650]);
subplot(2,1,1),plot(t,filtered_m,'b');
hold on;
title('m~(t)');
xlabel('time domain(s), Tb=25 msec');
ylabel('m~(t)');
xlim([0 8.5*tb]);
set(gca,'XTick',0:tb:8*tb);
set(gca,'XTickLabel',{'0','Tb','2Tb','3Tb','4Tb','5Tb','6Tb','7Tb','8Tb'});
hold on;
%FT of m~(t)
filtered_Mft=fftshift(fft(filtered_m,nfft+1));
mag_filtered_Mft=abs(filtered_Mft)/nfft;
subplot(2,1,2), plot(freq,mag_filtered_Mft,'r');
axis([-fs/2,fs/2,0,1]);
title('FT of m~(t)');
xlabel('frequency domain(Hz)');
ylabel('Magnitude');
% envelope detection by using a hilbert Tranform
hil_m=hilbert(r);
envelope_r = sqrt(hil_m.*conj(hil_m));
figure('Position',[0 300 650 650]);
plot(t,envelope_r);
hold on;
title('The Envelope Detection');
xlabel('time domain(s), Tb=25 msec');
ylabel('Envelope');
xlim([0 8.5*tb]);
set(gca,'XTick',0:tb:8*tb);
set(gca,'XTickLabel',{'0','Tb','2Tb','3Tb','4Tb','5Tb','6Tb','7Tb','8Tb'})
ylim([-4 4]);

키워드

매트랩,   통신시스템,   MATLAB,   Communication system,   DSB-SC,   아날로그,   변복조,   변조