信号处理新方法作业.ppt

《信号处理新方法作业.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、现代信号处理新方法班级：电子1201班姓名：刘悦学号：20123662TheCollegeofInformationScienceandEngineering题目一模拟随机试验计算圆周率pi题目二窗函数法设计FIR数字滤波器题目四数字滤波器过滤高斯白噪声主要内容TheCollegeofInformationScienceandEngineering题目一用计算机模拟蒲丰投针的实验，并估计随机事件发生的概率以计算圆周率pi。投掷n次相交次数为m平行线的距离为a针长度为l相交概率p=m/n=2l/(a*)=2l/(a*p)TheCollegeofInformatio

2、nScienceandEngineeringTheCollegeofInformationScienceandEngineeringmatlab实现方法a=1;%设置两条平行线之间的距离I=0.6;%投针的长度m=0;%针与平行线相交的次数n=10000000;%投掷次数x=unifrnd(0,a/2,1,n);%产生n个(0,a/2)之间均匀分布的随机数phi=unifrnd(0,pi,1,n);%产生n个(0,pi)之间均匀分布的随机数，这里pi是投针到最近的平行线的角度fori=1:nifx(i)

3、i(i))/2,则相交m=m+1;%相交次数加1endendfrequency=m/n;%计算相交的频率，即相交次数比总次数pi=2*I/(a*frequency)%从相交的频率求总的piTheCollegeofInformationScienceandEngineering试验次数：10000000试验次数：1000000试验次数：100000TheCollegeofInformationScienceandEngineering窗函数法滤波器特性图处理叠加的正弦信号,使信号通过滤波器后留下指定频率处理加噪声的语音信号，过滤噪声，得到清晰的语音信号题目二FIR

4、数字滤波器设计TheCollegeofInformationScienceandEngineering低通滤波器通带截止频率：200HZ，阻带起始频率：300HZ采样频率fs=2000HZ阶数m=66最小阻带衰减-50db（海明窗）滤波器指标TheCollegeofInformationScienceandEngineering设计过程FIR滤波器的窗函数法设计过程：是逼近的理想滤波器频率响应，是理想滤波器的单位脉冲响应TheCollegeofInformationScienceandEngineeringM=66;fp=200;fst=300;fc=250;f

5、s=2000;wc=2*pi*fc/fs;%将模拟滤波器的技术指标转换为数字滤波器技术指标alpha=(M-1)/2;n=[0:1:(M-1)];m=n-alpha;hd=[sin(wc*m)]./(pi*m);%hd为理想冲激相应w_ham=(hamming(M))';%汉明窗h=hd.*w_ham;%加窗[H,w]=freqz(h,1);%求频谱mag=abs(H);%求实部db=20*log(mag);%衰减TheCollegeofInformationScienceandEngineeringsubplot(221);stem(hd);title('理想

6、冲击响应');axis([0M-1-0.10.3]);ylabel('hd[n]');subplot(222);stem(w_ham);title('海明窗');axis([0M-101.1]);ylabel('w[n]');subplot(223);stem(h);title('实际冲击响应');axis([0M-1-0.10.3]);ylabel('h[n]');subplot(224);plot(w/(2*pi),db);title('衰减幅度');ylabel('H[w]');TheCollegeofInformationScienceandEngine

7、eringTheCollegeofInformationScienceandEngineering频率分别为：150HZ，300HZ，700HZ采样频率fs=2000HZ正弦信号TheCollegeofInformationScienceandEngineeringf1=150;f2=300;f3=700;%待滤波正弦信号频率fs=2000;%采样频率subplot(211);t=0:1/fs:0.25;%定义时间范围和步长s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);%滤波前信号plot(t,s);%滤波前的信号

8、图像xlabel('时间