Problem with a high pass filter

Question

Martin 2013-12-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/109023-problem-with-a-high-pass-filter

评论： Levi 2017-4-2

I would like to make a high pass filter so I have used fdatool for set up of it, which generated a function Hd, then I have imported everything to my source code but it doesn't work yet. I don't know why and can't to find and solve that problem.

function Hd = horni_propust_01
%HORNI_PROPUST_01 Returns a discrete-time filter object.
% MATLAB Code
% Generated by MATLAB(R) 8.2 and the Signal Processing Toolbox 6.20.
% Generated on: 07-Dec-2013 15:04:37
% Elliptic Highpass filter designed using FDESIGN.HIGHPASS.
% All frequency values are in Hz.
Fs = 48000;  % Sampling Frequency
N     = 10;    % Order
Fpass = 1000;  % Passband Frequency
Apass = 1;     % Passband Ripple (dB)
Astop = 80;    % Stopband Attenuation (dB)
% Construct an FDESIGN object and call its ELLIP method.
h  = fdesign.highpass('N,Fp,Ast,Ap', N, Fpass, Astop, Apass, Fs);
Hd = design(h, 'ellip');
% [EOF]
[y,Fs,nBits]=wavread('def.wav'); 
left=y(:,1);                        % Left channel 
right=y(:,2);                       % Right channel
hd = horni_propust_01;
fr = filter(hd,right);
zvuk = left - fr;
player = audioplayer(zvuk,Fs,nBits);
play(player);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Wayne King 2013-12-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/109023-problem-with-a-high-pass-filter#answer_117653

编辑：Wayne King 2013-12-7

在 MATLAB Online 中打开

Please specify what you mean by "...I have imported everything to my source code but it doesn't work yet"

What error message are you seeing?

You should just be able to execute the following in the workspace to filter one channel of your .wav file.

   Fs = 48000;  % Sampling Frequency
   N     = 10;    % Order
   Fpass = 1000;  % Passband Frequency
   Apass = 1;     % Passband Ripple (dB)
   Astop = 80;    % Stopband Attenuation (dB)
   h  = fdesign.highpass('N,Fp,Ast,Ap', N, Fpass, Astop, Apass, Fs);
   Hd = design(h, 'ellip');
   [y,Fs,nBits]=wavread('def.wav'); 
   left=y(:,1);                        % Left channel 
   right=y(:,2);                       % Right channel
   fr = filter(Hd,right);