Matlab仿真FM0解码程序

clear T=18; SR=100; %采样点时间间隔参数 FM0_code=[0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 2 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0]; % 标签编码 aaa=zeros(1000,1); %start for i=1:size(FM0_code,2) if FM0_code(i)==1 if aaa((size(aaa,1)))==0 aaa((size(aaa,1)+1):(size(aaa,1)+SR))=ones(SR,1); else if aaa((size(aaa,1)))==1 aaa((size(aaa,1)+1):(size(aaa,1)+SR))=zeros(SR,1); end end end if FM0_code(i)==0 if aaa((size(aaa,1)))==0 aaa((size(aaa,1)+1):(size(aaa,1)+SR/2))=ones(SR/2,1); aaa((size(aaa,1)+1):(size(aaa,1)+SR/2))=zeros(SR/2,1); else if aaa((size(aaa,1)))==1 aaa((size(aaa,1)+1):(size(aaa,1)+SR/2))=zeros(SR/2,1); aaa((size(aaa,1)+1):(size(aaa,1)+SR/2))=ones(SR/2,1); end end end if FM0_code(i)==2 if aaa((size(aaa,1)))==0 aaa((size(aaa,1)+1):(size(aaa,1)+SR))=zeros(SR,1); end if aaa((size(aaa,1)))==1 aaa((size(aaa,1)+1):(size(aaa,1)+SR))=ones(SR,1); end end end maxa=1; mina=0.8; for i=1:size(aaa,1) aaa(i)=mina+(maxa-mina)*aaa(i); Isig(i)=aaa(i)*sin(mod((i-1),T)*2*pi/T); Qsig(i)=aaa(i)*cos(mod((i-1),T)*2*pi/T); IQ(i)=sqrt(Isig(i)^2+Qsig(i)^2); end Y=Isig+j*Qsig; addpath 'E:\Program Files\Agilent\Signal Generators\WaveformDownloadAssistant' io = agt_newconnection('tcpip','192.168.0.98'); [status, status_description] = agt_sendcommand(io, 'SOURce:FREQuency 922500000'); [status, status_description] = agt_sendcommand(io, 'POWer 0'); [status, status_description] = agt_sendcommand(io,':source:rad:arb:state off'); if (status < 0) return; end [status, status_description] = agt_waveformload(io, Y, 'FM0', 80000000, 'play','no_normscale'); [status, status_description] = agt_sendcommand(io,':source:rad:arb:state on'); [status, status_description] = agt_sendcommand(io, 'OUTPut:STATe ON'); %agt_closeAllSessions;