MATLAB语言及应用选修课课程设计报告.doc

matlab语言及应用选修课课程设计报告一、问题重述：实作一：信号的产生、显示、谱分析、滤波1、理想信号的时域表示：x=asin（wt），a=444.128，w=100pi，0t1，设f=500hz； 2、产生信号并显示其波形、幅度谱、相位谱；3、在原始信号上加上-3db的白噪声重复上述过程；4、在上述-3db背景下采用matlab中的lms进行滤波，画出滤波后的信号和相应图谱，并对不同lms参数得到的结果进行分析、比对。实作二：语音信号的录制和处理（实作）1、录制一段自己的语音信号，并对录制的信号进行采样（界面上以打开音频文件体现）；2、回放语音信号，并画出采样后的语音信号的时域波形和频谱图；3、给原始语音信号加上-3db的白噪声，画出加噪声后的波形与频谱图；4、设计一带通滤波器对原始信号进行滤波，频带为1000hz至2000hz，并画出滤波器的响应和滤波后信号的时域波形及频谱，并对滤波后的信号进行对比，分析信号的变化；5、设计一个信号处理系统界面；（上述功能最好都能通过gui界面完成）6、按照要求完成课程设计报告。二、问题解答：实作一：源代码如下：% 原 始 信 号 %figure(1);n=0:500; a=444.128; fs=500; w0=100*pi; x=a*sin(w0*n/fs); subplot(3,1,1);plot(x);title(理想信号);fx=(fft(x); subplot(3,1,2);plot(fs*(0:250)/500,abs(fx(1:251);title(幅度谱);angx=angle(fx);subplot(3,1,3);plot(fs*(0:250)/500,angx(1:251) ; title (相位谱)% 加 噪 声 %x=awgn(x,-3,measured);figure(2);subplot(3,1,1);plot(x);title(加噪声后原始信号); fx=(fft(x); subplot(3,1,2);plot(fs*(0:250)/500,abs(fx(1:251);title(加噪声后幅度谱);angx=angle(fx);subplot(3,1,3);plot(fs*(0:250)/500,angx(1:251) ; title (加噪声后相位谱)% 滤 波 %x=x-mean(x);x=x./max(x);m=64;d = x(m:length(x) rand(1,m-1); mu = 0.00001; ha = adaptfilt.lms(m,mu);y,e = filter(ha,x,d); figure(3);subplot(3,1,1);plot(y); title(滤波后原始信号);fx=(fft(y);subplot(3,1,2);plot(fs*(0:250)/500,abs(fx(1:251);title(滤波后幅度谱);angx=angle(fx); subplot(3,1,3);plot(fs*(0:250)/500,angx(1:251) ;title (滤波后相位谱)% 结 束 %画图如下：实作二：源代码如下：% 原 始 信 号 %figure(1);n=0:500; a=500; fs=5000; w0=2000*2*pi; x=a*sin(w0*n/fs); subplot(3,1,1);plot(x);title(理想信号);fx=(fft(x); subplot(3,1,2);plot(fs*(0:250)/500,abs(fx(1:251);title(幅度谱);angx=angle(fx);subplot(3,1,3);plot(fs*(0:250)/500,angx(1:251) ; title (相位谱)% 加 噪 声 %x=awgn(x,-3,measured);figure(2);subplot(3,1,1);plot(x);title(加噪声后原始信号); fx=(fft(x); subplot(3,1,2);plot(fs*(0:250)/500,abs(fx(1:251);title(加噪声后幅度谱);angx=angle(fx);subplot(3,1,3);plot(fs*(0:250)/500,angx(1:251) ; title (加噪声后相位谱)% 滤 波 % x=x-mean(x);x=x./max(x);m=64;d = x(m:length(x) rand(1,m-1); mu = 0.00001; ha = adaptfilt.lms(m,mu);figure(3);n = 4; wn = 2*1000 2000./(fs);ftype = bandpass;b,a = butter(n,wn);y,e = filter(b,a,x);subplot(3,1,1);plot(y); title(滤波后原始信号);fx=(fft(y);subplot(3,1,2);plot(fs*(0:250)/500,abs(fx(1:251);title(滤波后幅度谱);angx=angle(fx); subplot(