dsp 프로젝트

소개글

dsp 프로젝트에 대한 보고서 자료입니다.

목차

1. spec

2.. 설계 목표

3. 설계 방향 및 방법

4. 설계 과정

5. 성능 평가

6. 고찰

본문내용

것이다.
⑥결정지어진 design spec을 매트랩 함수를 통해
구현하고 확인한다.
Ⅳ 설계 과정(using MATLAB)
①필터 스펙을 확인한다.
②스펙에 맞는 필터 구조를 선택한다.
③매트랩 함수를 통해 필터를 구현한다.
④비교를 위해 Reference filter의 주파수응답, 임펄스 응답,
위상, 양자화에러를 확인한다.
⑤design filter를 설계하고 값을 비교한후 부족한 부분은
튜닝을 통해 성능을 향상시켜나간다.
<소스코드>
① Reference filter
f=[700,1200]; %trasition band: 700~1200
a=[1,0]; %lowpass filter
[n0,f0,a0,w0]=remezord(f,a,[0.1, 0.01],10000);
%optimal n0를 계산
b0=remez(n0,f0,a0); %필터계수를 결정
[h0,w0]=freqz(b0);
figure(1) % #1 ideal filter, design filter plot
plot(f0/2,a0,w0/(2*pi),abs(h0)), grid;
legend('Ideal','firpm Design')
figure(2) % #2 필터 계수 stem
stem([0:n0],b0), grid;
figure(3) %#3 위상 plot
plot(w0/(2*pi),angle(h0)), grid;
%---------------양자화 에러체크--------%
bl=8;
m=1;
d1=abs(b0);
while fix(d1)>0
d1=abs(b)/(2^m);
m=m+1;
end
beq=fix(d1*2^bl);
beq=sign(b0).*beq.*2^(m-bl-1);
[H,w]=freqz(beq);
figure(4)
plot(w/(2*pi),abs(H)), grid;
%------------------------------------%
② Design filter
f=[1000,1200]; %trasition band: 1000~1200
a=[1,0];
[n0,f0,a0,w0]=remezord(f,a,[0.05, 0.005],10000);
% Stopband attenuation:46dB, Passband ripple:0.05dB
b0=remez(n0,f0,a0);
[h0,w0]=freqz(b0);
figure(1)
plot(f0/2,a0,w0/(2*pi),abs(h0)), grid;
legend('Ideal','firpm Design')
figure(2)
stem([0:n0],b0), grid;
figure(3)
plot(w0/(2*pi),angle(h0)), grid;
%---------------양자화 에러체크--------%
bl=8;
m=1;
d1=abs(b0);
while fix(d1)>0
d1=abs(b)/(2^m);
m=m+1;
end
beq=fix(d1*2^bl);
beq=sign(b0).*beq.*2^(m-bl-1);
[H,w]=freqz(beq);
figure(4) %8bit 양자화 plot
plot(f0/2,a0,w/(2*pi),abs(H)), grid;
legend('Ideal','firpm Design')
%----------------------------------%
Ⅴ 성능 평가
1. Reference spec
<크기응답> reference spec을 만족시킨다.
<필터계수>
<위상> 선형위상을 나타낸다.
<8bit 양자화 후> 양자화 에러가 발생하였다.
2. Design spec
<주파수 응답> design spec을 만족시킨다.
<필터계수> 차수가 크게 늘어났다.
<위상> 선형위상을 나타낸다.
<8bit 양자화 후> 양자화 에러가 해결되었다.
Ⅵ 고찰

키워드

dsp,   프로젝트,   선형위상,   양자화,   필터설계,   ripple,   matlab