五邑大学matlab实验报告

1.Conv_m.m:

function[y,ny]=conv_m(_,n_,h,nh)

nyb=n_(1)+nh(1);nye=n_(length(_))+nh(length(h));ny=[nyb:nye];

y=conv(_,h);

2.impseq.m:

function[_,n]=impseq(n0,n1,n2)

n=[n1:n2];_=[(n-n0)==0];

3.stepseq.m:

function[_,n]=stepseq(n0,n1,n2)

n=[n1:n2];_=[(n-n0)>=0];

2.38(1):

n=[-4:3];

_=2impseq(-3,-4,3)-impseq(-2,-4,3)+2impseq(0,-4,3)+4impseq(-1,-4,3);

subplot(2,2,1);stem(n,_);

_label('n');ylabel('_(n)');

2.38(2)

n=[0:12];

_=(0.8).^n.(stepseq(0,0,12)-stepseq(10,0,12));stem(n,_);

_label('n');ylabel('_(n)');

2.38(3)

n=[0:200];

_=5cos(0.04pin)+0.3randn(size(n));

subplot(2,1,1);plot(n,_);_label('n');ylabel('_(n)');

subplot(2,1,2);plot(n,_);_label('n');ylabel('_(n)');

2.38(4)

n=[-10:0.5:10];blta=-0.08+0.3j;

_=e_p(bltan);

subplot(2,2,1);stem(n,real(_));title('实部');_label('n');subplot(2,2,2);stem(n,imag(_));title('虚部');_label('n');subplot(2,2,3);stem(n,abs(_));title('幅度');_label('n');

subplot(2,2,4);stem(n,(180/pi)angle(_));title('相位');_label('n');

2.40

_=[1,2,3,4,5];n_=[0:4];

h=[1,-2,1,3];nh=[0:3];

[y,ny]=conv_m(_,n_,h,nh)

n=length(ny);_1=zeros(1,n);h1=zeros(1,n);

_1(find((ny>=min(n_))(ny<=ma_(n_))==1))=_;

h1(find((ny>=min(nh))(ny<=ma_(nh))==1))=h;

subplot(3,1,1);stem(ny,_1);_label('n');ylabel('_(n)');

subplot(3,1,2);stem(ny,h1);_label('n');ylabel('h(n)');

subplot(3,1,3);stem(ny,y);_label('n');ylabel('h(n)');

2.42

b=[1];a=[1,-1,0.9];

_=impseq(0,-20,100);n=[-20:100];

h=filter(b,a,_);

subplot(2,1,1);stem(n,h);title('冲激响应');

_label('n');ylabel('h(n)');

_=stepseq(0,-20,100);

s=filter(b,a,_);

subplot(2,1,2);stem(n,s);title('阶跃响应');

_label('n');ylabel('s(n)');

z=roots(a);magz=abs(z)

2.45(3)

n=0:7;

k=0:500;w=(pi/500)k;

_=(e_p(-j0.3pi/500)).^(n'k);

mag_=abs(_);ang_=angle(_);

subplot(2,1,1);plot(w/pi,mag_);title('幅度响应');grid;ylabel('幅度');_label('以pi为单位的频率');

subplot(2,1,2);plot(w/pi,ang_);title('相位响应');grid;ylabel('相位/pi');_label('以pi为单位的频率');

2.47(3)

b=[1,-1,0.5];a=[1,0.25,0.125];

[R,P,C]=residuez(b,a);

disp('R=');disp(R');

disp('P=');disp(P');

disp('C=');disp(C');

Dft.m:

function[_k]=dft(_n,N)

n=[0:1:N-1];

k=[0:1:N-1];

WN=e_p(-j2pi/N);

nk=n'k;

WNnk=WN.^nk;

_k=_nWNnk;

Idft.m:

function[_n]=idft(_k,N)

n=[0:1:N-1];

k=[0:1:N-1];

WN=e_p(-j2pi/N);

nk=n'k;

WNnk=WN.^(-nk);

_n=(_kWNnk)/N;

3.22(1)

_=[1,1,1,1];n=0:3;

k=-200:200;

w=(pi/100)k;

_=_(e_p(-jpi/100)).^(n'k);

mag_=abs(_);ang_=angle(_)180/pi;

figure(1)

subplot(2,1,1);plot(w/pi,mag_);

%a_is([_min_ma_yminyma_])

a_is([-1104]);grid;

title('DTFT的幅度');_label('以pi为单位的频率');ylabel('幅度');subplot(2,1,2);plot(w/pi,ang_);grid;

a_is([-11-200200]);title('DTFT的相位');

_label('以pi为单位的频率');ylabel('相位');

N=4;

_=dft(_,N);

mag_=abs(_);pha_=angle(_)180/pi

k=0:3;

figure(2)

subplot(2,1,1);stem(k,mag_);

title('DFT的幅度');_label('k');

subplot(2,1,2);stem(k,pha_);

title('DFT的相位');_label('k');

3.22(2)

_=[1,1,1,1,0,0,0,0];n=0:7;

k=-200:200;

w=(pi/100)k;

_=_(e_p(-jpi/100)).^(n'k);

mag_=abs(_);ang_=angle(_)180/pi;subplot(2,1,1);stem(k,mag_);title('DFT的幅度');_label('k');subplot(2,1,2);stem(k,pha_);title('DFT的相位');_label('k');

3.23

_=[1,1,1,1,0,0,0,0];n=0:7;

k=-200:200;

w=(pi/100)k;

_=_(e_p(-jpi/100)).^(n'k);

mag_=abs(_);ang_=angle(_)180/pi;figure(1)

subplot(2,1,1);plot(w/pi,mag_);N=8;

_=dft(_,N);

mag_=abs(_);pha_=angle(_)180/pik=0:7;

subplot(2,1,1);stem(k,mag_);title('DFT的幅度');_label('k');subplot(2,1,2);stem(k,pha_);title('DFT的相位');_label('k');

3.29

_=[1,1,1,1,0,0,0,0];n=0:7;

k=-200:200;

w=(pi/100)k;

_=_(e_p(-jpi/100)).^(n'k);

mag_=abs(_);ang_=angle(_)180/pi;figure(1)

subplot(2,1,1);plot(w/pi,mag_);N=8;

_=dft(_,N);

mag_=abs(_);pha_=angle(_)180/pik=0:7;

subplot(2,1,1);stem(k,mag_);title('DFT的幅度');_label('k');subplot(2,1,2);stem(k,pha_);title('DFT的相位');_label('k');