N=1024;
n=0:N-1;
fs=972800000;
t=n/fs;
sin = sin(2*pi*950000*t); %生成sin波形
sin_q =round(sin /max(abs(sin))*8192);
cos = cos(2*pi*950000*t);%生成sin波形
cos_i =round(cos /max(abs(cos))*8192);
fid = fopen('d:/sin_q .txt', 'wt');
fprintf(fid, '%16.0f\n', sin_q);
fclose(fid);
*****************************************************************************************************************
fs = 100000000; %100M 采样率
t = 0:1/fs:0.00000999; % 共1000个点
f1 = 100000; % 信号的频率是100k
signal1 = sin(2*pi*f1*t); %生成sin波形
plot(signal1, 'b'); %观察波形
% 由于正、余弦波形的值在[0,1]之间,需要量化成16bit,先将数值放大
y1 = signal1 * 32768; % 32 * 1024 = 32768
fid = fopen('d:/matlab_wave_data.txt', 'wt');
fprintf(fid, '%x\n', y1);
fclose(fid);
plot(y1, 'b');
****************************************************************************************************************************************
N=1024;
n=0:N-1;
cos_i=1000*cos(0.002*pi*n); %频率为100m左右取值1024个,多于一个周期
cos_i=round(cos_i/max(abs(cos_i))*32767); % abs() 求绝对值,max() 求最大值,round() 四舍五入
fid = fopen('d:/cos.txt','wt'); %在d盘目录下新建文件
fprintf(fid,'%16.0f\n',cos_i);%将滤波器系数写入文件,以16位浮点数格式保存。16表示16位,0表示没有小数,f表示浮点。
fclose(fid);%文件操作结束,否则,该文件处于只读状态,不能更改和删除
----------------------------------
N=1024;
n=0:N-1;
sin_q = sin(0.002*pi*n); %频率为100m
sin_q =round(sin_q /max(abs(sin_q ))*32767); % abs() 求绝对值,max() 求最大值,round() 四舍五入
fid = fopen('d:/sin.txt','wt'); %在d盘目录下新建文件
fprintf(fid,'%16.0f\n',sin_q);%将滤波器系数写入文件,以16位浮点数格式保存。16表示16位,0表示没有小数,f表示浮点。
fclose(fid);%文件操作结束,否则,该文件处于只读状态,不能更改和删除
**********************************************************************************************************************************************
打开notpad++,ctrl+f,选择替换,原文本: \n 或者 \r 。 替换为: ,
**********************************************************************************************************************************************
fs=100;
N=128; %采样频率和数据点数
n=0:N-1;
t=n/fs; %时间序列
x=0.5*sin(2*pi*15*t);%信号
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y); %求得Fourier变换后的振幅
f=n*fs/N; %频率序列
subplot(2,2,1),plot(f,mag); %绘出随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');
title('N=128');
grid on;
subplot(2,2,2);
plot(x);%(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('时间');
ylabel('振幅');
title('N=128');
grid on;