![仂亞舒仄仄舒 1. 亠舒仍亳亰舒亳 亠亠仆亳 仗亳仄亠舒 1 于 仂仄亠 m-舒亶仍舒
MATLAB
FS=1000; % 丼舒仂舒 亟亳从亠亳亰舒亳亳 FN=FS/2; % 丼舒仂舒
舒亶从于亳舒 N=73; % 仍亳仆舒 亳仍舒
beta=5.65; % 舒舒仄亠 仆亠舒于仆仂仄亠仆仂亳 仆从亳舒亶亰亠舒
fc1=125/FN; % 仂仄亳仂于舒仆仆亠 舒仂 亠亰舒 fc2=275/FN;
FC=[fcl fc2]; % 亠从仂 从舒亠于 舒仂
hn=fir1(N-l,FC,kaiser(N,beta)); % 仂仍亳
从仂亳亳亠仆 亳仍舒
[H,f]=freqz(hn,1,512,FS); % 亳仍亳 舒仂仆
舒舒从亠亳亳从
mag=20*loglO(abs(H));
plot(f,mag),
grid on
xlabel ('丼舒仂舒 ())
ylabel('仄仗仍亳亟仆舒 舒舒从亠亳亳从舒 (亟)')](https://image.slidesharecdn.com/pkrdspav-140214010401-phpapp02/85/P-kr-dsp-a-v-8-320.jpg)
![亳仄亠 2
仗仂仍亰 仂仗亳仄舒仍仆亶 仄亠仂亟, 于亳仍亳亠 从仂亳亳亠仆 亳 亳亰仂弍舒亰亳亠
舒仂仆 舒舒从亠亳亳从 仗仂仍仂仂于仂亞仂 丱足亳仍舒 仍亳仆亠亶仆仂亶 舒亰仂于仂亶
舒舒从亠亳亳从仂亶, 亟仂足于仍亠于仂ﾑ亠亞仂 仍亠亟ﾑ亳仄 仗亠亳亳从舒亳礆:
仗仂仍仂舒 仗仂仗从舒仆亳
1000足1500 ,
仗仂仍仂舒 仗亠亠仂亟舒
500 ,
亟仍亳仆舒 亳仍舒
41,
舒仂舒 亟亳从亠亳亰舒亳亳 10 000 .
仂亞舒仄仄舒 2. 亠舒仍亳亰舒亳 于亳仍亠仆亳 从仂亳亳亠仆仂于 仂仗亳仄舒仍仆仂亞仂
丱-亳仍舒 亳 于于仂亟舒 仆舒 从舒仆 舒仂仆仂亶 舒舒从亠亳亳从亳 于 仂仄亠 m舒亶仍舒 MATLAB
Fs=1000; % 丼舒仂舒 亟亳从亠亳亰舒亳亳 N=41; % 仍亳仆舒 亳仍舒
=[0 0 1 1 0 0]; % 亠仍舒亠仄舒 舒仄仗仍亳亟仆舒 舒舒从亠亳亳从舒
F=[0, 0.1, 0.2, 0.3, 0.4, 1]; % 舒 仗仂仍仂
b=remez(N-l, F, ); % 亳仍亳 从仂亳亳亠仆 亳仍舒
[, f]=freqz(b, 1, 512, Fs); % 亳仍亳 舒仂仆
舒舒从亠亳亳从
mag=20*logl0(abs()); % 亳仍舒 亳 仆舒亳仂于舒 亠亠
plot(f ,mag)
xlabel ('Frequency (Hz))
ylabel ('Magnitude (dB))](https://image.slidesharecdn.com/pkrdspav-140214010401-phpapp02/85/P-kr-dsp-a-v-10-320.jpg)
![亳仄亠 3
舒亰舒弍仂从舒 仗仂仂亞仂 亳仍舒 仆亳亢仆亳 舒仂 亳仗仂仍亰仂于舒仆亳亠仄 弍亳仍亳仆亠亶仆仂亞仂 z仗亠仂弍舒亰仂于舒仆亳 亳 MATLAB. 丐亠弍亠 亳仂于仂亶 丱-亳仍
舒舒从亠亳亳从仂亶 舒亠于仂舒, 亟仂于仍亠于仂ﾑ亳亶 仍亠亟ﾑ亳仄 仗亠亳亳从舒亳礆:
舒仂舒 亠亰舒
150 ,
舒仂舒 亟亳从亠亳亰舒亳亳
1,28 从,
仗仂磲仂从 亳仍舒, N
2.
亠亠仆亳亠. 弌仂仂于亠于ﾑ亳亶 m-舒亶仍 MATLAB 仗亳于亠亟亠仆 于 仗仂亞舒仄仄亠 3 .
1.
弌 仗仂仄仂 m-舒亶仍舒 仗仂仍舒亠仄 仍亠亟ﾑ亳亠 于亠从仂 从仂亳亳亠仆仂于 (b 亳
a) 仆仍亠亶 亳 仗仂仍ﾑ仂于 (z 亳 ) 丱-亳仍舒:
b= [0.0878, 0.1756, 0.0878]
a= [1.0000, -1.0048, 0.3561]
z =[-1,-1]
= [0.5024 賊 0.3220i]
k = 0.0878
仗仂仍亰 从仂亳亳亠仆, 亰舒仗亳于舒亠仄 仗亠亠亟舒仂仆 仆从亳:
H ( z) =
0, 0878(1 + 2 z 1 + z 2 )
1 1, 0308 z 1 + 0,3553 z 2](https://image.slidesharecdn.com/pkrdspav-140214010401-phpapp02/85/P-kr-dsp-a-v-14-320.jpg)
![仂亞舒仄仄舒 3. 亠亠仆亳亠 仗亳仄亠舒 9 于 仂仄亠 m-舒亶仍舒
MATLAB
% 仂仂亶 亳仍 仆亳亢仆亳 舒仂 %
N=2; % 仂磲仂从 亳仍舒
Fs=1280; % 丼舒仂舒 亟亳从亠亳亰舒亳亳
FN=Fs/2;
fc=150; % 丼舒仂舒 亠亰舒
Fc=fc/FN; % 仂仄亳仂于舒仆仆舒 舒仂舒 亠亰舒
[b, a]=butter (N,Fc); % 弌仂亰亟舒 亳 仂亳仂于舒
舒仆舒仍仂亞仂于亶 亳仍
[z, , k]=butter (N, Fc);
subplot (2,1,1) % 于亠亳 仆舒 从舒仆 舒仄仗仍亳亟仆仂舒仂仆 舒舒从亠亳亳从
[, f]=freqz(b, 舒, 512, Fs);
plot(f, abs(H))
x label ('Frequency (Hz)[1])
ylabel('Magnitude Response (dB)')
subplot(2,1,2) % 于亠亳 仆舒 从舒仆 亟亳舒亞舒仄仄
仆仍亠亶 亳 仗仂仍ﾑ仂于 zplane(b, 舒)
仄仗仍亳亟仆仂-舒仂仆舒 舒舒从亠亳亳从舒 亳
亟亳舒亞舒仄仄舒 仆仍亠亶 亳 仗仂仍ﾑ仂于](https://image.slidesharecdn.com/pkrdspav-140214010401-phpapp02/85/P-kr-dsp-a-v-15-320.jpg)
More Related Content
Viewers also liked (20)
Similar to P kr dsp a&v (20)
![[FEE] 4. Bipolar transistors](https://cdn.slidesharecdn.com/ss_thumbnails/4-210620131609-thumbnail.jpg?width=560&fit=bounds)
More from student_kai (20)
P kr dsp a&v
- 1. 丶亂 弌弌丐亂弌弌 丕弌丐丐 弌丕弌丐 丐亂 丕丼 乘弌丿 个弌弌亂 亊 弌 弌丕弌丐乘 丐丱丼弌 丕弌丐丐 亳仄. ..丐丕 .., ..弌丕, .., ..弌乘丕 丶个亊 丐 丕 弌 亠仂亟亳亠从亳亠 从舒亰舒仆亳 仗仂 于仗仂仍仆亠仆亳 从仂于仂亶 舒弍仂 亠仂亟亳亠从亳亠 从舒亰舒仆亳 仗亠亟仆舒亰仆舒亠仆 亟仍 弍舒从舒仍舒于仂于 亅丐, 仂弍舒ﾑ亳 仗仂 弍舒从舒仍舒于从仂亶 仗仂亞舒仄仄亠 仆舒仗舒于仍亠仆亳 210400 束舒亟亳仂亠仆亳从舒損, 仗仂亳仍 束亟亳仂于亳亰舒仍仆舒 亠仆亳从舒損 亂 2010
- 2. 仂于舒 舒弍仂舒 仗仂 从 束丶亳仂于舒 仂弍舒弍仂从舒 舒亟亳仂 亳 于亳亟亠仂亳亞仆舒仍仂于損 于仗仂仍仆磳 亟亠仆舒仄亳 舒亟亳仂亠仆亳亠从亳 仗亠亳舒仍仆仂亠亶, 仂弍舒ﾑ亳 仗仂 仗仂亳仍 仗仂亟亞仂仂于从亳 束亟亳仂于亳亰舒仍仆舒 亠仆亳从舒損 仆舒仗舒于仍亠仆亳 210400 束舒亟亳仂亠仆亳从舒損. 仍 于仗仂仍仆亠仆亳 从仂于仂亶 舒弍仂 仗仂 从 束丶亳仂于舒 仂弍舒弍仂从舒 舒亟亳仂 亳 于亳亟亠仂亳亞仆舒仍仂于損 亟亠仆 亟仂仍亢亠仆 仗亳仄亠仆亳 亰仆舒仆亳, 仗仂仍亠仆仆亠 仗亳 亳亰亠仆亳亳 从仂于 束舒亟亳仂亠仆亳亠从亳亠 亠仗亳 亳 亳亞仆舒仍損, 束弌亠仄仂亠仆亳从舒 舒仆舒仍仂亞仂于 仍亠从仂仆仆 仂亶于損, 束仆仂于 亠仍亠于亳亟亠仆亳晛, 束亅仍亠从仂舒从亳从舒損, 束丶亳仂于舒 仂弍舒弍仂从舒 舒亟亳仂 亳 于亳亟亠仂亳亞仆舒仍仂于損. 仂仄亠 仂亞仂, 亟仍 于仗仂仍仆亠仆亳 从仂于仂亶 舒弍仂 于 仆舒仂亳 仄亠仂亟亳亠从亳 从舒亰舒仆亳 仗亳于仂亟亳 弍亳弍仍亳仂亞舒亳亠从亳亶 仗亳仂从 亠从仂仄亠仆亟亠仄仂亶 仍亳亠舒.
- 3. 丶亂 丼 丕弌 丐亊 仂于舒 舒弍仂舒 于仗仂仍仆磳 亟仍 亰舒从亠仗仍亠仆亳 亳 舒亳亠仆亳 亠仂亠亳亠从亳 亰仆舒仆亳亶 仗仂 从 束丶亳仂于舒 仂弍舒弍仂从舒 舒亟亳仂 亳 于亳亟亠仂亳亞仆舒仍仂于損, 仗仂仍亠仆亳 仆舒于从仂于 亳仆亠亰舒 亳 舒亠舒 亳仂于 亳仍仂于. 亳 于仗仂仍仆亠仆亳亳 从仂于仂亶 舒弍仂 亟亠仆 亟仂仍亢亠仆: - 仆舒亳 舒仄仂仂亠仍仆仂亶 舒弍仂亠 亠仆亳亠从仂亶 仍亳亠舒仂亶 亳 仆仂仄舒亳于仆仄亳 亟仂从仄亠仆舒仄亳; - 仂亰仆舒从仂仄亳 仄亠仂亟舒仄亳 仄仂亟亠仍亳仂于舒仆亳 亳仂于 亳仍仂于. 仗仂仍仆亠仆亳亠 仗仂亠从舒 亟舒亠 仂仗 于仗仂仍仆亠仆亳 亳仆亢亠仆亠仆 舒亠仂于, 仗舒从亳从 仂舒于仍亠仆亳 亳 仂仂仄仍亠仆亳 仗仂仆亳亠仍仆仂亶 亰舒仗亳从亳, 于仗仂仍仆亠-仆亳 仆从亳仂仆舒仍仆 亠仄, 仗亠亟舒于仍亠仆亳 亳 亰舒亳 亠仆亳亠从亳 亠亠仆亳亶. 舒 于亠仄 从仂于仂亞仂 仗仂亠从亳仂于舒仆亳 亟亠仆 仗亳仂弍亠舒亠 仆舒于从亳 亟仍 于仗仂仍仆亠仆亳 亳仂亞仂于仂亶 舒弍仂 于 仆亳于亠亳亠亠 亟亳仗仍仂仄仆仂亞仂 仗仂亠从-舒 (舒弍仂).
- 4. 弌 丕弌 丐乘 舒亟舒舒 亟亠仆舒 - 仗仂亠从亳仂于舒 仂亟亳仆 亳亰 4- 亳仗仂于 亳仍仂于 (个丼, 个丼, 个 亳 个) 亟于 从仍舒仂于 丱 亳 丱, 仗仂 亠仄 仄亠仂亟亳从舒仄 (亟仍 丱 仄亠仂亟 于亰于亠亳于舒仆亳, 仄亠仂亟 舒仂仆仂亶 于弍仂从亳, 仂仗亳仄亳亰舒亳仂仆仆亠 仄亠仂亟; 亟仍 丱 仄亠仂亟 亳仆于舒亳舒仆仆仂亞仂 仗亠仂弍舒亰仂于舒仆亳 亳仄仗仍仆仂亶 舒舒从亠亳亳从亳, 弍亳仍亳仆亠亶仆仂亠 仗亠仂弍舒亰仂于舒仆亳亠, 舒亰仄亠亠仆亳亠 仆仍亠亶 亳 仗仂仍ﾑ仂于). 舒亠 仗仂亳亰于仂亟亳 于 仗仂亞舒仄仄仆仂仄 仗仂亟从亠 Matlab, 亳仄亠ﾑ亠仄 仂弍亳仆亠 于仂亰仄仂亢仆仂亳 亟仍 仄仂亟亠仍亳仂于舒仆亳 于亠于仂亰仄仂亢仆 亳仗仂于 丶个. 仆亟亳于亳亟舒仍仆仂亠 亰舒亟舒仆亳亠 仂仆仂于舒仆仂 仆舒 从仂仆从亠仆仂仄 仗亳仄亠仆亠仆亳亳 亳仂于仂亶 亳仍舒亳亳 于 舒亟亳仂 亳仍亳 于亳亟亠仂 仗亳仍仂亢亠仆亳.
- 5. 弌仂亟亠亢舒仆亳亠 仂亠舒 亳仍仆亶 仍亳 亰舒亟舒仆亳亠 仆舒 仗仂亠从亳仂于舒仆亳亠 仍亳亳仆亞 仗仂亞舒仄仄, 于仗亳舒仆仆亠 于 仆亠亞仂 舒于仂从仂亶 从仂仄仄亠仆舒亳亳 从 从仂仄舒仆亟舒仄 仗仂亞舒仄仄 亠亰仍舒 仗仂亠从亳仂于舒仆亳: - 舒弍仍亳 从仂亳亳亠仆舒 仗亠亠亟舒亳, - 亞舒亳从 丼丱 - 亞舒亳从 个丼丱 - 亞舒亳从 舒仗仂仍仂亢亠仆亳 仆仍亠亶 亳 仗仂仍ﾑ仂于 - 亞舒亳从 亠舒从亳亳 丶个 仆舒 亳仄仗仍仆仂亠 于仂亰亟亠亶于亳亠 - 从仆亠 亠仄 亳仆亠亰亳仂于舒仆仆 丶个 于于仂亟.
- 6. 乂 弌亊 乘 弌丼丐 1.舒亰舒弍仂从舒 丱-亳仍舒 仗仂仄仂 MATLAB 亠仂亟 于亰于亠亳于舒仆亳. 亅舒仗 于亳仍亠仆亳 从仂亳亳亠仆仂于 舒仆亟舒仆仂亞仂 舒仂仆仂-亳亰弍亳舒亠仍仆仂亞仂 丱-亳仍舒 仍亳仆亠亶仆仂亶 舒亰仂于仂亶 舒舒从亠亳亳从仂亶 仗仂仄仂 仄亠仂亟舒 于亰于亠亳于舒仆亳 仄仂亢仆仂 仗仂磲仂亳 仍亠亟ﾑ亳仄 仂弍舒亰仂仄. 1.舒亟舒 仆亠仂弍仂亟亳仄 舒仂仆 舒舒从亠亳亳从. 2.弍舒 于亠仂于 仆从亳 亳 仂亠仆亳 亳仍仂 从仂亳亳亠仆仂于 亳仍舒, N. 3.仂仍亳 亳亟亠舒仍仆 舒仂仆 舒舒从亠亳亳从, hD(n) (亠亠仆仆 亟仂 N 亰仆舒亠仆亳亶). 4.仂仍亳 N 从仂亳亳亠仆仂于 于亠仂于仂亶 仆从亳亳, (仗). 5.仂仍亳 从仂亳亳亠仆 丱-亳仍舒, 于仂亰亟亠亶于仂于舒于 仆舒 舒仂仆 舒舒从亠亳亳从 于亠仂于仂亶 仆从亳亠亶, h(n) = hD(n) (n). 仂仄 仍舒亠 亳仗仂仍亰亠 舒从仂亶 亳仆舒从亳: b = fir1(N-1, Fc, '亳仗 亳仍舒') b = fir1(N-1, Fc, window) b = fir1(N -1, Fc, '亳仗 亳仍舒', window) MATLAB 仗仂亟亟亠亢亳于舒亠 亳仗仂仍亰仂于舒仆亳亠 舒亰仍亳仆 于亠仂于 仆从亳亶, 于从仍ﾑ舒 仆从-亳亳 丱仄仄亳仆亞舒 (Hamming), 丱亠仆仆亳仆亞舒 (Hanning), 仗礆仂亞仂仍仆 (boxcar), 舒亶亰亠舒 (Kai-ser) 亳 丼亠弍亠于舒 (Chebyshev). w = boxcar(N) w = blackman(N) w = hamming(N) w = hanning(N) w = kaiser(N, beta)
- 7. 亳仄亠 1 亳仍亠仆亳亠 从仂亳亳亠仆仂于 丱-亳仍舒 亳仗仂仍亰仂于舒仆亳亠仄 于亠仂于仂亶 仆从亳亳 舒亶亰亠舒. 仗亠亟亠仍亳亠 从仂亳亳亠仆 亳 亳亰仂弍舒亰亳亠 舒仄仗仍亳亟仆仂-舒仂仆 舒-舒从亠亳亳从 仗仂仍仂仂于仂亞仂 丱-亳仍舒, 亳仗仂仍亰 于亠仂于 仆从亳 舒亶亰亠舒 亳 MAT-LAB. 个亳仍 亟仂仍亢亠仆 亟仂于仍亠于仂 仍亠亟ﾑ亳仄 仗亠亳亳从舒亳礆: 仗仂仍仂舒 仗仂仗从舒仆亳 150-250 , 亳亳仆舒 仗仂仍仂 仗亠亠仂亟舒 50 , 仆亠舒于仆仂仄亠仆仂 于 仗仂仍仂亠 仗仂仗从舒仆亳 0,1 亟, 亰舒舒仆亳亠 于 仗仂仍仂亠 仗仂亟舒于仍亠仆亳 60 亟, 舒仂舒 亟亳从亠亳亰舒亳亳 1 从. 仂亳亳亠仆 亳仍舒 丼丱 亳仍舒 n h(n) n h(n) n h(n) n h(n) 0 -0,0001 10 -0, 0005 20 0,0069 30 0,0094 -0,0004 11 -0, 0044 21 0,0189 31 0,0856 2 -0,0001 12 -0,0022 22 0,0029 32 0,0453 3 -0,0001 13 0,0069 23 0, 0044 33 -0,1665 4 -0,0007 14 0,0066 24 0,0188 34 -0,2066 5 0, 0005 15 -0,0016 25 -0,0125 35 0,0891 6 0, 0023 16 0, 0000 26 -0,0520 36 0, 2998 7 0, 0008 17 0,0022 27 -0,0165 8 -0,0017 18 -0,0117 28 0,0333 9 -0,0005 19 -0,0164 29 0,0104 1
- 8. 仂亞舒仄仄舒 1. 亠舒仍亳亰舒亳 亠亠仆亳 仗亳仄亠舒 1 于 仂仄亠 m-舒亶仍舒 MATLAB FS=1000; % 丼舒仂舒 亟亳从亠亳亰舒亳亳 FN=FS/2; % 丼舒仂舒 舒亶从于亳舒 N=73; % 仍亳仆舒 亳仍舒 beta=5.65; % 舒舒仄亠 仆亠舒于仆仂仄亠仆仂亳 仆从亳舒亶亰亠舒 fc1=125/FN; % 仂仄亳仂于舒仆仆亠 舒仂 亠亰舒 fc2=275/FN; FC=[fcl fc2]; % 亠从仂 从舒亠于 舒仂 hn=fir1(N-l,FC,kaiser(N,beta)); % 仂仍亳 从仂亳亳亠仆 亳仍舒 [H,f]=freqz(hn,1,512,FS); % 亳仍亳 舒仂仆 舒舒从亠亳亳从 mag=20*loglO(abs(H)); plot(f,mag), grid on xlabel ('丼舒仂舒 ()) ylabel('仄仗仍亳亟仆舒 舒舒从亠亳亳从舒 (亟)')
- 9. 仗亳仄亳亰舒亳仂仆仆亠 仄亠仂亟 弌亠亟于仂 Signal Processing Toolbox 于 MATLAB 仂亟亠亢亳 仆亠从仂仍从仂 仗仂亞舒仄仄 舒亰足舒弍仂从亳 亳 仆从亳亶 亟仍 仂亰亟舒仆亳 仂仗亳仄舒仍仆 丱足亳仍仂于 仆舒 仂仆仂于亠 舒仍亞仂亳仄仂于 舒从舒足舒从足 仍亳仍仍舒仆舒 亳 亠仄亠亰舒. 仆仂于仆仂亶 从仂仄舒仆亟仂亶 亟仍 于亳仍亠仆亳 从仂亳亳亠仆仂于 仗仂仄仂 仂仗亳仄舒仍仆仂亞仂 仄亠仂亟舒 磦仍磳 remez. 仂仄舒仆亟舒 仄仂亢亠 亳仗仂仍亰仂于舒 亟仍 舒亰舒弍仂从亳 仄仆仂亞仂仗仂仍仂仆 丱足亳仍仂于 仍亳仆亠亶仆仂亶 舒亰仂于仂亶 舒舒从亠亳亳从仂亶. 舒仆亟舒仆仂亶 仂仄亠 仂仆舒 亳仄亠亠 仍亠亟ﾑ亳亶 亳仆舒从亳: b = remez(N - 1,F,M), 亞亟亠 N 亟仍亳仆舒 亳仍舒, F 于亠从仂 仆仂仄亳仂于舒仆仆 亞舒仆亳仆 舒仂, 舒 于亠从仂 亢亠仍舒亠仄仂亶 舒仄仗仍亳亟仆仂亶 舒舒从亠亳亳从亳 亳仍舒 仆舒 亰舒亟舒仆仆 亞舒仆亳仆 舒仂舒. 舒足 仆亳仆亠 舒仂 仆仂仄亳仂于舒仆 仆舒 仗仂仍仂于亳仆 舒仂 亟亳从亠亳亰舒亳亳 亳 仍亠亢舒 于 亟亳舒仗舒亰仂仆亠 仂 0 亟仂 1 (舒仂舒 舒亶从于亳舒 仂仂于亠于亠 1).
- 10. 亳仄亠 2 仗仂仍亰 仂仗亳仄舒仍仆亶 仄亠仂亟, 于亳仍亳亠 从仂亳亳亠仆 亳 亳亰仂弍舒亰亳亠 舒仂仆 舒舒从亠亳亳从 仗仂仍仂仂于仂亞仂 丱足亳仍舒 仍亳仆亠亶仆仂亶 舒亰仂于仂亶 舒舒从亠亳亳从仂亶, 亟仂足于仍亠于仂ﾑ亠亞仂 仍亠亟ﾑ亳仄 仗亠亳亳从舒亳礆: 仗仂仍仂舒 仗仂仗从舒仆亳 1000足1500 , 仗仂仍仂舒 仗亠亠仂亟舒 500 , 亟仍亳仆舒 亳仍舒 41, 舒仂舒 亟亳从亠亳亰舒亳亳 10 000 . 仂亞舒仄仄舒 2. 亠舒仍亳亰舒亳 于亳仍亠仆亳 从仂亳亳亠仆仂于 仂仗亳仄舒仍仆仂亞仂 丱-亳仍舒 亳 于于仂亟舒 仆舒 从舒仆 舒仂仆仂亶 舒舒从亠亳亳从亳 于 仂仄亠 m舒亶仍舒 MATLAB Fs=1000; % 丼舒仂舒 亟亳从亠亳亰舒亳亳 N=41; % 仍亳仆舒 亳仍舒 =[0 0 1 1 0 0]; % 亠仍舒亠仄舒 舒仄仗仍亳亟仆舒 舒舒从亠亳亳从舒 F=[0, 0.1, 0.2, 0.3, 0.4, 1]; % 舒 仗仂仍仂 b=remez(N-l, F, ); % 亳仍亳 从仂亳亳亠仆 亳仍舒 [, f]=freqz(b, 1, 512, Fs); % 亳仍亳 舒仂仆 舒舒从亠亳亳从 mag=20*logl0(abs()); % 亳仍舒 亳 仆舒亳仂于舒 亠亠 plot(f ,mag) xlabel ('Frequency (Hz)) ylabel ('Magnitude (dB))
- 11. 仂亳亳亠仆 亳仍舒 丼丱 亳仍舒 n h(n) n h(n) 0 -0,0001 11 -0,0044 1 -0,0004 12 -0, 0022 2 -0,0001 13 0,0069 3 -0,0001 14 0,0066 4 -0,0007 15 -0,0016 5 0,0005 16 0,0000 6 0,0023 17 0,0022 7 0,0008 18 -0,0117 8 -0,0017 19 -0,0164 9 -0,0005 20 0,0069 10 -0,0005 21 0,0189
- 12. 亠仂亟 舒仂仆仂亶 于弍仂从亳. 仍 舒亰舒弍仂从亳 丱足亳仍仂于 仗仂亳亰于仂仍仆仄亳 舒仂仆仄亳 舒舒从亠亳亳从舒仄亳, 仗仂足亟仂弍仆 亳仍舒仄, 舒仄仂亠仆仆仄 仗亳 亳仗仂仍亰仂于舒仆亳亳 仄亠仂亟舒 舒仂仆仂亶 于弍仂从亳, 仗亠亟足仆舒亰仆舒亠仆舒 从仂仄舒仆亟舒 fir2. 弌亳仆舒从亳 舒仆亟舒仆仂亶 从仂仄舒仆亟: b = fir2(N-1,F,H) 仂仄舒仆亟舒 fir2 于亳仍磳 从仂亳亳亠仆 丱足 亳仍舒 亟仍亳仆 N. 亠从仂 F 亰舒亟舒亠 仆仂仄亳仂于舒仆仆亠 舒仂 于 亳仆亠于舒仍亠 仂 0 亟仂 1 (仗亳亠仄 舒仂, 从舒从 亳 舒仆亠亠, 仆仂仄亳仂于舒足仆 仆舒 仗仂仍仂于亳仆 舒仂 亟亳从亠亳亰舒亳亳). 亠从仂 H 仂仗亠亟亠仍磳 亢亠仍舒亠仄 舒仄仗仍亳亟仆 舒舒从亠亳亳从 于 仂从舒, 亰舒亟舒仆仆 于 F. 弍舒 于亠从仂舒 亟仂仍亢仆 亳仄亠 仂亟亳仆舒从仂于 亟仍亳仆.
- 13. 2. 舒亰舒弍仂从舒 丱-亳仍仂于 仗仂仄仂 MATLAB 弌亠亟于仂 MATLAB Signal Processing Toolbox 仗亠亟仍舒亞舒亠 仄仆仂亞仂 仗仂仍亠亰仆 仆从亳亶 亟仍 舒亰舒弍仂从亳 亳 舒仆舒仍亳亰舒 从仍舒亳亠从亳 丱-亳仍仂于 (仆舒仗亳仄亠, 舒亠于仂舒, 丼亠弍亠于舒 亳仗舒 I 亳 II, 仍仍亳仗亳亠从亳 亳仍仂于) 亟仍 亟舒仆仆仂亞仂 仆舒弍仂舒 仗亠亳亳从舒亳亶 (仆舒仗亳-仄亠, 亞舒仆亳仆亠 舒仂 仗仂仍仂 仗仂仗从舒仆亳 亳 仗仂亟舒于仍亠仆亳, 仆亠舒于仆仂仄亠仆仂 于 仗仂仍仂亠 仗仂仗从舒仆亳 亳 亰舒舒仆亳亠 于 仗仂仍仂亠 仗仂亟舒于仍亠仆亳). 舒仆仂亳, Toolbox 仗亠亟仍舒亞舒亠 仆从亳亳 仗亠仂弍舒亰仂于舒仆亳 从仍舒亳亠从亳 舒仆舒仍仂亞仂于 亳仍仂于 于 从于亳于舒仍亠仆仆亠 亳仍 亟亳从亠仆仂亞仂 于亠仄亠仆亳. 舒仗仂仄仆亳仄, 仂 从亳亳仆仄 舒仗仂仄 舒亰舒弍仂从亳 亳仂于仂亞仂 丱-亳仍舒 磦仍磳 舒亠 从仂亳亳亠仆仂于. 仍 从仍舒亳亠从亳 亳仂于 丱-亳仍仂于 仄仂亢仆仂 仍亠亟ﾑ亳仄 仂弍舒亰仂仄 仗仂亟仂亢亳 亟亠亶于亳 仆舒 仂仄 舒仗亠. 1.舒亟舒 亟亳从亠仆亶 亳仍. 2.仗亠亟亠仍亳 仗仂亟仂亟亳亶 舒仆舒仍仂亞仂于亶 亳仍 仆亳亢仆亳 舒仂, 仆舒仗亳仄亠, 舒舒从亠亳-亳从仂亶 舒亠于仂舒, 丼亠弍亠于舒 亳仗舒 I, 丼亠弍亠于舒 亳仗舒 II 亳仍亳 仍仍亳仗亳亠从仂亶 舒舒从-亠亳亳从仂亶. 3.亠仂弍舒亰仂于舒 舒仆舒仍仂亞仂于亶 亳仍-仗仂仂亳仗 于 亳仍 仆亳亢仆亳 舒仂, 于亠仆亳 舒-仂, 仗仂仍仂仂于仂亶 亳仍亳 亠亢亠从仂仆亶. 4.亠亠于亠亳 仗亠仂弍舒亰仂于舒仆仆亶 亳仍 于 从于亳于舒仍亠仆仆亶 亳仍 亟亳从亠仆仂亞仂 于亠仄亠仆亳 (仆舒仗亳仄亠, 亳仗仂仍亰 仄亠仂亟 亳仆于舒亳舒仆仆仂亞仂 仗亠仂弍舒亰仂于舒仆亳 亳仄仗仍仆仂亶 舒舒从亠亳-亳从亳 亳仍亳 弍亳仍亳仆亠亶仆仂亞仂 z-仗亠仂弍舒亰仂于舒仆亳).
- 14. 亳仄亠 3 舒亰舒弍仂从舒 仗仂仂亞仂 亳仍舒 仆亳亢仆亳 舒仂 亳仗仂仍亰仂于舒仆亳亠仄 弍亳仍亳仆亠亶仆仂亞仂 z仗亠仂弍舒亰仂于舒仆亳 亳 MATLAB. 丐亠弍亠 亳仂于仂亶 丱-亳仍 舒舒从亠亳亳从仂亶 舒亠于仂舒, 亟仂于仍亠于仂ﾑ亳亶 仍亠亟ﾑ亳仄 仗亠亳亳从舒亳礆: 舒仂舒 亠亰舒 150 , 舒仂舒 亟亳从亠亳亰舒亳亳 1,28 从, 仗仂磲仂从 亳仍舒, N 2. 亠亠仆亳亠. 弌仂仂于亠于ﾑ亳亶 m-舒亶仍 MATLAB 仗亳于亠亟亠仆 于 仗仂亞舒仄仄亠 3 . 1. 弌 仗仂仄仂 m-舒亶仍舒 仗仂仍舒亠仄 仍亠亟ﾑ亳亠 于亠从仂 从仂亳亳亠仆仂于 (b 亳 a) 仆仍亠亶 亳 仗仂仍ﾑ仂于 (z 亳 ) 丱-亳仍舒: b= [0.0878, 0.1756, 0.0878] a= [1.0000, -1.0048, 0.3561] z =[-1,-1] = [0.5024 賊 0.3220i] k = 0.0878 仗仂仍亰 从仂亳亳亠仆, 亰舒仗亳于舒亠仄 仗亠亠亟舒仂仆 仆从亳: H ( z) = 0, 0878(1 + 2 z 1 + z 2 ) 1 1, 0308 z 1 + 0,3553 z 2
- 15. 仂亞舒仄仄舒 3. 亠亠仆亳亠 仗亳仄亠舒 9 于 仂仄亠 m-舒亶仍舒 MATLAB % 仂仂亶 亳仍 仆亳亢仆亳 舒仂 % N=2; % 仂磲仂从 亳仍舒 Fs=1280; % 丼舒仂舒 亟亳从亠亳亰舒亳亳 FN=Fs/2; fc=150; % 丼舒仂舒 亠亰舒 Fc=fc/FN; % 仂仄亳仂于舒仆仆舒 舒仂舒 亠亰舒 [b, a]=butter (N,Fc); % 弌仂亰亟舒 亳 仂亳仂于舒 舒仆舒仍仂亞仂于亶 亳仍 [z, , k]=butter (N, Fc); subplot (2,1,1) % 于亠亳 仆舒 从舒仆 舒仄仗仍亳亟仆仂舒仂仆 舒舒从亠亳亳从 [, f]=freqz(b, 舒, 512, Fs); plot(f, abs(H)) x label ('Frequency (Hz)[1]) ylabel('Magnitude Response (dB)') subplot(2,1,2) % 于亠亳 仆舒 从舒仆 亟亳舒亞舒仄仄 仆仍亠亶 亳 仗仂仍ﾑ仂于 zplane(b, 舒) 仄仗仍亳亟仆仂-舒仂仆舒 舒舒从亠亳亳从舒 亳 亟亳舒亞舒仄仄舒 仆仍亠亶 亳 仗仂仍ﾑ仂于