1. Home
  2. Tài liệu khác
  3. Bài tập môn học kỹ thuật truyền dẫn

Bài tập môn học kỹ thuật truyền dẫn - pdf 15

Download miễn phí Bài tập môn học kỹ thuật truyền dẫn
Bài tập môn học: Kỹ Thuật truyền dẫn

Bài I: Biểu diễn tín hiệu của một số mã đường dây và vẽ phổ của chúng.
I.1. Mã nguồn chương trình Matlab:


1.Mã NRZ: file nrz_l.m

function y=nrz_l(x)
for i=1:length(x)
if x(i)==0
temp(i) = -1;
else
temp(i) = +1;
end;
end
y=temp;

2.Mã CMI: file cmi_1.m
function y = cmi(x);
polar = +1; %chon cuc tinh cua bit 1 dau tien la duong
for i=1:length(x)
if x(i)==1
temp(2*i-1) = polar;
temp(2*i) = polar;
polar = -polar; %xac dinh cuc tinh cho bit 1 tiep theo
else
temp(2*i-1) = -1;
temp(2*i) = +1;
end
end
y = temp;


Để tải bản Đầy Đủ của tài liệu, xin Trả lời bài viết này, Mods sẽ gửi Link download cho bạn sớm nhất qua hòm tin nhắn.
Ai cần download tài liệu gì mà không tìm thấy ở đây, thì đăng yêu cầu down tại đây nhé:
Nhận download tài liệu miễn phí

Tóm tắt nội dung tài liệu:

