%%??bpsk信號的產生
T?=?1;%采樣時間
f?=?1e5;%采樣率??fs=100000
t?=?0:1/f:（T-1/f）;?%?采樣點
n?=?length（t）;?%?采樣點數=100000

binary_code?=?‘1001100110‘;
len_code?=?length（binary_code）;??%得出二進制數的長度

%?一個碼元對應的脈沖長度???采樣點的個數和信號的長度一致
f_bpsk?=?1000;?????%?載波頻率的頻率為1000Hz，周期1ms
T_bpsk?=?1/f_bpsk;?%?載波頻率的周期
A_bpsk?=?10;???????%?載波頻率的幅值

t_mayuan?=?0:1/f:（T_bpsk-1/f）;???????????????%?一個碼元所持續的時間內的采樣點時刻
y_mayuan?=?A_bpsk*sin（2*pi*f_bpsk*t_mayuan）;?%?一個碼元的波形，用于表示0
t_bpsk?=?0:1/f:（T_bpsk*len_code-1/f）;????????%bpsk所持續的時間內的采樣點時刻

y_bpsk?=?[];???????%通過下面循環將其拼接在一起
for?i=1:len_code
????tmp?=?str2num（binary_code（i））;??%str2num?表示?Convert?string?to?number
????if?tmp?==?0
????????y_bpsk?=?[y_bpsky_mayuan];
????end
???if?tmp?==?1
????????y_bpsk?=?[y_bpsk-y_mayuan];
????end
end
%?subplot（221）;plot（t_bpsky_bpsk‘green‘）;??%繪制對應的bpsk信號
%?xlabel（‘t/s‘）;title（‘‘1001100110’的BPSK編碼信號‘）;
%?axis（[0?1e-2?-2*A_bpsk?2*A_bpsk]）;

%%?bpsk信號疊加噪聲
n=randn（size（y_bpsk））;??%產生隨機噪聲
%?subplot（331）;plot（size（y_bpsk）randn）;%繪制噪聲信號
y_bpsk=y_bpsk+n;%受燥后的信號
subplot（221）;plot（t_bpsky_bpsk‘green‘）;??%繪制對應的bpsk信號
xlabel（‘t/s‘）;title（‘‘1001100110’的BPSK編碼信號（受rand）‘）;
axis（[0?1e-2?-2*A_bpsk?2*A_bpsk]）;

%%?對BPSK信號做FFT變換
Y=fft（y_bpsk）;%對信號做傅里葉變換

%%??幅度譜
%信號采集的頻率?Fs?=?1000?Hz，采集的數據點個數?N?=?1000????信號長度N=1000和fs=1000采樣頻率要一一對應
L=length（t_bpsk）;???%disp（L）;??%L=1000
P1=2*abs（Y）/L;%幅度譜???要得到真實的振幅值的大小，只要將得到的變換后結果乘以2除以L即可。
P=P1（1:L）;%和下一行語句采樣點的個數對應
%?disp（P）;
f=（0:L-1）;
subplot（222）;plot（fP）;
xlabel（‘f/Hz‘）;title（‘BPSK信號的幅度譜‘）;

%%?相位圖
%?L=1000;

Pyy=[1:L];
for?i=1:L
????Pyy（i）=phase（Y（i））;?%計算相位
????Pyy（i）=Pyy（i）*180/pi;?%換算為角度
end;
F=（[1:L]-1）*f_bpsk/L;?%換算成實際的頻率值
subplot（223）;plot（F（1:L）Pyy（1:L））;???%顯示相位圖
xlabel（‘f/Hz‘）;title（‘相位-頻率曲線圖‘）;

%%??載波頻率
Fs?=?2500;%信號采集的頻率，Fs>2*f_bpsk采樣定理的要求
T=1/Fs;
L?=?2500;%采集的數據點個數????L和Fs必須相等??
t=（0:L-1）*T;
f_bpsk?=?1000;?%?載波頻率1000Hz，周期1ms
y_mayuan?=?10*sin（2*pi*f_bpsk*t）;
Y=fft（y_mayuan）;%對信號做傅里葉變換
P1=2*abs（Y）/L;%幅度譜
P2=P1（1:L/2）;%和下一行語句采樣點的個數對應
f=（0:L/2-1）;
subplot（224）;plot（fP2）;
xlabel（‘f/Hz‘）;title（‘高頻載波的幅度譜‘）;
fprintf（‘觀察點1\n‘）;

