?? param_tconv.m
字號:
function Bml= param_tconv(Bml1, Bml2, N, beta)% function Bml= param_tconv(Bml1, Bml2, N, beta)% This file is part of the TFPM toolbox v1.0 (c)% michael.jachan@tuwien.ac.at and underlies the GPL.% % Computes the twisted convolution Bml= Bml2#Bml1 of TFMA% parameters. The systems Bml1, Bml2 are causal! No factor N!% Noncommutative! %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;tfpm;MAR = 1;LAR = 1;MMA = 1;LMA = 1;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;% Compute Bml and B1Bml= [];for m= 0:M b= []; for l= -L:L bml= 0; for mp= 0:M2 for lp= -L2:L2 bml= bml + param_get(Bml2, mp , lp ) * ... param_get(Bml1, m-mp, l-lp) * ... exp(-j*2*pi/N*tfpm_phi(m, l, mp, lp, beta)); end; end; b= [b; bml]; end; Bml= [Bml b];end;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% TEST:e= randn(N, 1);u = tfarma_gen(e, 1, Bml1, beta);y = tfarma_gen(u, 1, Bml2, beta);y0= tfarma_gen(e, 1, Bml , beta);norm(y-y0)/norm(y)figure(1);clf;hold onplot(real(y), 'r')plot(real(y0), 'g')plot(real(y0-y), 'b')%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -