基于matlab的信號的采樣與恢復、采樣定理的仿真(已改無錯字)
2023-02-16 10:51:39 本頁面
【正文】 linspace(0,8*pi,length(x1))。 figure(1) subplot(211) plot(x1,cos(x1))。 cos1=cos(x1)。 n=0:(length(x1)1)。 subplot(212) plot(w,fft1(w,cos1,n))。程序運行結(jié)果如圖：程序中的fft1函數(shù)為FFT變換的定義，其定義程序如下： function result=fft1(w,hanshu,n) a=cell(1,length(w))。 for i=1:length(w)。 m=hanshu.*((exp(j*(i1)*pi/100)).^n)。 a{i}=sum(m)。 end for i=1:length(w) result(i)=a{i}。 end信號的采樣 分別對產(chǎn)生的連續(xù)余弦函數(shù)做6點、20點、40點、80點、160點和200點的等間隔采樣，并對采樣后的序列信號做FFT頻譜分析，其程序如下： n1=input(39。請輸入采樣點數(shù)n:39。)。 n=0:n11。 zb=size(n)。 figure(2) cosf=cos(4*pi*n/zb(2))。 subplot(211)。 stem(n,cosf,39。.39。)。 w=0:(pi/100):8*pi。 subplot(212) plot(w,fft1(w,cosf,n)) 程序運行結(jié)果如圖：N=6N=20N=40N=80N=160N=200低通濾波器設(shè)計設(shè)計巴特沃斯低通濾波器程序如下： [B,A]=butter(8,350/500)。 %設(shè)置低通濾波器參數(shù) [H,w]=freqz(B,A,512,2000)。 figure(3)。subplot(2,1,1)。 %繪制低通頻譜圖 plot(w*2000/(2*pi),abs(H))。 xlabel(39。Hz39。)。ylabel(39。頻率響應幅度39。)。 title(39。低通濾波器39。)。 grid。程序運行結(jié)果如下：信號的恢復采用上述低通濾波器，分別對N=6,20,40,80,160，200點的采樣信號進行恢復，其恢復程序如下： y=filter(B,A,cosf)。 figure(4)。subplot(2,1,1)。plot(y)。 xlabel(39。t39。)。ylabel(39。x(t)39。)。 title(39。連續(xù)信號39
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1