%%采样频率为256Hz，采样点数为4096个点；
%信号含有2V的直流分量， 频率为50Hz、相位为-30度、幅度为3V的交流信号， 频率为75Hz、相位为90度、幅度为1.5V的交流信号。 
%用数学表达式就是如下： S=2+3*cos(2*pi*50*t-pi*30/180)+1.5*cos(2*pi*75*t+pi*90/180)  
%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%对信号采样数据为4096点的处理
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;close all;
clear all;close;
fs=256;N=4096; %采样频率和数据点数
n=0:N-1;
t=n/fs; %时间序列
signal = 2 + 3*cos(2*pi*50*t-pi*30/180) + 1.5*cos(2*pi*75*t+pi*90/180); %信号
figure;
plot(signal);
y=fft(signal,N); %对信号进行快速Fourier变换
mag=abs(y); %求得Fourier变换后的振幅
f=n*fs/N; %频率序列
figure;
subplot(2,1,1),
plot(f,mag); %绘出随频率变化的振幅
xlabel('频率/Hz');  %%频率为：Fn=(n-1)*Fs/N
ylabel('振幅');
title('N=4096');
subplot(2,1,2);
plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('N=4096');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%将信号写入txt文件
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
image_data = signal*0;
for x = 1:N
    signal(x)= round(signal(x));
    image_data(x) = round(image_data(x));
end
fidr = fopen('sin_process_real_input.txt','wt');%将信号进行打印输??
fidi = fopen('sin_process_image_input.txt','wt');%将信号进行打印输??
fprintf(fidr,'%d\n',signal);
fprintf(fidi,'%d\n',image_data);
fclose(fidr);
fclose(fidi);

clc;close all;
clear all;close;
N = 1024;
sin_fidro = fopen('ip_fft_real_output_ver.txt','r');
sin_fidio = fopen('ip_fft_imag_output_ver.txt','r');
sin_fideo = fopen('ip_fft_exponent_output_ver.txt','r');
sin_yreal_out = fscanf(sin_fidro,'%d');
sin_yimag_out = fscanf(sin_fidio,'%d');
sin_yexp_out = fscanf(sin_fideo,'%d');
fclose(sin_fidro);
fclose(sin_fidio);
fclose(sin_fideo);
figure;
plot(sin_yreal_out);
title('sin_yreal_out');
figure;
plot(sin_yimag_out);
title('sin_yimag_out');
figure;
plot(sin_yexp_out);
title('sin_yexp_out');
for i=1:N*4-1
sin_real_out(i) = sin_yreal_out(i)/(2^(sin_yexp_out(i)));%%实部和缩放因子的处理
sin_imag_out(i) = sin_yimag_out(i)/(2^(sin_yexp_out(i)));%%虚部和缩放因子的处理
sin_mod_out(i) = sqrt(sin_real_out(i)^2 + sin_imag_out(i)^2);%计算幅值
end
figure;
plot(abs(sin_mod_out));
title('sin mod out');
figure;
plot(sin_real_out);
title('sin real out');
figure;
plot(sin_imag_out);
title('sin imag out');