dds.m

驱动信号源芯片的Matlab系统建模及频谱分析：

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%驱动信号源芯片的Matlab系统建模及频谱分析：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;
N=21;             %累加器的位数；
K=fix((2^N)*0.222);%频率控制字，即累加的步长；
M=11;             %截取累加器的高10位； 
D=9;              %8bit DAC
g=gcd(K,2^N);
pe=2^N/(g);
n=1:pe;
pp=pe+1;          %频谱分析点数
add_y=mod(n*K,2^N);  %累加器的输出表达式；
error=mod(n*K,2^(N-M));
rom_y=cos(2*pi*(add_y-error)/(2^N));     %相位到幅度映射的输出；
pat = - 1 + 1/ (2^ (D - 1) ):1/ (2^ (D - 1) ):1 - 1/ (2^ (D -1) ); %设定量化区间
codebook = - 1 + 1/ (2^D):1/ (2^(D - 1) ) :1 - 1/ (2^D); %设定量化码本值
[Pyy,wy]=periodogram(rom_y,[],'onesided',pp);
figure; %输出波形（未加低通滤波器）
plot(n,rom_y);
axis([0 2000 -1.2 1.2]); %%一个周期2*PI是3490个抽样点
% figure;%figure2 只有相位截断，没有量化 figure1
% psdplot(Pyy/max(Pyy),wy);%      归一化显示出图；
% rom_x=cos(2*pi*add_y/(2^N));
% [index,quants2] = quantiz(rom_x,pat,codebook);%量化
% [Pzz,wz]=periodogram(quants2,[],'onesided',pp);
% figure;%没有相位截取和ＤＡＣ的非线性等，ｆｉｇｕｒｅ2,只有量化
% psdplot(Pzz/max(Pzz),wz);%      归一化显示出图；
% [index,quants1] = quantiz(rom_y,pat,codebook);%量化
% [Pxx,wx]=periodogram(quants1,[],'onesided',pp);    %频谱分析；
% figure;       %ｆｉｇｕｒｅ3 相位截断加量化
% psdplot(Pxx/max(Pxx),wx);%     归一化显示出图；
% grid on