Matlab相位截斷誤差仿真

Matlab仿真誤差頻譜仿真:N=14 M=4l=8k=1程序：clear;Fs=1;%采樣頻率1MHzM=4;%截取累加器的高M位；N=14;%累加器的位數(shù)；1=8;k=1;K=2八1*(2*k+1);L=2人(N-M-l);pe=100000;n=1:pe;pp=pe+1;%總的采樣點數(shù)add_y=(n*K)/2八N;error=(round(n*(2*k+1))/L)/(2八M);NFFT=2八nextpow2(pp);%設(shè)定頻譜分析點數(shù)f=Fs/2*1inspace(0,1,NFFT/2+1);%只顯示0?Fs/2的頻率rom_y=sin(2*pi*error);%有相位截斷誤差Y=fft(rom_y,NFFT)/pp;%NFFT點的FFT變換va1ue_Y=abs(Y(1:NFFT/2+1));%求FFT變換值的絕對值DB_Y=20*1og10(va1ue_Y);%求功率的分貝p1ot(f,va1ue_Y);%繪制頻譜圖，只有相位截斷誤差axis([0100.5]);N=14M=4l=7k=4程序：clear;Fs=1;%采樣頻率1MHzM=4;%截取累加器的高M位；N=14;%累加器的位數(shù)；l=7;k=4;K=2八l*(2*k+1);%頻率控制字L=2人(N-M-l);pe=100000;n=1:pe;pp=pe+1;%總的采樣點數(shù)add_y=(n*K)/2八N;error=(round(n*(2*k+1))/L)/(2八M);NFFT=2八nextpow2(pp);%設(shè)定頻譜分析點數(shù)f=Fs/2*linspace(0,1,NFFT/2+1);%只顯示0?Fs/2的頻率rom_y=sin(2*pi*error);%有相位截斷誤差Y=fft(rom_y,NFFT)/pp;%NFFT點的FFT變換value_Y=abs(Y(1:NFFT/2+1));%求FFT變換值的絕對值DB_Y=20*log10(value_Y);%求功率的分貝plot(f,value_Y);%繪制頻譜圖，只有相位截斷誤差axis([0100.5]);%包括坐標軸范圍，axis([xminxmaxyminymax])N=14M=4l=9k=0程序：clear;Fs=1;%采樣頻率1MHzM=4;%截取累加器的高M位；N=14;%累加器的位數(shù)；l=9;k=0;%頻率控制字，即累加的步長；K=2八l*(2*k+1);L=2人(N-M-l);m=1:L-1pe=100000;n=1:pe;pp=pe+1;%總的采樣點數(shù)error=(2*m*pi)/2八(N-l);NFFT=2八nextpow2(pp);%設(shè)定頻譜分析點數(shù)f=Fs/2*linspace(0,1,NFFT/2+1);%只顯示0?Fs/2的頻率rom_y=sin(2*pi*error);%有相位截斷誤差Y=fft(rom_y,NFFT)/pp;%NFFT點的FFT變換value_Y=abs(Y(1:NFFT/2+1));%求FFT變換值的絕對值DB_Y=20*log10(value_Y);%求功率的分貝pl；t(f,value_Y);%繪制頻譜圖沒有幅度量化誤差，只有相位截斷誤差axis([0100.5]);N=14M=5 l=10 k=0程序：clear;Fs=1;%采樣頻率1MHzM=5;%截取累加器的高M位；N=14;%累加器的位數(shù)；l=10;k=0;%頻率控制字，即累加的步長；K=2八l*(2*k+1);L=2人(N-M-l);pe=100000;n=1:pe;pp=pe+1;%總的采樣點數(shù)add_y=(n*K)/2八N;NFFT=2八nextpow2(pp);%設(shè)定頻譜分析點數(shù)f=Fs/2*linspace(0,1,NFFT/2+1);%只顯示0?Fs/2的頻率rom_y=sin(2*pi*add_y);%有相位截斷誤差Y=fft(rom_y,NFFT)/pp;%NFFT點的FFT變換value_Y=abs(Y(1:NFFT/2+1));%求FFT變換值的絕對值DB_Y=20*log10(value_Y);%求功率的分貝plOt(f,value_Y);%繪制頻譜圖沒有幅度量化誤差，只有相位截斷誤差axis([0100.5]);N=14M=5l=10k=1程序：clear;Fs=1;%采樣頻率1MHzM=5;%截取累加器的高M位；N=14;%累加器的位數(shù)；l=10;k=1;%頻率控制字，即累加的步長；K=2八l*(2*k+1);L=2人(N-M-l);pe=100000;n=1:pe;pp=pe+1;%總的采樣點數(shù)add_y=(n*K)/2八N;NFFT=2八nextpow2(pp);%設(shè)定頻譜分析點數(shù)f=Fs/2*linspace(0,1,NFFT/2+1);%只顯示0?Fs/2的頻率rom_y=sin(2*pi*add_y);%有相位截斷誤差Y=fft(rom_y,NFFT)/pp;%NFFT點的FFT變換value_Y=abs(Y(1:NFFT/2+1));%求FFT變換值的絕對值DB_Y=20*log10(value_Y);%求功率的分貝plOt(f,value_Y);%繪制頻譜圖沒有幅度量化誤差，只有相位截斷誤差axis([0100.5]);時域波形仿真:N=14M=4k=1l=8程序：function[y,t]=dds_matlab(N,M,Fs,endtime)k=1;l=8;delta_F=Fs/2八N;disp(['ThefrequencyresolutionofthisDDSis:'...,num2str(delta_F)'Hz']);disp(['Outputfrequenciesare:0Hz,',num2str(delta_F),'Hz,'...num2str(2*delta_F)'Hz, ,',num2str(Fs/2),'Hz']);%ConstructtheSineWaveLUTstoreinarrayLUTt=2*pi*(0:2八N-1)/2八N;LUT=sin(t);t=0:1/Fs:endtime;S=length(t);n=1;y=zeros(1,S);IND=zeros(1,S);%InputFrequencyWordK=2八l*(2*k+1);index=0;whilen<=SIND(n)=index;index=index+K;index=mod(index,2八N-2八(N-M));n=n+1;endIND=IND+1;y=LUT(IND);plot(t,y);figure(1);axis([0endtime-1.21.2]);gridon

