lijinjia.m

来自「这是一个用MATLAB语言编写的程序」· M 代码 · 共 38 行

clear all 
clc

nd = 1024;
dec = floor(randn(1,nd)*20);
n = 1:nd;

%*****************理论上的fft输出***************
dec_fft = fft(dec);
dec_fft_real = real(dec_fft);
dec_fft_imag = imag(dec_fft);
figure(1)
stem(n,dec_fft_real,'*');
title('实部')
figure(2)
stem(n,dec_fft_imag,'o');
title('虚部')

%********************将DEC转换为对应的有符号的十进制数******

for k = 1:nd
    if dec(k)>=0
        a((k-1)*8+1:8*k) = dec2bin(dec(k),8);
        %fprintf('%d\n',a((k-1)*8+1) );
    else
        dec_xin = -1*dec(k);
        a((k-1)*8+1:8*k) = dec2bin(dec_xin,8);
        a((k-1)*8+1:8*k) = a((k-1)*8+1:8*k)+[1 0 0 0 0 0 0 0];
        %fprintf('%d%d%d%d%d%d%d%d\n',a((k-1)*8+1:8*k));
    end
end
a2 = [1:1024]';
a1 = (reshape(a,8,nd))'
a11 = [a2 a1];
%for k = 1:nd
%    fprintf('%s%s%s%s%s%e%f%d\n',a1(k,:));
%end