QAM.rar_qaM MATLAB_qam

%qam_display %qam分析的图形化显示 %参数 if parameter_display==1 disp(['调制使用的同相分量进制数M=',num2str(Mc),'个']); disp(['调制使用的正交分量进制数M=',num2str(Ms),'个']); disp(['调制使用的进制数M=',num2str(M),'个']); disp(['同相信号幅度距离dc=',num2str(dc),'v']); disp(['正交信号幅度距离ds=',num2str(ds),'v']); disp(['比特周期Tb=',num2str(Tb),'s']); disp(['比特速率Rb=',num2str(Rb),'Hz']); disp(['符号周期Ts=',num2str(Ts),'s']); disp(['符号速率Rs=',num2str(Rs),'波特']); disp(['脉冲持续时间L=',num2str(L),'*',num2str(Ts),'s']); disp(['一符号周期含有的载波数=',num2str(num_carriers_per_symbol),'个']); disp(['载波频率fc=',num2str(fc),'Hz']); disp(['脉冲类型=',num2str(gpulse_type),' 说明，1:矩形脉冲；2:升余弦脉冲；3:高斯脉冲']); disp(['仿真的符号数N=',num2str(N),'个']); disp(['仿真的比特数n=',num2str(bit_num),'个']); disp(['一个载波周期采样的点数=',num2str(num_sample_per_carrier),'点']); disp(['仿真结果时间采样间隔=',num2str(dt),'s']); end %波形 figure(1) stem(b);title('信源输出的比特信息流');xlabel('n');ylabel('b(n)');grid on;axis([0 0 -0.2 0.2]+axis);line(t,0); figure(2) subplot(211);stem(Amc);title('同相符号');xlabel('n');ylabel('Amc(n)');grid on;axis([1 1 1.2 1.2].*axis); subplot(212);stem(Ams);title('正交符号');xlabel('n');ylabel('Ams(n)');grid on;axis([1 1 1.2 1.2].*axis); figure(3) plot(t,g);title('成型脉冲g(t)');xlabel('t');ylabel('g(t)');grid on figure(4) plot(t,s);title('QAM调制信号波形s(t)');xlabel('t');ylabel('s(t)');grid on figure(5) subplot(311);plot(f,10*log10(Ps));title('QAM调制信号s(t)的功率谱密度');xlabel('f/Hz');ylabel('PSD(dB)/Watts/Hertz');u=axis;grid on; subplot(312);plot(f_c,10*log10(Psc));title('QAM调制同相信号sc(t)的功率谱密度');xlabel('f/Hz');ylabel('PSD(dB)/Watts/Hertz');axis(u);grid on; subplot(313);plot(f_s,10*log10(Pss));title('QAM调制正交信号sc(t)的功率谱密度');xlabel('f/Hz');ylabel('PSD(dB)/Watts/Hertz');axis(u);grid on; %星座图 figure(6) xc=[1:Mc];xs=[1:Ms]; yc=(2*xc-(Mc+1))*dc;ys=(2*xs-(Ms+1))*ds; for k=1:Mc plot(yc(k),ys,'ro');hold on end grid on; axis(1.2*axis); title([num2str(M),'QAM星座图'])