sdr.rar_DDC DUC_ddc_duc_interpolation

clear clc %% set constelation % QPSK - 1 % 8PSK - 2 % 16QAM - 3 % 64QAM - 4 set_constel = 3; switch set_constel case 1 dif_cod = 1; base = 4; f1 = [1 0 0 0]; case 2 dif_cod = 2; base = 8; f1 = [0 1 0 0]; case 3 dif_cod = 4; base = 16; f1 = [0 0 1 0]; otherwise dif_cod = 16; base = 64; f1 = [0 0 0 1]; end % Tsym - 縊徵蝓贛�� 縈褒褓靦 閱擇滓鞣 慚縴整皖態 Tsym = 1/(1 ); % Eb/N0 - 褕瘟羅璃� �/� � 蒙熱雒 AWGN EbN = 20; % Ph_offset - 糟麂順� 篩閎� 縊廉畫� � 蒙熱雒 Ph_offset = pi/3; % FRQ_mu - 闐羲褕� 愿慮獗璃� 鳩皓粵�� "撰" (mu) 頡滓鈇粳蝶� 熱 慚鋼暹蝓�錙� % 哮曄褊 溫� 縴整截蝴痢�� 滔螂瘟� 頑滇擂蒐 (愿慮獗璃� 頜蝙翮蝠摘 頡 % 縊皚縝彰畫贛衛� 頑蜻皚 � 番鞍摘擊蝗 {-1...1} � 闐羲褕蝗 FRQ_mu). FRQ_mu = 0.001*Tsym; % AMP_f - 鳩皓粵徹 愿慮獗璃� 闐羲褕� 閱擇滓鈇粳蝶� 縊廉畫� � � (愿慮獗璃� % 頜蝙翮蝠摘 頡 縊皚縝彰畫贛衛� 頑蜻皚 � 番鞍摘擊蝗 (0...0,05/Tsym - 5% 褕 % 縊徵蝓贛蝗 縈褒褓錚) � 闐羲褕蝗 FREQ_f AMP_f = 0; % 0.05/Tsym FREQ_f = 0; % 2*pi*200*Tsym %% %看蜘蝴�� 闐羲褕� 海� FreqDAC = 1200 ; %看蜘蝴�� 闐羲褕� ADC FreqADC = 192; % FREQ_del - 闐羲褕� 愿慮獗璃� 獵餃� 頑滇擂蒐 縊廉畫� � 蒙熱雒 FREQ_del = 0; % 2*pi*0*Tsym %押翰麴� 闐羲褕� � 頜�憐廈 394 慣� 螂擋錮蝶� 羲隆閤 434 呸� Freq_car = 394 ; %閃� 篦錦� 闐羲褕 �奏 N = 4 %1...16 Freq_car_step = N*0.465 ; %縝幌痘蝴痢磋� 繃鼻隧 RRC_up = 6; rolloff = 0.9; %閱順� dds 32呸� 閱擇羲曉疑� 獗翰嚼� Fdds1 = 30 ; Nup1 = Fdds1*(Tsym/RRC_up); %頡餉� RRC %second dds 1200 ? Fdds2 = Freq_car; Nup2 = FreqDAC/(Fdds1); %青擇滇縈擇錚頑甕� %N_up = FreqDAC/Fdds1; car_out = Freq_car_step+Freq_car; Nyq_zone = floor(car_out/FreqADC); rxNCO =car_out - Nyq_zone*FreqADC; % %看萄摟� 縊皚築 � 蜻縊皚築 % Ns = 512; % n=0:360/Ns:360-360/Ns; Bandpas =2/((1+rolloff)*Tsym); %MHz Filterband = Bandpas*2/FreqADC; %down conversion Ndownrx = FreqADC/RRC_up * 4; Nuprx = 4; %cic_ii_0_example_design_fir_comp_coeff(R, L, Fs, Fc) h = cic_ii_0_example_design_fir_comp_coeff(FreqADC/RRC_up, 31, FreqADC, RRC_up/2);%fir_compensator(FreqADC/RRC_up, 4, FreqADC, Bandpas ); %fir_compensator( R, N, Fs, Fc ) %0.625*1/Tsym B = 12; %% Coeffi. Bit-width hz = round(h*power(2,B-1)-1); %% Fixed point coefficients KcicD = N*20*log10( FreqADC/RRC_up) %up conversion h1 = cic_ii_0_example_design_fir_comp_coeff(Nup1, 31, RRC_up*1/Tsym, 0.6/Tsym);%fir_compensator(FreqADC/RRC_up, 4, FreqADC, Bandpas ); %fir_compensator( R, N, Fs, Fc ) %0.625*1/Tsym %h1 = fir_compensator(Nup1,4, RRC_up*1/Tsym, Bandpas/2 ); hz1 = round(h1*power(2,B-1)-1); %% Fixed point coefficients KcicI = N*20*log10( Nup1)