% To generate DM signals with different step sizes, we apply

the

% same signal g(t) as used in the PCM example.

% The MATLAB program ExDM.m applies three step sizes:

% △1=0.2,?△2=2△1, and?△3=4△1.

%?( file name: ExDM.m )

% Example of sampling, quantization, and zero-order hold

clear; clf;

td=0.002;?% original sampling

rate 500 Hz

t=[0:td:1.];?% time interval of 1 second

xsig=sin(2*pi*t)-sin(6*pi*t);?%

1Hz+3Hz sinusoids

Lsig=length(xsig);

ts=0.02;

Nfact=ts/td;

%?send the signal through a

16-level uniform quantizer

Delta=0.2;?% First select a small Delta=0.2 in

DM

s_DMout1=deltamod(xsig,Delta,td,ts);

%?obtained the DM signal

%?plot the original signal and

the DM signal in time domain

figure(1);

subplot(311); sfig1=plot(t,xsig,'r',t,s_DMout1(1:Lsig),'b');

set(sfig1,'Linewidth',2);

title(' Signal {\it g}({\it t}) and DM signal ')

xlabel('time (sec.)'); axis([0 1 -2.2 2.2]);

ylabel('△1=0.2');

%

%?Apply DM again by doubling

the Delta

Delta2=2*Delta;

s_DMout2=deltamod(xsig,Delta2,td,ts);

%?obtained the DM signal

%?plot the original signal and

the DM signal in time domain

subplot(312); sfig2=plot(t,xsig,'r',t,s_DMout2(1:Lsig),'b');

set(sfig2,'Linewidth',2);

title(' Signal {\it g}({\it t}) and DM signal with doubled

stepsize')

xlabel('time (sec.)'); axis([0 1 -2.2 2.2]);

ylabel('△2=2△1');

%

Delta3=2*Delta2;?% Double the DM Delta again

s_DMout3=deltamod(xsig,Delta3,td,ts);

%?obtained the DM signal

%?plot the original signal and

the DM signal in time domain

subplot(313); sfig3=plot(t,xsig,'r',t,s_DMout3(1:Lsig),'b');

set(sfig3,'Linewidth',2);

title(' Signal {\it g}({\it t}) and DM signal with quadrupled

stepsize')

xlabel('time (sec.)'); axis([0 1 -2.2 2.2]); ylabel('△3=4△1');

%

-------------------------------------------------------------

%?程序調(diào)用函數(shù)

%

-------------------------------------------------------------

% Delta Modulation

% We illustrate the practical effect of stepsize selection △

% in the design of DM encoder. The basic function to

% implement DM is given deltamod.m.

% ( file name: deltamod.m )

function s_DMout= deltamod(sig_in,Delta,td,ts)

%?Usage

%?s_DMout = deltamod(xsig,Delta,td,ts)

%?Delta - DM stepsize

%?sig_in - input signal

vector

%?td - original signal sampling

period of sig_in

%?ts - new sampling

period

% NOTE: td*fs must be a positive integer;

%?Function outputs:

%?s_DNout - DM sampled output

if (rem(ts/td,1)==0)

nfac=round(ts/td);

p_zoh=ones(1,nfac);

s_down=downsample(sig_in,nfac);

Num_it=length(s_down);

s_DMout(1)=-Delta/2;

for

k=2:Num_it

xvar=s_DMout(k-1);

s_DMout(k)=xvar+Delta*sign(s_down(k-1)-xvar);

end

s_DMout=kron(s_DMout,p_zoh);

else

warning('error! ts/td is not an integer!');

s_DMout=[];

end

end

%

-------------------------------------------------------------------------------------

% 參考文獻(xiàn)

% B.P.拉茲，丁峙. 現(xiàn)代數(shù)字與模擬通信系統(tǒng)[M]. 第4版. 北京：電子工業(yè)出版社，2011:368-371

總結(jié)

以上是生活随笔為你收集整理的增量调制原理matlab仿真,增量调制（DM  △M）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。