ms30_1.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 20 行

b=[0.0114747,0,-0.034424,0,0.034424,0,-0.0114747];
a=[1,0,-2.13779,0,-1.76935,0,-0.539758 ]; 
N=length(b)-1;
bx=abs(b); ax=abs(a);       %去掉符号
bmax=max(bx);amax=max(ax);
maxba=max(bmax,amax);       %求b和a中绝对值最大的数
b1=round(b./maxba*32767);   %将b系数归一化后转换为Q0格式
a1=round(a./maxba*32767);   %将a系数归一化后转换为Q0格式
%在当前目录下打开由文件名指定的文件,准备进行写操作
fid=fopen('iirdata.inc','w'); 
	fprintf(fid,'N    .set  %d\n',N);
	fprintf(fid,'table:    .sect    "table"\n');
	for i=0:N
   	fprintf(fid,'    .word		%d\n',b1(i+1));
	end   
	for i=0:N
   	fprintf(fid,'    .word		%d\n',a1(i+1));
	end
fclose(fid)                 %关闭文件