信号处理IIRFIRFFT算法滤波器

信号处理课程设计 题 1：自相关与FFT 得分 一、 （ 25 分） 已知某正弦波信号混有白噪声信号（正弦波信号频率=学号后两位/100， 其他参数任选）。对上述信号进行自相关和FFT 变换。 1、设定正选信号的频率为79/100=0.79HZ，抽样频率为100HZ； 2、设定N(0,0.25)高斯白噪声，及噪声功率为0.25W； 要求：1.给出Matlab 源程序及说明； 2.给出原始波形、自相关函数图形和FFT 变换后的图形； 3.根据结果分析信号当中混有的主要频率成分； 4.对比说明上述两种信号分析方法各有什么特点。 1.答：源程序及说明 fs=100; fc=0.79;%频率为学号/100 的带有白噪声的正弦信号 x=(0:1/fs:2); n=201; y1=sin(2*pi*fc*x); %原正弦信号，频率为10 a=0;b=0.5; %均值为a，方差为b^2 subplot(3,2,1); plot(x,y1,'r'); title('y=sin(20pi*x)'); ylabel('y'); xlabel('x/20pi'); grid; y2=a+b*rand(1,n); %均匀白噪声 subplot(3,2,2); plot(x,y2,'r'); title('N(0， 0.25)的均匀白噪声'); ylabel('y'); xlabel('x/20pi'); grid; y=y1+y2; %加入噪声之后的信号 subplot(3,2,3); plot(x,y,'r'); title('叠加了均匀白噪声的sinx'); ylabel('y'); xlabel('x/20pi'); grid; FY=fft(y); %傅里叶变换得出频谱函数 FY1=fftshift(FY); %频谱校正 f=(0:200)*fs/n-fs/2; subplot(3,2,4); plot(f,abs(FY1),'r'); title('函数频谱图'); ylabel('F(jw)'); xlabel('w'); grid; 2.答： Butterworth 数字滤波器的模拟滤波器的传递函数： 554329.967*10( )0.51250.13140.02080.0020.0001G ssssss 数字滤波器的转移函数为： 33132333435123450.0598*100.299*10*0.598*10*0.598*10*0.299*10*0.598*10*( )1 3.9845*6.4349*5.2536*0.299*0.0598*zzzzzH zzzzzz 根据以下运行结果： bs2 =[9.9565*10^-4,0,0,0]， as2 =[ 1.0000,0.0946,.0233,.0010]， bz2 =[ 0.0009,0.0027,0.0027,0.0009]， az2 =[ 1.0000, -2.7423,2.5786, -0.8291] Chebyshev 字滤波器的模拟滤波器的传递函数： 4329.9565*10( )0.09460.2330.001G ssss 数字滤波器的转移函数为： 1231230.00090.0027*0.0027*0.0009*( )12.7423*2.5786*0.8291*zzzH zzzz 3 .答：根据第三个图,-1 和 0.8Hz 为主要成分 4 .答：...