【正文】 據可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多 強大功能集成在一個易于使用的視窗環(huán)境中，為科學研究、工程設計以及必須進行有效數值計算的眾多科學領域提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設計語言（如 C、 Fortran）的編輯模式，代表了當今國際科學計算軟件的先進水平。 （ 4） 有利于治學態(tài)度的培養(yǎng)。 u=zeros(1,n)。)。)。dt=。 時域波形 ： ② 單位脈沖序列 ( ) ( )x n n?? 的實現程序： k1=3。 時域波形 ： （ 3）單邊衰減指數信號 ( ) ( )tf t e t? ??? 和實指數序列 ()nxn a? 的實現 ① 單邊衰減指數信號 ( ) ( )tf t e t? ??? 的實現程序： clear t1=1。單邊衰減指數信號 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è)出版社。當然主要原因還是自己不夠努力，讓自己的缺點一直保持為缺陷。 Matlab的應用領域極為廣泛，除數學計算和分析外，還被廣泛用于自動控制、系統(tǒng)仿真、數字信號處理、圖形圖像分析、數理統(tǒng)計、人工智能、虛擬現實技術、通信工程、金融系統(tǒng)等領域。Hz39。 plot(f,abs(y1(1:512)))。time n39。)。噪聲信號時域圖 39。 figure(2) freqz(x) %繪制原始語音信號的頻率響應圖 title(39。 濾波器頻率響應圖： （ 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。單位階躍函數 39。信號的位移 39。) 時域波形 ： ③ 實現兩個連續(xù)信號的移位，即 00( ) ( )f t t f t t??， 或 者 ，常數 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?? 的實現 ① 單位沖激信號 ( ) ( )f t t?? 的實現程序： clear。)。幅值 (f)39。 t=t1:dt:t0。 （ 2） 有利于邏輯思維的鍛煉 程序設計 是公認的、最能直接有效地訓練學生的創(chuàng)新思維，培養(yǎng)分析問題、解決問題能力的學科之一。 MATLAB 可以進行 矩陣 運算、繪制函數和數據、實現算法、創(chuàng)建用戶界面、連接其他編程語言的程序等，主要應用于工程計算、控制設計、信號 處理與通訊、 圖像處理 、 信號檢測 、金融建模設計與分析等領域。 課程 設計中，語句的語法和常量變量的定義都有嚴格的要求，有時輸了一個中文標點、打錯了一個字母，程序無法正常運行。 u3=ones(1,n3)。 時域波形 ： ② 單位階躍序列 ( ) ( )x n u 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)) 運行結果： 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ū)別 ： 從原始語音信號和加噪后的信號頻譜對比圖中我們可以看出：在原始信號的基礎上加噪后，在 04000Hz 的波形趨于平緩，而在 5500z 時的波形發(fā)生突變，即加了高頻信號的時候，