信号处理新方法作业

现代信号处理新方法

班级：电子1201班

姓名：刘悦

学号：20123662TheCollegeofInformationScienceandEngineering题目一模拟随机试验计算圆周率pi题目二窗函数法设计FIR数字滤波器题目四数字滤波器过滤高斯白噪声

主要内容TheCollegeofInformationScienceandEngineering题目一用计算机模拟蒲丰投针的实验，并估计随机事件发生的概率以计算圆周率pi。投掷n次相交次数为m平行线的距离为a针长度为l相交概率p=m/n=2l/(a*)=2l/(a*p)TheCollegeofInformationScienceandEngineeringTheCollegeofInformationScienceandEngineeringmatlab实现方法

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)<I*sin(phi(i))/2

%只要x小于I*sin(phi(i))/2,则相交m=m+1;

%相交次数加1endendfrequency=m/n;

%计算相交的频率，即相交次数比总次数pi=2*I/(a*frequency)

%从相交的频率求总的piTheCollegeofInformationScienceandEngineering试验次数：10000000试验次数：1000000试验次数：100000TheCollegeofInformationScienceandEngineering窗函数法滤波器特性图处理叠加的正弦信号,使信号通过滤波器后留下指定频率处理加噪声的语音信号，过滤噪声，得到清晰的语音信号题目二FIR数字滤波器设计

TheCollegeofInformationScienceandEngineering低通滤波器通带截止频率：200HZ，阻带起始频率：300HZ采样频率fs=2000HZ阶数m=66最小阻带衰减-50db（海明窗）滤波器指标TheCollegeofInformationScienceandEngineering设计过程FIR滤波器的窗函数法设计过程：是逼近的理想滤波器频率响应，是理想滤波器的单位脉冲响应TheCollegeofInformationScienceandEngineeringM=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);

%衰减TheCollegeofInformationScienceandEngineeringsubplot(221);stem(hd);title('理想冲击响应');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]');TheCollegeofInformationScienceandEngineeringTheCollegeofInformationScienceandEngineering频率分别为：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);

%滤波前的信号图像xlabel('时间/秒');ylabel('幅度');title('信号滤波前时域图');subplot(212);Fs=fft(s,512);

%将信号变换到频域AFs=abs(Fs);

%信号频域图的幅值f=(0:255)*fs/512;

%频率采样plot(f,AFs(1:256));

%滤波前的信号频域图xlabel('频率/赫兹');ylabel('幅度');title('信号滤波前频域图');TheCollegeofInformationScienceandEngineering即将被过滤频率的信号TheCollegeofInformationScienceandEngineeringsf=filter(h,1,s);%使用filter函数对信号进行滤波subplot(211);plot(t,sf);%滤波后的信号图像xlabel('时间/秒');ylabel('幅度');title('信号滤波后时域图');axis([0.20.25-22]);%限定图像坐标范围subplot(212);Fsf=fft(sf,512);%滤波后的信号频域图AFsf=abs(Fsf);%信号频域图的幅值f=(0:255)*fs/512;%频率采样plot(f,AFsf(1:256))%滤波后的信号频域图xlabel('频率/赫兹');ylabel('幅度');title('信号滤波后频域图');TheCollegeofInformationScienceandEngineering保留150HZ频率信号TheCollegeofInformationScienceandEngineering过滤语音信号录制一段语音信号，频率在300HZ~3.4kHZ范围内噪声是10kHZ的正弦信号滤波器的截止频率是4kHZTheCollegeofInformationScienceandEngineeringf4=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('加噪前音频语音波形的时域图');y=0.05*sin(2*pi*f4*t);Y=Y(:,1);%取Y矩阵的第一列Y1=y+Y';%给音频加噪声subplot(212);plot(t,Y1);xlabel('时间(t)');ylabel('幅度(Y)');title('加噪后音频语音波形的时域图');TheCollegeofInformationScienceandEngineeringFY=abs(fft(Y,L));%原始音频信号的傅里叶变换subplot(211)plot(f,FY);gridon;%加噪前音频的频域图xlabel('频率(f)');ylabel('幅度(FY)');title('加噪前音频波形的频域图');FY1=abs(fft(Y1,L));%加噪声的傅里叶变换subplot(212)plot(f,FY1);gridon;%加噪后音频的频域图xlabel('频率(f)');ylabel('幅度(FY1)');title('加噪声后音频波形的频域图');TheCollegeofInformationScienceandEngineeringTheCollegeofInformationScienceandEngineeringTheCollegeofInformationScienceandEngineeringM=66;fc=2000;fst=3000;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],1000,'whole');%求频谱mag=abs(H);%求实部db=20*log(mag);%衰减TheCollegeofInformationScienceandEngineeringsf=filter(h,1,Y1);%使用filter函数对信号进行滤波subplot(211);plot(t,sf);%滤波后的信号图像xlabel('时间/秒');ylabel('幅度');title('信号滤波后时域图');subplot(212);Fsf=fft(sf,L);%滤波后的信号频域图AFsf=abs(Fsf);%信号频域图的幅值plot(f,AFsf)%滤波后的信号频域图xlabel('频率/赫兹');ylabel('幅度');title('信号滤波后频域图');sound(sf,fs);wavwrite(sf,fs,'F:\hao');wavwrite(Y1,fs,'F:\noise');TheCollegeofInformationScienceandEngineeringTheCollegeofInformationScienceandEngineering题目四编写MATLAB程序，计算在的加性高斯白噪声环境下，幅度为0.3，频率为10Hz的正弦信号通过50Hz低通滤波器的情况。

编写功率谱为的加性高斯白噪声叠加到正弦信号上加噪信号通过滤波器TheCollegeofInformationScienceandEngineering实现高斯白噪声环境下的正弦信号fs=1000;fc=10;t=(0:1/fs:1);y1=0.3*sin(2*pi*fc*t);y2=wgn(1,length(y1),5*10^-6,'dBm');%生成高斯白噪声y=y1+y2;figure(1);subplot(211);plot(y2);subplot(212);plot(y);gridon;TheCollegeofInformationScienceandEngineeringTheCollegeofInformationScienceandEngineering50HZ低通滤波器滤波