%%?求3dB帶寬?
%找到頻譜的最大值
disp（P）;
max=0;??%將max的值初始化
N=length（P）;??%??求出P的長度為1000
%?fprintf（‘觀察點2\n‘）;
%P（11）??檢測P（1，i）的值
for?i=1:N
????if?P（1i）>=max
????????max=P（1i）;
????else?
????????max?=?max;
????end
end
fprintf（‘幅頻的最大值為?%g\n‘max）;
%求出最大值對應的頻率
global?fh;
for?i=1:N
????if?P（1i）==max
????????fh=i;
????end
end
fprintf（‘最大值對應的頻率為?%gHz\n‘fh）;

%求出下降3dB對應的幅度大小
global?m;
for?i=1:N
????if?P（1i）>=0.7*max?&??P（1i）<=0.8*max
????????m=P（1i）;
????end
end
fprintf（‘幅頻的0.707*最大值為?%g\n‘m）;
%求出下降3dB對應的幅度所對應的頻率
global?fl;
for?i=1:N
????if?P（1i）==m
????????fl=i;
????end
end
fprintf（‘下降3dB對應的幅度所對應的頻率為?%gHz\n‘fl）;

%求出3dB帶寬
global?B;
B=fh-fl;
fprintf（‘3dB帶寬為?%gHz\n‘B）;

%%?碼元速率的計算
Rb=B/2;
fprintf（‘碼元速率為?%g波特\n‘Rb）;

%%?信號的功率?????計算公式為?pow=sum（xn.^2）/（length（x）/Fs），xn為信號序列，Fs為采樣率
%?disp（y_bpsk）;??%顯示y_bpsk序列的各點的值
%disp（leng

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2020-05-25?10:26??PSK調制\
?????目錄???????????0??2020-05-26?18:07??PSK調制\BPSK\
?????文件?????????551??2020-05-25?10:05??PSK調制\BPSK\DTDT_bianhuan.m
?????文件????????3761??2020-05-25?10:19??PSK調制\BPSK\bpsk_rand.m
?????文件????????5181??2020-05-25?09:42??PSK調制\BPSK\bpsk_tiaozhi_jietiao.m
?????文件????????3930??2020-05-25?10:16??PSK調制\BPSK\compute_Rb.m
?????文件????????3129??2020-05-25?10:22??PSK調制\BPSK\compute_Rb_way2.m
?????文件????????2994??2020-05-26?17:17??PSK調制\BPSK\erjingzhi_tiaozhi.m
?????文件?????????861??2020-05-25?09:28??PSK調制\BPSK\fft_FFT.m
?????文件?????????644??2020-05-14?09:10??PSK調制\BPSK\generate_bpsk.m
?????文件????????2955??2020-05-15?09:39??PSK調制\BPSK\sample_5.m
?????文件????????1028??2020-05-25?10:00??PSK調制\BPSK\test_2.m
?????文件????????3697??2020-05-18?13:27??PSK調制\BPSK\源代碼1.txt
?????文件????????3126??2020-05-18?08:33??PSK調制\BPSK\碼元速率的估計.txt
?????目錄???????????0??2020-05-27?16:12??PSK調制\MQAM\
?????文件????????1466??2020-05-27?16:10??PSK調制\MQAM\MQAM_error.m
?????文件????????3359??2020-05-27?15:56??PSK調制\MQAM\QAM16.m
?????文件????????1775??2020-05-27?16:00??PSK調制\MQAM\QAM_16_tiaozhi_jietiao.m
?????文件?????????230??2020-05-27?16:04??PSK調制\MQAM\sample_1.m
?????目錄???????????0??2020-06-09?10:10??PSK調制\MSK\
?????文件???????14181??2020-06-07?07:20??PSK調制\MSK\1.docx
?????文件????????3490??2020-06-10?22:33??PSK調制\MSK\MSK_MSK_1.m
?????文件????????2458??2020-06-09?10:06??PSK調制\MSK\MSK_MSK_2.m
?????文件????????2548??2020-06-01?12:49??PSK調制\MSK\MSK_signal.m
?????文件????????2930??2020-06-09?09:43??PSK調制\MSK\MSK_signal_1.m
?????文件????????3156??2020-06-05?11:03??PSK調制\MSK\MSK_signal_2.m
?????文件????????9391??2020-06-02?16:13??PSK調制\MSK\Untitled2.m
?????文件????????2861??2020-05-30?15:10??PSK調制\MSK\msk.m
?????文件????????1679??2020-05-31?13:23??PSK調制\MSK\msk_pingpu.m
?????文件????????2434??2020-06-01?10:20??PSK調制\MSK\msk_pingpu_1.m
?????文件?????????111??2020-06-01?19:47??PSK調制\MSK\test.m
............此處省略17個文件信息