clc;
clear;
%%?比特序列產生
N?=?20;?????????????%設置比特序列的長度
x?=?randi（21N）-1;
%x?=?[1?1?0?1?0?0?1?1?1?1?0?1?0?0?1?1?0?1?1?1]
I_4PSK?=?zeros（1N/2）;
Q_4PSK?=?zeros（1N/2）;
I_QPSK?=?zeros（1N/2）;
Q_QPSK?=?zeros（1N/2）;

%%?4PSK編碼
for?m=1:N/2
????if（（x（2*m-1）==0）&&（x（2*m）==0））
????????I_4PSK（m）?=?1;
????????Q_4PSK（m）?=?0;
????end
????if（（x（2*m-1）==1）&&（x（2*m）==0））
????????I_4PSK（m）?=?0;
????????Q_4PSK（m）?=?1;
????end
????if（（x（2*m-1）==1）&&（x（2*m）==1））
????????I_4PSK（m）?=?-1;
????????Q_4PSK（m）?=?0;
????end
????if（（x（2*m-1）==0）&&（x（2*m）==1））
????????I_4PSK（m）?=?0;
????????Q_4PSK（m）?=?-1;
????end
end

%%?QPSK編碼
for?m=1:N/2
????if（（x（2*m-1）==0）&&（x（2*m）==0））
????????I_QPSK（m）?=?1;
????????Q_QPSK（m）?=?1;
????end
????if（（x（2*m-1）==1）&&（x（2*m）==0））
????????I_QPSK（m）?=?-1;
????????Q_QPSK（m）?=?1;
????end
????if（（x（2*m-1）==1）&&（x（2*m）==1））
????????I_QPSK（m）?=?-1;
????????Q_QPSK（m）?=?-1;
????end
????if（（x（2*m-1）==0）&&（x（2*m）==1））
????????I_QPSK（m）?=?1;
????????Q_QPSK（m）?=?-1;
????end
end


%%?低通濾波器參數設置
t?=?-10:0.1:10;
L?=?length（t）;
h?=?sinc（t）;
y_4PSK_I?=?zeros（N/2L）;
y_4PSK_Q?=?zeros（N/2L）;
Y_4PSK_I?=?zeros（1L）;
Y_4PSK_Q?=?zeros（1L）;

y_QPSK_I?=?zeros（N/2L）;
y_QPSK_Q?=?zeros（N/2L）;
Y_QPSK_I?=?zeros（1L）;
Y_QPSK_Q?=?zeros（1L）;

%%?4PSK濾波過程展示
for?i?=?1:1:N/2
????y_4PSK_I（i1:L）?=?I_4PSK（i）*sinc（t-i）;
????y_4PSK_Q（i1:L）?=?Q_4PSK（i）*sinc（t-i）;
end
for?i?=?1:1:N/2
????Y_4PSK_I（1:L）?=?Y_4PSK_I（1:L）?+?y_4PSK_I（i1:L）;
????Y_4PSK_Q（1:L）?=?Y_4PSK_Q（1:L）?+?y_4PSK_Q（i1:L）;
end

%%?4PSK濾波過程展示
for?i?=?1:1:N/2
????y_QPSK_I（i1:L）?=?I_QPSK（i）*sinc（t-i）;
????y_QPSK_Q（i1:L）?=?Q_QPSK（i）*sinc（t-i）;
end
for?i?=?1:1:N/2
????Y_QPSK_I（1:L）?=?Y_QPSK_I（1:L）?+?y_QPSK_I（i1:L）;
????Y_QPSK_Q（1:L）?=?Y_QPSK_Q（1:L）?+?y_QPSK_Q（i1:L）;
end


%%?平移疊加法?4PSK?I路?實部
figure（1）
subplot（211）
for?i?=?1:1:N/2
????plot（ty_4PSK_I（i1:L））;
????hold?on;
end
axis（[0?10?-1?1]）
title（‘4PSK實部?疊加前波形‘）;
subplot（212）
plot（tY_4PSK_I）
axis（[0?10?-2?2]）
title（‘4PSK實部?疊加后波形‘）;

%%?平移疊加法?4PSK?Q路?虛部
figure（2）
subplot（211）
for?i?=?1:1:N/2
????plot（ty_4PSK_Q（i1:L））;
????hold?on;
end
axis（[0?10?-1?1]）
title（‘4PSK虛部?疊加前波形‘）;
subplot（212）
plot（tY_4PSK_Q）
axis（[0?10?-2?2]）
title（‘4PSK虛部?疊加后波形‘）;

%%?平移疊加法?QPSK?I路?實部
figure（3）
subplot（211）
for?i?=?1:1:N/2
????plot（ty_QPSK_I（i1:L））;
????hold?on;
end
axis（[0?10?-1?1]）
title（‘QPSK實部?疊加前波形‘）;
subplot（212）
plot（tY_QPSK_I）
axis（[0?10?-2?2]）
title（‘QPSK實部?疊加后波形‘）;

%%?平移疊加法?QPSK?Q路?虛部
figure（4）
subplot（211）
for?i?=?1:1:N/2
????plot（ty_QPSK_Q（i1:L））;
????hold?on;
end
axis（[0?10?-1?1]）
title（‘QPSK虛部?疊加前波形‘）;
subplot（212）
plot（tY_QPSK_Q）
axis（[0?10?-2?2]）
title（‘QPSK虛部?疊加后波形‘）;

%%?卷積法?4PSK?實部
figure（5）
for?i?=?1:1:N/2
????a?=?zeros（1N/2）;
????a（i）?=?I_4PSK（i）
????a?=?upsample（a10）;
????y?=?conv（ah）;
????plot（y‘k‘）;
????hold?on;
end
I_4PSK?=?upsample（I_4PSK10）;
y?=?conv（I_4PSKh）;
plot（y‘r‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????141975??2019-01-26?11:03??4PSK和QPSK調制及成型濾波sinc\4PSK和QPSK調制及成型濾波sinc.docx
?????文件????????3975??2019-01-26?10:52??4PSK和QPSK調制及成型濾波sinc\Sim_NO2.m
?????目錄???????????0??2019-01-26?11:03??4PSK和QPSK調制及成型濾波sinc\