matlab (한국어)

MATLAB 도구 상자의 freqz에있는 유한 시퀀스의 DTFT가 DTFT를 얻는 방법이라고 생각합니다. 순서의. freqz는 다음 주파수 응답을 계산할 수 있습니다.

H (z) = (Num) / (Den)

다음과 같이 유한 시퀀스 x(n)의 z 변환을 쉽게 계산할 수 있습니다.

H (z) = x (0) z ^ 0 + x (1) z ^ 1 + …

위의 식에서 Den는 1이라는 것을 알고 있습니다.

상기 표현 : freeqz(num,den,n)는 n 지점에서 단계 응답을 제공합니다. x는 x (n)의 벡터가되므로

[x1freqz, x1freqzw]=freqz(x,1,3000,"whole"); 

반드시 DTFT를 제공해야합니다.

1) 그렇습니까 (위 문장) 옳은? 다항식을 이동하면 어떻게 될까요 ?? 왜?

두 번째 방법은 다음과 같이 DTFT 공식을 완전히 계산하는 것입니다.

[X, W]=me_dtft(x1",pi,3000); figure title("my") % plot(W/pi,20*log10(abs(X))); plot(W/pi,abs(X)) ax = gca; % ax.YLim = [-40 70]; xlabel("Normalized Frequency (\times\pi rad/sample)") ylabel("Magnitude (dB)") function [X, w]=me_dtft(x,whalfrange, nsample) w= linspace(-whalfrange,whalfrange,nsample); t=0:1:size(x,2)-1; X=zeros(1,size(w,2)); for i=1:1:size(w,2) X(i)=x*exp(-t*1i*w(i))"; end end 

2) 혼란 스럽습니다. 위 코드의 매개 변수 t, 중요합니까?

3)이 구현이 올바른가요? 그 이유는 무엇입니까?

그림에서 불협화음이있을 것 같습니다. 우리에게 잘못된 것을 말하고 있습니다! 이 변환은 순수 사인파 (그림의 코드)에서 가져옵니다. 오른쪽에서 fft, freqz 매너는 위에서 설명한대로 상단과 (왼쪽) DTFT를 설명했습니다.

“Jason R” “의 설명 후 편집 : 좋습니다. 또한 대수 스케일을 제거했습니다. 헷갈 리기 때문입니다. 그 후 직관적으로 thy are 비슷하게, 다음 이미지에서 볼 수 있지만 정확히 동일하지 않은 이유는 무엇입니까 (로그 스케일로 마지막 이미지 참조?)?

freqz :

[x1freqz, x1freqzw]=freqz(fliplr(XX"),1,3000,"whole"); figure title("freqz") % plot((x1freqzw/pi)-1,20*log10(abs(fftshift(x1freqz)))) plot((x1freqzw/pi)-1,abs(fftshift(x1freqz))) ax = gca; % ax.YLim = [-40 70]; ax.XTick = -1:.5:2; xlabel("Normalized Frequency (\times\pi rad/sample)") ylabel("Magnitude") 

사인 샘플 :

Fs=1000; Ts=1/Fs; time=0:Ts:1; Freqs=500; Xs=zeros(length(Freqs),length(time)); for i=1:length(Freqs) Xs(i,:)= cos(2*pi*Freqs(i)*time); end XX=Xs; XX=XX./ max(abs(XX)); figure;plot(time, XX); axis(([0 time(end) -1 1])); xlabel("Time (sec)"); ylabel("Singal Amp."); title("A sample audio signal"); sound(XX,Fs) 

댓글

• 사용자가 작성한 플롯을 생성하는 완전한 스크립트를 게시해야합니다.
• 한 가지 차이점은 코드에서 DTFT를 -fs / 2에서 + fs / 2로 플로팅하는 반면 FFT는 0에서 fs로 비교한다는 것입니다. " fft shitt ") 모든 코드를 보지 않고 두 플롯의 차이는 더 많은 점을 포함하거나 배정 밀도 부동 소수점을 사용하여 줄일 수있는 반올림 오류라고 생각합니다. DTFT는 샘플링 속도의 정수 하위 배수 인 주파수로 사인파를 생성하지 않는 한 FFT가 샘플이 될 Sinc 함수입니다. 또한 FFT를 제로 패딩하여 DTFT 샘플을 얻을 수 있습니다.
• @DanBoschen (fft 그림, 1000은 fs / 2를 의미합니까?)
• 1000이 없으면 fs, 파형 is at fs / 2
• @DanBoschen 왜 그것들 (내 DTFT와 freqz, 내 말은)이 정확히 같지 않은지 설명해 주시겠습니까? fft를 채우면 sinc가 나타 납니까?