?? nonwssus_rs.m
字號:
function RS= nonWSSUS_RS(Lalpha)% function RS= nonWSSUS_RS(Lalpha)% 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;tfpmalpha= 1/2;N= 256;J= 5;B= J+1;offset= 200;M= 4;L= 3;filename= sprintf('chpar%02d%04d%02d%02d-%04d', J, N, M, L, offset)load(filename)BML= B0L;Lalpha= vtfarma_weyl(AML, BML, N, B, 1/2);for j= 0:J for jp= 0:J mesh(abs(reshape(Lalpha(j+1, jp+1, :, :), N, N))) drawnow end;end;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Dimensions[J, J, N, N]= size(Lalpha);J= J-1;RS= zeros(J+1, N, J+1, N);for tau= 0:J for Dtau= 0:tau rs= zeros(N); for taup= 0:J L1= reshape(Lalpha(tau+1 , taup+1, :, :), N, N); L2= reshape(Lalpha(tau-Dtau+1, taup+1, :, :), N, N); rs= rs + L1.*conj(L2); end; rs1= fftshift(fft(rs), 1); RS(tau+1, :, Dtau+1, :)= reshape(rs1.', 1, N, 1, N); end;end;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%RRSS= fftshift(ifft(fftshift(RS, 4), [], 4), 4);% PLOT LSF @ k= 0k= 0;for n= 0:N-1 rrss= fftshift(fft(RRSS(:, :, :, n+1), N, 3), 3); n figure(1);mesh(real( rrss(:, :, k+N/2+1))) figure(2);tf_show(abs(rrss(:, :, k+N/2+1)))% tf_show(rrss(:, :, k+N/2+1)) drawnowend;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -