signaltest.m

this base on DSP a iir design

%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('信号的幅度谱（离散频谱）');