untitled2.m

来自「matlab的调制解调」· M 代码 · 共 40 行

% 通过这个程序可以学到很多东西：对于qpsk调制来说，收发双方的
% 相位同步是很重要的；

% 试验一下这个程序；
% close all; clear all;  clc;
N=1000;
x=Jin_Zhi_2(N);

fc=512*1; fs=2048*1;   symbol_t=0.001;   theta=pi/4;
% 注意这里的时间是一个四进制码元下的持续时间；
n=ceil(symbol_t*fs);
t=(0:n-1)/fs;
[x_mod]=qpsk_mod_ZiBian(x,fc,fs,symbol_t,theta,0);

power_bit=ones(1,n)*ones(1,n)'/n;

ber=[];
snr_plot=-20:20;

for snr=-20:20
% 高斯白噪声功率；
sigma_2=power_bit/(10^(snr/10));
sigma=sqrt(sigma_2);
randn('state',0);
x_noise=sigma*randn(size(x_mod));
x_mod_noise=x_mod+x_noise;
[y_demod_bit]=...
    qpsk_demod_ZiBian(x_mod_noise,fc,fs,symbol_t,theta);
xx=find(x-y_demod_bit);
ber_temp=length(xx)/N;
ber=[ber,ber_temp];
end;


figure; semilogy(snr_plot,ber,'r*-'); grid ; 
xlabel('信噪比(dB)'); ylabel('误比特率');
title('误比特率曲线图');

pause(10);
figure;plot(peaks);  grid;