householder变换qr分解matlab_【基础教程】Matlab实现傅里叶变换
傅立葉變換
傅立葉變換是一種常見的分析方法,傅立葉變換將滿足一定條件的函數(shù)表示為一些函數(shù)的加權(quán)和(或者積分)。可以分為四個類別:?
1. 非周期連續(xù)性信號?
對應(yīng)于傅里葉變換,頻域連續(xù)非周期?
2. 周期性連續(xù)性信號?
對應(yīng)于傅立葉級數(shù),頻域離散非周期?
3. 非周期離散信號?
對應(yīng)于DTFT(離散時間傅立葉變換),頻域連續(xù)周期?
4. 周期性離散信號?
對應(yīng)于DFT(離散時間傅立葉變換),頻域離散周期
傅立葉級數(shù)
首先從傅立葉級數(shù)開始分析,傅立葉級數(shù)是將一個信號在一組正交基上進(jìn)行分解的體現(xiàn)。
連續(xù)時間傅立葉變換
離散時間傅立葉變換
離散時間傅立葉變換在頻域上是連續(xù)的,但由于計(jì)算機(jī)無法表示無限長的時間片段,已經(jīng)無法表示全部頻率,一般取一定頻域的分量。
二維傅立葉變換
傅立葉變換實(shí)現(xiàn)
只有離散傅里葉變換才可以實(shí)現(xiàn),在MATLAB中實(shí)現(xiàn)有fft,fft2進(jìn)行傅里葉變換,同樣可以手動進(jìn)行變換。
一維傅立葉變換
% xn是信號,n是坐標(biāo),N是點(diǎn)數(shù)% N =8;% n = [0:1:N-1];% xn = 0.5.^n; % 指數(shù)信號function [] = DFTusefft(xn,n,N) figure(1); Xk=fft(xn,N); % 傅立葉變換 subplot(211); stem(n,xn); title('原信號'); subplot(212); stem(n,abs(Xk)); title('FFT變換')endDTFT?由于DTFT的頻域是連續(xù)的而且是無窮的,當(dāng)我們選擇的最高頻域足夠高時,可以基本代表信號特征,可以進(jìn)行編程。
function [] = DFT(xn,n,N) Xk = zeros(1,N); for k=1:N sn =0.0; for i=1:N sn = sn+xn(i)*exp(-j*2*pi*i*k/N); end Xk(k) = sn; end figure(2); subplot(211); stem(n,xn); title('原信號'); subplot(212); stem(n,abs(Xk)); title('DFT')end二維傅立葉變換
function [] = imageDFT2fft() I=imread('rice.png');????????I=im2double(I); [x,y] = size(I); Ax = ones(x,y); ans = ones(x,y); com = 0+1i; % 對每一列進(jìn)行DFT for m=1:y Ax(:,m) = fft(I(:,m)); end % 對每一行進(jìn)行DFT for k=1:x ans(k,:) = fft(Ax(k,:)); end F=fftshift(ans); F= abs(F); F=log(F+1); figure(7); imshow(F,[]);end往期回顧>>>>>>
基于柵格地圖——遺傳算法的機(jī)器人最優(yōu)路徑規(guī)劃【模式識別】Matlab指紋識別【基礎(chǔ)教程】matlab繪制疫情地圖matlab自動識別銀行卡號【優(yōu)化問題】遺傳算法求最值無參考圖像質(zhì)量評價之基于多特征的增強(qiáng)圖像質(zhì)量評價總結(jié)
以上是生活随笔為你收集整理的householder变换qr分解matlab_【基础教程】Matlab实现傅里叶变换的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: python递归算法_python递归算
- 下一篇: matlab 定义一个有自变量的方程_常