newsubsampling315.m

是信号处理方面欠采样的源程序 ,是很不错的程序,可以这方面的学者一些参考价值

  clear all;
  clc;
  close all;
  k=zeros(1,18);
%   f01_temp=zeros(1,21);
  errork=zeros(1,18);
  errorf=zeros(1,18);
  sumerrork=zeros(1,18);
  sumerrorf=zeros(1,18);
  del=5*10^-6;
  del2=8*10^-6;
  unit=1e6;
  fs=100*unit;
  f1=150*unit;Bw=600*unit;T=20*10^-6;k1=Bw/T;
  t1=0:1/fs:T;
  lt1=length(t1);
  nn1=round(del/T*lt1);
  nn2=round(del2/T*lt1);
  ad1=fs;%100*unit;
%    for a=-7:10;
%     for b=1:100;
         %********估计k值*********%
          SNR=10;
          j=sqrt(-1);
          s1=exp(j*2*pi*(f1*t1+k1*t1.^2/2) +j*0.6*pi);
          s1=awgn(s1,SNR,'measured');
          r1=s1(nn1+1:lt1).*conj(s1(1:(lt1-nn1)));
          r2=s1(nn2+1:lt1).*conj(s1(1:(lt1-nn2)));
          fr=abs(fft(r1));
          fr2=abs(fft(r2));
%           figure
%           plot(fr);
%           figure
%           plot(fr2);
          [marg,ipeak]=max(fr);
          [marg2,ipeak2]=max(fr2);
          [f11_1,a11]=pufeng_xihua_fushu(fr,fs,ipeak);
          [f12_1,a12]=pufeng_xihua_fushu(fr2,fs,ipeak2);
           for ii=1:2
              for jj=1:2
                  kk1=f11_1+fs*ii;
                  kk2=f12_1+fs*jj;
                  compare1=kk2/kk1;
                  if(abs(compare1-8/5)<0.05)
                      kt1_1=kk1;
                      kt2_1=kk2;
                      break;
                  end
              end
          end
          k0_1=kt1_1/del;
%            k00_1=kt2_1/del2;
         x2t=exp(j*2*pi*(f1*t1+(k1-k0_1)*t1.^2/2));
         fr3=abs(fft(x2t));
         [marg3,ipeak3]=max(fr3);
         [f12_3,a13]=pufeng_xihua_fushu(fr3,fs,ipeak3);
         f0test=fs+f12_3;
         errof=abs(f0test-f1)/fs
%          errork(a+8)=errork(a+8)+(k0_1-k1)^2;
%          errorf(a+8)=errorf(a+8)+(f0test-f1)^2;
%       end
%          sumerrork(a+8)=sqrt(errork(a+8)/100);
%          sumerrorf(a+8)=sqrt(errorf(a+8)/100);
%      end
%        xs=linspace(-7,10,18);
%        figure   
%        plot(xs,sumerrork);
%        xlabel('SNR/dB');
%        ylabel('RMSE(k)/(Hz/s)');   
%        figure
%        plot(xs,sumerrorf);
%        xlabel('SNR/dB');
%        ylabel('RMSE(f0)/Hz');