基于matcom的仿真信号采集与频谱分析

1，核心代码

Mm y;//采样信号Mm t;//时间序列Mm n;//采样序列Mm f;//频率序列int N;//采样点数double fo;//信号频率double fs;//采样频率fo=3;//3Hzfs=2.56*fo;N=1024;n=linspace(0,N-1,N);f=n*fs/N;t=n/fs;y=sin(2*pi*fo*t)+0.5*cos(2*pi*fo/3*t); //两种信号频率分别为fo,fo/3Mm m_line=plot((CL(f),abs(fft(y))*2.0/N));//FFT变换并绘图title((CL(TM("fft"))));xlabel((CL(TM("f(Hz)"))));ylabel((CL(TM("y(A)"))));set(h,(CL(TM("Box")),TM("on")));set(h,(CL(TM("Color")),TM("white")));set(m_line,(CL(TM("Color")),TM("red")));grid();drawnow();

2，运行结果