【正文】 據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多 強大功能集成在一個易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計以及必須進行有效數(shù)值計算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設(shè)計語言（如 C、 Fortran）的編輯模式，代表了當(dāng)今國際科學(xué)計算軟件的先進水平。 （ 4） 有利于治學(xué)態(tài)度的培養(yǎng)。 u=zeros(1,n)。)。)。dt=。 時域波形 ： ② 單位脈沖序列 ( ) ( )x n n?? 的實現(xiàn)程序： k1=3。 時域波形 ： （ 3）單邊衰減指數(shù)信號 ( ) ( )tf t e t? ??? 和實指數(shù)序列 ()nxn a? 的實現(xiàn) ① 單邊衰減指數(shù)信號 ( ) ( )tf t e t? ??? 的實現(xiàn)程序： clear t1=1。單邊衰減指數(shù)信號 39。(k)39。 xlabel(39。 for i=1:n u(i)=b*u(i)*(t(i)t0)。 plot(t,y,39。 f1=stepfun(t1,0)。ylabel(39。y(t)39。f2 與 f2 的卷積 39。 ylabel(39。 ylabel(39。)。 xlabel(39。 f=f*p。 title(39。,h)。 plot(t,y)。 信號 f(t)的傅里葉變換如圖所示： ii 試畫出信號 f(t)= 32 ()3 tet?? 的波形及其幅頻特性曲線。 MATLAB程序 ： R= 。) 。 t=5:r:5。f(t)39。 xlabel(39。 %f(t)的傅里葉變換 F1=abs(F)。 hold on。相位 (度 )39。 %f2(t)傅里葉變換 F1=real(F1)。 hold on。 syms t。 w2=40000。k39。 [Z,P,K]=buttap(N)。) 。 %傅里葉變換 y2=fftshift(y1)。) s=length(c1)。%噪聲幅值 d=[Au*sin(2*pi*5500*t)]39。)。 xlabel(39。)。我用此次機會，把自己把知識在實踐中的到了運用，把理論的抽象具體化，使 我把知識和實踐得到很好的融合。電子工業(yè)出版社。當(dāng)然主要原因還是自己不夠努力，讓自己的缺點一直保持為缺陷。 Matlab的應(yīng)用領(lǐng)域極為廣泛，除數(shù)學(xué)計算和分析外，還被廣泛用于自動控制、系統(tǒng)仿真、數(shù)字信號處理、圖形圖像分析、數(shù)理統(tǒng)計、人工智能、虛擬現(xiàn)實技術(shù)、通信工程、金融系統(tǒng)等領(lǐng)域。Hz39。 plot(f,abs(y1(1:512)))。time n39。)。噪聲信號時域圖 39。 figure(2) freqz(x) %繪制原始語音信號的頻率響應(yīng)圖 title(39。 濾波器頻率響應(yīng)圖： （ 2） MATLAB 語音信號的采集與處理 MATLAB對原始 語 音 信號進行分析，畫出它的時域波形和頻譜 。 Fs=10000。 magH=abs(H)。 subplot(3,1,3)。 t=2:R:2。頻譜 F1(jw)39。 F1=f1*exp(j*t39。 xlabel(39。f(t1)39。 w=k*W/N。F(jw)39。grid xlabel(39。) 。 %頻譜圖幅值左右對稱 subplot(2,1,1) 。 MATLAB程序 R=。 subplot(2,1,1)。 sys=tf(b,a)。)。f1(t)39。 f1=sin(k1)。 x=cos(t)。 y=lsim(sys,x,t)。 plot(t,y)。 plot(t,y)。axis([5,5,0,5000])。 title(39。單位階躍函數(shù) 39。信號的位移 39。) 時域波形 ： ③ 實現(xiàn)兩個連續(xù)信號的移位，即 00( ) ( )f t t f t t??， 或 者 ，常數(shù) 0 0t? 程序： clear all。 b=3。 for i=1:n u(i)=b*u(i)*(t(i)t0)。 title(39。ylabel(39。)。)。 時域波形 ： （ 2）單位沖激信號 ( ) ( )f t t?? 和單位脈沖序列 ( ) ( )x n n?? 的實現(xiàn) ① 單位沖激信號 ( ) ( )f t t?? 的實現(xiàn)程序： clear。)。幅值 (f)39。 t=t1:dt:t0。 （ 2） 有利于邏輯思維的鍛煉 程序設(shè)計 是公認的、最能直接有效地訓(xùn)練學(xué)生的創(chuàng)新思維，培養(yǎng)分析問題、解決問題能力的學(xué)科之一。 MATLAB 可以進行 矩陣 運算、繪制函數(shù)和數(shù)據(jù)、實現(xiàn)算法、創(chuàng)建用戶界面、連接其他編程語言的程序等，主要應(yīng)用于工程計算、控制設(shè)計、信號 處理與通訊、 圖像處理 、 信號檢測 、金融建模設(shè)計與分析等領(lǐng)域。 課程 設(shè)計中，語句的語法和常量變量的定義都有嚴格的要求，有時輸了一個中文標點、打錯了一個字母，程序無法正常運行。 u3=ones(1,n3)。 時域波形 ： ② 單位階躍序列 ( ) ( )x n u n? 的實現(xiàn)程序： k0=0。 hold off。 t=t1:dt:t2。 k2=6。t2=10。)。)。時間 (t)39。 end y=sin(2*pi*t)。39。 f2=stepfun(t1,1/a)stepfun(t1,1/a)。f2(t)39。)。)。y(t)39。y(t)39。ylabel(39。時間 (t)39。 k0=k1(1)+k2(1)。f2(t)39。 title(39。 xlabel(39。 MATLAB 程序 : syms t 。 %頻譜線寬度 t=2 :R :2 。ylabel(39。 N=200。)。w39。 P1=angle(F)。 xlabel(39。)。 F2=real(F2)。 xlabel(39。 f=fourier(exp(3*t)*Heaviside(t)) 運行結(jié)果： f =1/(3+i*w) (2) 試畫出信號 f(t)= 3 ()tet?? ， f(t4)以及信號 f(t) 4jte? 的頻譜圖。 Ap=3。)。 [A,B,C,D]=zp2ss(Z,P,K)。 t=0:1/22050:(length(x)1)/22050。 %對頻譜圖進行平移 f=0:fs/n:fs*(n1)/n。 %噪聲信號的長度 xc=fft(c1,s)。 x2=x1+d。 subplot(2,1,2)。Hz39。 圖 加噪語音信號的時域波形 圖 原始語音信號和加噪后的信號頻譜對比 與原始信號對比，區(qū)別 ： 從原始語音信號和加噪后的信號頻譜對比圖中我們可以看出：在原始信號的基礎(chǔ)上加噪后，在 04000Hz 的波形趨于平緩，而在 5500z 時的波形發(fā)生突變，即加了高頻信號的時候，