Subversion Repositories svnkaklik

Rev

Rev 640 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log

Rev 640 Rev 641
Line 374... Line 374...
374
            delayr[1] = n;
374
            delayr[1] = n;
375
            r = correlationr[n];
375
            r = correlationr[n];
376
        }
376
        }
377
    }
377
    }
378
 
378
 
-
 
379
//spocitejj frekvencni spektrum pro levy kanal
379
    for(i=delayl[1]; i < delayl[1] + chirp_size; i++) inchirp[i-delayl[1]] = L_signal[i];
380
    for(i=delayl[1]; i < delayl[1] + chirp_size; i++) inchirp[i-delayl[1]] = L_signal[i];
380
    fftw_execute(fft_plan_chirp);
381
    fftw_execute(fft_plan_chirp);
381
    for(i=0; i < frequency_bins; i++) lecho_spect[i] = sqrt(outchirp[i][0] * outchirp[i][0] + outchirp[i][1] * outchirp[i][1]);
382
    for(i=0; i < frequency_bins; i++) lecho_spect[i] = sqrt(outchirp[i][0] * outchirp[i][0] + outchirp[i][1] * outchirp[i][1]);
382
 
383
 
-
 
384
 
-
 
385
// napln pole daty z praveho kanalu a spocitej frekvencni spektrum
383
    for(i=delayr[1]; i < delayr[1] + chirp_size; i++) inchirp[i-delayr[1]] = R_signal[i];
386
    for(i=delayr[1]; i < delayr[1] + chirp_size; i++) inchirp[i-delayr[1]] = R_signal[i];
384
    fftw_execute(fft_plan_chirp);
387
    fftw_execute(fft_plan_chirp);
385
    for(i=0; i < frequency_bins; i++) recho_spect[i] = sqrt(outchirp[i][0] * outchirp[i][0] + outchirp[i][1] * outchirp[i][1]);
388
    for(i=0; i < frequency_bins; i++) recho_spect[i] = sqrt(outchirp[i][0] * outchirp[i][0] + outchirp[i][1] * outchirp[i][1]);
386
 
389
 
387
    printf("Writing output files\n");
390
    printf("Writing output files\n");