%writen by Ziwei Chen 
clc;close all;clear all;
f1=500;f2=2000;
%計算采樣周期Ts及采樣點數N,對原時間信號進行截斷，選取Tp時間進行觀察
fs=10000;Ts=1/fs;
N=128;                                 %N=fs/F
Tp=N*Ts                               %最小記錄時間（對信號的觀察時間）
fo=fopen('testdata.dat','w');
n=0:N-1;
%對模擬信號進行采樣，得到序列信號x(n)
xn=5*sin(2*pi*f1*n*Ts)+2.5*cos(2*pi*f2*n*Ts);
xn=xn.*4096;
fprintf(fo,'  .word %d\r\n',fix(xn)); 
fclose(fo);
%對序列信號x(n)進行譜分析，得到信號的頻譜
h1=figure;
k=[0:N-1];
XK=fft(xn,N);
set(h1,'name','頻譜分析的結果（頻譜信號）');
plot(fs*k/N,abs(XK));grid on;
set(gca,'XTickMode','manual','XTick',[0,500,1000,1500,2000]);
xlabel('頻率f  (單位：Hz)');ylabel('幅度|Z(j2\pif)|');
title('信號的幅度譜（離散頻譜）');