%%%%調(diào)頻信號的累加器實(shí)現(xiàn)
%%%%Grape
clear all; 
clc;    		%清除所有變量，清屏
M  =  2^12;  		%采樣的點(diǎn)數(shù)
a(1  :  1  :  M)  =  0; %頻率累加器的數(shù)值
b(1 : 1 : M) = 0;  	%相位累加器的數(shù)值
y(1 : 1 : M) = 0;  	%輸出的波形數(shù)據(jù)
K = 400;  		%調(diào)頻斜率
Kc = 10;  		%初始頻率控制字
N = 10;   		%幅度量化位數(shù)
L = 24;   		%相位累加器位數(shù)

%下面一段實(shí)現(xiàn)頻率累加器
a(1) = 0.5 * K;  	%初始頻率步進(jìn)量
for i = 2 : 1 : M
a(i) = a(i - 1) + K; 
end

%下面一段實(shí)現(xiàn)相位累加器
b(1) = Kc + a(1);  	%相位初始值
for i = 2 : 1 : M
b(i) = b(i - 1) + (Kc + a(i)); 
end

%下面一段實(shí)現(xiàn)了查找表ROM以及進(jìn)行幅度量化
for i = 1 : 1 : M
y(i) = floor(2^(N) * cos(2 * pi/(2^L) * b(i))); 
end

%下面一段畫出相應(yīng)的圖形
figure(1); plot(y); axis([0 M, - 2^N - 100 2^N + 100]); 
figure(2); freqz(y)

%下面一段實(shí)現(xiàn)相位累加器
b(1) = Kc + a(1);  	%相位初始值
for i = 2 : 1 : M
b(i) = b(i - 1) + (Kc + a(i)); 
end

%下面一段實(shí)現(xiàn)了查找表ROM以及進(jìn)行幅度量化
for i = 1 : 1 : M
y(i) = floor(2^(N) * cos(2 * pi/(2^L) * b(i))); 
end

%下面一段畫出相應(yīng)的圖形
figure(1); plot(y); axis([0 M, - 2^N - 100 2^N + 100]); 
figure(2); freqz(y)