Bài tập môn học: Kỹ Thuật truyền dẫn
Bài I: Biểu diễn tín hiệu của một số mã đường dây và vẽ phổ của chúng.
I.1. Mã nguồn chương trình Matlab:
1.Mã NRZ: file nrz_l.m
function y=nrz_l(x)
for i=1:length(x)
if x(i)==0
temp(i) = -1;
else
temp(i) = +1;
end;
end
y=temp;
2.Mã CMI: file cmi_1.m
function y = cmi(x);
polar = +1; %chon cuc tinh cua bit 1 dau tien la duong
for i=1:length(x)
if x(i)==1
temp(2*i-1) = polar;
temp(2*i) = polar;
polar = -polar; %xac dinh cuc tinh cho bit 1 tiep theo
else
temp(2*i-1) = -1;
temp(2*i) = +1;
end
end
y = temp;
3.Mã HDB3: File hdb3.m
function y = hdb3(x)
next_polar = +1;
b=0; %so bit B giua 2 bit V
i=1;
while i<= length(x)
if x(i)==1
temp(i)=next_polar;
next_polar = -next_polar;
b=b+1;
i=i+1;
else
if (i<=length(x)-3)&(x(i+1)==0) & (x(i+2)==0) & (x(i+3)==0) %truong hop 4 bit 0 lien tiep
if mod(b,2)==1
temp(i) = 0; %bit 0
temp(i+1) = 0; %bit 0
temp(i+2) = 0; %bit 0
temp(i+3) = -next_polar;%bit V
next_polar = -next_polar;
b=0;
i=i+4;
else
temp(i) = next_polar; %bit B
temp(i+1) = 0; %bit 0
temp(i+2) = 0; %bit 0
temp(i+3) = next_polar; %bit V
next_polar = -next_polar;
b=0;
i=i+4;
end
else
temp(i)=0; %truong hop 3 bit ke tiep khac 000
i=i+1;
end
end
end
y=temp;
4. Mã AMI: File ami.m
function y=ami(x)
b=+1;
for i=1: length(x)
if x(i)==1
temp(i)=b;
b=-b;
else
temp(i)=0;
end
end
y=temp;
5. Mã Manchester: File manchester.m
function y= manchester(x)
for i=1 : length(x)
if x(i)==1
temp(2*i-1) = +1;
temp(2*i) = -1;
else
temp(2*i-1) = -1;
temp(2*i) = +1;
end
end
y=temp;
6. Mã Polar RZ: File polar_RZ
function y = polar_RZ(x)
for i=1:length(x)
if x(i)==1
temp(2*i-1) = 1;
temp(2*i) = 0;
else
temp(2*i-1) = -1;
temp(2*i) = 0;
end;
end
y=temp;
7. Mã Unipolar RZ: File unipolar_RZ.m
function y=unipolar_RZ(x)
for i=1:length(x)
if x(i)==1
temp(2*i-1) =1;
temp(2*i) =0;
else
temp(2*i-1) =0;
temp(2*i) =0;
end
end
y=temp;
I.2.Chương trình vẽ và phân tích phổ của tín hiệu ở dạng mã đường dây
File maduongday.m
clear; %Xoa cac bien hien co
bits = 30; %so bit se khao sat
t = 1:30;
data2 = round(rand(1, bits));
%Ma hoa chuoi du lieu goc%
%hdb3 la mot ham duoc dinh nghia o file hdb3.m
%ami la mot ham duoc dinh nghia o file ami.m
%unipolar_RZ la mot ham duoc dinh nghia o file unipolar_RZ.m
disp('chon dang ma duong day');
disp('1: ma NRZ');
disp('2: ma HDB3');
disp('3: ma AMI');
disp('4: ma NRZ_L');
disp('----');
disp('10: RZ');
disp('20: ma polar_RZ');
disp('30: ma Manchester');
disp('40: ma CMI');
linecode = input('chon dang ma se phan tich: ');
switch linecode
%%Truong hop cac ma don bit day du
case {1}
data3 = data2;
TenDoThi='NRZ';
case {2}
data3 = hdb3(data2);
TenDoThi='HDB3';
case {3}
data3=ami(data2);
TenDoThi='AMI';
case {4}
data3=nrz_l(data2);
TenDoThi='NRZ-L';
%Truong hop cac ma co ve khong hay nua bit
case {10}
data3=RZ(data2);
data2=reshape([data2; data2],1,length(data2)*2);
t=0.5:0.5:30;
TenDoThi='RZ';
case {20}
data3=polar_RZ(data2);
data2=reshape([data2; data2],1,length(data2)*2);
TenDoThi='polar_RZ';
t=0.5:0.5:30;
case {30}
data3=manchester(data2);
data2=reshape([data2; data2],1,length(data2)*2);
TenDoThi='Manchester';
t=0.5:0.5:30;
case {40}
data3=cmi(data2);
data2=reshape([data2; data2],1,length(data2)*2);
TenDoThi='CMI';
t=0.5:0.5:30;
otherwise
error('nhap gia tri sai');
end
%Do thi thoi gian cua du lieu goc
subplot(3, 1, 1);
stairs(t,data2)
ylabel('Bien do');
title(TenDoThi);
axis([0 bits+1 -1.2 1.2]);
grid off;
%Do thi thoi gian cua du lieu da ma hoa
subplot(3, 1, 2);
stairs(t,data3);
ylabel(['Bien do']);
axis([0 bits+1 -1.2 1.2]);
grid off;
%Enconding using NRZ-L%
%Mo rong khong gian lay mau%
sampleRate = 2;
if mod(linecode,10)==0
Y1 = reshape(ones(sampleRate, 1)*data3, 1, sampleRate*bits*2);
else
Y1 = reshape(ones(sampleRate, 1)*data3, 1, sampleRate*bits);
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%tinh FFT va ve do thi pho cong suat%
N=sampleRate*bits;
if mod(linecode,10) == 0
FFT_Y1 = abs(fft(Y1, N*2))/N/2;
FFT_Y1 = fftshift(FFT_Y1);
subplot(3, 1, 3);
F = [-N:N-1]./N;
plot(F, FFT_Y1);
xlabel('frequency/F_s voi F_s = 4/T_b= 4*f_C_L_K' )
else
FFT_Y1 = abs(fft(Y1, N))/N;
FFT_Y1 = fftshift(FFT_Y1);
F = [-N/2:N/2-1]./N*2;
subplot(3, 1, 3);
axis([0 1 0 1]);
plot(F, FFT_Y1);
xlabel('frequency/F_s voi F_s = 2/T_b= 2*f_C_L_K' )
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
shg;
%Ket thuc chuong trinh
Một số kết quả cho các trường hợp:
Bài II: Biểu diễn tín hiệu điều chế số
II.1.Mã nguồn chương trình
%ve tin hieu dieu che so bang goc va phan tich pho cua cac tin hieu nay
%cac loai dieu che duoc minh hoa o doan chuong trinh nay bao gom
% 1: ask
% 2: fsk
% 3: psk
% 4: qask
%voi so muc tin hieu la 8
clear;
%thiet laptop cac thong so co ban
M=8; %moi muc ung voi 3 bit tin hieu lien tiep
Fc=20;
Fd=10;
Fs=50;
% lua chon dang dieu che so%
disp('chon mot trong cac phuong thuc dieu che sau');
disp('1: ask (8 muc)');
disp('2: fsk (8 muc)');
disp('3: psk (8 muc)');
disp('4: qask (8 muc)');
cach_dieu_che = input('chon phuong thuc dieu che ');
%Tao day tin hieu ban dau la mot chuoi bit ngau nhien
symbols=5;
bits = symbols*3; % so luong bit trong day tin hieu se dieu che
% 3 bit tin hieu tao ra mot muc tuong ung
data= round(rand(1,bits)); % tao chuoi bit ban dau
%ve do thi tin hieu goc ban dau%
subplot(3,1,1);
stairs([1:bits],data);
axis([1 bits -0.2 1.2] );
subplot(2,1,2)
%tinh toan va chuyen 3 bit lien tiep thanh cac muc tin hieu tuong ung
for i=1:symbols
temp(i)=data(3*i-2)*4+data(3*i-1)*2+data(3*i);
end;
data1= temp;
subplot(3,1,2);
stairs([0:symbols-1],data1);
axis([0 symbols -0.2 M] );
%mo rong khong gian lay mau%
data1= reshape(ones(Fd,1)*data1, 1, Fd*symbols);
%ve do thi tin hieu sau khi dieu che%
switch cach_dieu_che
case {1}
subplot(3,1,3);
[y,t]=dmod(data1,Fc,Fd,Fs,'ask',M);
plot(t,y);
case {2}
subplot(3,1,3)
[y,t]=dmod(data1,Fc,Fd,Fs,'fsk',M);
plot(t,y);
case {3}
subplot(3,1,3)
[y,t]=dmod(data1,Fc,Fd,Fs,'psk',M);
plot(t,y);
case {4}
subplot(3,1,3)
[y,t]=dmod(data1,Fc,Fd,Fs,'qask',M);
plot(t,y);
otherwise
error('nhap sai gia tri yeu cau');
end
%%%%%%%%%xem ket qua duoi dang do thi%%%%%%%%%%%%
shg
II.2 Kết quả mô phỏng thu được
ASK
QASK
PSK
FSK
...
Yêu cầu Download

Tài liệu, ebook tham khảo khác

Học thêm

Music ♫

Copyright: Tài liệu đại học © DMCA.com Protection Status