matlab 怎么运行不出来 ，错误在哪里？

来源：百度知道 编辑：UC知道 时间：2024/09/19 16:03:05
[H,w ] = freqz ( b, a, 1024, 'whole') ;
% b和a分别为H ( z)分子和分母多项式系数(对F IR而言, b = a, a = 1).
H = (H (1: 501) ) ;
w = (w (1: 501) ) ;
mag = abs(H) ;
db = 203 log10 ( (mag + ep s) /max(mag) ) ;
pha = angle (H) ;
grd = grpdelay( b, a,w) ;
IDEAL_LP. M
function hd = ideal_lp (wc,M)
%理想低通滤波器; hd为0到M - 1之间的理想脉冲响应; wc为截止频率;M为理想滤波器的长度
alpha = (M - 1) /2;
n = 0:M - 1;
m = n - alpha + ep s;
hd = sin (wc3 m). / (p i3 m) ;
ACTUAL_LP. m
wp = 0. 103 pi;
ws = 0. 203 pi;
tr_width =ws - wp;
M = ceil(6. 63 p i/ tr_width) + 1;
disp ( [ ' 滤波器的长度为', num2 str (M) ] ) ;
n = 0:M - 1;
wc = (ws +wp) /2;
276
hd = IDEAL_LP (wc,M) ;
w_han = ( hanning (M) ) ;
h = hd. 3 w_han;
[ db,mag, pha, grd,w ] = FREQZ_M ( h, [ 1 ] ) ;
delta_w = 23 pi/1000;
Rp = - (min ( db (1: 1: wp /delta_w + 1) ) ) ;
disp ( ['实际通带波动为', num2 str (Rp) ] ) ;
As = - round

试试吧！
[H,w ] = freqz ( b, a, 1024, 'whole') ;
% b和a分别为H ( z)分子和分母多项式系数(对F IR而言, b = a, a = 1).
H = (H (1: 501) ) ;
w = (w (1: 501) ) ;
mag = abs(H) ;
db = 203*log10((mag + eps) /max(mag));
pha = angle (H) ;
grd = grpdelay( b, a,w) ;
IDEAL_LP. M
functionname hd = ideal_lp (wc,M)
%理想低通滤波器; hd为0到M - 1之间的理想脉冲响应; wc为截止频率;M为理想滤波器的长度
alpha = (M - 1) /2;
n = 0:M - 1;
m = n - alpha + eps;
hd = sin (wc3m)./(pi3m) ;
ACTUAL_LP. m
wp = 0.103*pi;
ws = 0.203*pi;
tr_width =ws - wp;
M = ceil(6.63*pi/tr_width) + 1;
disp ( [ ' 滤波器的长度为', num2 str (M) ] ) ;
n = 0:M - 1;
wc = (ws +wp) /2;
276
hd = IDEAL_LP (wc,M) ;
w_han = ( hanning (M) ) ;
h = hd.*3*w_han;
[ db,mag, pha, grd,w ] = FREQZ_M ( h, [ 1 ] ) ;
delta_w = 23*pi/1000;
Rp = - (min ( db (1: 1: wp /delta_w + 1) ) ) ;
disp ( ['实际通带波动为', num2 str