资源描述：
《混沌通信中QCSK调制matlab代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、clearall;x=randsrc(20,1,[0:1]);%产生二进制随机码stairs(x);axis([0,20,-0.1,1.1]);title('二进制随机序列');clcclearcloseall%q=99;%k=[1:99];%x(k)=sin(k*pi/q);x(1)=0.;fork=1:99;x(k+1)=4*x(k)*(1-x(k));endplot(x);legend('混沌信号x');gridon;%加网格clcclearcloseall%q=99;%k=[1:99];%x(k)=sin(k*pi/q);x(1)=0.;fork=1:99;x(k+1)=
2、4*x(k)*(1-x(k));endy=hilbert(x);figure(1)plot(imag(y));legend('希尔伯特变换y');gridonclcclearcloseall%q=99;%k=[1:99];%x(k)=sin(k*pi/q);x(1)=0.;%x的初植fork=1:99;x(k+1)=4*x(k)*(1-x(k));endy=hilbert(x);%x的希尔伯特变换figure(1)plot(imag(y));gridonlegend('加密后的信号ms');clcclearcloseall%q=99;%k=[1:99];%x(k)=sin(k*p
3、i/q);x(1)=0.;fork=1:99;x(k+1)=4*x(k)*(1-x(k));endy=hilbert(x);%希尔伯特变换figure(1)plot(imag(y));gridonlegend('加密后的信号ms');y2=AWGN(imag(y),0.8,1);%imag(y)为已调信号，0.8为信噪比，1为信号功率figure(2)plot(y2);gridonlegend('加噪声后的调制信号y2');bs=1;Eb=1;Zb=bs.*Eb;%解调出的信号stem(Zb,'r');title('解调出的信号Zb')clc;x=0:18;%信噪比的取植范围xS
4、NR=10.^(x/20);%化成dB形式fori=1:length(x);k=2;y1(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+k*((1/xSNR(i))^2))));%求k=2的系统理论的误码率endsemilogy(x,y1);%绘出x和y1的图fori=1:length(x)holdon;plot(x(i),y1(i),'g*')endgtext('QCSK,K=2');%在光标的位置放置给定的文字fori=1:length(x);k=16;y2(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+k*((1/xSNR(i))^2))
5、));%求k=16的系统理论的误码率endsemilogy(x,y2);%绘出x和y2的图fori=1:length(x)holdon;plot(x(i),y2(i),'r.')endgtext('QCSK,K=16');fori=1:length(x);k=64;y3(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+k*((1/xSNR(i))^2))));%求k=64的系统理论的误码率endsemilogy(x,y3);%绘出x和y3的图fori=1:length(x)holdon;plot(x(i),y3(i),'b+')endgtext('QCSK,K=6
6、4');gridonxlabel('Eb/No[dB]');ylabel('BER');title('QCSK的系统误码率');clc;x=0:18;%信噪比的取植范围xSNR=10.^(x/20);%化成dB形式fori=1:length(x);k=2;y1(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+k*((1/xSNR(i))^2))));%求QCSK，k=2的系统理论的误码率endsemilogy(x,y1);fori=1:length(x)holdon;plot(x(i),y1(i),'g*')endfori=1:length(x);k=2;y2(i
7、)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+2*k*((1/xSNR(i))^2))));%求DCSK，k=2的系统理论的误码率endsemilogy(x,y2);fori=1:length(x)holdon;plot(x(i),y2(i),'r.')endgtext('绿色的*表示QCSKK=2;红色的.表示DCSKK=2');gridonxlabel('Eb/No[dB]');ylabel('BER');title('K=2的系统误码率');clc;x=0:
