主程序：.doc

主程序： function pb=main(SNR) if nargin1 SNR=5; end frame_num=2000; input=randn(1,frame_num*6); for i=1:frame_num*6 if input(i)0 input(i)=1; else input(i)=0; end end Tb=0.01; fc=2000; frame_length=6; %figure(1) %subplot(2,1,1); %stem(input); %xlabel(n); %title(input bit stream); %input1=diff_qam(input); size1=size(input); bit_out=zeros(1,size1(2)); for i=1:size1(2)/frame_length m=(i-1)*frame_length+1; y=qam_module(input(1,m:(m+5)),Tb,i);%调制部分 y=raylfading_channel(y);%瑞利衰落信道 y=awgn(y,SNR); %高斯白噪声信道 y=equalization(y);%频域均衡 bit_out(1,[((i-1)*frame_length+1):(i*frame_length)])=qam_demodule(y,Tb,fc,i);%解调部分 end %bit_out=int_qam(bit_out1); %figure(1) %subplot(2,1,2); %stem(bit_out); %xlabel(n); %title(output bit stream); error=abs(bit_out-input); error_num=sum(error); pb=sum(error)/(frame_num*6); %误码率 fprintf(\n); fprintf(************** Eb/No = %5.2f dB ***************\n, SNR); fprintf(***********************************************\n); fprintf(Frame size = %d\n, frame_num); fprintf(%d bits transmitted, %d bits in error.\n,frame_num*6, error_num); fprintf(Bit Error Rate %1.4e.\n, pb); fprintf(***********************************************\n\n); 64QAM调制程序： function [yout]=bandpass_trans(x,Tb,sequence_num,fc) if nargin4 fc=2000 if nargin3 sequence_num=1 end if nargin2 Tb=0.01 end if nargin1 x=[0 0 0 0 0 0] end Ts=Tb/100 [x1(1)]=xor(0,x(1)) [x1(2)]=xor(x1(1),x(2)) [x1(3)]=xor(x1(2),x(3)) [y1(1)]=xor(0,x(4)) [y1(2)]=xor(y1(1),x(5)) [y1(3)]=xor(y1(2),x(6)) xm=-4+(x1(1)*4+x1(2)*2+x1(3)) if xm=0 xm=xm+1 end ym=-4+(y1(1)*4+y1(2)*2+y1(3)) if ym=0 ym=ym+1 end t=(sequence_num-1)*Tb:Ts:sequence_num*Tb-Ts yout=xm*cos(2*pi*fc*t)+ym*sin(2*pi*fc*t) figure(2) plot(t,yout) hold on grid end 瑞利衰落信道程序： function[yout]=raylfading_channel(signal) if nargin1 signal=ones(1,200) end size1=size(signal) length=size1(2) noise=raylrnd(1,1,length) yout=noise.*signal x=su