数字信号处理课程设计毕业设计word格式内容摘要:

%输出连续谱图 grid on title(39。 原始音乐信号频域波型 39。 ) xlabel(39。 Frequency/Hz39。 ) ylabel(39。 Magnitude39。 ) n=0:N1。 x1=cos(n*pi*)。 %进行调制的余弦信号 figure subplot(2,2,1)。 plot(x1) %输出余弦信号的连续波形 title(39。 余弦信号时域波型 39。 ) xlabel(39。 Time39。 ) ylabel(39。 Magnitude39。 ) F2=fft(x1,N)。 %对余弦信号进行 N点 fft变换 subplot(2,2,2)。 plot(w,abs(F2)) %输出余弦信号的连续谱图 grid on title(39。 余弦信号频域波型 39。 ) xlabel(39。 Frequency/Hz39。 ) ylabel(39。 Magnitude39。 ) X=x.*x139。 subplot(2,2,3)。 plot(X) %输出调制后音乐信号的连续波形 F3=fft(X,N)。 grid on title(39。 高频调制后音乐信号时域波型 39。 ) xlabel(39。 Time39。 ) ylabel(39。 Magnitude39。 ) subplot(2,2,4)。 plot(w,abs(F3)) %输出调制后音乐信号的连续谱图 grid on title(39。 高频调制后音乐信号频域波型 39。 ) xlabel(39。 Frequency/Hz39。 ) ylabel(39。 Magnitude39。 ) X1=X.*x139。 %对调制后的音乐信号进行解调 F4=fft(X1,N)。 %对解调后的音乐信号进行 N点 fft变换 figure subplot(2,1,1)。 plot(X1)。 %输出解调后的音乐信号的连续波形 grid on title(39。 解制后音乐信号时域波型 39。 ) xlabel(39。 Time39。 ) ylabel(39。 Magnitude39。 ) subplot(2,1,2)。 plot(w,abs(F4))。 %输出解调后音乐信号的连续谱图 grid on title(39。 解制后音乐信号频域波型 39。 ) xlabel(39。 Frequency/Hz39。 ) ylabel(39。 Magnitude39。 ) Rp=1。 Rs=15。 %Rp 为通带内最大波 纹, Rs 为阻带内最小衰减 [N1,Wc]=buttord(,Rp,Rs)。 %N1为滤波器的阶数, Wc为巴特沃斯滤波器的截止频率 [B,A]=butter(N1,Wc)。 %B为巴特沃斯滤波器的分子, A为巴特沃斯滤波器的分母 [H,W]=freqz(B,A)。 %H为数字滤波器的频率响应 figure subplot(1,1,1)。 plot(W/pi,abs(H))。 %输出滤波器的频率响应图 grid on title(39。 数字巴特沃斯滤波器 39。 ) xlabel(39。 Digital Frequency/pi39。 ) ylabel(39。 Magnitude39。 ) y=filter(B,A,X1)。 %对解调后的音乐信号进行滤波 % sound(y,fs)。 %播放滤波后的音乐信号 figure subplot(2,1,1) plot(y)。 %输出用巴特 沃斯滤波后的连续波形 grid on title(39。 巴特沃斯滤波后音乐信号时域波型 39。 ) xlabel(39。 Time39。 ) ylabel(39。 Magnitude39。 ) y=fft(y,N)。 subplot(2,1,2) plot(w,abs(y))。 %输出巴特沃斯滤波后的连续谱图 grid on title(39。 滤波后音乐信号频域波型 39。 ) xlabel(39。 Frequency/Hz39。 ) ylabel(39。 Magnitude39。 ) N2=33。 wc=*pi。 %N为理想低通滤波器的阶数, Wc为滤波器的截止频率 hd=ideal(N2,wc)。 %调用理想低通滤波器函数 w2=boxcar(N2)。 %产生矩形窗函数 w3=blackman(N2)。 %产生布莱克曼窗函数 h1=hd.*w239。 %加窗设计矩形窗滤波器 h2=hd.*w339。 %加窗设计布莱克曼窗滤波器 M=length(x)。 fh1=fft(h1,M)。 db1=20*log10(abs(fh1(1)./(abs(fh1)+eps)))。 %矩形窗的幅度分贝表示 fh2=fft(h2,M)。 db2=20*log10(abs(fh2(1)./(abs(fh2)+eps)))。 %布莱克曼窗的幅度分贝表示 figure w2=2/M*[0:M1]。 %w为连续频谱的数字角频率横坐标 plot(w2,db1)。 %输出矩形窗的分贝图 grid on title(39。 矩形窗滤波器 39。 ) xlabel(39。 Frequency39。 ) ylabel(39。 DB39。 ) y1=conv(X1,h1)。 %用矩形窗对解调后的音乐信号进行滤波 sound(y1,fs) %播放滤波后的音乐信号 figure subplot(2,1,1)。 plot(y1)。 %输 出滤波之后的音乐信号的连续波型 grid on title(39。 矩形窗滤波后的音乐信号时域波形 39。 ) xlabel(39。 Time39。 ) ylabel(39。 Magtitude39。 ) y1=fft(y1,M)。 subplot(2,1,2)。 plot(w2,abs(y1))。 %输出滤波之后音乐信号的连续谱图 grid on title(39。 矩形窗滤波后的音乐信号频域波形 39。 ) xlabel(39。 Frequency39。 ) ylabel(39。 Magtitude39。 ) figure plot(w2,db2) %输出布莱克曼窗的分贝图 grid on title(39。 布莱克曼窗滤波器 39。 ) xlabel(39。 Frequency39。 ) ylabel(39。 DB39。 ) y2=conv(X1,h2)。 %用布莱克曼窗对解调后的音乐信号进行滤波 sound(y2,fs) %播放滤波之后的音乐信号 figure subplot(2,1,1)。 plot(y2)。 %输出用布莱克曼 窗滤波之后的连续波形 grid on title(39。 布莱克曼窗滤波后的音乐信号时域波形 39。 ) xlabel(39。 Time39。 ) ylabel(39。 Magtitude39。 ) y2=fft(y2,M)。 subplot(2,1,2)。 plot(w2,abs(y2))。 %输出用布莱克曼窗滤波之后的连续谱图 grid on title(39。 布莱克曼窗滤波后的音乐信号频域波形 39。 ) xlabel(39。 Frequency39。 ) ylabel(39。 Magtitude39。 ) 3. 输出波形 0 0 . 5 1 1 . 5 2 2 . 5x 1 051 0 . 8 0 . 6 0 . 4 0 . 200 . 20 . 40 . 60 . 81原始音乐信号时域波型T i m eMagnitude 0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2050010001500202025003000原始音乐信号频域波型Fr e q u e n c y / H zMagnitude 0 1 2 3x 1 051 0 . 500 . 51余弦信号时域波型T i m eMagnitude0 0 . 5 1 1 . 5 2051015x 1 04余弦信号频域波型Fr e q u e n c y / H zMagnitude0 1 2 3x 1 051 0 . 500 . 51高频调制后音乐信号时域波型T i m eMagnitude0 0 . 5 1 1 . 5 2050010001500高频调制后音乐信号频域波型Fr e q u e n c y / H zMagnitude 0 0 . 5 1 1 . 5 2 2 . 5x 1 051 0 . 500 . 51解制后音乐信号时域波型T i m eMagnitude0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2050010001500解制后音乐信号频域波型Fr e q u e n c y / H zMagnitude 0 0 . 1 0 . 2 0 . 3 0 . 4 0 . 5 0 . 6 0 . 7 0 . 8 0 . 9 100 . 20 . 40 . 60 . 811 . 21 . 4数字巴特沃斯滤波器D i g i t a l Fr e q u e n c y / p iMagnitude 0 0 . 5 1 1 . 5 2 2 . 5x 1 05 0 . 500 . 51巴特沃斯滤波后音乐信号时域波型T i m eMagnitude0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2050010001500滤波后音乐信号频域波型Fr e q u e n c y / H zMagnitude 0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2 1 8 0 1 6 0 1 4 0 1。
阅读剩余 0%
本站所有文章资讯、展示的图片素材等内容均为注册用户上传(部分报媒/平媒内容转载自网络合作媒体),仅供学习参考。 用户通过本站上传、发布的任何内容的知识产权归属用户或原始著作权人所有。如有侵犯您的版权,请联系我们反馈本站将在三个工作日内改正。