N=14M=8k=4l=7程序：function[y,t]=dds_matlab(N,M,Fs,endtime)k=4;l=7;delta_F=Fs/2八N;disp(['ThefrequencyresolutionofthisDDSis:'...,num2str(delta_F)'Hz']);disp(['Outputfrequenciesare:0Hz,',num2str(delta_F),'Hz,'...num2str(2*delta_F)'Hz, ,',num2str(Fs/2),'Hz']);%ConstructtheSineWaveLUTstoreinarrayLUTt=2*pi*(0:2八N-1)/2八N;LUT=sin(t);

:■I'0f000l乍乍口樣■［擋HI一理口<<UOPT?([2,I2,I-9UIT4PU90])STX^?(I)1(Ay4)40yd-(QNI)im=A?I+QNI=aNIpus/T+U=U-((W-N)vS-NvSyxepuT)poui=xepuT/H+X9pUT=X9pUT/X9pUT=(U)QNIS=〉uSfTqM10=xepuT-(1+習(xí)￥乙)^TvS=MpjOM z^ndui%?(S7I)SOJ9Z=QNI1(SyI)sojsz=A?I=u-(3)UKU9[=S/9UIT4PU9:S^/T:0=4function[y,t]=dds_matlab(N,M,Fs,endtime)k=0;l=9;delta_F=Fs/2八N;disp(['ThefrequencyresolutionofthisDDSis:'...,num2str(delta_F)'Hz']);disp(['Outputfrequenciesare:0Hz,',num2str(delta_F),'Hz,'...num2str(2*delta_F)'Hz, ,',num2str(Fs/2),'Hz']);%ConstructtheSineWaveLUTstoreinarrayLUTt=2*pi*(0:2八N-1)/2八N;LUT=sin(t);t=0:1/Fs:endtime;S=length(t);n=1;y=zeros(1,S);IND=zeros(1,S);%InputFrequencyWordK=2八l*(2*k+1);index=0;whilen<=SIND(n)=index;index=index+K;index=mod(index,2八N-2八(N-M));n=n+1;endIND=IND+1;y=LUT(IND);plot(t,y);figure(1);axis([0endtime-1.21.2]);gridon

N=14M=15k=0 l=10程序：function[y,t]=dds_matlab(N,M,Fs,endtime)k=0;l=10;delta_F=Fs/2八N;disp(['ThefrequencyresolutionofthisDDSis:'...,num2str(delta_F)'Hz']);disp(['Outputfrequenciesare:0Hz,',num2str(delta_F),'Hz,'...num2str(2*delta_F)'Hz, ,',num2str(Fs/2),'Hz']);%ConstructtheSineWaveLUTstoreinarrayLUTt=2*pi*(0:2八N-1)/2八N;LUT=sin(t);t=0:1/Fs:endtime;

ZHOOS■■^HIOSSI'0fsHSE0I90'0fzH0：3IEsaiDuanbajj日式日。

zHS￡0I90'0：ETSdd￡邛11ZHOOSLQT%(Wl迫也)哩p警WPP<<_H■0口土4%每悠七V）翌F1& （M）aw cd苜工①y輯向尊回爵警UOPT?([2,I2,I-9UIT4PU90])STX^?(I)9Jh6tj/(Ay4)40yd-(QNI)im=A?I+QNI=aNIpus/T+U=U/(NvSyxepuT)poui=xepuT

/H+X9pUT=X9pUT

/X9pUT=(U)QNIS=〉uSfTqM/0=xepuT-(1+習(xí)￥乙)^TvS=MpjOMAouenbej^z^ndui%?(S7I)SOJ9Z=QNI1(SyI)sojsz=A?I=u-(3)UKU9[=SN=14M=5k=1l=10程序：function[y,t]=dds_matlab(N,M,Fs,endtime)k=1;l=10;delta_F=Fs/2八N;disp(['ThefrequencyresolutionofthisDDSis:'...,num2str(delta_F)'Hz']);disp(['Outputfrequenciesare:0Hz,',num2str(delta_F),'Hz,'...num2str(2*delta_F)'Hz, ,',num2str(Fs/2),'Hz']);%ConstructtheSineWaveLUTstoreinarrayLUTt=2*pi*(0:2八N-1)/2八N;LUT=sin(t);t=0:1/Fs:endtime;S=length(t);n=1;y=zeros(1,S);IND=zeros(1,S);%InputFrequencyWordK=2八l*(2*k+1);index=0;whi