哪位好人 能把下面的Matlab程序转化成C++或者C语言啊 毕业求教  坐等




N = (M*Fs)/df;
N = pow2 (nextpow2(N)); % windowing record for FFT

X = fft(x,N);   % fft of the truncated (or zero padded) time series
X = fftshift(X);% shift components of fft
Xc = conj(X);                 % precompute the complex conjugate vector

S = zeros (N,N);              % size of the Spectral Correlation Density matrix
f = zeros (N,N);              % size of the frequency matrix;
alfa = zeros (N,N);           % size of the cycle frequency matrix
F = Fs/(2*N);                 % precompute constants -  F = Fs/(2*N);     
G = Fs/N;                     % precompute constants -  G = Fs/N;  
m = -M/2+1:M/2;               % set frequency smoothing window index

for k = 1:N                                % fix k
    % computes vectors of f and alfa,
    % store frequency and cycle frequency data for given k.
   
    k1 = 1:N;
    f(k,k1) = F*(k+k1-1) - Fs/2;          % Computes f values and shift them to center in zero (f = (K+L)/2N) [1]
    alfa(k,k1) = G*(k-k1 + N-1) - Fs;       % Computes alfa values and shift them to center in zero (alfa = (K-L)/N) [1]
   
    for k1 = 1:N %fix k1 = J
        %calculate X(K+m) & conj (X(J+m)) for arguments of X(1:N) only
        B = max(1-k, 1-k1);          % Largest min of 1 <= (K+m)| (J+m) <= N
        A = min (N-k, N-k1);         % Smallest max of 1 <= (K+m)| (J+m) <= N
        n = m((m<=A) & (m>=B));   %fix the index out of range problem by
                                                   % truncating the window
        if isempty(n)
            S(k,k1) = 0;
        else
            p = k+n; q = k1+n;
            Y = X(p).*Xc(q);
            S(k,k1) = sum(Y);
        end
    end
end

其中的参数均假设已经给定 希望大家给点建议啊~~~~~~~~~~~~~

大神 你懂matlab吧

额 这个我用matlab自己编的好吧 只是现在导师要转化成C++ 时间催的紧好吧