www.pudn.com > tfarma10.rar > param_tconv_b2.m, change:2004-01-19,size:1601b

function B2= param_tconv_b2(Bml1, Bml2, N, beta) % function B2= param_tconv_b2(Bml1, Bml2, N, beta) % This file is part of the TFPM toolbox v1.0 (c) % michael.jachan@tuwien.ac.at and underlies the GPL. % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clear;tfpm; MAR = 7; LAR = 5; MMA = 5; LMA = 3; N = 128; re_im= 'i'; mo_no= 'n'; tfpm_file_gen; %------------- beta= -1/2; Bml1= Aml; Bml2= Bml; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Dimensions: [M1, L1]= param_dim(Bml1); [M2, L2]= param_dim(Bml2); M= M1+M2; L= L1+L2; B2= zeros((M+1)*(2*L+1), (M1+1)*(2*L1+1)); % Compute B2 for m= 0:M for l= -L:L for mp= 0:M1 for lp= -L1:L1 B2((m*(2*L+1))+l+L+1, mp*(2*L1+1)+lp+L1+1)= ... param_get(Bml2, m-mp, l-lp) * exp(j*2*pi/N*tfpm_phi(m, l, mp, lp, -beta)); end; end; end; end; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Test B2: norm(param_destack_ml(B2*param_stack_ml(Bml1), M, L)-param_tconv(Bml1, Bml2, N, beta)) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%