资源描述
混沌通信中 QCSK 调制 matlab 代 码;allclear x=randsrc(20,1,0:1);%产生二进制随机码stairs(x);axis(0,20,-0.1,1.1);title(二进制随机序列);二讲扎轴祐:闻clcclearclose all % q=99;%k=1:99;%x(k)=sin(k*pi/q); x(1)=0.212345;for k = 1:99;x(k+1)=4 * x(k) * (1 - x(k);plot(x);legend(混沌信号x); grid on;%加网格clcclearclose all % q=99;%k=1:99;%x(k)=sin(k*pi/q);x(1)=0.212345;for k=1:99;x(k+1)=4*x(k)*(1-x(k);y=hilbert(x);figure(1) plot(imag(y);legend(希尔伯特变换y);grid onclcclearclose all% q=99;%k=1:99;%x(k)=sin(k*pi/q);x(1)=0.212345;%x 的初植for k=1:99;x(k+1)=4*x(k)*(1-x(k); end y=hilbert(x);%x 的希尔伯特变换 figure(1) plot(imag(y);grid onlegend(加密后的信号ms);clcclearclose all% q=99;%k=1:99; %x(k)=sin(k*pi/q); x(1)=0.212345;for k=1:99; x(k+1)=4*x(k)*(1-x(k);endy=hil bert(x);%希尔伯特变换 figure(1)plot(imag(y);grid onlegend(加密后的信号ms);y2=AWGN(imag(y),0.8,1);%imag(y) 为已调信号,0.8 为信噪比,1 为信号功率 figure(2)plot(y2);grid ony2);加噪声后的调制信号 legend(bs=1;Eb=1;Zb=bs.*Eb;%解调出的信号stem(Zb,r);title(解调出的信号Zb)clc;x=0:18;%信噪比的取植范围xSNR=10(x/20);%化成 dB 形式for i=1:length(x);k=2;y1(i)=1/2*erfc(1/(sqrt(4*1/xSNR (i)+k*(1/xSNR(i)A2);%求 k=2 的系统理论的误码率endsemilogy(x,y1);% 绘出 x 和 y1 的图for i=1:length(x)hold on;plot(x(i),y1(i),g*)end gtext(QCSK,K=2);%在光标的位置放置给定的文字for i=1:length(x);k=16;y2(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+k*(1/xSNR(i)A2);%求k=16 的系统理论的误码率endsemilogy(x,y2);% 绘出 x 和 y2 的图 for i=1:length(x)hold on; plot(x(i),y2(i),r.)end gtext(QCSK,K=16);for i=1:length(x);k=64;y3(i)=1/2*erfc(1/(sqrt(4*1/xSNR (i)+k*(1/xSNR( 的系统理论的误码率 k=64 求 i)A2);%endsemilogy(x,y3);%绘出 x 和 y3 的图 for i=1:length(x)hold on;plot(x(i),y3(i),b+) end gtext(QCSK,K=64);grid on xlabel(Eb/NodB); ylabel(BER); title(QCSK 的系统误码率);clc;x=0:18;%信噪比的取植范围 xSNR=10(x/20);%化成 dB 形式for i=1:length(x);k=2;y1(i)=1/2*erfc(1/(sqrt(4*1/xSNR (i)+k*(1/xSNR(i)A2);%求QCSK,k=2 的系统理论的误码率end semilogy(x,y1);for i=1:length(x)hold on;plot(x(i),y1(i),g*)endfor i=1:length(x); k=2;y2(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+2*k*(1/xSNR(i)入2);%求DCSK,k=2 的系统理论的误码率 endsemilogy(x,y2);for i=1:length(x)hold on;plot(x(i),y2(i),r.) endgtext(f绿色的*表示QCSK K=2 ;红色的. 表示 DCSK K=2);grid on xlabel(Eb/NodB); ylabel(BER););的系统误码率 title(K=2 clc;x=0:18;%信噪比的取植范围 xSNR=10(x/20);%化成 dB 形式for i=1:length(x);k=16;y1(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+k*(1/xSNR(i)A2);%求QCSK,k=16 的系统理论的误码率 endsemilogy(x,y1);for i=1:length(x)hold on;plot(x(i),y1(i),g*) end gtext(QCSK,K=16);for i=1:length(x); k=16;y2(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+2*k*(1/xSNR(i)入2);%求DCSK,k=16 的系统理论的误码率end semilogy(x,y2);for i=1:length(x)hold on; plot(x(i),y2(i),r.) endgtext(DCSK,K=16);grid on xlabel(Eb/NodB); ylabel(BER););时的系统误码率 title (K=16clc;x=0:18;%信噪比的取植范围 xSNR=10(x/20);%化成 dB 形式for i=1:length(x); k=64;y1(i)=1/2*erfc(1/(sqrt(4*1/xSNR (i)+k*(1/xSNR(i)A2);%求QCSK,k=64 的系统理论的误码率end semilogy(x,y1);for i=1:length(x)hold on;plot(x(i),y1(i),g*)endgtext(QCSK K=64);for i=1:length(x); k=64;y2(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+2*k*(1/xSNR(i)入2);%求DCSK,k=64 的系统理论的误码率 endsemilogy(x,y2);for i=1:length(x)hold on;plot(x(i),y2(i),r.) endgtext(DCSK K=64); grid on xlabel(Eb/NodB); ylabel(BER););的系统误码率 title(K=64
展开阅读全文