『壹』 matlab 如何調用時頻分析工具箱
找到工具箱是指toolbox嗎 ?之後怎麼辦呢?麻煩說具體一點,要是有圖就更好了。新手報到
『貳』 求wigner ville 分布的程序,要fortran or c寫的。已附上matlab程序,轉換也成
if (nargin .eq. 0) then
error('At least one parameter required')
end if
xrow = size(x,1)
xcol = size(x,2)
if (nargin .eq. 1) then
t=1:xrow
N=xrow
trace=0
elseif (nargin .eq. 2) then
N=xrow
trace=0
elseif (nargin .eq. 3) then
trace = 0
end if
if (N<0) then
error('N must be greater than zero')
end if
trow = size(t,1)
tcol = size(t,2)
if ((xcol.eq.0).or.(xcol>2)) then
error('X must have one or two columns')
elseif (trow.neq.1) then
error('T must only have one row')
elseif (2^nextpow2(N).neq.N) then
print*,'For a faster computation, N should be a power of two'
end if
tfr= 0
if (trace .neq. 0) then
print*,'Wigner-Ville distribution'
end if
do icol=1,tcol,
ti= t(icol)
taumax=min([ti-1,xrow-ti,round(N/2)-1])
tau=-taumax:taumax
indices= rem(N+tau,N)+1
tfr(indices,icol) = x(ti+tau,1) * conj(x(ti-tau,xcol))
tau=round(N/2)
if (ti<=xrow-tau)&(ti>=tau+1) then
tfr(tau+1,icol) = 0.5 * (x(ti+tau,1) * conj(x(ti-tau,xcol)) + ...
x(ti-tau,1) * conj(x(ti+tau,xcol)))
end if
if (trace .neq. 0) then
disprog(icol,tcol,10)
end if
end do
tfr= fft(tfr)
if (xcol.eq.1) then
tfr=real(tfr)
end if
if (nargout.eq.0) then
tfrqview(tfr,x,t,'tfrwv')
elseif (nargout.eq.3) then
f=(0.5*(0:N-1)/N)'
end if
說明:還有幾個MATLAB自帶的函數我沒有改,比如error, tfrqview,基本都是顯示輸出的函數,你要自己改一下。
『叄』 matlab調用不了信號處理工具箱的函數
說明你的matlab工具箱里沒有這個函數,從別的可以運行的matlab中的toolbox中查找這個函數文件(m文件),然後復制到你的toolbox里,再嘗試應該就沒問題了
『肆』 matlab在信號系統的常用函數
MATLAB信號處理常用函數
【一】、波形產生
函數名功能
sawtooth產生鋸齒波或三角波
Sinc產生sinc或函數sin(pi*t)/(pi*t)
Square產生方波
Diric產生Dirichlet或周期sinc函數
【二】、濾波器分析和實現
函數名功能
Abs求絕對值(幅值)
Freqs模擬濾波器頻率響應
Angle求相角
Freqspace頻率響應中的頻率間隔
Conv求卷積
Freqz數字濾波器頻率響應
Fftfilt重疊相加法FFT濾波器實現
Grpdelay平均濾波器延遲(群延遲)
Filter直接濾波器實現
Impz數字濾波器的沖激響應
Filtfilt零相位數字濾波
Zplane離散系統零極點圖
FiltieFilter函數初始條件選擇
【三】、線性系統變換
函數名功能
Convmtx卷積矩陣
Ss2tf變系統狀態空間形式為傳遞函數形式
Ploy2rc從多項式系數中計算反射系數
Ss2zp變系統狀態空間形式為零極點增益形式
Rc2ploy從反射系數中計算多項式系數
Tf2ss變系統傳遞函數形式為狀態空間形式
ResiezZ變換部分分式展開或留數計算
Tf2zp變系統傳遞函數形式為零極點增益形式
Sos2ss變系統二階分割形式為狀態空間形式
Zp2sos變系統零極點形式為二階分割形式
Sos2zp變系統二階分割形式為零極點增益形式
Zp2tf變系統零極點增益形式為傳遞函數形式
Ss2sos變系統狀態空間形式為二階分割形式
【四】、IIR濾波器設計
BesselfBessel(貝塞爾)模擬濾波器設計
Cheby2Chebyshev(切比雪夫)II型模擬濾波器設計
ButterButterworth(巴特沃思)模擬濾波器設計
Ellip橢圓模擬濾波器設計
Cheby1Chebyshev(切比雪夫)I型模擬濾波器設計
Yulewalk遞歸數字濾波器設計
【五】、IIR濾波器階選擇
ButtordButterworth(巴特沃思)濾波器階的選擇
Cheb2ordChebyshev(切比雪夫)II型濾波器階的選擇
EhebordChebyshev(切比雪夫)I型濾波器階的選擇
Clipord橢圓濾波器設計階的選擇模擬原型濾波器設計
BesselapBessel模擬低通濾波器原型
Cheb2apChebyshev(切比雪夫)II型低通濾波器原型
ButtapButterworth(巴特沃思)模擬低通濾波器原型
Ellipap橢圓模擬低通濾波器原型
Cheb1apChebyshev(切比雪夫)I型低通濾波器原型
【六】、頻率變換
Lp2bp低通到帶通模擬濾波器轉換
Lp2bs低通到帶阻模擬濾波器變換
Lp2hp低通到高通模擬濾波器變換
Lp2lp低通到低通模擬濾波器轉換
【七】、濾波器離散化
Blinear雙線性變換
Impinvar沖激響應不變法
【八】、FIR濾波器設計
Fir1基於窗函數的FIR濾波器設計—標准響應
Intfilt內插FIR濾波器設計
Fir2基於窗函數的FIR濾波器設計—任意響應
RemezFirls最小二乘FIR濾波器設計
RemezordParks-McCellan最優FIR濾波器j階估計
【九】、窗函數
Boxcar矩形窗
HanningHanning(漢寧)窗
Triang三角窗
BlackmanBlackman(布萊克曼)窗
BartlettBartlett(巴特得特)窗
ChebwinChebyshev(切比雪夫)窗
HammingHamming(漢明)窗
KaiserKaiser(凱澤)窗
【十】、變換
Ctz線性調頻Z變換
Fft一維快速傅里葉變換
Dct離散餘弦變換
Ifft一維快速傅里葉逆變換
Idct逆離散餘弦變換
Fftshift重新排列fft的輸出
Dftmtx離散傅里葉變換矩陣
HilbertHilbert(希爾伯特)變換
【十一】、統計信號處理
Cov協方差矩陣
Psd信號功率譜密度(PSD)估計
Xcov互協方差函數估計
Tfe從輸入輸出中估計傳遞函數
Corrcoef相關系數矩陣
Periodogram採用周期圖法估計功率譜密度
Xcoor互相關系數估計
Pwelch採用Welch方法估計功率譜密度
Cohere相關函數平方幅值估計
Rand生成均勻分布的隨機數
Csd互譜密度估計
Randn生成正態分布的隨機數
【十二】、自適應濾波器部分
Adaptfilt.lms最小均方(LMS)自適應演算法
Adaptfilt.rls遞推最小二乘(RLS)自適應演算法
Adaptfilt.nlms歸一化最小均方(NLMS)自適應演算法
【十三】、時頻分析與小波變換部分
Spectrogram短時傅里葉變換
Idwt單級離散一維小波逆變換
Waveinfo介紹小波工具箱中所有小波的信息
Wavedec多級離散一維小波分解
Cwt連續一維小波變換
Appcoef一維小波變換近似系數
Dwt單級離散一維小波變換
Detcoef一維小波變換細節系數
【十四】、二維信號處理
Conv2二維卷積
Xcorr2二維互相關參數
Fft2二維快讀傅里葉變換
Dwt2單級離散二維小波變換
Ifft2二維逆快速傅里葉變換
Idwt2單級離散二維小波逆變換
Filter2二維數字濾波器
Waverec2多級離散二維小波分解
『伍』 matlab 時頻分析工具箱怎麼使用
將文件夾拷貝到MATLAB下work目錄中。
在MATLAB中file菜單下的setpath中添加整個work目錄後保存。
『陸』 我要用Matlab做一個離散信號的Wigner分布,調用了tfrwv函數,但做出來的wigner分布不對,求一個示常式序。
關於時頻工具箱,我倒沒用過。我做數據一般都是用FFT,然後設計濾波演算法。
關於你的問題,我查了一下,這邊有一個步驟,供你參考下。
詳細步驟:
1、x=sig;(也就是你要處理的信號)
2、如果信號是實信號,則求其解析信號:y=hilbert(x);
3、求其WVD,[tfr,t,f]=tfrwv(y);
4、作圖:contour(t,f,tfr),其中f為歸一化頻率,可以乘上信號的采樣頻率的一般還原,看其頻率分布,另外tfr ,最好寫成abs(tfr),做出的圖好看。