驱动信号源芯片的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;%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
