목차
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 코드
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 코드
본문내용
data 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]);
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]);
추천자료
전파통신시스템
[매트랩]FM 통신 시스템
[통신시스템]매트랩을 이용한 DSB-SC(Double Side Band-Suppressed Carrier)변복조 구현
[광섬유][광통신][광섬유 정의][광섬유 연혁][광섬유 이점][광통신 구현 전제조건][광통신시...
통신시스템 프로젝트 - MATLAB을 이용한 DSBSC 변복조 구현
통신시스템 DMB 수신보드 제어 및 통신 최종설계보고서
1. 컴퓨터의 각종 기능설명, 2. 이동통신시스템의 세대별 정리(1~5세대), 3. 바이오인식의 특...
QAM 방식의 통신시스템 MATLAB 코딩 M-file
디지털 통신시스템의 구성
소개글