bmod2.m

交全面的给出了好几种调制解调方案的matlab代码,又ASK,BPSK,FSK,QAM等的源代码 对通信领域的研究人员很有帮助

function []=bmod2(x,y,Kmod,fs,fb,fc,y1);
%
% BBI 2000;绘制图形
n=length(y);
yy=fft(y)/n;yy=abs(yy(1:fix(n/2)))*2;
I=find(yy<4e-7);yy(I)=4e-7;yy=20*log10(yy)+10;
yy1=fft(y1)/n;yy1=abs(yy1(1:fix(n/2)))*2;
I=find(yy1<4e-7);yy1(I)=4e-7;yy1=20*log10(yy1)+10;
f1=fs/n;f=0:f1:(length(yy)-1)*f1;
if Kmod==3|Kmod==4;
    subplot('position',[0.05 0.11 0.75 0.3439]);
    plot(f,yy1,f,yy,'r');grid on;
    h=get(gca,'chidren');set(h(2),'color',[0 .99 .99]);
    title('Spectrum(dBm)'); xlabel('f(MHz)'); zoom xon;
elseif Kmod<3;
    subplot('postition',[0.05 0.11 0.45 0.3439]);
    plot(f,yy1,f,yy,'r');grid on;
    h=get(gca,'children');set(h(2),'color',[0 .99 .99]);
    title('Spectrum(dBm)'); xlabel('f(MHz)'); zoom xon;
    ph=phdet(y1,fs/fc);
    subplot('position',[0.52 0.11 0.3 0.3439]);
    H=polar(ph,ones(size(ph)),'ro');
    set(H,'markers',8,'linnew',2);xlabel('Phase');
elseif Kmod>4;
    subplot('position',[0.05 0.11 0.45 0.3439]);
    plot(f,yy1,f,yy,'r');grid on;
    h=get(gca,'children');set(h(2),'color',[0 .99 .99]);
    title('Spectrum(dBm)');xlabel('f(MHz)');zoom xon;
    subplot('position',[0.5 0.11 0.344 0.344]);
    constel(y1,fs,fb,fc);xlabel('Constellation');
end;zoom xon;