QPSK调制与解调在MATLAB平台上的实现

程序没有错，不是直接运行的，你需要建立三个.m文件，分别把两个Function和最上面的程序粘贴上去。我修改了一下你只需要建立两个.m文件然后运行第一个文件就可以得到结备辩差果。

SNRindB1=0:2:10

SNRindB2=0:0.1:10

for i=1:length(SNRindB1)

[pb,ps]=cm_sm32(SNRindB1(i))

smld_bit_err_prb(i)=pb

smld_symbol_err_prb(i)=ps

end

for i=1:length(SNRindB2)

SNR=exp(SNRindB2(i)*log(10)/10)

theo_err_prb(i)=(1/2)*erfc(sqrt(2*SNR)/sqrt(2))

end

title('QPSK误码率分析灶斗')

semilogy(SNRindB1,smld_bit_err_prb,'*')

axis([0 10 10e-8 1])

hold on

% semilogy(SNRindB1,smld_symbol_err_prb,'仿皮o')

semilogy(SNRindB2,theo_err_prb)

legend('仿真比特误码率','理论比特误码率')

hold off

另一个.m文件

function [pb,ps]=cm_sm32(SNRindB)

N=10000

E=1

SNR=10^(SNRindB/10)

sgma=sqrt(E/SNR)/2

s00=[1 0]

s01=[0 1]

s11=[-1 0]

s10=[0 -1]

for i=1:N

temp=rand

if (temp<0.25)

dsource1(i)=0

dsource2(i)=0

elseif (temp<0.5)

dsource1(i)=0

dsource2(i)=1

elseif (temp<0.75)

dsource1(i)=1

dsource2(i)=0

else

dsource1(i)=1

dsource2(i)=1

end

end

numofsymbolerror=0

numofbiterror=0

for i=1:N

n=sgma*randn(size(s00))

if((dsource1(i)==0)&(dsource2(i)==0))

r=s00+n

elseif((dsource1(i)==0)&(dsource2(i)==1))

r=s01+n

elseif((dsource1(i)==1)&(dsource2(i)==0))

r=s10+n

else

r=s11+n

end

c00=dot(r,s00)

c01=dot(r,s01)

c10=dot(r,s10)

c11=dot(r,s11)

c_max=max([c00 c01 c10 c11])

if (c00==c_max)

decis1=0decis2=0

elseif(c01==c_max)

decis1=0decis2=1

elseif(c10==c_max)

decis1=1decis2=0

else

decis1=1decis2=1

end

symbolerror=0

if(decis1~=dsource1(i))

numofbiterror=numofbiterror+1

symbolerror=1

end

if(decis2~=dsource2(i))

numofbiterror=numofbiterror+1

symbolerror=1

end

if(symbolerror==1)

numofsymbolerror=numofsymbolerror+1

end

end

ps=numofsymbolerror/N

pb=numofbiterror/(2*N)

记得要保存为cm_sm32.m,然后运行前面那个.m文件就可以了。

QPSK调制程序的testbench程序如下：

********************************************************************************

`timescale 1ns/1ns //单位时间，时间精度

module qpsk_tb//qpsk调制的testbench

regclk

regrst

regx

wire y

qpsk qpsk(.clk(clk),.rst(rst),.x(x),.y(y)) //将tb输入输出与主程序qpsk关联起来

always #25 clk=~clk //#25表示25个时间单位,

// 这句表示clk信号的周期是50（*1ns）

initial

begin

clk=0x=0rst=0//赋初值

#10rst=1x=1 /宴埋判/间隔10ns，将复位信号置为1

#10x=0

#400x=1 //时钟周期是50ns，晌改这里的信号表示的波形为

#400x=0 //1111111100000000111100001100

#200x=1

#200x=0

#100x=1

#100x=0

end

endmodule

QPSK解调程序的testbench程序如下：

********************************************************************************

`timescale 1ns/1ns //单位时间，时间精度

module qpsk_detb //qpsk解调的testbench

regclk

regrst

regx

wire y

qpsk_de qpsk_de(.clk(clk),.rst(rst),.x(x),.y(y))//将液茄tb输入输出与主程序qpsk_de关联起来

always #25 clk=~clk //#25表示25个时间单位

initial

begin

clk=0x=0rst=0

#10rst=1x=1#10x=0 //时钟周期是50ns，这里的信号表示的波形为

#400x=1#400x=0 //1111111100000000

#200x=1#200x=0 //11110000

#200x=0#200x=1 //00001111

#100x=0#200x=1#100x=0 //00111100

#100x=1#200x=0#100x=1 //11000011

#100x=1#100x=0 //1100

end

endmodule

---