pie_code.asv

pie edcode编码 程序设计

clc;
clear all;

fs=0.08*10^6;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
Tari=25*10^(-6);

count0=0;
count1=0;
for j=1:20
    x=round(rand);
    if x==0

%  %fs*Tari=1250;   fs=50*10^6;   Tari=25*10^(-6);
%         for t=1 : 1250/2                                     %fs*Tari=1250
%             Q(1,count0*1250+count1*1250*2+t)=1;
%         end
%         for t=1250/2+1 : 1250
%             Q(1,count0*1250+count1*1250*2+t)=0;
%         end
%         count0=count0+1;
%     else
%         for t=1 : 1250*3/2
%             Q(1,count0*1250+count1*1250*2+t)=1;
%         end
%         for t=1250*3/2+1 : 1250*2
%             Q(1,count0*1250+count1*1250*2+t)=0;
%         end
%             count1=count1+1;
%     end
%     
    
%  %fs*Tari=100;   fs=4*10^6;   Tari=25*10^(-6);
%         for t=1 : 100/2                                     %fs*Tari=100
%             Q(1,count0*100+count1*100*2+t)=1;
%         end
%         for t=100/2+1 : 100
%             Q(1,count0*100+count1*100*2+t)=0;
%         end
%         count0=count0+1;
%     else
%         for t=1 : 100*3/2
%             Q(1,count0*100+count1*100*2+t)=1;
%         end
%         for t=100*3/2+1 : 100*2
%             Q(1,count0*100+count1*100*2+t)=0;
%         end
%             count1=count1+1;
%     end    
    


%  %fs*Tari=20;   fs=0.8*10^6;   Tari=25*10^(-6);
            for t=1 : 2/2                                     %fs*Tari=100
            Q(1,count0*2+count1*2*2+t)=1;
        end
        for t=2/2+1 : 2
            Q(1,count0*2+count1*2*2+t)=0;
        end
        count0=count0+1;
        
% % % % % % % % % % % % % % % % % % % % % % % % % % %          
                        for t0=1 : 20/2                                     %fs*Tari=100
            Q0(1,count00*20+count10*20*2+t0)=1;
        end
        for t=20/2+1 : 20
            Q0(1,count00*20+count10*20*2+t0)=0;
        end
        count00=count00+1;
% % % % % % % % % % % % % % % % % % % % % % % % % % %         
        
    else
        for t=1 : 2*3/2
            Q(1,count0*2+count1*2*2+t)=1;
        end
        for t=2*3/2+1 : 2*2
            Q(1,count0*2+count1*2*2+t)=0;
        end
            count1=count1+1;

% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %    
        for t=1 : 20*3/2
            Q0(1,count00*20+count10*20*2+t0)=1;
        end
        for t=20*3/2+1 : 20*2
            Q0(1,count00*20+count10*20*2+t)=0;
        end
            count10=count10+1;
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %
% % % 
end
end
Q1=upsample(Q,8);
Y=firrcos(150,40000,0,640000);%
Y=conv(Y,Q1);
% Y=rcosflt(Q,80000,640000,'fir/normal',0,48);


% Y1=fir1(150,0.1);
% Y1=conv(Y1,Q);


figure(1)
subplot(221)
plot(Q);
subplot(222)
plot(20*log10(abs(fft(Q))));

subplot(223)
plot(Y);
subplot(224)
plot(20*log10(abs(fft(Y))));
% 
% subplot(325)
% plot(Y1);
% subplot(326)
% plot(20*log10(abs(fft(Y1))));
figure(2)
subplot(321)
plot(Q);
subplot(322)
plot(20*log10(abs(fft(Q))));

subplot(323)
plot(Q1);
subplot(324)
plot(20*log10(abs(fft(Q1))));

subplot(325)
plot(Q0);
subplot(326)
plot(20*log10(abs(fft(Q0))));