Hệ thống viễn thông 1(semilar)

Điều chế ASK, PSK, FSK
cái này chỉ là điều chế ASK.

function ASK_FSK_PSK(msglen)
%msglen= so bit truyen

%%===============================================
n=msglen;
b=randint(1,n);
f1=1;f2=2;
t=0:1/30:1-1/30;
%ASK
sa1=sin(2*pi*f1*t);
E1=sum(sa1.^2);
sa1=sa1/sqrt(E1); %unit energy
sa0=0*sin(2*pi*f1*t);

%MODULATION
ask=[];
for i=1:n
if b(i)==1
ask=[ask sa1];

else
ask=[ask sa0];

end
end
figure(1)
subplot(411)
stairs(0:10,[b(1:10) b(10)],'linewidth',1.5)
axis([0 10 -0.5 1.5])
title('Message Bits');grid on
subplot(412)
tb=0:1/30:10-1/30;
plot(tb, ask(1:10*30),'b','linewidth',1.5)
title('ASK Modulation');grid on

%AWGN
for snr=0:20
askn=awgn(ask,snr);


%DETECTION
A=[];F=[];P=[];
for i=1:n
%ASK Detection
if sum(sa1.*askn(1+30*(i-1):30*i))>0.5
A=[A 1];
else
A=[A 0];
end

end

%BER
errA=0;errF=0; errP=0;
for i=1:n
if A(i)==b(i)
errA=errA;
else
errA=errA+1;
end

end
BER_A(snr+1)=errA/n;

end

figure(2)
subplot(411)
stairs(0:10,[b(1:10) b(10)],'linewidth',1.5)
axis([0 10 -0.5 1.5]);grid on
title('Received signal after AWGN Channel')
subplot(412)
tb=0:1/30:10-1/30;
plot(tb, askn(1:10*30),'b','linewidth',1.5)
title('Received ASK signal');grid on

figure(3)
semilogy(0:20,BER_A, 'b','linewidth',2)
title('BER Vs SNR')
grid on;
hold on

xlabel('Eo/No(dB)')
ylabel('BER')
hold off
legend('ASK');

PSK: nè
function ASK_FSK_PSK(msglen)

n=msglen;
b=randint(1,n);
f1=1;f2=2;
t=0:1/30:1-1/30;

%PSK
sp0=-sin(2*pi*f1*t)/sqrt(E1);
sp1=sin(2*pi*f1*t)/sqrt(E1);

%MODULATION
psk=[];
for i=1:n
if b(i)==1

psk=[psk sp1];

else

psk=[psk sp0];

end
end
figure(1)
subplot(411)
stairs(0:10,[b(1:10) b(10)],'linewidth',1.5)
axis([0 10 -0.5 1.5])
title('Bits gui di');grid on

subplot(414)
plot(tb, psk(1:10*30),'k','linewidth',1.5)
title('Điều chế PSK');grid on
xlabel('Time');ylabel('Amplitude')
%AWGN
for snr=0:20

pskn=awgn(psk,snr);


%DETECTION
A=[];F=[];P=[];
for i=1:n

%PSK Detection
if sum(sp1.*pskn(1+30*(i-1):30*i))>0
P=[P 1];
else
P=[P 0];
end
end

%BER
errA=0;errF=0; errP=0;
for i=1:n

if P(i)==b(i)
errP=errP;
else
errP=errP+1;
end
end

BER_P(snr+1)=errP/n;
end

figure(2)


subplot(414)
plot(tb, pskn(1:10*30),'k','linewidth',1.5)
title('Tín hiệu nhận');grid on
figure(3)
semilogy(0:20,BER_A, 'b','linewidth',2)
title('BER Vs SNR')
grid on;
hold on

semilogy(0:20,BER_P, 'k','linewidth',2)
xlabel('Eo/No(dB)')
ylabel('BER')
hold off
legend('PSK');



có gì everybody cho ý kiến nha!

Bạn chỉ đưa lên 1 đoạn code mà ko giải thích j thì thực tình thì mình ko hiểu j cả. Bạn có thể nói qua về phần mềm soạn thảo code đó, giải thích 1 số đoạn code cơ bản và đưa 1 sơ đồ nguyên lý ra cho mình xem với nhe ^^