OFDM.rar_OFDM系统仿真_hatpy7_ofdm_ofdm流图

%综合实验二 N=256; %子载波数,FFT取256 NUM=2; %每个子载波传递的信号数，设为2 NG=1/16*N; %CP长度，复制16个到前面形成CP Xn=randi([0,1],N,NUM); %Xn为比特流（0和1） 数字信号 Sn_1=Xn(1:N,1);Sn_2=Xn(1:N,2); %由于抽头数较多，超过前缀长度，因此考虑两个相邻OFDM符号,但是主要考察的是Sn_2信号 scatterplot(Sn_2);title('输入信号星座图'); %画出输入星座图 sn_1=(Sn_1-0.5).*2;sn_2=(Sn_2-0.5).*2; %BPSK调制，将0,1序列转化为-1，1序列 s_t1=ifft(sn_1);s_t2=ifft(sn_2); SER=zeros(1,5); %存放误码率 for H_NUM=[10,20,64] %抽头数为10,20,64 %%%%%%%%%%%功率均匀分布的抽头%%%%%%%%%% figure; for i=0:5:20 %信噪比 s_T_add=[s_t1(N-NG+1:N);s_t1;s_t2(N-NG+1:N);s_t2].'; %插入保护循环前缀 h=randn(1,H_NUM); %抽头为正态分布的10个数 s_h=conv(s_T_add,h); %通过多径信道；加抽头的作用相当于循环卷积 s_h_noise=awgn(s_h,i,'measured'); %加0-20dB噪声 s_out_t=s_h_noise(N+2*NG+1:2*N+2*NG); %去掉前缀和后缀,只保留中间部分 s_out=fft(s_out_t); H=fft(h,N); S_out=s_out./H; %相当于功率归一化 S_OUT=S_out/2+0.5; %解调 S_OUT=S_OUT'; %串/并行转化 x_out=round(real(S_OUT),3); subplot(2,3,i/5+1); plot(S_OUT,'.'); axis([-2 2 -2 2]); title(['加',num2str(i),'dB噪声|抽头数',num2str(H_NUM),'（功率均匀）对应输出信号星座图']); X_out=1*(real(x_out)>0.5); %就近判断输出信号，靠近0判定为0，靠近1判定为1，以0.5为标准 SER(i/5+1)=symerr(X_out,Sn_2)./N*100; end subplot(2,3,6); plot(0:5:20,SER); xlabel('信噪比/dB');ylabel('误码率/%'); title(['多径信道数为',num2str(H_NUM),'|抽头功率均匀分布的误码率图']); %%%%%%%%%功率指数分布的抽头%%%%%%% expo=power(1:H_NUM,0.5); %指数分布系数设为0.5 figure; SER_2=zeros(1,5); for i=0:5:20 %信噪比 s_T_add=[s_t1(N-NG+1:N);s_t1;s_t2(N-NG+1:N);s_t2].'; %插入保护循环前缀 h=randn(1,H_NUM)./expo; %抽头为正态分布的10个数 s_h=conv(s_T_add,h); %通过多径信道；加抽头的作用相当于循环卷积 s_h_noise=awgn(s_h,i,'measured'); %加0-20dB噪声 s_out_t=s_h_noise(N+2*NG+1:2*N+2*NG); %去掉前缀和后缀,只保留中间部分 s_out=fft(s_out_t); H=fft(h,N); S_out=s_out./H; %相当于功率归一化 S_OUT=S_out/2+0.5; %解调 S_OUT=S_OUT'; %串/并行转化 x_out=round(real(S_OUT),3); subplot(2,3,i/5+1); plot(S_OUT,'.'); axis([-2 2 -2 2]); title(['加',num2str(i),'dB噪声|抽头数',num2str(H_NUM),'功率指数分布对应输出信号星座图']); X_out=1*(real(x_out)>0.5); %就近判断输出信号，靠近0判定为0，靠近1判定为1，以0.5为标准 SER_2(i/5+1)=symerr(X_out,Sn_2)./N*100; end subplot(2,3,6); plot(0:5:20,SER_2); xlabel('信噪比/dB');ylabel('误码率/%'); title(['多径信道数为',num2str(H_NUM),'|抽头功率指数分布的误码率图']); end