FIR低通濾波器+matlab編程+濾波前后圖形

Matlab實(shí)現(xiàn)振動(dòng)信號(hào)低通濾波附件txt中的數(shù)字是一個(gè)實(shí)測(cè)振動(dòng)信號(hào)，采樣頻率為5000Hz，試設(shè)計(jì)一個(gè)長(zhǎng)度為M=32的FIR低通濾波器，截止頻率為600Hz，用此濾波器對(duì)此信號(hào)進(jìn)行濾波。要求：（1） 計(jì)算數(shù)字截止頻率；（2） 給出濾波器系數(shù)；（3） 繪出原信號(hào)波形；（4） 繪出濾波后的信號(hào)波形；解答過(guò)程：第一部分：數(shù)字截止頻率的計(jì)算數(shù)字截止頻率等于截止頻率除以采樣頻率的一半，即=600/5000/2=0.24第二部分：濾波器系數(shù)的確定在matlab中輸入如下程序，即可得到濾波器系數(shù)：n=32Wn=0.24b=fir1(n,Wn)得到的濾波器系數(shù)b為Columns 1 through 9 -0.0008 -0.0018 -0.0024 -0.0014 0.0021 0.0075 0.0110 0.0077 -0.0054Columns 10 through 18 -0.0242 -0.0374 -0.0299 0.0087 0.0756 0.1537 0.2166 0.2407 0.2166 Columns 19 through 27 0.1537 0.0756 0.0087 -0.0299 -0.0374 -0.0242 -0.0054 0.0077 0.0110 Columns 28 through 33 0.0075 0.0021 -0.0014 -0.0024 -0.0018 -0.0008第三部分：原信號(hào)波形將附件4中的dat文件利用識(shí)別軟件讀取其中的數(shù)據(jù)，共1024個(gè)點(diǎn)，存在TXT文檔中，取名bv.txt，并復(fù)制到matlab的work文件夾。在matlab中編寫如下程序：x0=load(zhendong.txt); %找到信號(hào)數(shù)據(jù)地址并加載數(shù)據(jù)。t=0:1/5000:1023/5000; %將數(shù)據(jù)的1024個(gè)點(diǎn)對(duì)應(yīng)時(shí)間加載figure(1);plot(t,x0);xlabel(t/s);ylabel(幅值);運(yùn)行之后就得到如下波形，即振動(dòng)信號(hào)的原始波形圖：第四部分：濾波后的波形圖在matlab窗口中輸入如下程序，即可得到濾波后的波形圖：x0=load(zhendong.txt);t=0:1/5000:1023/5000;figure(1);plot(t,x0);xlabel(t/s);ylabel(幅值);fs=5000;n=32;Wn=0.24;b=fir1(n,Wn); %求濾波系數(shù)y0=filter(b,1,x0); %用matlab自帶的filter函數(shù)進(jìn)行濾波figure(2);p