混沌通信中QCSK调制matlab代码

混沌通信中 QCSK调制 matlab 代码clear all;x=randsrc(20,1,[0:1]); %产生二进制随机码stairs(x); axis([0,20,-0.1,1.1]); title('二进制随机序列'); clc clear close 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)); end plot(x); legend('混沌信号 x'); grid on;% 加网格clc clear close 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)); end y=hilbert(x); figure(1) plot(imag(y)); legend('希尔伯特变换 y'); grid on clc clear close 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 on legend('加密后的信号 ms'); clc clear close 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)); end y=hilbert(x);% 希尔伯特变换figure(1) plot(imag(y)); grid on legend('加密后的信号 ms'); y2=AWGN(imag(y),0.8,1);%imag(y)为 已 调 信号， 0.8 为信噪比， 1 为信号功率figure(2) plot(y2); grid on legend('加噪声后的调制信号y2'); 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))^2))));% 求 k=2 的系统理论的误码率end semilogy(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))^2))));% 求 k=16 的系统理论的误码率end semilogy(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(i))^2))));% 求 k=64 的系统理论的误码率end semilogy(x,y3);% 绘出 x 和 y3 的图for i=1:length(x) hold on; plot(x(i),y3(i...