信号处理新方法作业分析.ppt

现代信号处理新方法 班级：电子1201班 姓名：刘悦 学号 题目一 用计算机模拟蒲丰投针的实验，并估计随机事件发生的概率以计算圆周率pi。 投掷n次 相交次数为m 平行线的距离为a 针长度为l 相交概率p=m/n=2l/(a* ) =2l/(a*p) matlab实现方法 a=1; %设置两条平行线之间的距离 I=0.6; %投针的长度 m=0; %针与平行线相交的次数 n %投掷次数 x=unifrnd(0,a/2,1,n); %产生n个(0,a/2)之间均匀分布的随机数 phi=unifrnd(0,pi,1,n); %产生n个(0,pi)之间均匀分布的随机数，这里pi是投针到最近的平行线的角度 for i=1:n if x(i)I*sin(phi(i))/2 %只要x小于I*sin(phi(i))/2,则相交 m=m+1; %相交次数加1 end end frequency=m/n; %计算相交的频率，即相交次数比总次数 pi=2*I/(a*frequency) %从相交的频率求总的pi 题目二 FIR数字滤波器设计 窗函数法 滤波器特性图 处理叠加的正弦信号,使信号通过滤波器后留下指定频率 处理加噪声的语音信号，过滤噪声，得到清晰的语音信号 低通滤波器 通带截止频率：200HZ，阻带起始频率：300HZ 采样频率fs=2000HZ 阶数m=66 最小阻带衰减-50db（海明窗） 设计过程 FIR滤波器的窗函数法设计过程： M=66; fp=200; fst=300; fc=250; fs=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); %衰减 subplot(221);stem(hd); title(理想冲击响应); axis([0 M-1 -0.1 0.3]); ylabel(hd[n]); subplot(222);stem(w_ham);title(海明窗); axis([0 M-1 0 1.1]); ylabel(w[n]); subplot(223);stem(h); title(实际冲击响应); axis([0 M-1 -0.1 0.3]); ylabel(h[n]); subplot(224); plot(w/(2*pi),db); title(衰减幅度);ylabel(H[w]); 频率分别为：150HZ，300HZ，700HZ 采样频率fs=2000HZ 过滤语音信号 录制一段语音信号，频率在300HZ~3.4kHZ范围内 噪声是10kHZ的正弦信号 滤波器的截止频率是4kHZ f4=10000;%所加噪声正弦函数的频率 [Y,fs,bits]=wavread(F:\b.wav);%利用wavread产生音频的函数及采样频率 L=length(Y); t=0:1/fs:(L-1)/fs;%定义时间的范围及步长 f=0:fs/L:fs*(L-1)/L; subplot(211); plot(t,Y);xlabel(时间(t));ylabel(幅度(Y)); title(加噪