fm0.m

使用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;