繼續(xù)講解！如果不調(diào)用通信工具箱，那么之前的程序如何改寫才能完成相同的功能呢？先來看看MATLAB里面強(qiáng)大的工具箱！使用方法見上一篇文章中的程序。

我們重點關(guān)注調(diào)制和解調(diào)的相關(guān)函數(shù)！

有這么多！

再看看Modulation方面有哪些函數(shù)！

太多知識了，

簡直就是應(yīng)有盡有！

真的是學(xué)習(xí)通信基礎(chǔ)知識的好幫手！

我在通信仿真的系列文章里面會詳細(xì)介紹！

看函數(shù)的語言就知道如何使用了！

學(xué)習(xí)仿真，help必須經(jīng)常查閱！

%%%% File: BPSK_theory_sim2.m %%%%

%%% date: 2022-09-10 author: 算法工匠 %%%

%%% 程序說明

% BPSK調(diào)制解調(diào)理論仿真。

% 可以幫助大家消化通信模型！

% 通過函數(shù)來理解通信模型，這也是最簡單的學(xué)習(xí)路徑！

% 不使用工具箱！

%%% 仿真環(huán)境

% 版本：2021a

clear;close all;format long;

%%****************** 程序主體 ***********************%%

EsNo = 10;

% 蒙特卡洛仿真

% 請查找資料，了解什么是蒙特卡洛仿真？

err_num = 0;

for counter = 1:100

txData = randi([0 1],500,1); % Generate data 信源

% 調(diào)制方式 PSK

modSig = 2*txData - 1 ; % Modulate 發(fā)射機(jī)

rxSig = awgn(modSig,EsNo); % Pass through AWGN 傳輸

rxData = (sign(rxSig)+1)/2; % Demodulate 接收機(jī)

err_num = sum(abs(rxData-txData)) + err_num;

biterr_rate = err_num/(counter*length(txData));

end

fprintf('Error rate = %f Number of errors = %d ', biterr_rate, err_num)

% 請掌握這個函數(shù)！

figure

plot(rxSig-modSig,'*-')

title('噪聲')

審核編輯：劉清