clear;
N = 50000;
Fs = 50e3;
[f,y] = a_weighting(N,Fs);
plot(f,y)
fid = fopen('aweighting.inc','w');
for j =1:1:length(y)
fprintf(fid,'%d\t%6.2f\n',f(j),y(j));
end
fid = fopen('aweighting.inc1','w');
for j =1:2:length(y)
fprintf(fid,'%d\t%6.2f\n',f(j),y(j));
end
fid = fopen('aweighting.inc2','w');
for j =1:5:length(y)
fprintf(fid,'%d\t%6.2f\n',f(j),y(j));
end