【正文】 plot(f1,abs(fftshift(N1)))。xlabel(39。頻率39。)。ylabel(39。幅值39。)。axis([500,500,0,])。title(39。抽樣正常時的信號頻譜39。)。axis([500,500,])。t0=10。 %信號持續(xù)的時間ts2=。 %不滿足抽樣條件的抽樣間隔fs2=1/ts2。t2=[t0/2:ts2:t0/2]。 %定義不滿足抽樣條件的時間序列x2=sin(200*t2)。 m2=x2./(200*t2+eps)。 w2=t0/(2*ts2)+1。 m2(w2)=1。 %修正 t=0 時的信號值m2=m2.*m2。 %定義信號[M2,mn2,df2]=fft_seq(m2,ts2,df)。%對不滿足抽樣條件的信號進行傅立葉變換M2=M2/fs2。N2=[M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2]。f2=[7*df2*length(mn2):df2:6*df2*length(mn2)df2]fs2/2。figure(3) subplot(2,1,1)。 stem(t2,m2)。xlabel(39。時間39。)。ylabel(39。幅值39。)。title(39。抽樣失真(fs=100Hz)時的信號波形 39。)。axis([,0,1])。subplot(2,1,2)plot(f2,abs(fftshift(N2)))。xlabel(39。頻率39。)。ylabel(39。幅值39。)。axis([500,500,0,])。title(39。抽樣失真時的信號頻譜39。)。axis([500,500,])。function [M,m,df]=fft_seq(m,ts,df)fs=1/ts。if nargin==2 n1=0else n1=fs/dfendn2=length(m)。n=2^(max(nextpow2(n1),nextpow2(n2)))。M=fft(m,n)。m=[m,zeros(1,nn2)]。df=fs/n程序 2：function average()t=[0::4*pi]。山東輕工業(yè)學院 2022 屆本科生畢業(yè)設(shè)計（論文）27y=sin(t)。w=jylh(y,1,64)。subplot(2,1,1)。plot(t,y)。xlabel(39。時間39。)。ylabel(39。幅度39。)。axis([0,4*pi,])。title(39。原始信號 39。)。subplot(2,1,2)。plot(t,w)。xlabel(39。時間39。)。ylabel(39。幅度39。)。axis([0,4*pi,])。title(39。均勻量化后的信號39。)。function h=jylh(f,V,L)n=length(f)。t=2*V/L。p=zeros(1,L+1)。for i=1:L+1,p(i)=V+(i1)*t。endfor i=1:n if f(i)V,h(i)=V。end if f(i)=V,h(i)=V。end flag=0。 for j=2:L/2+1 if(flag==0) if(f(i)p(j)) h(i)=p(j1)。 flag=1。 end。 end。 end。 for j=L/2+2:L+1 if(flag==0) if(f(i)p(j)) h(i)=p(j)。 flag=1。山東輕工業(yè)學院 2022 屆本科生畢業(yè)設(shè)計（論文）28 end end endend nq=V^2/(3*L^2)。 程序 3：function a_quantize()t=0::。y=sin(8000*pi*t)。figure subplot(2,1,1)plot(t,y)axis([0 ])xlabel(39。時間39。) ylabel(39。幅度39。)title(39。原始信號 39。)z=a_pcm(y,)。subplot(2,1,2)plot(t,z)axis([0 ])xlabel(39。時間39。) ylabel(39。幅度39。)title(39。A 律量化后的信號39。)function y=a_pcm(x,a)t=1/a。for i=1:length(x) if x(i)=0 if(x(i)=t) y(i)=(a*x(i))/(1+log(a))。 else y(i)=(1+log(a*x(i)))/(1+log(a))。 end else if(x(i)=t)山東輕工業(yè)學院 2022 屆本科生畢業(yè)設(shè)計（論文）29 y(i)=(a*x(i))/(1+log(a))。 else y(i)=(1+log(a*x(i)))/(1+log(a))。 end endend程序 4：function a_13code()t=0::。y=sin(8000*pi*t)z=line13(y)c=pcmcode(z)function y=line13(x)x=x/max(x)。z=sign(x)。x=abs(x)。for i=1:length(x) if((x(i)=0)amp。(x(i)1/64)) y(i)=16*x(i)。 else if((x(i)=1/64)amp。(x(i)1/32)) y(i)=8*x(i)+1/8。 else if((x(i)=1/32)amp。(x(i)1/16)) y(i)=4*x(i)+2/8。 else if((x(i)=1/16)amp。(x(i)1/8)) y(i)=2*x(i)+3/8。 else if((x(i)=1/8)amp。(x(i)1/4)) y(i)=x(i)+4/8。 else if((x(i)=1/4)amp。(x(i)1/2)) y(i)=1/2*x(i)+5/8。山東輕工業(yè)學院 2022 屆本科生畢業(yè)設(shè)計（論文）30 else if((x(i)=1/2)amp。(x(i)=1)) y(i)=1/4*x(i)+6/8。 end end end end end end endendy=z.*y。function f=pcmcode(y)f=zeros(length(y),8)。z=sign(y)。y=y.*128。y=fix(y)。y=abs(y)。for i=1:length(y) if(y(i)==128) y(i)=。 endendfor i=1:length(y) for j=6:1:0 f(i,8j)=fix(y(i)/2^j)。 y(i)=mod(y(i),(2^j))。 endendfor i=1:length(y)。 if(z(i)==1) f(i,1)=0。 else f(i,1)=1。 end山東輕工業(yè)學院 2022 屆本科生畢業(yè)設(shè)計（論文）31end程序運行結(jié)果： a_13codey = Columns 1 through 7 0 Columns 8 through 11 z = Columns 1 through 7 0 Columns 8 through 11 c = 1 0 0 0 0 0 0 0 0 1 1 1 0 0 1 1 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 0 1 1 1 0 0 1 1 0 0 0 0 0 0 0 0 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 0 0 0 0 0 0 0程序 5：Ts=1e3。t=0:Ts:20*Ts。x=sin(2*pi*50*t)+*sin(2*pi*150*t)。山東輕工業(yè)學院 2022 屆本科生畢業(yè)設(shè)計（論文）32delta=。D(1+length(t))=0。for k=1:length(t) e(k)=x(k)D(k)。 e_q(k)=delta*(2*(e(k)=0)1)。 D(k+1)=e_q(k)+D(k)。 codeout(k)=(e_q(k)0)。endsubplot(3,1,1)。plot(t,x,39。o39。)。axis([0 20*Ts,2 2])。hold on。subplot(3,1,2)。stairs(t,codeout)。axis([0 20*Ts,2 2])。Dr(1+length(t))=0。for k=1:length(t) eq(k)=delta*(2*codeout(k)1)。 xr(k)=eq(k)+Dr(k)。 Dr(k+1)=xr(k)。endsubplot(3,1,3)。stairs(t,xr)。hold on。subplot(3,1,3)。plot(t,x