|
dan@0
|
1 function [df3,A_sys,b_sys] = ddm_fft_3(win, win_der, mf_ders, sig, N, N_fft, fs)
|
|
dan@0
|
2 R = 3;
|
|
dan@0
|
3 Q = R; % nr of unknows is the same as the number of model funcs
|
|
dan@0
|
4 [A_sys b_sys] = ddm_lin_sys_fft(Q, R, win, win_der, mf_ders, sig, N, N_fft, fs);
|
|
dan@0
|
5 % A_sys1 = cat(2,shiftdim(A(1:N_fft-2,1),-2),shiftdim(A(1:N_fft-2,2),-2),shiftdim(A(1:N_fft-2,3),-2));
|
|
dan@0
|
6 % A_sys2 = cat(2,shiftdim(A(2:N_fft-1,1),-2),shiftdim(A(2:N_fft-1,2),-2),shiftdim(A(2:N_fft-1,3),-2));
|
|
dan@0
|
7 % A_sys3 = cat(2,shiftdim(A(3:N_fft ,1),-2),shiftdim(A(3:N_fft ,2),-2),shiftdim(A(3:N_fft ,3),-2));
|
|
dan@0
|
8 % A_sys = cat(1,A_sys1,A_sys2,A_sys3);
|
|
dan@0
|
9 % b_sys = cat(1,shiftdim(b(1:N_fft-2).',-1),shiftdim(b(2:N_fft-1).', -1),shiftdim(b(3:N_fft ).', -1));
|
|
dan@0
|
10 df3 = lin_solve_dgr_3(A_sys,b_sys,N_fft-2); %hardcoded degree 3 solver (fast)
|
|
dan@0
|
11
|
|
dan@0
|
12 end |
