信號與系統(tǒng)課程設(shè)計(jì)(存儲版)
【正文】 后的這個序列可以看做對離散的門函數(shù)加窗得到的，所以補(bǔ)零后的頻譜接近于離散的 Sa函數(shù)；若補(bǔ)零點(diǎn)數(shù)遠(yuǎn)小于窗的時(shí)域長度， 補(bǔ) 零 之后的這個序列可以近似看做對原序列（ ( ) 1( )g n n? ?? ? ? ?）加窗得到的，所以補(bǔ) 零 后的頻譜近似于原序列的頻譜。) % 標(biāo)注橫坐標(biāo) ylabel(39。 start(AI) % 開始采樣 try [data1,time]=getdata(AI)。TriggerDelayUnits39。)。TriggerType39。% Add channels Add one channel to AI. % Configure property values Assign values to the basic setup properties, and %create the variables blocksize and Fs, which are used for subsequent analysis. %The actual sampling rate is retrieved since it may be set by the engine to a %value that differs from the specified value. set(AI,39。)。) 0 1 2 3 4 5 6 7 800 . 511 . 522 . 533 . 544 . 55對序列插零并進(jìn)行快速傅里葉變換 K ? G（K） ? 題目二： 方法一： 用計(jì)算機(jī)的聲音編輯工具錄制一段語音信號，生成 .wav文件，得到語音信號 。 對加窗后所得的 5 點(diǎn)序列，每隔一個點(diǎn)插一個零，源程序及運(yùn)行所得圖像如下： g=[1 0 1 0 1 0 1 0 1] G=fft(g,9) %進(jìn)行 9 點(diǎn)快速傅里葉變換 n=0:8 G(n+1)=G stem(0:8,abs(G),39。) axis([0 100 0 100]) title(39。對 100 點(diǎn)的序列補(bǔ)上 10 個 0 并進(jìn)行 110 點(diǎn)快速傅里葉變換 39。\fontsize{14} \bfK \rightarrow39。) 0 10 20 30 40 50 60 70 80 90 1000102030405060708090100對 g ( n ) 加時(shí)域?yàn)?100 點(diǎn)的窗并進(jìn)行 100 點(diǎn)的快速傅里葉變換 K ? G（K） ? ② 加窗后的序列進(jìn)行補(bǔ)零處理 ，分兩種情況進(jìn)行補(bǔ)零。 G=fft(g,100) %進(jìn)行 100 點(diǎn)的 快速 傅里葉變換 stem(0:99,abs(G),39。fill39。對加窗得到的信號進(jìn)行 插零 ， 可得到插零后的頻譜。（提示：可先做理論推導(dǎo)，然后將理論推導(dǎo)和編程得到的結(jié)果結(jié)合起來分析） 題目二： 在 Matlab 環(huán)境中，利用編程方法對語音信號進(jìn)行采集和譜分析。 系統(tǒng)功能 對于題目一，先對幅度為 1的直流信號進(jìn)行均勻采樣形成序列，然后對得到的序列進(jìn)行加窗、補(bǔ)零和插零處理，觀察信號頻譜結(jié)構(gòu)的變化情況。 計(jì)算 題目一： 設(shè) g（ t） =1 的傅里葉變換為 G（ jΩ） . 由公式： ? ? ? ? 可得： G（ jΩ） =2π δ（ Ω） . 設(shè) g（ n） =1 的 N 點(diǎn)離散傅里葉 為 G（ K） . 由公式： 可得： G（ K） =2π . 題目二：設(shè)語音信號為 x（ t），則其序列 x（ n）的離散傅 里 葉變換記為 X（ K） . 由離散傅里葉變換的定義公式： . 可以得到語音信號的頻譜。) xlabel(39。對 g(n)加時(shí)域?yàn)?100 點(diǎn)的窗并進(jìn)行 100 點(diǎn)的快速傅里葉變換 39。 end for n=5:10
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號-1