瀹為獙5FIR鏁板瓧婊ゆ尝鍣ㄨ�捐�′笌杞�浠跺疄鐜版眹鎬?- 鐧惧害鏂囧簱

棣栭〉

搴撳簱AI

骞垮満

鍘嗗彶

鏀惰棌

浜戠洏

鐭ヨ瘑搴?/div>

淇℃伅闄?/p>

 

 

 

 

14

鐢典俊锛堝笀鑼冿級

 

 

 

 

瀹為獙浜旓細

FIR

鏁板瓧婊ゆ尝鍣ㄨ�捐�′笌杞�浠跺疄鐜�

 

涓�銆佸疄楠屾寚瀵?/p>

 

1

锛庡疄楠岀洰鐨?/p>

 

锛?/p>

1

锛夋帉鎻＄敤绐楀嚱鏁版硶璁捐��

FIR

鏁板瓧婊ゆ尝鍣ㄧ殑鍘熺悊鍜屾柟娉曘�?/p>

 

锛?/p>

2

锛夋帉鎻＄敤绛夋尝绾规渶浣抽�艰繎娉曡�捐�?/p>

FIR

鏁板瓧婊ゆ尝鍣ㄧ殑鍘熺悊鍜屾柟娉曘�?/p>

 

锛?/p>

3

锛夋帉鎻?/p>

FIR

婊ゆ尝鍣ㄧ殑蹇�閫熷嵎绉�瀹炵幇鍘熺悊銆?/p>

 

锛?/p>

4

锛夊�︿細璋冪�?/p>

MATLAB

鍑芥暟璁捐�′笌瀹炵�?/p>

FIR

婊ゆ尝鍣ㄣ�?/p>

 

2

锛?/p>

 

瀹為獙鍐呭�瑰強姝ラ�?/p>

 

锛?/p>

1

锛夎�ょ湡澶嶄範绗�涓冪珷涓�鐢ㄧ獥鍑芥暟娉曞拰绛夋尝绾规渶浣抽�艰繎娉曡�捐�?/p>

FIR

鏁板瓧婊ゆ尝鍣ㄧ殑鍘熺悊锛?/p>

 

锛?/p>

2

锛夎皟鐢ㄤ俊鍙蜂骇鐢熷嚱鏁?/p>

xtg

浜х敓鍏锋湁鍔犳�у櫔澹扮殑淇″彿

xt

锛屽苟鑷�鍔ㄦ樉绀�

xt

鍙婂叾棰戣氨锛屽�傚�?/p>

1

鎵�绀?/p>

锛?/p>

 

鍥?/p>

1 

鍏锋湁鍔犳�у櫔澹扮殑淇″彿

x(t)

鍙婂叾棰戣氨濡傚浘

 

绋嬪簭浠ｇ爜锛?/p>

锛堜俊鍙蜂骇鐢熷嚱鏁?/p>

xtg

绋嬪簭娓呭崟锛?/p>

 

function

 xt=xtg(N)

 

%脢碌脩茅脦氓脨脜潞脜x(t)虏煤脡煤,虏垄脧脭脢戮脨脜潞脜碌脛路霉脝碌脤脴脨脭脟煤脧脽

 

%xt=xtg(N) 虏煤脡煤脪禄赂枚鲁陇露脠脦陋N,脫脨录脫脨脭赂脽脝碌脭毛脡霉碌脛碌楼脝碌碌梅路霉脨脜潞脜xt,虏脡脩霉脝碌脗脢Fs=1000Hz

 

%脭脴虏篓脝碌脗脢fc=Fs/10=100Hz,碌梅脰脝脮媒脧脪虏篓脝碌脗脢f0=fc/10=10Hz.

 

N=1000;Fs=1000;T=1/Fs;Tp=N*T;

 

t=0:T:(N-1)*T;

 

fc=Fs/10;f0=fc/10; 

%脭脴虏篓脝碌脗脢fc=Fs/10拢卢碌楼脝碌碌梅脰脝脨脜潞脜脝碌脗脢脦陋f0=Fc/10;

 

mt=cos(2*pi*f0*t);  

%虏煤脡煤碌楼脝碌脮媒脧脪虏篓碌梅脰脝脨脜潞脜mt拢卢脝碌脗脢脦陋f0

 

ct=cos(2*pi*fc*t);  

%虏煤脡煤脭脴虏篓脮媒脧脪虏篓脨脜潞脜ct拢卢脝碌脗脢脦陋fc

 

xt=mt.*ct;          

%脧脿鲁脣虏煤脡煤碌楼脝碌碌梅脰脝脨脜潞脜xt

 

nt=2*rand(1,N)-1;   

%虏煤脡煤脣忙禄煤脭毛脡霉nt

 

%=======脡猫录脝赂脽脥篓脗脣虏篓脝梅hn,脫脙脫脷脗脣鲁媒脭毛脡霉nt脰脨碌脛碌脥脝碌鲁脡路脰,脡煤鲁脡赂脽脥篓脭毛脡霉=======

 

fp=150; fs=200;Rp=0.1;As=70;     

% 脗脣虏篓脝梅脰赂卤锚

 

fb=[fp,fs];m=[0,1];                

% 录脝脣茫remezord潞炉脢媒脣霉脨猫虏脦脢媒f,m,dev

 

dev=[10^(-As/20),(10^(Rp/20)-1)/(10^(Rp/20)+1)];

 

[n,fo,mo,W]=remezord(fb,m,dev,Fs);   

% 脠路露篓remez潞炉脢媒脣霉脨猫虏脦脢媒

 

hn=remez(n,fo,mo,W);     

% 碌梅脫脙remez潞炉脢媒陆酶脨脨脡猫录脝,脫脙脫脷脗脣鲁媒脭毛脡霉nt脰脨碌脛碌脥脝碌鲁脡路脰

 

yt=filter(hn,1,10*nt);      

%脗脣鲁媒脣忙禄煤脭毛脡霉脰脨碌脥脝碌鲁脡路脰拢卢脡煤鲁脡赂脽脥篓脭毛脡霉yt

 

%================================================================

 

xt=xt+yt;           

%脭毛脡霉录脫脨脜潞脜

 

fst=fft(xt,N);k=0:N-1;f=k/Tp;

 

subplot(3,1,1);plot(t,xt);grid;xlabel(

't/s'

);ylabel(

'x(t)'

);