목차
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 신호의 전송 대역폭을 계산해 보시오.
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
%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
추천자료
Adaptive Modulation Matlab Source Code
FIR Filter using MATLAB
IIR Filter using MATLAB
가스 터빈 / 추진기관 2 성능 곡선 그래프 및 매트랩 코드
매틀랩 ISI 분석
MATLAB 매틀랩 ode solver 비교
MATLAB ode solver의 특징
생산 경영 - 향후 수요 예측 과 그에 대한 MATLAB 프로그래밍
고체역학 보의 처짐 (이론,Matlab)
디지털신호처리 설계프로젝트 #1 - 신호처리 프로젝트 matlab을 이용한 음성신호 분석
MATLAB 입 출력 함수
matlab 사용법 - MATLAB 을 이용한 제어시스템 해석 및 설계
MATLAB에 의한 ANALOG 신호처리
matlab 레포트
소개글