5-5.m

《DSP嵌入式系统开发典型案例》此书的源代码

clear all;                    	     				   	%清寄存器值
clf;                              				   	%清屏
N=256;                            			 	   	%数据点数
fs=10000;                            			        	%采样频率
dt=1/fs;                                		    		%采样时间间隔
tt=0.:dt:dt*(N-1);
dis=(rand(size(tt)));                              		    	%随机噪声
for k=1:N;                    
   f1=3000;                         			 		%信号1频率
   f2=1000;                         			 		%信号2频率
   y(k)=sin(2*pi*f1*k*dt)+sin(2*pi*f2*k*dt)+dis;  %产生合成信号
end
lp=3500;                                   	 		%截止频率
wn1=2*lp/fs;                              	 			%函数的参数
[z1,p1,k1] = CHEBY1(3,0.5,wn1);            	 			%滤波器的极零点表示
[b1,a1] = CHEBY1(3,0.5,wn1);              	 			%滤波器的传递函数表示
b1=b1/(8*1.0711);                    		 			%将参数按比例缩小
a1=a1/(8*1.0711);                             		%为DSP程序做准备，MATLAB中不需要

yy1=filter(b1,a1,y);                   	      			%滤波

y=fft(y,N);                         	      			%将原始信号做FFT变换
pyy=y.*conj(y);                     	      			%做功率谱分析
f=(1:(N/2－1));
for i=1:N/2-1;
   f(i)=f(i)*fs/N;
end
figure(1);

plot(f,pyy(2:N/2))                              			%画出原始信号的频谱图

y=fft(yy1,N);                      	       			%将滤波后数据做功率谱分析
pyy=y.*conj(y);
f=(1:(N/2－1));
for i=1:N/2－1;
   f(i)=f(i)*fs/N;
end
figure(2);
plot(f,pyy(2:N/2));                               			%画出滤波后信号的频谱图