赞
踩
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
% qpsk调制,两天线误码率计算
clc
clear
N =10000; % 符号数
Es_N0_dB = [-3:15];
ipHat = zeros(1,N);
h1=chann(N,1);
h2=chann(N,1);
ip = (2*(rand(1,N)>0.5)-1) + j*(2*(rand(1,N)>0.5)-1); %
s = (1/sqrt(2))*ip; % normalization of energy to 1 归一化能量为一
n = 1/sqrt(2)*[randn(1,N) + j*randn(1,N)]; % 白噪声, 0dB
for ii = 1:length(Es_N0_dB)
for i=1:N
%%%%%
w1=conj(h1(i))/abs(h1(i))^2;
w2=conj(h2(i))/abs(h2(i))^2;
y(i)=(w1*h1(i)+w2*h2(i))/2.7*s(i)+10^(-Es_N0_dB(ii)/18)*n(i);
end;
%%%%%%%
% 解调
y_re = real(y); % 实部
y_im = imag(y); % 虚部
ipHat(find(y_re < 0 & y_im < 0)) = -1 + -1*j;
ipHat(find(y_re >= 0 & y_im > 0)) = 1 + 1*j;
ipHat(find(y_re < 0 & y_im >= 0)) = -1 + 1*j;
ipHat(find(y_re >= 0 & y_im < 0)) = 1 - 1*j;
nErr(ii) = size(find([ip- ipHat]),2) % 计算错误字符数
end
simSer_QPSK = nErr/N;
% figure
semilogy(Es_N0_dB,simSer_QPSK,'r+-');%红色
hold on
axis([-3 15 10^-5 1])
grid on
% legend( '2天线');
xlabel('Es/No(dB)')
ylabel('SER')
title('理想信道波束形成误符号率(QPSK调制)')
该程序是理想信道条件下的波束形成仿真
请问如何加入时延t使程序变为时延信道下波束形成仿真
紧急 谢谢啦
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。