matlab生成正弦波，fft观测频谱

L=128;
f1=100;
fs=400;
t=0:1/fs:(1/fs)*(L-1);
si=sin(2*pi*f1*t);
f=fft(si,L);
f=20*log(abs(f))/log(10);
ft=[0:(fs/L):fs/2];
f=f(1:length(ft));
subplot(3,2,1);plot(t,si);
xlabel('时间（s）');
ylabel('幅度（v）');title('时域信号波形');%%%%%%
subplot(3,2,2);plot(ft,f);
xlabel('频率（hz）');
ylabel('功率（dbw）');title('信号频谱图');

L=128;
f1=100;
fs=12800;
t=0:1/fs:(1/fs)*(L-1);
si=sin(2*pi*f1*t);
f=fft(si,L);
f=20*log(abs(f))/log(10);
ft=[0:(fs/L):fs/2];
f=f(1:length(ft));
subplot(3,2,3);plot(t,si);
xlabel('时间（s）');
ylabel('幅度（v）');title('时域信号波形');%%%%%%
 subplot(3,2,4);plot(ft,f);
xlabel('频率（hz）');
ylabel('功率（dbw）');title('信号频谱图');

 

L=10240;
f1=1000000;
fs=1024000000;
t=0:1/fs:(1/fs)*(L-1);
si=sin(2*pi*f1*t);
f=fft(si,L);
f=20*log(abs(f))/log(10);
ft=[0:(fs/L):fs/2];
f=f(1:length(ft));
subplot(3,2,5);plot(t,si);
xlabel('时间（s）');
ylabel('幅度（v）');title('时域信号波形');%%%%%%
subplot(3,2,6);plot(ft,f);
xlabel('频率（hz）');
ylabel('功率（dbw）');title('信号频谱图');

sin_q =round(si /max(abs(si))*8192);

fid = fopen('d:/sin_q .txt', 'wt');
fprintf(fid, '%16.0f\n', sin_q);
fclose(fid);