2023年新版通信原理实验报告.pdf

中 聋）夫 普 信 息 科 学 与 工 程 学 院 课 程 设 计 报 告 书 课 题：通 信 原 理 班 级：学 号：姓 名：指 导 教 师：郭 丽 梅-0-0 年 十 月目 录 令 目 录 一.-.-3，*令 一、课 程 设 计 目 的 一 一.一 一 一 一 一 Q 1令 二、课 程 设 计 题 目.一 一 一 一-.-2 J令 三、课 程 设 计 实 验 原 理-3令 四、建 立 模 型 描 述 一 一 一.一 一 6令 五、模 块 功 能 分 析 及 源 程 序 代 码 一-6令 六、倜 试 过 程 及 结 论-22令 七、调 试 分 析 一.-3 4令 八、基 于 MATLAB 的 M ASK、M F SK、M PSK调 制 拓 展(M=4;M=8)-3 4令 九、心 得 体 会 一-.一-A n 十、参 考 文 献-一.41摘 要：通 信 原 理 是 通 信 工 程 专 业 相 称 重 要 的 学 科，对 日 后 就 业 和 科 研 有 重 大 的 意 义,通 过 M A T L A B,我 们 可 以 清 楚 地 理 解 通 信 原 理 中 难 以 理 解 的 一 面，对 理 论 的 知 识 加 以 深 化。关 键 字：MATLAB 通 信 原 理 G U I 序 列 频 谱 相 位 相 干 非 相 干 一、课 程 设 计 题 目 1）、应 用 MATLAB编 制 信 号 生 成 程 序，并 对 信 号 进 行 HDB 3 编 码 和 译 码。2）、课 程 设 计 需 要 运 用 MAT L A B编 程 实 现 2ASK,2 FSK,2PSK,2DPSK调 制 解 调 过 程，并 且 输 出 其 源 码，调 制 后 码 元 以 及 解 调 后 码 元 的 波 形。二、课 程 设 计 目 的 通 信 原 理 是 一 门 以 算 法 为 核 心,理 论 和 实 践 性 较 强 的 学 科。是 电 子 信 息 工 程、通 信 工 程 专 业、电 子 信 息 科 学 与 技 术 专 业 的 一 门 重 要 的 专 业 技 术 基 础 课。通 信 原 理 课 程 设 计 是 在 学 习 完 通 信 原 理 的 相 关 理 论 后，进 行 的 综 合 性 训 练 课 程，其 目 的 是：1.使 学 生 进 一 步 巩 固 通 信 原 理 的 基 本 概 念、理 论、分 析 方 法 和 实 现 方 法；2.增 强 学 生 应 用 M a 11 ab语 言 编 写 数 字 信 号 解 决 的 应 用 程 序 及 分 析、解 决 实 际 问 题 的 能 力；3.全 面 复 习 课 程 所 学 理 论 知 识，巩 固 所 学 知 识 重 点 和 难 点，将 理 论 与 实 践 很 好 地 结 合 起 来；4.提 高 综 合 运 用 所 学 知 识 独 立 分 析 和 解 决 问 题 的 能 力；5.纯 熟 使 用 一 种 高 级 语 言 进 行 编 程 实 现。三、课 程 设 计 实 验 原 理 1.HDB3编 码 解 码 原 理 HDB3码：三 阶 高 密 度 双 极 性 码。HDB3码 与 二 进 制 序 列 的 关 系：(1)二 进 制 信 号 序 列 中 的“0”码 在 HDB3码 中 仍 编 为“0”码,二 进 制 信 号 中“1”码，在 HDB3码 中 应 交 替 地 成+1和-1码，但 序 列 中 出 现 四 个 连“0”码 时 应 按 特 殊 规 律 编 码(引 入 传 号 交 替 反 转 码 的“破 坏 点”V 码)；(2)二 进 制 序 列 中 四 个 连“0”按 以 下 规 则 编 码:信 码 中 出 现 四 个 连“0”码 时，要 将 这 四 个 连“0”码 用 0 0 0V 或 B00V取 代 节 来 代 替(B和 V 也 是“1”码，可 正、可 负)。这 两 个 取 代 节 选 取 原 则 是，使 任 意 两 个 相 邻 v 脉 冲 间 的 传 号 数 为 奇 数 时 选 用 000 V 取 代 节，偶 数 时 则 选 用 B0 0 V 取 代 节。2.二 进 制 数 字 调 制 技 术 原 理 数 字 信 号 的 传 输 方 式 分 为 基 带 传 输 和 带 通 传 输，在 实 际 应 用 中，大 多 数 信 道 具 有 带 通 特 性 而 不 能 直 接 传 输 基 带 信 号。为 了 使 数 字 信 号 在 带 通 信 道 中 传 输，必 须 使 用 数 字 基 带 信 号 对 载 波 进 行 调 制，以 使 信 号 与 信 道 的 特 性 相 匹 配。这 种 用 数 字 基 带 信 号 控 制 载 波,把 数 字 基 带 信 号 变 换 为 数 字 带 通 信 号 的 过 程 称 为 数 字 调 制。通 常 使 用 键 控 法 来 实 现 数 字 调 制，比 如 对 载 波 的 振 幅、频 率 和 相 位 进 行 键 控。(1)2ASK信 号 的 产 生 方 法 通 常 有 两 种：模 拟 调 制 和 键 控 法。解 调 有 相 干解 调 和 非 相 干 解 调。P=1时 f(t)=A co sw t;p=0时 f(t)=0;其 功 率 谱 密 度 是 基 带 信 号 功 率 谱 的 线 性 搬 移。解 调 原 理 图 定 时 脉 冲 图 1 2ASK非 相 干 解 调 接 受 系 统 图 2 2 ASK相 干 解 调 接 受 系 统(2)一 个 2FSK信 号 可 以 当 作 是 两 个 不 同 载 波 的 2ASK信 号 的 叠 加。其 解 调 和 解 调 方 法 和 A S K差 不 多。2FSK信 号 的 频 谱 可 以 当 作 是 f 1和 f 2 的 两 个 2A SK频 谱 的 组 合。解 调 原 理 图图 3 2 F S K 相 干 解 调 接 受 系 统 图 4 2 FSK非 相 干 解 调 接 受 系 统(3)2PSK以 载 波 的 相 位 变 化 作 为 参 考 基 准 的，当 基 带 信 号 为 0 时 相 位 相 对 于 初 始 相 位 为 0,当 基 带 信 号 为 1时 相 对 于 初 始 相 位 为 180。(4)2DPSK调 制 原 理 方 框 图 如 下 图：t图 5 间 接 法 信 号 调 制 器 原 理 方 框 图 2DPSK信 号 的 解 调，重 要 有 两 种 方 法，即 相 位 比 较 法 和 相 干 解 调 法。相 干 解 调 法 原 理 方 框 图 如 下 图:图 6 相 干 解 调 法 原 理 方 框 图 四、建 立 模 型 描 述 用 MA TLA B实 现 二 进 制 振 幅 键 控（2ASK）的 调 制 和 解 调 在 这 里 用 MATLAB函 数 来 对 这 个 二 进 制 振 幅 键 控（2 A S K）来 实 现 调 制 与 解 调 的 仿 真。二 进 制 振 幅 键 信 号 可 以 表 达 完 毕 一 个 单 极 性 矩 形 脉 冲 序 列 与 一 个 正 弦 型 载 波 的 乘 积。通 常 它 的 调 制 方 法 有 两 种，即 模 拟 幅 度 调 制 方 法 和 键 控 方 法，在 MATLAB里 我 们 采 用 模 拟 幅 度 调 制 的 方 法,解 调 采 用 相 干 解 调（包 络 检 波 法）的 方 式。我 们 用 SOURCE函 数 来 产 生 一 个 原 始 二 进 制 基 带 信 号，即 一 个 单 矩 形 脉 冲 序 列。以 a s k M。d u 函 数 来 进 行 模 拟 幅 度 调 制,得 到 一 个 已 调 2ASK信 号,并 用 此 函 数 进 行 此 2ASK信 号 的 频 谱 分 析。然 后 用 guss i an函 数 加 入 加 性 高 斯 白 噪 声，再 用 demoAS K函 数 进 行 想 干 解 调 并 分 别 输 出 各 点 的 输 出 波 形，最 后 通 过 抽 样 判 决 后 得 出 输 出 波 形。同 时 我 们 用 C h e e kR atePe函 数 来 得 出 误 码 率，最 后 运 营 主 函 数 ASK_ma in可 以 看 到 各 种 波 形。用 M ATLAB函 数 来 对 2 F S K、2PSK、2 D P S K 实 现 调 制 与 解 调 的 仿 真 与 2ASK的 调 制 与 解 调 过 程 的 仿 真 类 似。五、模 块 功 能 分 析 及 源 程 序 代 码 5.1 HDB3码 的 编 码 和 译 码 gl o ba 1 xx n=x;y n=xn；n u m=0;%计 数 器 初 始 化 for k=1:1 e n gth(xn)if xn(k)=l%1的 计 数 器 num=num+1;if num/2=fix(n um/2)/奇 数 个 一 时 输 出 一 1,进 行 记 性 交 替 yn(k)=l;e Iseyn(k)=-l；endendend%HDB3编 码 num=0;y h=yn;sign=0；v=zeros(1,1 ength(y n);B=ze r o s(1,length(y n)；%B 脉 冲 位 置 记 录 for k=l:len g t h(yn)if y n(k)=0num=num+1；%连 0个 数 计 数 if num=4%假 如 4连 0n u m=O；%计 数 器 清 零 yh(k)=l*y h(k-4)*最 后 一 个 0改 变 与 前 一 个 非 零 符 号 同 极 性 v(k)=y h(k);if yh(k)=s ign%假 如 当 前 V与 前 一 个 极 性 相 同 y h(k)=-l*yh(k);%则 当 前 V 符 号 极 性 反 转，yh(k-3)=yh(k);%添 加 B脉 冲 B(k-3)=yh(k)；%8脉 冲 位 置 v(k)=y h(k);%V脉 冲 位 置 y h(k+1:leng t h(yn)=1*yh(k+1:len g th(y n);endsign二 y h(k)；%计 算 前 一 个 V符 号 的 极 性 en de Isen uni=O;%当 前 输 入 为 0，则 计 数 器 清 零 ende ndre=xn，,yn，,yh),v，,B，;%HDB3解 码 i np u t=yh;d e c od e=in p ut;sign=0;为 极 性 标 志 初 始 化 for k=1:le n gth(yh)if input(k)=0if sign=yh(k)%假 如 当 前 码 与 前 一 个 非 零 码 的 极 性 相 同 d e code(k-3:k)=0 0 0 0;endsign=i npu t(k);%极 性 标 志 endenddec o de=a b s(d ecode)；su b plot(3,1,1);stairs(xn)；ax i s(l 1 e n gth(xn)-2 2);gr i d;ylabel(xn);title。HDB3码 的 编 码 前 图 形);s ubplo t(3,1,2);s t airs(y h)；axis(1 length(xn)-2 2);grid;yla b el(HDB 3 codes);title(,H D B 3 码 的 编 码 后(解 码 前)图 形)；subpl o t(3,1,3)；sta i rs(dec ode);axis(1 leng t h(xn)-2 2);gr i d;y labe 1(decoded x n);ti t leC HDB3码 的 解 码 后 图 形);c o del=y h;d ecode 1=d e c o de；5.2 ASK的 调 制 与 解 调 5.2.1 source 函 数 fu n ct i on s e nd S ign a 1=source(n,N)s e ndS i g nal=r and i n t(1,n)bi t=;for i=l:length(s e ndSi g nal)if s e ndSign a 1(i)=0b i t 1=zer o s(1,N);e 1 seb i tl=ones(1,N);e n dbit=b i t,b i t 1;endfigu r e(1)plot(l:len g th(b i t),bit),title(发 送 端 二 进 制 波 形),grid on；a x i s(0,N*lengt h(send S ignal),-2,2);end5.2.2 a sk M o d u函 数 源 程 序 fun c t ion t r an s mi t ted Sig n al=askModu(sign al,bitRate,f c,N)%signal=l 0 1 0 1 0 0 1；%bi t Rate=l 0 00 0 00;%fc=1000 0 0 0;%N=3 2；t=linsp a c e(0,1/bit R a te,N);c=s i n(2*pi*t*f c);t ran sm i t tec!Signal=;for i=l:length(sig nal)t r a nsmi t t e d S i g nal=t r ansmi t ted S i g n a 1,s i gna 1(i)*c;endfigu r e(2)plo t(1:leng t h(t ran s mi t tedSigna 1),tran s mitte d S i g n a 1);tit 1 e(A S K调 制 波 形);grid on;figur e(3)m=0:le n gth(tran s mitt e d Signa 1)-1;F=f ft(tran s mitt e dS i g n al)；pl o t(m,a bs(re a 1(F),ti t le(ASK仿 真 频 谱 分 析);gr i d o n;%f i gure(4)%plot(m,imag(F)；title。ASK fre q uency-domai n an a ly s is imag,)；%grid o n;End5.2.3 gus s i a n 函 数 源 程 序 f u nctio n s i g n a l=gussian(t r ansm i t tedSignal,noise)signal=sqrt(2)*t r an s m ittedSignal;sig n a l=aw g n(sign a 1,no i se);f i g u r e(5)plot(1：le n gth(s ignal),sig n a l);ti t 1 e(z包 含 噪 声 的 波 形)，gr i d on;en d5.2.4 CheckRatepe 函 数 源 程 序 func t i o n PeW r ong=CheckR a t e Pe(signalL signa 1 2,s)right s=0;wron g s=0;for ki=l:s 2if(s ig n all(ki)=signal2(ki)rights=rig h t s+1；elsew r o n g s=w r ongs+1;e ndendPeWron g=wr ongs/(wrongs+rig hts);End5.2.5 de mo ASK函 数 源 程 序 f unc t ion bits t re a m=demo ASK(receive d Si gnal,bi tRa t e,fc,n,N)load numsigna 1 1=r eceiv e d Sig n al;s ign a 12=abs(signal 1);%整 流 si g nal3=filte r(numl,1,s ign a 12);%LPF,包 络 检 波”IN=f i x(1 e ngth(n uml)/2);%延 迟 时 间 b i t s tr e am=口;LL=fc/bit Rate*N;i=IN+LL/2;whi 1 e(i=0.5 1;i=i+LL;endfigure(6)subplot(3,1,1);plot(l:leng th(sign a 11),s i gnal 1);ti t le(接 受 端 波 形(包 含 噪 声)；g r id on;s u bplot(3,1,2);p 1 o t(1：len g th(sign a 12),s i g nal2)；t i t le(，整 流 之 后 的 波 形)；g r i d on;subplot(3,1,3);p lot(l:len gth(sign a 1 3),s i gna 1 3);tit 1 e(LPF 滤 波 后 的 包 络 波 形)；g ridon;bit=;fo r i=l:1 e n g th(b i t s tream)if b i ts t r e am(i)=0b it 1=zeros(l,N);elseb i t 1=one s(1,N);e n db i t=bit,bit 1;en df igu r e(7)plot(bit),ti t leC 接 受 端 二 进 制 波 形),g r id on;axis(0,N*1 e n g th(bi t stream),-2.5,2.5);end5.2.6 ASK_m a i n 函 数 源 程 序 clos e al 1c 1 e a r al 1%(ti=0；fpeask=；st a rtn=-6;endn=18；for ti=star t n:en d nn=100 0;%n=16;fc=1 0 00000;%fc=bitR a te fc/b i t Ra t e为 每 个 包 含 sin周 期 个 数 b i tRate=1000000；N=5 0；%n o ise=ti;n o i se=l 0;s i g n a 1=source(n,N);transmit t e d S i g nal=as k Mo d u(sign a 1,b itRa t e,f c,N);signa ll=g uss ian(tran s mi t t e dS i g n al,n oise);c o n f i gue S i g n al=demoASK(s ig n all,b itRa t e,fc,n,N);%(conf i gue S i gn a 1；P=C h e ckRa t eP e(s i g n al,co n f i gueS i g n al,n)f p e as k=fpea s k,P；endfigur e(8);sem i logy(startn：l e n g t h(fpeas k)+st a r tn-l,f peask)；grid on;title(ASK误 码 率);xla b el C r/d B);ylabe 1 C PeA S K);sa v e PeRate.mat fpeask%5.3 FSK的 调 制 与 解 调 5.3.1 s o u r ce 函 数 f u n c ti o n sendSi g nal=s o u r c e(n,N)s e n dSignal=rand i n t(1,n)b i t=;for i=l:len g t h(send S i g na 1)i f s en d Si g nal(i)=0bit 1=ze r o s(1,N);e l s eb i t 1=one s(1,N)；endbit=bit,bitl；e ndf igure(1)pl o t(bit),ti t le C 发 送 端 二 进 制 波 形),g ri da x i s(0,N*length(sen d Si g n a 1),-2.5,2.5);en d5.3.2 F s k M o d u函 数 源 程 序 f un c t io n transmitt e dSignal=fskM o du(sig n al,b i t Rt=1 i nspace(0,1/bitRa t e,N);cl=sin(2*p i*t*fl);on；a t e,fl,f2,N)c 2=s in(2*pi*t*f2);tr a n sm i t t edSignal=;fo r i=l:le n gth(s i gnal)if sig n al(i)=ltransmi t t e dSi g na 1=transmi t ted Sig nal,c 1;el s et r a n smitted S i g n a l=t r a n s m i ttedSign a 1,c2;endendfigu r e(2)plot(1:1 e n gth(t ran s mi t t e dSignal),t ran smit t ed S i g nal);t itleC FSK调 制 波 形)；grid on;figu r e(3)m=0:1 ength(t r an s mitted S i gn a 1)-1;F=ff t(transmi t tedS i g nal);P 1 ot(m,ab s(re a 1(F),t it 1 e(FSK仿 真 频 谱 分 析)；gri d o n;end5.3.3 g u s s ia n函 数 源 程 序 fun c t i on sig n al=g u ssian(t ransmitted S i g nal,n o i s e)signal=s q r t(2)*t ransmitt e dSi g na 1;s ignal=aw g n(signa 1,no i s e);f i gure(4)P 1 ot(l:leng t h(signal),s i gnal),t i t le(z 包 含 噪 声 的 波 形)；grid o n;end5.3.4 C h ec k Rate p e 函 数 源 程 序 function P eWr o ng=Check R a tePe(signal 1,signa 1 2,s)r i g h ts=0；wr o ngs=0；for k i=1:s-2i f(s i gn a 1 1(ki)=s i gnal2(ki)r i g h ts=rig h t s+1;elsewr o ngs=w r on g s+1；endendPeWron g=wr o n g s/(w r ongs+r i g h t s);en d5.3.5 demoFSK函 数 源 程 序 function b i tstr e am=demoFSK(re c e i ve d Signa 1,bi t Ra t e,f 1,f2,N)1 o ad FSK n umsign a 11=rec e iv e dSig n al;sign a 1 2=f i Iter(g a o ton g,1,si g n all);为 通 过 H P F,得 到 高 频 分 量 s ign a 1 3=a b s(signa12);/整 流 signal3=f ilter(lowpass,1,s ig n al 3);%通 过 LPF,形 成 包 络 bit s tream=;I N 1=f i x(len g t h(lowp a s s)/2)+f i x(leng t h(g a o ton g)/2);%延 迟 时 间 bits t r e a ml=;LL=N;%每 个 bi l的 抽 样 点 数 i=IN 1+LL/2；w h i le(i=0.5;i=i+LL;en db i t s treamlfigu r e(5)sub plot(3,1,1);pl o t(1：1 e ng t h(signall),signall)接 受 端 波 形(包 含 噪 声)；g r id on;s u bpl o t(3,1,2)；pl o t(l:leng t h(s ign a 12),signa 1 2);title。通 过 H PF得 到 的 高 频 分 量 波 形);g rid on；s u bplot(3,1,3);P lot(1:leng t h(s i goal 3),signal 3)；titl e C 通 过 LPF后 的 包 络 波 形，);g ri d on;sign a 14=fi 1 ter(da i t o n g,1,s i gnall)；先 通 过 BP F,得 到 低 频 分 量 si g na 1 5=a b s(sig n a 14);%整 流 s i gnal5=fi 1 t er(1 owpa s s,1,s ignal5);%通 过 LPF,形 成 包 络 IN2=fix(1 en g t h(1 owp a ss)/2)+f i x(len g t h(da i tong)/2);%延 迟 时 间 b itstr e am2=口；LL=N;%每 个 b i t的 抽 样 点 数 i=IN2+LL/2;w h i l e(i=0.5;i=i+LL;e n dbitstream2fig u re(6)s ubplot(3,1,1);p 1 o t(l:len g th(s i gnall),s i gnall);ti t 1 e(接 受 端 波 形(包 含 噪 声)；gridon;s u b p lot(3,1,2);p lot(1：lengt h(si g nal4),s i gna14)；t i t l e。通 过 B PF得 的 低 频 分 量 波 形)；grid o n;sub p 1 ot(3,1,3);p 1 ot(l:len g th(s ignal5),s i g nal5);titl e(通 过 LPF后 的 包 络 波 形);gri d on；for i=1：m i n(leng t h(bitst r e am 1),le n gth(bit s t ream 2)%判 决 i f(b i tstreaml(i)bits t r e am2(i)bi t stream(i)=l;elsebitstream(i)=0;e n de n db i tst reamb i t二 口；外 接 受 端 波 形 fo r i=1：len g t h(bi t stre a m)if b i tstream(i)=0bit 1=z e ros(1,N);elsebitl=one s(1,N)；en dbit=bit,b i t 1;endfi g u r e(7)p 1 ot(bi t),title(接 受 端 波 形(解 调 后 波 形)，g r i d o n;axis(0,N*leng t h(bitstream),-2.5,2.5);en d5.3.6 FSK_m a i n 函 数 源 程 序 close a 1 1c 1 e ar all%(ti=0；fpe f sk=;startn=-6;en d n=1 8；f o r t i=s t ar t n:endnn=1 000;%n=16；f 1=180000 0 0;f 2=6 0 00000;bi t Ra t e=l 0 00000；N=5 0;%n o i s e=t i;n o i s e=10;s i gnal=F SKs o u r ce(n,N);t r a nsmit t e dSignal=fskModu(s i gna 1,b i tRate,f 1,f 2,N);sig n a 1 1=F SKgussian(transm it t edSig n a 1,noise);configueSigna 1 二 d e moFSK(s i gnal 1,bit Rat e,f 1,f 2,N);%(con fig ueSigna 1;P=CheckRa t eP e(sig n al,c onfigueS i g n al,n)fpe fsk=fpe fsk,P;endfi gu r e(8);s em i log y(star t n:1 en g th(fpef s k)+s ta r tn-1,f p e fsk);grid o n;title Bit Er r o r Rat e Of FSK);x 1 ab e 1(r/dB1);y 1 ab e IC PeFSK*);1 o ad PeR a tes ave P e R a t e.mat fp e fs k fp e a s k%5.4 PSK的 调 制 与 解 调 5.4.1 source 函 数fun c t ion s endSign a 1=s ource(n,N)sendSignal=randint(1,n)bit=;f o r i=l:1 en g th(se n dSigna 1)i f sen d Sign a 1(i)=0b it 1=zer o s(1,N);e 1 sebit 1=one s(1,N);endbi t=bit,bitl;endfigure(1)P 1 o t(bi t),t i tie(发 送 端 二 进 制 波 形)，g ri d on;ax i s(0,N*1 e ngth(s e ndSi g nal),-2.5,2.5);e n d5.4.2 PSK Mod u 函 数 源 程 序 fun c t io n tra n s mit t edS i gn a 1=b p s kMo d u(s ignal,b i t R a te,fc,N)t=linspac e(0,1/bi t R a t e,N);c l=si n(2*pi*t*fc)；c 2=s in(2*pi*t*f c+pi);t ran s mit t edSignal=1;for i=l:1 ength(s i gnal)if sig n a 1(i)=1t ransm i t t edSignal=t ransmit t e dSignal,cl;el s et r a n smittedSign a l=tr a nsm i tte d Signal,c 2；ende n dfigure(2)plot(1:1 e ng t h(tra n s m i t t edSignal),transm i tt e d S ignal)；t itle(BPSK)；g ri d o n;figur e(3)m=0:length(trans m i tte d Sig n al)-1;F=f f t(t r ansmittedS i gna 1);plo t(m,ab s(real(F),title。2PSK仿 真 频 谱 分 析)；g r i d on;e nd5.4.3 gu s s i a n 函 数 源 程 序 function si g n a 1=BP S Kgussi an(tra n smitted S i g n al,noi s e)sig n a 1=s qrt(2)*t ransmitte d S i gnal;s i gnal=a w g n(s ig n al,noise);fig u re(4)plot(1:length(s i g n a l),s i g n al),g r i d on;title(包 含 噪 声 的 波 形)end5.4.4 Ch e ckR a t e pe 函 数 源 程 序 fu n ct i on PeW r ong=B PSKC h eckRat e Pe(sign a 11,si g nal2,s)r ig h ts=0;w r o n g s=0;for k i=1：s-2if(sign a 11(ki)二 二 sign a 12(ki)rig h ts=r ights+1;elsewrongs=wro n gs+1;e n dendPeWro n g=w r o n gs/(wro n gs+r i g h t s);End5.4.5 dem oPSK函 数 源 程 序 functi o n b i tst r eam=demoBPSK(rec e i v ed S i gna 1,b i t Ra t e,f c,n,N)load n u m%读 取 num中 存 储 的 低 通 滤 波 用 的 数 据 s i g n a 1 1=r eceiv e dSign a 1;t=linspac e(0,1/bitR a te,N);c=si n(2*p i*t*fc)；signal二 口；for i=1:nsignal=signal,c；endsi g na12=sig n all.*si g nal;为 乘 同 频 同 相 s insigna 1 3=f i Iter(numl,1,signal 2)；%L P F,包 络 检 波 IN=f i x(1 ength(numl)/2)；%延 迟 时 间 bi t s tre a m=1；L L=f c/bitRate*N；i=IN+L L/2;while(i=0;i=i+LL;e n df i gu r e(5)subplot(3,1,1);p 1 ot(l:le n gth(signal 1),s i gna 1 1);title(接 受 端 波 形(包 含 噪 声);grid o n；subplot(3,1,2);plot(1：length(s i gnal2),signa 1 2)；丫 12 6。1(相 干 解 调)；t it 1 e(z乘 同 频 同 相 s i n之 后 的 波 形);grid o n;s u b p lot(3,1,3)；p 1 ot(1：1 e n gth(signal3),si g n al3);title。LPF滤 波 后 的 包 络 波 形)；grid on;b it=；for i=1:le n gth(bits t r e a m)i f bitstream(i)=Ob i t 1=zer o s(1,N);els eb itl=ones(1,N);e ndb it=b i t,bit 1;e ndfigure(6)p l o t(bi t)；t i tl e c 接 受 端 二 进 波 形(解 调 后 波 形)；gr i d on;axis(0,N*1 e ng t h(bit s t re a m),-2.5,2.5);e nd5.4.6 P S K _ m a i n 函 数 源 程 序 c 1 os e a 11c 1 e ar all%(ti=0;f pe p sk=口；star t n=-30;e ndn=6;for t i=start n:end nn=l 0 00；%n=1 6;fc=100 0 0 0 0；bit Rat e=10 0 0000;N=5 0；%n oi s e=ti;noi s e=1 0;sign a l=BPSKso u r ce(n,N);tra n smit t e dSig n al=bps k M o d u(si g n a 1,bitR ate,fc,N);s i gna 1 1=BPSK g us s ian(t ran s mittedSi g n al,noise);c onfigueSignal=dem o BPSK(s i g nail,b it Rate,f c,n,N)；%(configue S ignal;P=C h eckRatePe(s i g na.l,conf i gu e Sig n a 1,n)f p e p sk=fp e p sk,P；endfigu re(8);s e mi 1 o g y(sta r tn:1 e n g t h(f p eps k)