fir.m

来自「基于matlab的FIR滤波器实验的源程序及实验结果图形」· M 代码 · 共 29 行

f11=50;    %/Hz
f12=350;   %/Hz
fs=1000;   %/采样Hz
N=100      %数据个数
T=1/fs;    %采样周期
n=0:N;

x11=sin(2*pi*f11*n*T);
x12=0.7*sin(2*pi*f12*n*T);
x_base=(x11+x12);

%待滤波信号波形
figure(1)
plot(x_base)

%待滤波信号频谱
figure(2)
yff=abs(fft(x_base))
df=n*(fs/N)
plot(df,yff)

xout=x_base/max(x_base);%归一化

xto_ccs=round(32768*xout)

fid=fopen('input.dat','w');%打开文件
fprintf(fid,'1651 1 0 0 0\n');%输出文件头
fprintf(fid,'%d\n',xto_ccs);%输出
fclose(fid);