Bài tập lớn - Kỹ thuật truyền dẫn

Bài tập lớn truyền dẫnBài tập môn học: Kỹ Thuật truyền dẫnBà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;...
Nội dung trích xuất từ tài liệu:
Bài tập lớn - Kỹ thuật truyền dẫnBài tập lớn truyền dẫn Bài tập môn học: Kỹ Thuật truyền dẫnBà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.mfunction y=nrz_l(x)for i=1:length(x)if x(i)==0temp(i) = -1;elsetemp(i) = +1;end;endy=temp; 2.Mã CMI: file cmi_1.mfunction y = cmi(x);polar = +1; %chon cuc tinh cua bit 1 dau tien la duongfor 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; endendy = temp; -1-Bài tập lớn truyền dẫn3.Mã HDB3: File hdb3.mfunction y = hdb3(x)next_polar = +1;b=0; %so bit B giua 2 bit Vi=1;while iBài tập lớn truyền dẫnendy=temp;4. Mã AMI: File ami.mfunction y=ami(x)b=+1;for i=1: length(x) if x(i)==1 temp(i)=b; b=-b; else temp(i)=0; endendy=temp;5. Mã Manchester: File manchester.mfunction 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; endendy=temp;6. Mã Polar RZ: File polar_RZfunction 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; -3-Bài tập lớn truyền dẫnendy=temp;7. Mã Unipolar RZ: File unipolar_RZ.mfunction 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; endendy=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.mclear; %Xoa cac bien hien cobits = 30; %so bit se khao satt = 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.mdisp(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); -4-Bài tập lớn truyền dẫndisp(30: ma Manchester);disp(40: ma CMI);linecode = input(chon dang ma se phan tich: );switch linecode %%Truong hop cac ma don bit day ducase {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 bitcase {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); -5-Bài tập lớn truyền dẫn TenDoThi=CMI; t=0.5:0.5:30;otherwise error(nhap gia tri sai);end%Do thi thoi gian cua du lieu gocsubplot(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 hoasubplot(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); -6-Bài tập lớn truyền dẫn 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 trinhMột số kết quả cho các trường hợp: -7-Bài tập lớn truyền dẫn -8-Bài tập lớn truyền dẫn -9-Bài tập lớn truyền dẫn - 10 -Bài tập lớn truyền dẫnBà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 ba ...