matlab微分與積分ppt課件-資料下載頁(yè)
2025-05-05 18:17本頁(yè)面
【正文】 6 例 給定數(shù)學(xué)函數(shù) x(t)=12sin(2π 10t+π/4)+5cos(2π 40t) 取 N=128，試對(duì) t從 0~1秒采樣，用 fft作快速傅立葉變換，繪制相應(yīng)的振幅 頻率圖。 在 0~1秒時(shí)間范圍內(nèi)采樣 128點(diǎn)，從而可以確定采樣周期和采樣頻率。由于離散傅立葉變換時(shí)的下標(biāo)應(yīng)是從 0到 N1，故在實(shí)際應(yīng)用時(shí)下標(biāo)應(yīng)該前移 1。又考慮到對(duì)離散傅立葉變換來(lái)說(shuō)，其振幅 | F(k)|是關(guān)于 N/2對(duì)稱(chēng)的，故只須使 k從 0到 N/2即可。 MATLAB @SDU 17 程序如下： N=128。 % 采樣點(diǎn)數(shù) T=1。 % 采樣時(shí)間終點(diǎn) t=linspace(0,T,N)。 % 給出 N個(gè)采樣時(shí)間 ti(I=1:N) x=12*sin(2*pi*10*t+pi/4)+5*cos(2*pi*40*t)。 % 求各采樣點(diǎn)樣本值 x dt=t(2)t(1)。 % 采樣周期 f=1/dt。 % 采樣頻率 (Hz) X=fft(x)。 % 計(jì)算 x的快速傅立葉變換 X F=X(1:N/2+1)。 % F(k)=X(k)(k=1:N/2+1) f=f*(0:N/2)/N。 % 使頻率軸 f從零開(kāi)始 plot(f,abs(F),39。*39。) % 繪制振幅 頻率圖 xlabel(39。Frequency39。)。 ylabel(39。|F(k)|39。)
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1