Digital Signal Processing using matlab ED2/Ingle,Proakis/3장 연습문제 solution

소개글

Digital Signal Processing using matlab ED2/Ingle,Proakis/3장 연습문제 solution에 대한 보고서 자료입니다.

본문내용

|');
>> title('Magnitude Response of h(n)=(0.9)^{|n|}');
>> subplot(2,1,2); plot(w/pi,phaH);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase Response of h(n)=(0.9)^{|n|}');
3.11 (3)
>> [h1,n1]=stepseq(0,0,40);
>> [h2,n2]=stepseq(40,0,40);
>> [h3,n3]=sigadd(h1,n1,-h2,n2);
>> n=n3;
>> h=sinc(0.2*n).*h3;
>> w=[-300:300]*pi/300;
>> H=dtft(h,n,w);
>> magH=abs(H); phaH=angle(H)*180/pi;
>> subplot(2,1,1); plot(w/pi,magH);
>> xlabel('w/pi'); ylabel('|H|');
>> title('Magnitude Response of h(n)');
>> subplot(2,1,2); plot(w/pi,phaH);
>> axis([-1 1 -180 180]);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase Response of h(n)');
3.11 (5)
>> w=[-300:300]*pi/300;
>>H=0.5*0.75*ones(size(w))./(1.25-cos(w-(0.1*pi)))+0.5*0.75*ones(size(w))./(1.25-cos(w+(0.1*pi)));
>> magH=abs(H); phaH=angle(H)*180/pi;
>> subplot(2,1,1); plot(w/pi,magH);
>> xlabel('w/pi'); ylabel('|H|');
>> title('Magnitude Response');
>> subplot(2,1,2); plot(w/pi,phaH);
>> axis([-1 1 -180 180]);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase Response');
3.14 (2)
>> n=[0:40];
>> alpha=20; wc=0.5*pi;
>> fc=wc/(2*pi); h=2*fc*sinc(2*fc*(n-alpha));
>> Hs=stem(n,h);
>> xlabel('n'); ylabel('h(n)');
>> title('Truncated Impulse Response h(n)');
3.14 (3)
>> n=[0:40];
>> alpha=20; wc=0.5*pi;
>> fc=wc/(2*pi); h=2*fc*sinc(2*fc*(n-alpha));
>> K=500; w=[-K:K]*pi/K;
>> H=dtft(h,n,w);
>> magH=abs(H); phaH=angle(H);
>> H_d=zeros(1,length(w));
>> H_d(K/2+1:3*K/2+1)=exp(-j*alpha*w(K/2+1:3*K/2+1));
>> magH_d=abs(H_d); phaH_d=angle(H_d);
>> subplot(2,2,1); plot(w/pi,magH);
>> xlabel('w/pi'); ylabel('|H|');
>> title('Magnitude of H(e^{jw})');
>> subplot(2,2,2); plot(w/pi,phaH*180/pi);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase of H(e^{jw})');
>> subplot(2,2,3); plot(w/pi,magH_d);
>> axis([-1 1 0 1.2]);
>> xlabel('w/pi'); ylabel('|H_d|');
>> title('Magnitude of H_d(e^{jw})');
>> subplot(2,2,4); plot(w/pi,phaH_d*180/pi);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase of H_d(e^{jw})');
3.17 (1)
>> w=[-300:300]*pi/300;
>> a=[1]; b=[0.2 0.2 0.2 0.2 0.2];
>> [H]=freqresp(b,a,w);
>> magH=abs(H); phaH=angle(H)*180/pi;
>> subplot(2,1,1); plot(w/pi,magH);
>> xlabel('w/pi'); ylabel('|H|');
>> title('Magnitude Response');
>> subplot(2,1,2); plot(w/pi,phaH);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase Response');
3.17 (3)
>> w=[-300:300]*pi/300;
>> a=[1 -0.95 0.9025]; b=[1 -1 1];
>> [H]=freqresp(b,a,w);
>> magH=abs(H); phaH=angle(H)*180/pi;
>> subplot(2,1,1); plot(w/pi,magH);
>> xlabel('w/pi'); ylabel('|H|');
>> title('Magnitude Response');
>> subplot(2,1,2); plot(w/pi,phaH);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase Response');
3.17 (5)
>> w=[-300:300]*pi/300;
>> l=[0:5]; a=0.5.^l; b=[1];
>> H=freqresp(b,a,w);
>> magH=abs(H); phaH=angle(H)*180/pi;
>> subplot(2,1,1); plot(w/pi,magH);
>> xlabel('w/pi'); ylabel('|H|');
>> title('Magnitude Response');
>> subplot(2,1,2); plot(w/pi,phaH);
>> axis([-1 1 -180 180]);
>> xlabel('w/pi'); ylabel('Degrees');
>> title('Phase Response');

키워드

DSP ,   Digital Signal Processing using matlab ED2,   Ingle,   Proakis,   matlab