matlab課程設(shè)計(jì)---matlab語(yǔ)音信號(hào)的采集與處理(編輯修改稿)
2025-07-09 05:22 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 y=lsim(sys,x,t,zi)。 plot(t,y)。 xlabel(39。時(shí)間 (t)39。)。ylabel(39。y(t)39。)。title(39。系統(tǒng)的全響應(yīng) 39。)。 系統(tǒng)的全響應(yīng)如圖所示 ： ③ 已知某 LTI系統(tǒng)的激勵(lì)為 f1=sintε (t)，單位沖激響應(yīng)為 h(t)=te2tε (t)， 試給出系統(tǒng)零狀態(tài)響應(yīng) yf(t)的數(shù)學(xué)表達(dá)式。 MATLAB程序： p=。 k1=0:p:10。 f1=sin(k1)。 k2=k1。 f2=1*exp(2*k1)。 f=conv(f1,f2)。 f=f*p。 k0=k1(1)+k2(1)。 k3=length(f1)+length(f2)2。 k=k0:p:(k3*p+k0)。 subplot(2,2,1)。 plot(k1,f1)。 title(39。f1(t)39。)。 xlabel(39。t39。),ylabel(39。f1(t)39。)。 subplot(2,2,2)。 plot(k2,f2)。 title(39。f2(t)39。)。 xlabel(39。t39。),ylabel(39。f2(t)39。)。 subplot(2,2,3)。 plot(k,f)。 h=get(gca,39。position39。)。 h(3)=*h(3)。 set(gca,39。position39。,h)。 title(39。f(t)=f1(t)*f2(t)39。)。 xlabel(39。t39。),ylabel(39。f(t)39。)。 系統(tǒng)的零狀態(tài)響應(yīng)如圖所示 ④ 程序設(shè)計(jì)實(shí)驗(yàn) 計(jì)算下述系統(tǒng)在指數(shù)函數(shù)激勵(lì)下的零狀態(tài)響應(yīng)。 H(s)= ssssssssss 23456 234 ????? ???? MATLAB程序 ： a=[1 463 12131 0]。b=[ 576 19080]。 sys=tf(b,a)。 t=0::10。 x=exp(t)。y=lsim(sys,x,t)。 plot(t,y)。 xlabel(39。時(shí)間 (t)39。)。ylabel(39。y(t)39。)。title(39。零狀態(tài)響應(yīng) 39。) 系統(tǒng)的響應(yīng)如圖所示： （ 2）連續(xù) LTI 系統(tǒng)的頻域分析 ① 編程實(shí)現(xiàn)信號(hào)的傅里葉變換和傅里葉逆變換 I 傅里葉變換。 i 已知連續(xù)時(shí)間信號(hào) f(t)=e2|t|，通過(guò)程序完成信號(hào) f(t)的傅里葉變換 MATLAB程序 ： syms t x=exp(2*abs(t))。 F=fourier(x)。 subplot(2,1,1)。 ezplot(x)。 subplot(2,1,2) 。 ezplot(F) 。 信號(hào) f(t)的傅里葉變換如圖所示： ii 試畫出信號(hào) f(t)= 32 ()3 tet?? 的波形及其幅頻特性曲線。 MATLAB 程序 : syms t 。 f=2/3*exp(3*t)*sym(39。heaviside(t)39。) 。 F=fourier(f) 。 subplot(2,1,1) 。 ezplot(f) 。 subplot(2,1,2) 。 ezplot(abs(F)) 。 信號(hào) f(t)= 32 ()3 tet?? 的波形及其幅頻特性曲線如圖所示 II 傅里葉逆變換。 已知 f(jω )= 211??，求信號(hào) F(jω )的逆傅里葉變換。 MATLAB程序 R=。 t=2:R:2。 syms t w ifourier(1/(1+w^2),t) 運(yùn)行結(jié)果 ： ans =1/2*exp(t)*heaviside(t)+1/2*exp(t)*heaviside(t) III 傅里葉變換數(shù)值計(jì)算。 已知門函數(shù) f(t)=g2(t)=ε (t+1)ε (t1)， 試采用數(shù)值計(jì)算方法確定信號(hào)的傅里葉變換 F(jω )。 MATLAB程序 ： R= 。 %頻譜線寬度 t=2 :R :2 。 f=stepfun(t,1)stepfun(t,1) 。 W1=2*pi*5 。 %頻率寬度 N=500 。 %采樣數(shù)為 N k=0 :N 。 W=k*W1/N 。 %W 頻率正半周的采樣點(diǎn) F=f*exp(j*t39。*W)*R 。 %求 F（ jw） F=real(F) 。 %求 F（ jw）的實(shí)部 W=[fliplr(W),W(2 :501)] 。 %頻譜圖頻率寬度左右對(duì)稱 F=[fliplr(F),F(2 :501)] 。 %頻譜圖幅值左右對(duì)稱 subplot(2,1,1) 。plot(t,f) 。 xlabel(39。t39。) 。ylabel(39。f(t)39。) 。axis([2,2,2]) 。 title(39。f(t)= ε(t+1)ε(t1)39。) 。 subplot(2,1,2) 。plot(W,F) 。 title(39。f(t)傅氏變換 39。) 。 信號(hào)的傅里葉變換如圖所示： IV 傅里葉變換的時(shí)移特性 分別繪出信號(hào) f(t)= 21 ()2 tet?? 信號(hào) f(t1)的頻譜圖， 并觀察信號(hào)時(shí)移對(duì)信號(hào)頻譜的影響。 i f(t)= 21 ()2 tet??的頻譜。 MATLAB程序 ： r=。 t=5:r:5。 N=200。 W=2*pi。 k=N:N。 w=k*W/N。 f1=1/2*exp(2*t).*stepfun(t,0)。 F=r*f1*exp(j*t39。*w)。 F1=abs(F)。 P1=angle(F)。 subplot(3,1,1)。 plot(t,f1)。grid xlabel(39。t39。)。 ylabel(39。f(t)39。)。 title(39。f(t)39。)。 subplot(3,1,2)。 plot(w,F1)。 xlabel(39。w39。)。 grid。 ylabel(39。F(jw)39。)。 subplot(3,1,3) plot(w,P1*180/pi)。 grid。 xlabel(39。w39。)。 ylabel(39。相位 （ 度 ） 39。)。 信號(hào) f(t)= 21 ()2 tet?? 傅里葉變換的時(shí)頻特性如圖所示： ii f(t1)的頻譜。 MATLAB程序 ： r=。 t=5:r:5。 N=200。 W=2*pi。 k=N:N。 w=k*W/N。 f1=1/2*exp(2*(t1)).*stepfun(t,1)。 %f(t) F=r*f1*exp(j*t39。*w)。 %f(t)的傅里葉變換 F1=abs(
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1