應(yīng)用matlab對圖像信號進行頻譜分析及濾波-文庫吧
2025-06-08 17:21 本頁面
【正文】 灰度圖像；對于任務(wù)5，利用快速傅里葉變換算法（FFT）對各階段數(shù)據(jù)分別進行頻譜分析，并將它們的頻譜繪制在同一張圖上作為對比。Matlab代碼clc。close all。clear all。%%圖像的讀取以及轉(zhuǎn)換x=imread(39。39。)。%讀取jpg圖像x1=rgb2gray(x)。%生成M*N的灰度圖像矩陣[M,N]=size(x1)。%求圖像規(guī)模%%生成原始序列并求頻率響應(yīng)x2=im2double(x1)。x3=zeros(1,M*N)。%初始化for i=1:M for j=1:N x3(N*(i1)+j)=x2(i,j)。 endend %將M*N維矩陣變成1維矩陣fs=1000。%掃描頻率1kHzT=1/fs。 %掃描時間間隔L=length(x3)。%計算序列長度n=0:L1。Xk3=fft(x3)。%快速傅里葉變換Xf3=fftshift(Xk3)。%中心對稱變換f=(n/L1/2)*fs。%橫坐標變幻%%生成帶有噪聲的序列并求頻率響應(yīng)fz1=356。fz2=383。fz3=427。%三個噪聲頻率xz=*sin(2*pi*fz1*n*T)+*sin(2*pi*fz2*n*T)+*sin(2*pi*fz3*n*T)。%噪聲序列x4=x3+xz。%加入噪聲的序列x5=zeros(M,N)。for i=1:M for j=1:N x5(i,j)=x4(N*(i1)+j)。 endend %一維變M*N矩陣Xk4=fft(x4)。Xf4=fftshift(Xk4)。%%設(shè)計IIR濾波器并分析相關(guān)指標wp=250*2/fs。ws=300*2/fs。Rp=3。Rs=20。%通帶截止頻率250Hz,阻帶截止頻率300Hz,通帶衰減3dB,阻帶衰減20dB[Nm,Wc]=buttord(wp,ws,Rp,Rs)。%計算滿足指標的最小
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1