dsp鍚勭�嶇▼搴� - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

DSP

棰樺簱璇﹁В

 

 

 

1

宸茬煡

3

闃舵き鍦?/p>

IIR

鏁板瓧浣庨�氭护娉㈠櫒鐨勬�ц兘鎸囨爣涓猴細

閫氬甫鎴�姝㈤�戠巼

0.4

蟺

锛?/p>

閫氬甫娉㈢汗涓?/p>

0.6dB

锛?/p>

鏈�灏忛樆甯﹁“鍑忎负

32dB

銆傝�捐�′竴涓?/p>

6

闃跺叏閫氭护娉㈠櫒瀵瑰叾閫氬甫鐨勭兢寤舵椂杩涜�屽潎琛°�傜粯鍒朵綆閫氭护

娉㈠櫒鍜岀骇鑱旀护娉㈠櫒鐨勭兢寤舵椂銆?/p>

 

%Progranm 1 

% Group-delay equalization of an IIR filter. 

% 

[n,d] = ellip(3,0.6,32,0.4); 

[GdH,w] = grpdelay(n,d,512); 

plot(w/pi,GdH); grid 

xlabel('\omega/\pi'); ylabel('Group delay, samples'); 

title('Original Filter'); 

 

F = 0:0.001:0.4; 

g = grpdelay(n,d,F,2); 

 

 

% Equalize the passband 

Gd = max(g)-g; 

% Design the allpass delay equalizer 

[num,den,tau] = iirgrpdelay(6, F, [0 0.4], Gd); 

 

%

璁捐�″叚闃剁殑鍏ㄩ�氭护娉㈠櫒

 

[b,a]=iirgrpdelay(6,F,[0 0.4],Gd);  

 

 

He1=dfilt.df2(b,a);

 

He=dfilt.df2(n,d);

 

He_all=dfilt.cascade(He,He1);

 

grpdelay(He_all) 

 

 

 

 

%DFILT:Digital Filter Implementation. 

[GdA,w] = grpdelay(num,den,512); 

figure(2); 

plot(w/pi,GdH+GdA); grid 

xlabel('\omega/\pi');ylabel('Group delay, samples'); 

title('Group Delay Equalized Filter'); 

 

2

璁捐�″反鐗规矁鍏规ā鎷熶綆閫氭护娉㈠櫒锛屽叾婊ゆ尝鍣ㄧ殑闃舵暟鍜?/p>

3-dB

鎴�姝㈤�戠巼鐢遍敭鐩樿緭鍏ワ紝绋嬪簭鑳芥牴

鎹�杈撳叆鐨勫弬鏁帮紝缁樺埗婊ゆ尝鍣ㄧ殑澧炵泭鍝嶅簲銆?/p>

 

% Program 2 

% Program to Design Butterworth Analog Lowpass Filter 

% 

% Type in the filter order and passband edge frequency

璁捐�¤竟缂橀�戠巼

 

N = input('Type in filter order = '); 

Wn = input('3-dB cutoff angular frequency = '); 

% Determine the transfer function 

[num,den] = butter(N,Wn,'s'); 

 

 

 

 

 

%

璁＄畻婊ゆ尝鍣ㄧ殑鍒嗗瓙鍜屽垎姣?/p>

 

% 

 

Compute and plot the frequency response