1、我們使用開發的傅立葉變換程序代碼如下:
syms x
f = exp(-2*x^2); %our
(f,[-2,2]) % plot of our
FT = (f) %
將其寫入到我們的程序模塊中。
2、我們運行上面的傅立葉變換程序代碼傅里葉變換代碼實現,將得出運行結果:FT = (2^(1/2)*pi^(1/2)*exp(-w^2/8))/2。
3、如果我們需要更高級的顯示,我們修改上述代碼何海即可,如使用(FT)作傅里葉變換折愛往泥線圖。
1.離散序列的傅里葉變換DTFT( Time )
代碼:
1 N=8; %原離散信號有8點
2 n=[0:1:N-1] %原信號是1行8列的矩陣
3 xn=0.5.^n; %構建原始信號,為指數信號
4
5 w=[-800:1:800]*4*pi/800; %頻域共-800----+800 的長度(本應是無窮傅里葉變換代碼實現,高頻分量很少,故省去)
6 X=xn*exp(-j*(n'*w)); %求dtft變換,采用原始定義的方法,對復指數分量求和而得
7 (311)
8 stem(n,xn);
9 title('原始信號(指數信號)');
10 (312);
11 plot(w/pi,abs(X));
12 title('DTFT變換')
關鍵詞: