%產生一個函數的128個q15數據并把它寫在一個x.h的文件里面供dsp讀取用,qformart表示q值
t=0:0.001:0.128
f=20;
qformat=15;
y1=0.5*sin(2*pi*f*t)+0.5*sin(2*pi*100*t);%y1=0.5*sin(2*pi*f*t).*sin(2*pi*300*t)是自己想要產生q格式的函數,根據自己情況改寫
y=floor(y1*(2^qformat));
subplot(211);
plot(t,y);
bsize=length(y1);  %求yi的長度

for i=1:bsize                               % 數據處理防止它超過規定的范圍
   if y(i)==32768
      y(i)=32767;
   end
end
fid = fopen('x.h','w');%把產生的數據寫入指定的文件
fprintf(fid,'int x[]={');
for i=1:128
  fprintf(fid,'\t');
  fprintf(fid,'%d,',y(i));    
if(mod(i,10)==0)
      fprintf(fid,'\\');
      fprintf(fid,'\n');
end
end
fseek(fid,-1,0);
fprintf(fid,'};\n');
fclose(fid);
b=fft(y1,1024);
m=abs(b);
n=1000/1024*(0:511);
subplot(212);
plot(n,m(1:512)); title('magnitude spectrum');
xlabel('frequency (hz)');