混沌通信中QCSK調制matlab代碼

%產生二進制隨機碼clearall;%產生二進制隨機碼axis([0,20,-0.1,1.1]);title('二進制隨機序列')；clcclearcloseall%q=99;%k=[1:99];%x(k)=sin(k*pi/q);x(1)=0?212345;fork=1:99;x(k+1)=4*x(k)*(1-x(k));endplot(x);legend('混沌信號x');gridon;%加網(wǎng)格clcclearcloseall%q=99;%k=[1:99];%x(k)=sin(k*pi/q);x(1)=0?212345;fork=1:99;x(k+1)=4*x(k)*(1-x(k));endy=hilbert(x);figure⑴plot(imag(y));legend('希爾伯特變換y');gridonclcclearcloseall%q=99;

%k=[1:99];%x(k)=sin(k*pi/q);x(1)=0?212345;%x的初植fork=1:99;x(k+1)=4*x(k)*(1-x(k));endy=hilbert(x);%x的希爾伯特變換figure(1)plot(imag(y));gridonlegend(咖密后的信號ms');□.8□.5Cl.4□-0.2-0.4-0.6-□.81020304050GO708090100clcclearcloseall%q=99;%k=[1:99];%x(k)=sin(k*pi/q);x(1)=0?212345;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);bs=1;Eb=1;Zb=bs?*Eb;%解調岀的信號stem(Zb,'r');title('解調出的信號Zb')1.9.B.7.E.5.4.3.2□.1n解調出的信號孔IlliII11Illi11I1□ 0.2 0.4 0.6 0.8 1 12 1.4 16 1.8 2clc;x=0:18;%信噪比的取植范圍xSNR=10.A(x/20);%化成dB形式fori=1:length(x);k=2;y1(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+k*((1/xSNR(i))A2))));%求k=2的系統(tǒng)理論的誤碼率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))A2))));%求k=16的系統(tǒng)理論的誤碼率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))A2))));%求k=64的系統(tǒng)理論的誤碼率endsemilogy(x,y3);%繪出x和y3的圖fori=1:length(x)holdon;plot(x(i),y3(i),'b+')endgtext('QCSK,K=64');gridonxlabel('Eb/No[dB]');ylabel('BER');title('QCSK的系統(tǒng)誤碼率');

clc;EJx=0:18;%信噪比的取植范圍xSNR=10.A(x/20);%化成dB形式fori=1:length(x);EJk=2;semilogy(x,y1);fori=1:length(x)holdon;plot(x(i),y1(i),'g*')endfori=1:length(x);k=2;y2(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+2*k*((1/xSNR(i))A2))));%求DCSK,k=2的系統(tǒng)理論的誤碼endsemilogy(x,y2);fori=1:length(x)holdon;plot(x(i),y2(i),'r?')endgtext('綠色的*表示QCSKK=2;紅色的?表示DCSKK=2');gridonxlabel('Eb/No[dB]');clc;x=0:18;%信噪比的取植范圍xSNR=10.A(x/20);%化成dB形式fori=1:length(x);k=16;EJsemilogy(x,y1);fori=1:length(x)holdon;plot(x(i),y1(i),'g*')endgtext('QCSK,K=16');fori=1:length(x);k=16;y2(i)=1/2*erfc(1/(sqrt(4*1/xSNR(i)+2*k*((1/xSN率endsemilogy(x,y2);fori=1:length(x)holdon;plot(x(i),y2(i),'r.')endgtext('DCSK,K=16');gridonxlabel('Eb/No[dB]');clc;x=0:18;%信噪比的取植范圍xSNR=104(x/20);%化成dB形式fori=1:length(x);k=64;semilogy(x,y1);fori=1:length(x)holdon;plot(x(i),y1(i),'g*')endgtext('QCSKK=64');fori=1:length(x);k=64;y2(i)=1/2*erfc(1/(sqrt(4*1/