[MATLAB]매트랩을 이용한 FM신호 구현과 확인

소개글

[MATLAB]매트랩을 이용한 FM신호 구현과 확인에 대한 보고서 자료입니다.

목차

MATLAB 실습 2

FM 통신 시스템
1. message 신호에 대하여 m(t), |M(f)|의 출력을 보이시오. 신호의 BW(fm)를 계산해 보시오.
2. B 지점의 출력을 보이시오. △f(frequency deviation), β(modulation index)값을 계산해 보시오. Narrowband FM 인가?
3. C 지점의 출력과 스펙트럼을 보이시오. Carrier Amplitude는 1로 하시오.
4. D 지점의 출력을 보이시오.
5. 복조된 신호의 출력과 스펙트럼을 보이시오.
6. 위 FM 신호의 전송 대역폭을 계산해 보시오.

본문내용

843],'EdgeColor',[0 0 0]);
%s(t) at output C
ac=1;
s=ac*cos(2*pi*fc*t+2*pi*int_m);
figure('Position',[0 300 650 650]);
subplot(2,1,1), plot(t,s,'b');
hold on;
title('output at C');
xlabel('time domain(s), Tb=0.1sec');
ylabel('s(t)');
xlim([0 6*tb]);
set(gca,'XTick',0:tb:6*tb,'xgrid','on','ygrid','on');
set(gca,'XTickLabel',{'0','0.1','0.2','0.3','0.4','0.5','0.6'});
ylim([-1.1 1.1]);
%FT of s(t)
Sft=fftshift(fft(s,nfft+1));
mag_Sft=abs(Sft);
max_Sft=max(mag_Sft);
subplot(2,1,2), plot(freq,mag_Sft,'r');
axis([-2000,2000,0,max_Sft]);
title('FT of s(t)');
xlabel('frequency domain(Hz)');
ylabel('Magnitude');
text(1250,1250,{['-25KHz ~ 25KHz까지 출력시'],['스펙트럼을 정확히'],['알아보기 힘들어'],['-2000Hz ~ 2000Hz까지'],['출력하였습니다.'],['(주파수가 높아질수록'],['0에 가까워짐)']},'FontSize',8,'BackgroundColor',[0.8314 0.8157 0.7843],'EdgeColor',[0 0 0]);
%Phase Demodulator
hil_s=hilbert(s);
angle_s=angle(hil_s);
angle_s=angle_s-2*pi*fc*t;
demod=unwrap(angle_s)/(2*pi);
figure('Position',[0 300 650 650]);
plot(t,demod,'b');
hold on;
title('output at D');
xlabel('time domain(s), Tb=0.1sec');
ylabel('output at D');
xlim([0 6*tb]);
set(gca,'XTick',0:tb:6*tb,'xgrid','on','ygrid','on');
set(gca,'XTickLabel',{'0','0.1','0.2','0.3','0.4','0.5','0.6'});
%Differentiator
received_m=diff(demod)/(kf*ts);
figure('Position',[0 300 650 650]);
subplot(2,1,1), plot(0:ts:0.6-2*ts,received_m,'b');
hold on;
title('Received Message Signal');
xlabel('time domain(s), Tb=0.1sec');
ylabel('Received Message Signal');
xlim([0 6*tb]);
set(gca,'XTick',0:tb:6*tb,'xgrid','on','ygrid','on');
set(gca,'XTickLabel',{'0','0.1','0.2','0.3','0.4','0.5','0.6'});
ylim([-13 13]);
%FT of the received message signal
nfft=length(received_m);
Rft=fftshift(fft(m,nfft+1));
mag_Rft=abs(Rft);
max_Rft=max(mag_Rft);
fm=1/tb;
freq=-fs/2:fs/nfft:fs/2; % 주파수 범위 설정
subplot(2,1,2), plot(freq,mag_Rft,'r');
axis([-200,200,0,max_Rft]);
title('FT of received message signal');
xlabel('frequency domain(Hz)');
ylabel('Magnitude');
text(-180,80000,{['-25KHz ~ 25KHz까지 출력시'],['스펙트럼을 정확히'],['알아보기 힘들어'],['-200Hz ~ 200Hz까지'],['출력하였습니다.'],['(주파수가 높아질수록'],['0에 가까워짐)']},'FontSize',9,'BackgroundColor',[0.8314 0.8157 0.7843],'EdgeColor',[0 0 0]);
%calculate Bt
bt=2*f_dev + 2*fm

키워드

매트랩,   MATLAB,   통신시스템,   FM,   아날로그,   Analog,   전자전기,   에프엠