Ես ուզում եմ մոդելավորել omega k ալգորիթմը՝ կենտրոնացնելու սինթետիկ բացվածքի ռադարի հումքի տվյալները՝ հիմնված Քամինգի գրքի վրա՝ Սինթետիկ բացվածքի ռադարի տվյալների թվային մշակումը: Սկզբում ես մոդելավորեցի կետային թիրախային չմշակված տվյալները stripmap ռեժիմում և արեցի այն ամենը, ինչ նշված է գրքում: Բայց իմ թիրախը կենտրոնացած չէ: Համոզվելու համար, որ իմ չմշակված տվյալները իսկապես պատրաստված են, ես դրանք կենտրոնացրել եմ սովորական RDA ալգորիթմի միջոցով, և իմ կետային թիրախը կենտրոնացել է իրական դիրքում, ինչը նշանակում է, որ իմ չմշակված տվյալների մոդելավորման ռեժիմը Ok է:
Ահա իմ matlab կոդը օմեգա k ալգորիթմի համար.
%% __________________________________________________________________________
fr = linspace(-fs/2,fs/2,nfftr);
faz = linspace(-PRF/2,PRF/2,nffta);
fr_prime = sqrt((f0+fr).^2-(c*faz'/(2*vp)).^2)-f0;
Rref = rs(ceil(Ns/2));
theta_ref = 4*pi*Rref/c*(fr_prime+f0)+pi*fr.^2/kr;
%2D FFT
S_raw = fftshift(fft2(s_raw,nffta,nfftr));
%RFM
S_BC = S_raw.*exp(1j*theta_ref);
for idx = 1:Na
S_int(idx,:) = interp1(fr_prime(idx,:)+f0,S_BC(idx,:),fr+f0,'pchip');
end
S_c = S_int.*exp(-1j*4*pi*fr*Rref/c);
s_c = ifft2(S_c,Na,Nr);
%% __________________________________________________________________________
այս ծածկագրում.
f0: կենտրոնական հաճախականությունkr. Ծլվլոցների մակարդակը միջակայքումfs. Նմուշառման հաճախականությունը տիրույթումvp: հարթակի արագությունrs՝ միջակայքի զանգված (ձևավորել մոտ տիրույթից մինչև հեռավոր միջակայք)Rref. Հղման տիրույթ (լսում եմ, որ այն որպես միջին տիրույթի բջիջ)Ns՝ տիրույթի բջիջների քանակըNa. նմուշների քանակը ազիմուտումs_c: Կենտրոնացված պատկեր
երեք թիրախներ տեղակայված են [10 , Ns/2 , Ns-10] հեռավորության վրա և Na/2 ազիմուտում: ահա իմ արդյունքները.
Ժամանակի տիրույթում զանգվածային սեղմումից հետո տվյալներ
Տվյալներ ժամանակի տիրույթում stolt interpolation-ից հետո
Ես ուսումնասիրեցի մի քանի ինտերպոլացիայի մեթոդներ, ինչպիսիք են sinc interp-ը, linear interp-ը, pchip-ը և այլն, բայց դրանցից ոչ մեկն ինձ մոտ չաշխատեց: Ես գնահատում եմ բոլորին, ովքեր կարող էին օգնել ինձ և ասել, թե որն է իմ սխալը... շնորհակալություն...
