信号作业1.doc

数据采集与信号处理 1).已知信号x(n)=0.0*COS(2*3.14*SF*n/FS)+50.0*COS(4*3.14*SF*n/FS) ，n=0,1,2 ……N-1 式中： N-已知信号x(n)=0.0*COS(2*3.14*SF*n/FS) +50.0*COS(4*3.14*SF*n/FS) ，n=0,1,2 ……N-1 式中： N---采样点数，必须为2x SF---信号频率SF=10 Hz FS---采样频率 其FFT变换结果X(k)可用下面提供的FFT子程序求出，计算功率谱的公式为： G(k)=2(XR(k)2+XI(k)2)/N, k=0,1,2 ……N/2-1 式中： XR(k)--- X(k)的实部 XI(k)--- X(k)的虚部 程序如下： FS=300; SF=10; N=1024; n=0:N-1; t=n/FS; x=100.0*cos(2*3.14*SF*t)+50.0*cos(4*3.14*SF*t); subplot(211); plot(t,x); xlabel(t); ylabel(y); title(x=100.0*cos(2*3.14*SF*t)+50.0*cos(4*3.14*SF*t)时域波形); grid; y=fft(x,N); mag=abs(y); f=(0:leng