EEPW论坛

%显示原始信号
plot(S);
title('原始信号');

%显示原始的FFT模值结果

figure;
Y = fft(S,N);            %做FFT变换
Ayy = (abs(Y));          %abs--求绝对值----取模
plot(Ayy(1:N));          %显示原始的FFT模值结果
title('FFT 模值');

%幅度-频率曲线图-显示换算后的FFT模值结果--实际值

figure;
Ayy = Ayy / (N / 2);                    %换算成实际的幅度
%Ayy(1) = Ayy(1) / 2;                   %直流分量的幅值
F = ([1:N] - 1) * Fs / N;               %换算成实际的频率值
plot(F(1:N/2), Ayy(1:N/2));             %显示换算后的FFT模值结果
title('幅度-频率曲线图');

%相位-频率曲线图

figure;
Pyy = [1 : N/2];
for i = 1 : N/2
Pyy(i) = phase(Y(i));                %计算相位
Pyy(i) = Pyy(i) * 180 / pi;          %换算为角度
end;
plot(F(1 : N/2), Pyy(1 : N/2));      %显示相位图
title('相位-频率曲线图');

%数据导出部分
xout=S/max(S);                                 %归一化--(x-min)/(max-min)
xto_ccs=round(32767*xout)                      %Q15格式

fid=fopen('input_moudle.dat','w');             %创建input文件，写操作
fprintf(fid,'1651 1 0 0 0\n');                 %输出文件头
fprintf(fid,'%d\n',xto_ccs);                   %格式化输入数据到文件
fclose(fid);

                                             编译后用MATLAB看到的图形

                                                              生成的波形

                                                                  频谱图

                             fopen 格式参数说明